research · note · in motion
What SBP operators are actually doing at a boundary
A field engineer's read of summation-by-parts: where the dual norm earns its name, and where it doesn't.
The dual norm, where it earns its name
SBP operators give you a discrete summation-by-parts identity. The catch: they only give you that identity in a specific norm. When the scheme is set up properly that norm is also the energy norm of the problem, and a SAT-style boundary penalty term then gives you discrete energy stability for free. When the scheme is not set up properly, the same operator silently forfeits its order at the boundary, and the convergence rate quietly drops.
Where the thread stands
The MEng project this note grew out of is finished — SBP operators, SAT closures, and generalised Gregory quadrature studied on model PDEs, wrapped up at UTIAS in 2025. It closed as a project; there is no manuscript in the works. But the boundary question didn't stop being interesting because the degree ended. It now runs through SBPLab, a private C++ bench for SBP and SBP–SAT operators with high-order time integration, built for constructing operators and checking them against canonical problems where the exact answer is known — the entry case is 2-D diffusion of a single Fourier mode, whose discrete decay rate has an exact eigenvalue to answer to. Same relationship as coursework to CFDLab: the course ends, the questions keep going.
Next
SBPLab's ladder is the usual one: more operator families, more canonical problems with exact answers, and boundary closures checked against the convergence rates they advertise. A private bench with no deadline — the point is keeping the boundary instincts sharp.