Finite Volume Subcell Element Methods

• FVSE is a hybrid method that embeds a finite-volume subgrid within high-order DG/DGSEM frameworks to preserve accuracy and deliver robust performance.
• It blends element-local high-order fluxes with low-order finite-volume updates to stabilize shocks, discontinuities, and other nonlinear instabilities.
• The approach enables local activation of subcell corrections while ensuring conservation, entropy control, and positivity across complex geometries.

Searching arXiv for relevant papers on finite volume subcell element and closely related DG/FV subcell methods. Finite Volume Subcell Element (FVSE) denotes a class of hybrid high-order/low-order discretization strategies in which a high-order element-local representation—most commonly a discontinuous Galerkin (DG) or DG spectral element method (DGSEM)—is coupled to a finite-volume (FV) evolution on a finer subgrid of subcells inside each element. In this framework, the element is treated as a container for a conservative subcell finite-volume dynamics, either continuously through flux blending or intermittently through troubled-cell replacement. The defining objective is to preserve the high-order accuracy and subcell resolution of DG in smooth regions while recovering the robustness of FV discretizations near shocks, positivity threats, and other nonlinear instabilities. Although the exact label “Finite Volume Subcell Element” is not used uniformly across the literature, closely aligned formulations include local subcell monolithic DG/FV schemes on unstructured grids (Vilar, 2024), subcell limiting strategies for DGSEM on Legendre–Gauss–Lobatto (LGL) nodes (Rueda-Ramírez et al., 2022), a posteriori MOOD-type subcell limiters for ADER-DG and $P_NP_M$ schemes (Dumbser et al., 2014, Gaburro et al., 2020), positivity-preserving DGSEM/FV blending (Rueda-Ramírez et al., 2021), hybrid artificial-viscosity/FVSE shock capturing in aerospike simulations (Pyle et al., 14 Aug 2025), and h-adaptive FV subcell shock capturing on heterogeneous curvilinear meshes (Schwarz et al., 30 Apr 2026).

1. Conceptual definition and scope

FVSE is best understood as an embedded finite-volume representation inside a higher-order element method. The core idea appears in several closely related forms.

In the local subcell monolithic DG/FV formulation, the DG polynomial on each cell $\omega_c$ is reinterpreted through subcell mean values on a subdivision

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$

with subcell means

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$

The DG evolution is then rewritten as a conservative finite-volume-like update on this subgrid, so that the scheme is conservative at the subcell level (Vilar, 2024). This provides a direct FVSE-type interpretation: the “element” remains the DG cell, but the actual conservative update is expressed on internal finite-volume subcells.

In DGSEM-based work, the same structural idea is realized on LGL nodes. The nodal values $u_{ij}$ are interpreted as subcell-centered values on a local tensor-product subcell grid, and both the high-order DGSEM operator and the low-order subcell FV operator are written in the same conservative flux-difference form,

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$

thereby enabling direct convex blending (Rueda-Ramírez et al., 2022). This suggests that FVSE is not merely a fallback discretization, but a structural reinterpretation of the element-local DG operator itself.

A different but related interpretation arises in a posteriori subcell limiting. In those approaches, the primary solution remains a high-order DG or $P_NP_M$ polynomial on the main mesh, but troubled elements are scattered onto a finer subgrid, recomputed there by a robust FV scheme, and then reconstructed back to the polynomial space (Dumbser et al., 2014, Gaburro et al., 2020). In that setting, FVSE functions as a temporary finite-volume substructure embedded inside the element.

A plausible implication is that “FVSE” is best regarded as a family of subcell finite-volume embeddings rather than a single canonical algorithm.

2. Algebraic structure: from DG elements to subcell finite volumes

A defining feature of FVSE-type methods is that the high-order element discretization is made compatible with a conservative subcell formulation.

For the unstructured-grid monolithic DG/FV method, the starting point is the DG polynomial

$$u_h^c=\sum_{m=1}^{N_k} u_m^c(t)\,\sigma_m^c(\mathbf{x}), \qquad N_k=\frac{(k+1)(k+2)}{2}\ \text{in 2D},$$

together with the weak DG equations

$$\int_{\omega_c} \partial_t u_h^c\,\psi\,dV = \int_{\omega_c} F(u_h^c)\,\nabla\psi\,dV - \int_{\partial\omega_c} \psi\, F_n\, dS, \qquad \forall \psi\in \mathbb{P}^k(\omega_c).$$

After quadrature, this becomes

$M_c U_c = \Phi_c.$

The crucial step is then the construction of DG reconstructed fluxes $\omega_c$0, special subcell-face fluxes that make the finite-volume subgrid evolution algebraically equivalent to the original DG residual (Vilar, 2024). On faces touching neighboring cells, the reconstructed flux is simply the DG numerical flux; on interior subcell faces, it is obtained from

$\omega_c$1

This gives the unique set of interior subcell fluxes that reproduces the DG method exactly, independent of the particular subcell partition, as long as $\omega_c$2 is invertible (Vilar, 2024).

In the DGSEM setting, compatibility is achieved through a flux-differencing representation. The high-order DGSEM operator and the low-order subcell FV operator share the same stencil structure, with the high-order “fluxes” constructed recursively from volume terms, for example

$\omega_c$3

Because both operators are written on the same LGL-based subcell layout, convex mixing can be performed without interpolation to a separate mesh (Rueda-Ramírez et al., 2022).

On curvilinear meshes, compatibility additionally requires carefully designed metric terms. The subcell metric expressions

$\omega_c$4

satisfy the discrete subcell metric identities

$\omega_c$5

which guarantees local conservation and free-stream preservation (Rueda-Ramírez et al., 2022).

This body of work establishes the central algebraic point: FVSE is not a heuristic patch applied after discretization, but a conservative subcell reformulation or compatible embedding of the element-local operator.

3. Flux blending, switching, and limiting strategies

The most characteristic operational mechanism in FVSE-type methods is the local combination of a high-order DG-like update and a more robust low-order FV update.

In the monolithic DG/FV scheme, each face of each subcell carries two fluxes: a high-order DG reconstructed flux and a first-order FV flux

$\omega_c$6

The actual numerical flux is a convex blend

$\omega_c$7

Here $\omega_c$8 corresponds to the pure high-order DG reconstructed flux, $\omega_c$9 to the pure first-order FV flux, and intermediate values to a local blend (Vilar, 2024).

In DGSEM formulations, the same idea appears at either the element scale or the subcell-interface scale. Element-wise blending uses a single coefficient,

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$0

whereas subcell-wise blending blends the interface fluxes directly,

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$1

Subcell-wise blending is more localized and usually less dissipative, while element-wise blending is simpler and is the case where stronger entropy-dissipation results are currently available (Rueda-Ramírez et al., 2022).

Troubled-cell activation is implemented in several distinct ways across the literature.

A priori activation uses indicators or sensors before the update. In DGSEM subcell limiting, one option is a modal-energy detector based on

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$2

followed by a smooth mapping to $\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$3 (Rueda-Ramírez et al., 2022). In the heterogeneous-mesh DGSEM/FV scheme, switching is driven by

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$4

with hysteresis between lower and upper thresholds (Schwarz et al., 30 Apr 2026).

A posteriori activation computes a candidate high-order solution first and then checks admissibility. The ADER-DG MOOD limiter computes an unlimited candidate $\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$5, applies PAD and NAD criteria, and only if the cell fails is the solution scattered to subcells and recomputed by a robust FV method (Dumbser et al., 2014). The corresponding $\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$6 formulation follows the same logic with physical admissibility and a relaxed discrete maximum principle (RDMP) on subcell averages (Gaburro et al., 2020).

A plausible implication is that the distinction between “FVSE by blending” and “FVSE by replacement” is operational rather than conceptual: both rely on the same embedded subcell finite-volume representation.

4. Admissibility, positivity, and entropy control

FVSE-type methods are used because high-order DG discretizations can generate spurious oscillations, negative density or pressure, or code crashes in under-resolved or shock-dominated regimes. The subcell finite-volume component provides a robust admissible state toward which the solution can be pulled.

In the positivity-preserving DGSEM limiter for the Euler equations, the blended semi-discrete operator is

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$7

The limiter enforces

$\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$8

relative to a safe FV stage state obtained with $\omega_c = \bigcup_{m=1}^{N_s} S_m^c,$9 (Rueda-Ramírez et al., 2021). The paper uses $u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$0 in the numerical tests, and the correction is applied at the end of each RK stage. Density correction is algebraic, whereas pressure correction requires solving

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$1

with Newton’s method (Rueda-Ramírez et al., 2021). This is a strictly local, conservative fallback toward a positivity-safe subcell FV update.

The MOOD-style ADER-DG limiter also combines physical admissibility detection with a relaxed discrete maximum principle. For Euler, PAD requires

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$2

while NAD is imposed through

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$3

with

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$4

(Dumbser et al., 2014). The $u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$5 extension uses RDMP on subcell averages with

$u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$6

where $u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$7 and $u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$8 (Gaburro et al., 2020).

Entropy control has been analyzed most explicitly in the monolithic subcell DG/FV framework. Three levels are distinguished (Vilar, 2024).

1. Discrete subcell entropy inequality for any entropy gives the strongest guarantee, but forces the method toward first order.
2. Semi-discrete subcell entropy stability for one entropy improves accuracy to about second order, but does not preserve arbitrary high order.
3. Semi-discrete cell entropy stability for one entropy enforces entropy inequality only on the whole cell $u_m^c=\frac{1}{|S_m^c|}\int_{S_m^c} u_h^c\,dV.$9, leading to a continuous knapsack problem

$u_{ij}$0

solved greedily, and preserving high-order accuracy in smooth regions because

$u_{ij}$1

so the blended scheme reduces to pure DG up to

$u_{ij}$2

The paper also notes that this high-order-preserving entropy option requires modified FV fluxes evaluated on entropy-variable states and may destroy other properties such as positivity (Vilar, 2024).

The literature therefore treats positivity, invariant-domain preservation, maximum principles, and entropy stability as related but nonidentical admissibility targets. One common misconception is that “entropy stable” automatically means positivity preserving; the cited work explicitly separates these properties (Rueda-Ramírez et al., 2021, Vilar, 2024).

5. Subcell resolution, time stepping, and reconstruction

A recurrent design problem in FVSE-type schemes is how many subcells to place inside each element. The dominant choice in the a posteriori subcell-limiting literature is $u_{ij}$3 subcells per spatial dimension.

In the ADER-DG subcell limiter, troubled cells are scattered onto $u_{ij}$4 Cartesian subcells, with

$u_{ij}$5

This is described as optimal because it matches the maximum admissible time step of the subcell FV scheme to the maximum admissible time step of the DG scheme: $u_{ij}$6 for DG and

$u_{ij}$7

for the subcell FV method (Dumbser et al., 2014). The same rationale appears in the $u_{ij}$8 limiter, where the subgrid contains

$u_{ij}$9

subcells, chosen precisely so that the subcell FV scheme can run with the same timestep already admissible for the main discretization (Gaburro et al., 2020).

The heterogeneous-mesh DGSEM/FV scheme allows an arbitrary subcell resolution satisfying

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$0

but adopts

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$1

in the reported work (Schwarz et al., 30 Apr 2026). The paper states that at least $$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$2 subcells are needed for robust and accurate shock capturing, and the upper bound is used to avoid a stronger CFL restriction than DG.

Projection and reconstruction between the high-order element representation and the subcell FV state are equally central. In ADER-DG, scattering is performed by exact subcell averaging,

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$3

defining the projection operator $$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$4, and reconstruction is defined so that the reconstructed polynomial matches the subcell averages and preserves the whole-cell integral. When $$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$5, the reconstruction is well posed and

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$6

(Dumbser et al., 2014).

In the $$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$7 setting, projection and reconstruction satisfy

$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$8

and troubled cells may be recomputed with either second-order MUSCL-Hancock TVD FV or ADER-WENO FV on the subgrid (Gaburro et al., 2020). In the heterogeneous-mesh method, DG-to-FV and FV-to-DG transfer operators are written explicitly with conservation constraints, using least squares for reconstruction (Schwarz et al., 30 Apr 2026).

This suggests that subcell resolution is not only a geometric refinement device; it is also an algebraic mechanism for retaining the information content of the element-local polynomial while enabling robust low-order evolution.

6. Variants, applications, and practical regimes

The FVSE concept spans several numerical roles, from shock capturing to invariant-domain enforcement and hybrid stabilization.

Representative variants

Variant	Operational mode	Representative source
Monolithic local DG/FV on subcells	Continuous convex blending of DG reconstructed fluxes and FV fluxes on each subcell face	(Vilar, 2024)
DGSEM subcell limiting on LGL nodes	Element-wise or subcell-wise convex blending with compatible FV subcell operators	(Rueda-Ramírez et al., 2022)
A posteriori MOOD-type subcell limiter	Troubled-cell recomputation on subcells, then reconstruction back to DG/$$J_{ij} \dot{u}_{ij} = \frac{1}{\omega_i}\left( f_{(i-1,i)j} - f_{(i,i+1)j} \right) + \frac{1}{\omega_j}\left( f_{i(j-1,j)} - f_{i(j,j+1)} \right),$$9	(Dumbser et al., 2014, Gaburro et al., 2020)
Positivity-preserving DGSEM/FV blend	Element-local a posteriori blending toward a positivity-safe FV state	(Rueda-Ramírez et al., 2021)
Hybrid AV+FVSE shock capturing	AV as primary stabilizer, FVSE as localized fallback “parachute”	(Pyle et al., 14 Aug 2025)
h-adaptive mixed-element FV subcells	DGSEM/FV switching on hexahedra, prisms, tetrahedra, and pyramids	(Schwarz et al., 30 Apr 2026)

In practical shock capturing, the aerospike-nozzle solver uses a hybrid artificial-viscosity/FVSE strategy. There, artificial viscosity is the primary stabilizer and FVSE is invoked only where AV cannot sufficiently smooth the solution gradients. The FVSE residual is blended element-wise,

$P_NP_M$0

with a low-order finite-volume residual based on Lax–Friedrichs fluxes between subcells associated with DG nodal points (Pyle et al., 14 Aug 2025). The paper explicitly describes FVSE as a “parachute” in the context of numerical stability. It also emphasizes that broad use of FVSE is highly dissipative: a constant $P_NP_M$1 everywhere suppresses major wake dynamics, whereas the preferred “CM5” configuration $P_NP_M$2 uses FVSE only in extreme-gradient regions and captures the unsteady wake while preserving stability (Pyle et al., 14 Aug 2025).

In the positivity-preserving DGSEM limiter, the FV subcell method is not primarily a shock-capturing device, but a positivity-safe correction mechanism for density and pressure in under-resolved vortex-dominated and shock-dominated Euler flows (Rueda-Ramírez et al., 2021).

In the heterogeneous-mesh DGSEM/FV method, the important extension is topological generality. The mesh may contain hexahedra, prisms, tetrahedra, and pyramids, with non-hexahedral elements handled via collapsed coordinate transformations. The FV operator is a second-order TVD finite-volume scheme with generalized minmod or Barth–Jesperson limiting, and the paper verifies conservation, spatial convergence, and shock-capturing capability on curved mixed meshes, including a NACA 0012 airfoil case using a mesh with 1168 prisms and 12032 hexahedra (Schwarz et al., 30 Apr 2026).

Across these applications, the same practical pattern recurs: DG or DGSEM provides the high-fidelity baseline in smooth flow, while FVSE is activated locally to stabilize shocks, strong gradients, positivity threats, or difficult startup transients.

7. Relation to neighboring methods and major interpretive issues

FVSE belongs to the broader family of DG/FV hybridizations, but it is distinct from several nearby ideas.

It is not simply classical artificial viscosity. Artificial viscosity adds local diffusion inside the element and can degrade sub-element resolution; the cited mixed-mesh DGSEM/FV work contrasts this explicitly with subcell limiting, which preserves sub-element resolution better by redistributing information onto a fine internal FV grid (Schwarz et al., 30 Apr 2026). The aerospike solver likewise treats AV and FVSE as complementary rather than interchangeable: AV is computationally efficient and less intrusive in smooth regions, whereas FVSE is a more dissipative fallback (Pyle et al., 14 Aug 2025).

It is not standard finite volume. In all cited formulations, the method is not globally first-order FV. Rather, it retains a high-order DG or DGSEM representation almost everywhere and uses the subcell FV structure locally, either continuously or intermittently (Rueda-Ramírez et al., 2022, Rueda-Ramírez et al., 2021).

It is not identical to classical finite volume element methods in the older FEM sense. The literature summarized here uses “finite volume subcell element” to denote a high-order element method with an internal finite-volume substructure, not a separate finite-volume discretization on arbitrary control volumes. This distinction is especially clear in the DGSEM-based work, where the FV scheme is explicitly co-located with the DG nodal grid (Rueda-Ramírez et al., 2021, Rueda-Ramírez et al., 2022).

Several interpretive issues recur in the literature.

One is the tradeoff between localization and theory. Subcell-wise blending is more localized and less dissipative than element-wise blending, but a full entropy-stability proof for subcell-wise blending is still missing in the DGSEM framework (Rueda-Ramírez et al., 2022). Another is the tradeoff between entropy control and high-order accuracy. The monolithic unstructured-grid analysis shows that the strongest entropy condition drives the method essentially to first order, whereas cell-level semi-discrete entropy stability for one entropy can preserve high order but at increased complexity and possible loss of positivity (Vilar, 2024).

A further issue concerns whether FVSE should be viewed primarily as a limiter or as a reformulation. The unstructured-grid monolithic scheme and DGSEM flux-differencing formulations support the latter view, because the high-order operator itself is recast as a subcell finite-volume update (Vilar, 2024, Rueda-Ramírez et al., 2022). The MOOD-type and $P_NP_M$3 approaches emphasize the former, because subcells are activated only in troubled regions (Dumbser et al., 2014, Gaburro et al., 2020). A plausible synthesis is that FVSE is both: an algebraic reinterpretation of element-local high-order methods and a practical limiting architecture built on that reinterpretation.

In current usage, FVSE therefore names a robust and technically flexible paradigm for embedding finite-volume subcell dynamics inside high-order element methods. Its enduring significance lies in precisely this dual capacity: preserving DG-level accuracy and subcell resolution where the solution is smooth, while restoring finite-volume robustness where nonlinear dynamics would otherwise make the high-order scheme fail.