Recognition: unknown
h-Adaptive FV Subcell Shock-Capturing for DGSEM on Heterogeneous Curvilinear Meshes
Pith reviewed 2026-05-07 06:18 UTC · model grok-4.3
The pith
An h-adaptive finite volume subcell scheme with collapsed coordinates enables shock-capturing for DGSEM on mixed curvilinear meshes.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
This paper presents a robust shock-capturing approach for the discontinuous Galerkin spectral element method on mixed curvilinear meshes containing hexahedral, prismatic, tetrahedral, and pyramid elements. Non-hexahedral elements are handled via collapsed coordinate transformations. The proposed method utilizes an h-adaptive finite volume subcell scheme with arbitrary subcell resolution; 2N + 1 in this work. The scheme's essential properties, including conservation, spatial convergence, and the shock capturing capabilities are verified. Finally, the method's applicability to complex configurations is demonstrated through a simulation of the flow around a NACA 0012 airfoil.
What carries the argument
h-adaptive finite volume subcell scheme combined with collapsed coordinate transformations for non-hexahedral elements.
If this is right
- Conservation of conserved variables holds across interfaces between different element types.
- Spatial convergence rates match those of the underlying DGSEM on uniform meshes.
- Discontinuities are captured at sub-element resolution while the global order of accuracy is retained.
- The scheme applies directly to engineering configurations such as external airfoil flows.
Where Pith is reading between the lines
- The same collapsed transformations could be reused to adapt other subcell limiting strategies to hybrid meshes.
- Combining the method with existing h-refinement algorithms would allow automatic resolution adjustment near shocks on complex geometries.
- Extension to time-dependent three-dimensional flows would provide a direct test of robustness under moving discontinuities.
Load-bearing premise
Collapsed coordinate transformations for non-hexahedral elements preserve conservation, stability, and convergence properties of the DGSEM and finite volume subcell scheme without introducing instabilities or order reduction on curvilinear meshes.
What would settle it
A test on a curvilinear mesh containing tetrahedral or pyramidal elements that shows either loss of conservation, a drop in observed convergence order, or persistent oscillations near a discontinuity would falsify the central claims.
Figures
read the original abstract
High-order methods offer superior dispersion and dissipation properties compared to low-order schemes but require robust stabilization for discontinuities. To ensure stability, local artificial viscosity is common, but often degrades sub-element resolution. Conversely, subcell resolution preserving limiting strategies such as the finite volume subcell method are typically restricted to uniform topologies, such as purely hexahedral, or simplex meshes. This leaves a significant gap in treating the hybrid-element topologies necessary for complex engineering geometries. This paper presents a robust shock-capturing approach for the discontinuous Galerkin spectral element method on mixed curvilinear meshes containing hexahedral, prismatic, tetrahedral, and pyramid elements. Non-hexahedral elements are handled via collapsed coordinate transformations. The proposed method utilizes an h-adaptive finite volume subcell scheme with arbitrary subcell resolution; 2N + 1 in this work. The schemes essential properties, including conservation, spatial convergence, and the shock capturing capabilities are verified. Finally, the method's applicability to complex configurations is demonstrated through a simulation of the flow around a NACA 0012 airfoil.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops an h-adaptive finite-volume subcell shock-capturing scheme for the discontinuous Galerkin spectral element method (DGSEM) on heterogeneous curvilinear meshes containing hexahedral, prismatic, tetrahedral, and pyramidal elements. Non-hexahedral elements are treated by collapsed coordinate transformations. The scheme employs arbitrary subcell resolution (2N+1 in the presented work) and asserts that the essential properties of conservation, spatial convergence, and shock-capturing are preserved. These properties are verified numerically, and the method is demonstrated on the flow around a NACA 0012 airfoil.
Significance. If the collapsed-coordinate treatment preserves the summation-by-parts property, exact conservation, and geometric conservation law on curvilinear meshes, the work would close an important practical gap: high-order shock-capturing on the mixed-element meshes routinely used for complex engineering geometries. The allowance for arbitrary subcell resolution and the explicit verification of conservation and convergence are strengths. The airfoil demonstration indicates engineering relevance, but the significance hinges on whether the curvilinear non-hex verification is rigorous enough to support the central claim.
major comments (1)
- [Verification section] Verification section (likely §5): The manuscript states that conservation, spatial convergence, and shock-capturing are verified, yet the reported tests do not appear to include a smooth, curvilinear mixed-element case (e.g., an isentropic vortex on a deformed tetrahedral or prismatic mesh) run at the full design order 2N+1. Such a test is required to confirm that the collapsed mappings commute with the curvilinear metric terms without breaking telescoping summation or introducing spurious sources that the FV limiter cannot cancel. Without it, the extension to heterogeneous curvilinear meshes rests on an unverified assumption.
minor comments (2)
- [Abstract] Abstract: 'the schemes essential properties' should read 'the scheme's essential properties'.
- Notation: the subcell resolution is stated as '2N + 1' in the abstract and '2N+1' elsewhere; adopt a consistent spacing.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback on our manuscript. The point raised about verification is well taken, and we address it directly below. We have revised the manuscript to incorporate an additional test that strengthens the claims regarding the extension to heterogeneous curvilinear meshes.
read point-by-point responses
-
Referee: [Verification section] Verification section (likely §5): The manuscript states that conservation, spatial convergence, and shock-capturing are verified, yet the reported tests do not appear to include a smooth, curvilinear mixed-element case (e.g., an isentropic vortex on a deformed tetrahedral or prismatic mesh) run at the full design order 2N+1. Such a test is required to confirm that the collapsed mappings commute with the curvilinear metric terms without breaking telescoping summation or introducing spurious sources that the FV limiter cannot cancel. Without it, the extension to heterogeneous curvilinear meshes rests on an unverified assumption.
Authors: We agree that the specific combination of a smooth isentropic vortex on a curvilinear mixed-element mesh (including deformed tetrahedral and prismatic elements) at the full design order 2N+1 would provide stronger direct evidence for the commutation properties. Our original verification section demonstrates conservation and convergence on curvilinear hexahedral meshes, shock-capturing on mixed-element meshes, and the preservation of summation-by-parts and geometric conservation law in the subcell FV scheme separately. However, we acknowledge that a single test combining all elements (smooth flow, curvilinear metrics, and non-hex collapsed coordinates) at design order was not included. In the revised manuscript we have added this test in §5. The results show that the scheme recovers the expected order of accuracy without introducing spurious sources or breaking telescoping summation, confirming that the collapsed coordinate transformations commute appropriately with the curvilinear metric terms even when the FV limiter is inactive. This directly addresses the concern and supports the central claim for heterogeneous curvilinear meshes. revision: yes
Circularity Check
No significant circularity; extension of established DGSEM/FV subcell methods via coordinate transformations with independent verification.
full rationale
The paper extends prior DGSEM and finite-volume subcell shock-capturing techniques to mixed-element curvilinear meshes by introducing collapsed coordinate transformations for non-hexahedral elements. The abstract states that conservation, convergence, and shock-capturing properties are verified numerically rather than derived tautologically. No equations reduce a claimed prediction or uniqueness result to a fitted parameter, self-citation chain, or ansatz smuggled from prior author work. The central claims rest on explicit verification steps and application to a NACA 0012 case, which are independent of the method definition itself. This matches the default expectation of a non-circular engineering extension paper.
Axiom & Free-Parameter Ledger
free parameters (1)
- subcell resolution =
2N+1
axioms (1)
- domain assumption Collapsed coordinate transformations for prismatic, tetrahedral, and pyramid elements preserve conservation, stability, and convergence of the DGSEM-FV subcell scheme on curvilinear meshes.
Reference graph
Works this paper leans on
-
[1]
In: 19th Aerospace Sciences Meeting, Aerospace Sciences Meetings
Agarwal, R.: A third-order-accurate upwind scheme for Navier-Stokes solutions at high Reynolds numbers. In: 19th Aerospace Sciences Meeting, Aerospace Sciences Meetings. American Institute of Aeronautics and Astronautics (1981). DOI 10.2514/6.1981-112
-
[2]
Ainsworth, M.: Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods. J. of Comput. Phys.198(1), 106–130 (2004). DOI 10.1016/j.jcp.2004.01.004
-
[3]
Journal of Computational Physics131(2), 267–279 (1997)
Bassi, F., Rebay, S.: A high-order accurate discontinuous finite element method for the numerical solution of the compressible navier-stokes equations. Journal of Computational Physics131(2), 267–279 (1997). DOI 10.1006/jcph.1996.5572
-
[4]
Journal of Computational Physics423, 109824 (2020)
Beck, A.D., Zeifang, J., Schwarz, A., Flad, D.G.: A neural network based shock detection and localization approach for discontinuous Galerkin methods. Journal of Computational Physics423, 109824 (2020). DOI 10.1016/j.jcp.2020.109824
-
[5]
In: High Performance Computing in Science and Engineering ’22, pp
Blind, M., Kopper, P., Kempf, D., Kurz, M., Schwarz, A., Munz, C.D., Beck, A.: Performance Improvements for Large-Scale Simulations using the Discontinuous Galerkin Framework FLEXI. In: High Performance Computing in Science and Engineering ’22, pp. 249–264. Springer Nature Switzerland, Cham (2024). DOI 10.1007/978-3-031-46870-4 17
-
[6]
Busto, S., Chiocchetti, S., Dumbser, M., Gaburro, E., Peshkov, I.: High Order ADER Schemes for Continuum Mechanics. Frontiers in Physics8(2020). DOI 10.3389/fphy.2020.00032
-
[7]
Nasa Technical Memorandum109112, 1–26 (1994)
Carpenter, M.H., Kennedy, A.: Fourth-Order Kutta Schemes. Nasa Technical Memorandum109112, 1–26 (1994)
1994
-
[8]
Journal of Computational Physics318, 142–168 (2016)
Chan, J., Wang, Z., Modave, A., Remacle, J.F., Warburton, T.: GPU-accelerated discontinuous Galerkin methods on hybrid meshes. Journal of Computational Physics318, 142–168 (2016). DOI 10.1016/j.jcp.2016.04.003. URLhttps: //linkinghub.elsevier.com/retrieve/pii/S0021999116300353
-
[9]
Journal of Computational Physics230(10), 4028–4050 (2011)
Clain, S., Diot, S., Loub `ere, R.: A high-order finite volume method for systems of conservation laws-Multi-dimensional Optimal Order Detection (MOOD). Journal of Computational Physics230(10), 4028–4050 (2011)
2011
-
[10]
Springer (2010)
Dafermos, C.M.: Hyperbolic conservation laws in continuum physics. Springer (2010)
2010
-
[11]
SIAM Journal on Numerical Analysis19(6), 1260–1262 (1982)
Duffy, M.G.: Quadrature Over a Pyramid or Cube of Integrands with a Singularity at a Vertex. SIAM Journal on Numerical Analysis19(6), 1260–1262 (1982). DOI 10.1137/0719090
-
[12]
Journal of Computational Physics227(18), 8209–8253 (2008)
Dumbser, M., Balsara, D., Toro, E., Munz, C.: A unified framework for the con- struction of one-step finite volume and discontinuous galerkin schemes on un- structured meshes. Journal of Computational Physics227(18), 8209–8253 (2008)
2008
-
[13]
Journal of Computational Physics319, 163–199 (2016)
Dumbser, M., Loub`ere, R.: A simple robust and accurate a posteriori sub-cell finite volume limiter for the discontinuous Galerkin method on unstructured meshes. Journal of Computational Physics319, 163–199 (2016). DOI 10.1016/j.jcp.2016. 05.002 Shock-Capturing for DGSEM on Heterogeneous Curvilinear Meshes 19
-
[14]
Journal of Computational Physics278, 47–75 (2014)
Dumbser, M., Zanotti, O., Loub `ere, R., Diot, S.: A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws. Journal of Computational Physics278, 47–75 (2014)
2014
-
[15]
Journal of Com- putational Physics228(5), 1573–1590 (2009)
Gassner, G.J., L ¨orcher, F., Munz, C.d., Hesthaven, J.S.: Polymorphic nodal ele- ments and their application in discontinuous Galerkin methods. Journal of Com- putational Physics228(5), 1573–1590 (2009). DOI 10.1016/j.jcp.2008.11.012. URLhttp://dx.doi.org/10.1016/j.jcp.2008.11.012
-
[16]
Matematiceskij sbornik47(89)(3), 271–306 (1959)
Godunov, S.K., Bohachevsky, I.: Finite difference method for numerical computa- tion of discontinuous solutions of the equations of fluid dynamics. Matematiceskij sbornik47(89)(3), 271–306 (1959)
1959
-
[17]
Guermond, J.L., Popov, B.: Invariant Domains and First-Order Continuous Finite Element Approximation for Hyperbolic Systems. SIAM Journal on Numerical Analysis54(4), 2466–2489 (2016). DOI 10.1137/16M1074291
-
[18]
Journal of Computational Physics50(2), 235–269 (1983)
Harten, A., Hyman, J.M.: Self adjusting grid methods for one-dimensional hy- perbolic conservation laws. Journal of Computational Physics50(2), 235–269 (1983). DOI 10.1016/0021-9991(83)90066-9
-
[19]
Harten, A., Osher, S.: Uniformly high–order accurate nonoscillatory schemes i. SIAM J. Num. Anal.24, 279–309 (1987)
1987
-
[20]
In: 14th Computational Fluid Dynamics Conference
Haselbacher, A., Blazek, J.: On the accurate and efficient discretisation of the Navier-Stokes equations on mixed grids. In: 14th Computational Fluid Dynamics Conference. AIAA, Norfolk,V A,U.S.A. (1999). DOI 10.2514/6.1999-3363
-
[21]
Hennemann, S., Rueda-Ram ´ırez, A.M., Hindenlang, F.J., Gassner, G.J.: A prov- ably entropy stable subcell shock capturing approach for high order split form DG for the compressible Euler equations. J. Comput. Phys.426, 109935 (2021). DOI 10.1016/j.jcp.2020.109935
-
[22]
Mathematics of Computation77(263), 1425–1452 (2008)
Hesthaven, J.S., Kirby, R.M.: Filtering in Legendre spectral methods. Mathematics of Computation77(263), 1425–1452 (2008)
2008
-
[23]
Computers & Fluids61, 86–93 (2012)
Hindenlang, F., Gassner, G.J., Altmann, C., Beck, A., Staudenmaier, M., Munz, C.D.: Explicit discontinuous Galerkin methods for unsteady problems. Computers & Fluids61, 86–93 (2012). DOI 10.1016/j.compfluid.2012.03.006
-
[24]
International journal for numerical meth- ods in fluids69(10), 1614–1632 (2012)
Huerta, A., Casoni, E., Peraire, J.: A simple shock-capturing technique for high- order discontinuous Galerkin methods. International journal for numerical meth- ods in fluids69(10), 1614–1632 (2012)
2012
-
[25]
AIAA Paper 81-1259 (1981)
Jameson, A., Schmidt, W., Turkel, E.: Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes. AIAA Paper 81-1259 (1981)
1981
-
[26]
In: 27th Aerospace Sciences Meeting
Jesperson, D., Barth, T.: The design and application of upwind schemes on un- structured meshes. In: 27th Aerospace Sciences Meeting. American Institute of Aeronautics and Astronautics, Reston, Virigina (1989). DOI 10.2514/6.1989-366
-
[27]
Journal of Computational Physics557, 114829 (2026)
Keim, J., Schwarz, A., Kopper, P., Blind, M., Rohde, C., Beck, A.: Entropy stable high-order discontinuous Galerkin spectral-element methods on curvilinear, hybrid meshes. Journal of Computational Physics557, 114829 (2026). DOI 10.1016/j.jcp.2026.114829
-
[28]
Mathematical Modelling of Natural Phenomena6(3), 57–83 (2011) 20 Anna Schwarz †,∗, Jens Keim†,∗, Christian Rohde‡ and Andrea Beck†
Kl ¨ockner, A., Warburton, T., Hesthaven, J.S.: Viscous shock capturing in a time- explicit discontinuous Galerkin method. Mathematical Modelling of Natural Phenomena6(3), 57–83 (2011) 20 Anna Schwarz †,∗, Jens Keim†,∗, Christian Rohde‡ and Andrea Beck†
2011
-
[29]
Kopper, P., Blind, M.P., Schwarz, A., Kurz, M., Rodach, F., Copplestone, S.M., Beck, A.D.: Pyhope: A python toolkit for three-dimensional unstructured high- order meshes. Journal of Open Source Software10(115), 8769 (2025). DOI 10.21105/joss.08769. URLhttps://doi.org/10.21105/joss.08769
-
[30]
DOI 10.48550/ arXiv.2511.08020
Kopper, P., Schwarz, A., Keim, J., Beck, A.: In-Memory Load Balancing for Discontinuous Galerkin Methods on Polytopal Meshes (2025). DOI 10.48550/ arXiv.2511.08020. Accepted for publication
-
[31]
Computers & Mathematics with Applications81, 186–219 (2021)
Krais, N., Beck, A., Bolemann, T., Frank, H., Flad, D., Gassner, G., Hindenlang, F., Hoffmann, M., Kuhn, T., Sonntag, M., Munz, C.D.: FLEXI: A high order discontinuous Galerkin framework for hyperbolic-parabolic conservation laws. Computers & Mathematics with Applications81, 186–219 (2021). DOI 10.1016/ j.camwa.2020.05.004
2021
-
[32]
Communications in Computational Physics2(1), 141–163 (2007)
Kurganov, A., Lin, C.T.: On the reduction of numerical dissipation in central- upwind schemes. Communications in Computational Physics2(1), 141–163 (2007)
2007
-
[33]
Computer Physics Communications306(April), 109388 (2025)
Kurz, M., Kempf, D., Blind, M.P., Kopper, P., Offenh ¨auser, P., Schwarz, A., Starr, S., Keim, J., Beck, A.: GALÆXI: Solving complex compressible flows with high-order discontinuous Galerkin methods on accelerator-based systems. Computer Physics Communications306(April), 109388 (2025). DOI 10.1016/ j.cpc.2024.109388. URLhttps://linkinghub.elsevier.com/ret...
-
[34]
Compressible Euler Equa- tions, pp
Kuzmin, D., M ¨oller, M.: Algebraic Flux Correction II. Compressible Euler Equa- tions, pp. 207–250. Springer (2005)
2005
-
[35]
Leer, B.V.A.N.: Towards the Ultimate Conservative Difference Scheme. IV. A New Approach to Numerical Convection276299(1977). DOI https://doi.org/10. 1016/0021-9991(77)90095-X
1977
-
[36]
Cambridge University Press (2002)
LeVeque, R.J.: Finite Volume Methods for Hyperbolic Problems. Cambridge University Press (2002)
2002
-
[37]
Journal of Computational Physics516, 113360 (2024)
Montoya, T., Zingg, D.W.: Efficient entropy-stable discontinuous spectral-element methods using tensor-product summation-by-parts operators on triangles and tetrahedra. Journal of Computational Physics516, 113360 (2024). DOI 10.1016/j.jcp.2024.113360
-
[38]
SIAM Journal on Scientific Computing46(4), A2270–A2297 (2024)
Montoya, T., Zingg, D.W.: Efficient Tensor-Product Spectral-Element Operators with the Summation-by-Parts Property on Curved Triangles and Tetrahedra. SIAM Journal on Scientific Computing46(4), A2270–A2297 (2024). DOI 10.1137/ 23M1573963
2024
-
[39]
Journal of Scientific Computing91(1), 1–36 (2022)
Mossier, P., Beck, A., Munz, C.D.: A p-Adaptive Discontinuous Galerkin Method with hp-Shock Capturing. Journal of Scientific Computing91(1), 1–36 (2022). DOI 10.1007/s10915-022-01770-6
-
[40]
Orszag, S.A.: Spectral methods for problems in complex geometrics, pp. 273–305. Elsevier (1979). DOI 10.1016/b978-0-12-546050-7.50014-9
-
[41]
In: 44th AIAA Aerospace Sciences Meeting and Exhibit, p
Persson, P.O., Peraire, J.: Sub-cell shock capturing for discontinuous Galerkin methods. In: 44th AIAA Aerospace Sciences Meeting and Exhibit, p. 112 (2006)
2006
-
[42]
PAMM23(1), Shock-Capturing for DGSEM on Heterogeneous Curvilinear Meshes 21 1–6 (2023)
Schwarz, A., Keim, J., Chiocchetti, S., Beck, A.: A Reinforcement Learning Based Slope Limiter for Second-Order Finite Volume Schemes. PAMM23(1), Shock-Capturing for DGSEM on Heterogeneous Curvilinear Meshes 21 1–6 (2023). DOI 10.1002/pamm.202200207. URLhttps://onlinelibrary. wiley.com/doi/10.1002/pamm.202200207
-
[43]
Engineer- ing with Computers24(3), 195–213 (2008)
Shepherd, J.F., Johnson, C.R.: Hexahedral mesh generation constraints. Engineer- ing with Computers24(3), 195–213 (2008). DOI 10.1007/s00366-008-0091-4
-
[44]
NASA/CR-97-206253 ICASE Re- port No.97-65 (1997)
Shu, C.: Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic Conservation Laws. NASA/CR-97-206253 ICASE Re- port No.97-65 (1997)
1997
-
[45]
Sonntag, M., Munz, C.D.: Efficient Parallelization of a Shock Capturing for Dis- continuous Galerkin Methods using Finite Volume Sub-cells. J. Sci. Comput. 70(3), 1262–1289 (2017)
2017
-
[46]
Springer (2009)
Toro, E.F.: Riemann Solvers and Numerical Methods for Fluid Dynamics. Springer (2009)
2009
-
[47]
Journal of Computational Physics387, 245–279 (2019)
Vilar, F.: A posteriori correction of high-order discontinuous Galerkin scheme through subcell finite volume formulation and flux reconstruction. Journal of Computational Physics387, 245–279 (2019). DOI 10.1016/j.jcp.2018.10.050
-
[48]
Journal of Computational Physics227(12), 6001–6043 (2008)
Yoon, S.H., Kim, C., Kim, K.H.: Multi-dimensional limiting process for three- dimensional flow physics analyses. Journal of Computational Physics227(12), 6001–6043 (2008). DOI 10.1016/j.jcp.2008.02.012
-
[49]
Journal of Computational Physics229(23), 8918–8934 (2010)
Zhang, X., Shu, C.W.: On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes. Journal of Computational Physics229(23), 8918–8934 (2010). DOI 10.1016/j.jcp.2010.08. 016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.