Shifted Boundary Method (SBM)

• Shifted Boundary Method (SBM) is an unfitted finite element technique that shifts boundary conditions from curved physical boundaries to nearby surrogate boundaries.
• It employs high-order polynomial corrections to restore optimal convergence rates and simplify simulations without the need for curved meshes.
• SBM integrates seamlessly with existing finite element and finite volume codes, extending its applicability to multiphysics problems, compressible flows, and shock simulations.

The Shifted Boundary Method (SBM) is an unfitted finite element and finite volume technique designed to treat boundary value problems for partial differential equations (PDEs) on complex geometries. SBM "shifts" boundary data from the true, generally curved domain boundary to a nearby surrogate boundary that conforms to an underlying mesh, typically linear, Cartesian, or octree. Boundary conditions at mesh-aligned surrogate boundaries are corrected using high-order polynomial extrapolation based on the solution within the mesh cell, thereby eliminating the need for explicit high-order Taylor expansions or integration over cut cells. This approach delivers high-order accuracy without requiring curved meshes and is compatible with existing finite element and finite volume codes for both continuous and discontinuous Galerkin discretizations. The method is applicable to compressible and incompressible flows, diffusion, elasticity, and multiphysics, including problems with shocks.

1. Geometric Formulation and Closest-Point Maps

SBM operates by replacing the true physical domain $\Omega$ with a surrogate domain $\Omega_h$ composed of mesh cells fully interior to $\Omega$. The surrogate boundary $\Gamma_h = \partial \Omega_h$ generally does not coincide with the true boundary $\Gamma = \partial \Omega$; for each quadrature point $\hat{x} \in \Gamma_h$, the signed distance $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$ is computed, where positive values denote interior points. The outward unit normal $n(x)$ is measured at the closest point $x = \mathcal{R}(\hat{x})$ on $\Gamma$, with the closest-point map $\Omega_h$0 defined as $\Omega_h$1. Maintaining $\Omega_h$2 and $\Omega_h$3 is essential for consistency.

2. Polynomial Correction of Boundary Conditions

For high-order consistency, SBM introduces a correction based on extrapolating the cell-local solution polynomial to the true boundary. On each boundary cell $\Omega_h$4, the degree-$\Omega_h$5 polynomial approximation is

$\Omega_h$6

where $\Omega_h$7 is a basis for $\Omega_h$8.

Rather than explicitly computing derivatives for a Taylor expansion,

$\Omega_h$9

SBM relies on polynomial extrapolation:

$\Omega$0

with $\Omega$1, $\Omega$2. This correction is exact for degree-$\Omega$3 polynomials, requiring only one additional local evaluation per face quadrature point.

3. Numerical Boundary Condition Enforcement

SBM replaces standard discrete boundary states—such as ghost states in finite volume and DG methods—with shifted, corrected states at surrogate boundary points:

• Far-field (Dirichlet): $\Omega$4, applied at each quadrature point $\Omega$5.
• Characteristic far-field: $\Omega$6 is constructed from prescribed Riemann invariants, evaluated at the projected $\Omega$7.
• Slip-wall: For $\Omega$8 on $\Omega$9, the surrogate normal $\Gamma_h = \partial \Omega_h$0 at $\Gamma_h = \partial \Omega_h$1 is decomposed into directions relative to $\Gamma_h = \partial \Omega_h$2; only the normal component is corrected:

  $\Gamma_h = \partial \Omega_h$3

The ghost state $\Gamma_h = \partial \Omega_h$4 combines the density and tangential velocity from $\Gamma_h = \partial \Omega_h$5, the corrected normal velocity $\Gamma_h = \partial \Omega_h$6, and consistent enthalpy values.

The corrected ghost state $\Gamma_h = \partial \Omega_h$7 enters the numerical flux $\Gamma_h = \partial \Omega_h$8 in the boundary DG/FV computations.

4. Algorithmic Implementation

Implementing SBM on a linear mesh with precomputed $\Gamma_h = \partial \Omega_h$9 and $\Gamma = \partial \Omega$0 entails:

1. Identifying all boundary faces and quadrature points $\Gamma = \partial \Omega$1 on $\Gamma = \partial \Omega$2.
2. For each $\Gamma = \partial \Omega$3, compute projected $\Gamma = \partial \Omega$4.
3. Evaluate the cell polynomial $\Gamma = \partial \Omega$5 at $\Gamma = \partial \Omega$6.
4. Compute the correction $\Gamma = \partial \Omega$7.
5. Form the corrected boundary/ghost state $\Gamma = \partial \Omega$8.
6. Use $\Gamma = \partial \Omega$9 in the boundary flux operator.

The additional computational expense is strictly one local polynomial evaluation per boundary quadrature point; required data structures need only store $\hat{x} \in \Gamma_h$0 and $\hat{x} \in \Gamma_h$1 at the quadrature points.

5. Theoretical Analysis and Convergence Rates

Formal analysis (Ciallella et al., 2022) shows that with $\hat{x} \in \Gamma_h$2 boundary approximation (i.e., distance between $\hat{x} \in \Gamma_h$3 and $\hat{x} \in \Gamma_h$4 scales as $\hat{x} \in \Gamma_h$5) and degree-$\hat{x} \in \Gamma_h$6 accuracy in the interior, the SBM polynomial shift restores optimal convergence:

$\hat{x} \in \Gamma_h$7

Numerical rates for Dirichlet BCs on 2D curved domains:

• DG–P2 without SBM: rate $\hat{x} \in \Gamma_h$8, with SBM: rate $\hat{x} \in \Gamma_h$9.
• DG–P3 without SBM: rate $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$0, with SBM: rate $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$1.

In analogous 3D tests (spheres), DG–P3 with SBM achieves $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$2 L₂ convergence.

6. Extension to Shocked Flows

SBM boundary correction applies unchanged in time-dependent ADER-DG frameworks with MOOD limiters for flows with shocks. Boundary cells flagged as troubled and downgraded to subcell FV still use the shifted Dirichlet or slip-wall state for fluxes. For slip-wall treatments, a Rusanov-type penalty term can be added for stability:

$d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$3

without compromising accuracy.

7. Practical Significance and Applicability

The SBM polynomial correction method is compatible with standard finite element and finite volume codes, demanding only minimal extension for data storage of $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$4 and $d(\hat{x}) = \mathrm{dist}(\hat{x}, \Gamma)$5, and a negligible additional computational cost for boundary point evaluation. By eliminating the requirement for curved meshes and explicit high-order Taylor expansions, SBM streamlines high-order simulations of compressible flows in domains with curved boundaries. It yields high-order boundary consistency and optimal convergence for all types of Euler boundary conditions—Dirichlet, characteristic far-field, and slip-wall. This is especially beneficial for problems involving curvature, shocks, and complex geometries (Ciallella et al., 2022).

The SBM has been adopted and validated for a broad class of PDEs, including Poisson, elasticity, Navier–Stokes, and hyperbolic conservation laws. Its formulation supports efficient matrix-free and parallel implementations, robust scaling under geometric complexity, and seamless extension to immersed and adaptive mesh frameworks.

---

Key reference:

"Shifted boundary polynomial corrections for compressible flows: high order on curved domains using linear meshes" (Ciallella et al., 2022)