pith. sign in

arxiv: 2605.04660 · v1 · submitted 2026-05-06 · 🧮 math.NA · cs.NA

A third-order multi-moment cell-centered Lagrangian scheme for hydrodynamics with an accurate 2D nodal solver

Xiaoteng Zhang , Xun Wang , Zhijun Shen , Chao Yang This is my paper

Pith reviewed 2026-05-08 15:57 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Lagrangian hydrodynamicsmulti-moment finite volumenodal Riemann solvercell-centered schemehigh-order accuracycompressible flows2D hydrodynamicsmesh movement
0
0 comments X

The pith

A new Lagrangian scheme for 2D compressible hydrodynamics reaches third-order accuracy by combining multi-moment cell averages with a nodal Riemann solver that enforces pressure-velocity compatibility and nodal conservation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a cell-centered Lagrangian method called LMCV for simulating two-dimensional compressible flows. It merges the multi-moment constrained finite volume method, which evolves both point values at cell vertices and cell-averaged quantities, with a 2D nodal Riemann solver built on EUCCLHYD principles. New jump conditions accurately connect each cell's surface pressure to its nodal velocity, while balance conditions maintain nodal conservation and stabilize velocities as the mesh deforms with the fluid. This pairing preserves high-order accuracy together with conservation and robustness properties. The solver is presented as a natural high-order extension of HLLC and HLLC-2D methods, and numerical experiments confirm its performance on standard problems.

Core claim

By bridging the multi-moment constrained finite volume method with a 2D nodal Riemann solver, the LMCV scheme maintains high-order accuracy while inheriting the conservation and robust properties of the nodal Riemann solver. The new jump condition provides a high-accurate formulation linking the surface pressure of each cell to its nodal velocity, while the balance condition ensures nodal conservation and stabilizes the velocity field without losing accuracy. This nodal solver can be regarded as a natural high-order extension of the HLLC and HLLC-2D solvers.

What carries the argument

The 2D nodal Riemann solver with new jump and balance conditions, which links surface pressure of each cell to nodal velocity and ensures nodal conservation without accuracy loss during mesh movement.

If this is right

  • Rigorous numerical conservation is achieved by simultaneously evolving point-values at cell vertices and volume-integrated averages.
  • Compatibility between mesh movement and numerical fluxes is accomplished through the nodal solver formulations.
  • The solver acts as a direct high-order extension of existing HLLC and HLLC-2D methods in two dimensions.
  • Accuracy and robustness are demonstrated across a variety of numerical experiments on standard hydrodynamics problems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same jump-and-balance structure could be generalized to three dimensions to handle more complex mesh motions.
  • The method may allow stable long-time simulations of flows with large deformations without needing frequent remeshing.
  • Coupling the nodal solver with other high-order reconstruction techniques could further reduce computational cost while retaining the conservation guarantees.

Load-bearing premise

The new jump and balance conditions in the 2D nodal solver deliver high-accurate pressure-to-velocity linking and nodal conservation without accuracy loss when the mesh moves.

What would settle it

A convergence study on a smooth 2D flow with strong mesh deformation that shows the scheme dropping below third-order accuracy or violating global conservation would disprove the central claim.

Figures

Figures reproduced from arXiv: 2605.04660 by Chao Yang, Xiaoteng Zhang, Xun Wang, Zhijun Shen.

Figure 1
Figure 1. Figure 1: The transformation between reference cell view at source ↗
Figure 2
Figure 2. Figure 2: gradient reconstruction on cell surface x i r x i r+1 3.3. Evolution of VIA moment and vertices Similar to the EUCCLHYD [24], (5-7) is discretized as m i view at source ↗
Figure 3
Figure 3. Figure 3: General cells. By means of quadrature rules, the cell boundary integral is usually discretized as d dt Z ωi U¯ dω + 1 mi X q∈Qi lqF¯ qnq = 0, (30) where Q i is the integral node set of cell ω i , nq and lq are the unit normal vector and weight length at node q, F¯ q is the numerical flux computed at q with certain Riemann solver. Moreover, the shape control points of ω i should be advected through d dt xq … view at source ↗
Figure 4
Figure 4. Figure 4: Control points with neighboring cells. In order to give consistent numerical fluxes and velocities at the cell vertices without losing conservation or other properties, various multidimensional Riemann solvers [9, 24, 25] have been developed for Lagrangian schemes, which directly handle the calculation of nodal velocities on straight-edge meshes. However, complex￾ities arise when addressing curvilinear mes… view at source ↗
Figure 5
Figure 5. Figure 5: Surface pressures on ω i ∩ω j defined in different schemes. (a) EUCCLHYD [24] defines two different pressures for each half surface. (b) Schemes [47, 30] using both 2D nodal solver and 1D solver at each surface. The nodal solver is applied at both ends, introducing two pressures for each end. The 1D solver is used to provide both the pressure and the normal velocity at the middle point. The nonlinearity of… view at source ↗
Figure 6
Figure 6. Figure 6: Notations used in nodal solver. Substitute Eqs.(39) into Eqs.(38), we can get v˜ i r by solving the linear system X e∈E(x i r ) Le (αe+ + αe−)  v˜ i r · ne − Vee  ne = 0, (40) where Vee = αe+v i r · ne + αe−v j k · ne + Pe+ − Pe− αe+ + αe− (41) which could be regarded as a weighted average of normal velocity given by classical 1D acoustic Riemann solver for half surface e, more discussion can be found in… view at source ↗
Figure 7
Figure 7. Figure 7: Notations around surface f. velocity fields v˜ are acquired in different way. v˜ for our nodal solver is solved following Section 4.4, while the nodal velocity of HLLC-2D solver is acquired as Section 4.2. Equivalence of Lagrangian flux (13-15) and (43) can be checked that 1 2  v˜ i r + v˜ i r+1  · n i r,r+1 = Z 1 0 u˜(ζ)dζ, 1 2  Pei r,r+ 1 2 + Pei r+ 1 2 ,r+1  = Z 1 0 Pei L (ζ)dζ, 1 2  Pei r,r+ 1 2 v… view at source ↗
Figure 8
Figure 8. Figure 8: The close-up view of pressure contours with mesh us view at source ↗
Figure 9
Figure 9. Figure 9: The global view of pressure contours with mesh usin view at source ↗
Figure 10
Figure 10. Figure 10: The global view (top) and close-up view (below) of view at source ↗
Figure 11
Figure 11. Figure 11: The scatter plots of density at t = 0.2 for the Sod problem with EUCCLHYD (left) and LMCV (right). The black solid line represents the exact density distribution. 24 view at source ↗
Figure 12
Figure 12. Figure 12: The meshes and density contours for the Sedov blas view at source ↗
Figure 13
Figure 13. Figure 13: The scatter plots of density at t = 1 for the Sedov blast problem with EUCCLHYD (left) and LMCV (right). The black solid line represents the exact density distribution from [40]. 7.6. Noh problem Similar to Sedov blast problem, Noh problem [32] is also often used to test the robustness of Lagrangian schemes. It contains a shock caused by the convergence of a uniform gas towards the origin. More specifical… view at source ↗
Figure 14-16
Figure 14-16. Figure 14-16: From comparison in Fig. 16, LMCV gives more accur view at source ↗
Figure 14
Figure 14. Figure 14: The meshes and density contours for Noh problem at view at source ↗
Figure 15
Figure 15. Figure 15: The close-up of the meshes behind the shock at view at source ↗
Figure 16
Figure 16. Figure 16: The scatter plots of density at t = 0.6 for Noh problem with EUCCLHYD (left) and LMCV (right). 28 view at source ↗
Figure 17
Figure 17. Figure 17: The initial mesh for Saltzman problem. 7.8. Triple point problem Finally we consider a three-state 2D Riemann problem, i.e. triple point problem [23]. The computational domain Ω = [0, 7] × [0, 3] is surrounded by rigid walls and split into three regions, and the regions and states is initialized as Ω 1 = [0, 1] × [0, 3] : ρ = 1, v = 0, P = 1, Ω 2 = [1, 7] × [0, 1.5] : ρ = 1, v = 0, P = 0.1, Ω 3 = [1, 7] ×… view at source ↗
Figure 18
Figure 18. Figure 18: The density contours with mesh for Saltzman probl view at source ↗
Figure 19
Figure 19. Figure 19: The scatter plots of density for Saltzman problem view at source ↗
Figure 20
Figure 20. Figure 20: The density contours with mesh given by EUCCLHYD ( view at source ↗
Figure 21
Figure 21. Figure 21: New notations around xk,l Lemma 1 (Geometric properties of Mh). For a set of uniformly refined quadrilateral mesh {Mh}h∈R defined above, we have 1. Lk+ 1 2 ,l , Lk,l+ 1 2 = O(h), 2. |Lk+ 1 2 ,l+1 − Lk+ 1 2 ,l |, |Lk+1,l+ 1 2 − Lk,l+ 1 2 | = O(h 2 ), 3. |nk+ 1 2 ,l+1 − nk+ 1 2 ,l |, |nk+1,l+ 1 2 − nk,l+ 1 2 | = O(h), Proof of Lemma 1. By definition of Mh, it can be easily checked for one dimension that Lk+… view at source ↗
read the original abstract

This paper presents a novel high-order cell-centered Lagrangian scheme for 2D compressible hydrodynamics by bridging the multi-moment constrained finite volume method (MCV) [16, 51, 52] with a nodal Riemann solver. This scheme (denoted by LMCV) not only maintains high-order accuracy as MCV but also inherits the conservation and robust properties of the nodal Riemann solver. On the one hand, the MCV employs and evolves both the point-values (PV) at cell vertexes and the volume-integrated averages (VIA) on computational mesh, which ensures the rigorous numerical conservation and establishes an adequate foundation for the computation of Lagrangian fluxes with high accuracy. On the other hand, we developed a 2D Riemann solver based on EUCCLHYD [24], it takes fully advantage of numerical formulations from high-order scheme and accomplishes the compatibility between the mesh movement and numerical fluxes. The main new features of the solver are the introduction of a new set of jump and balance conditions. The jump condition provides a high-accurate formulation linking the surface pressure of each cell to its nodal velocity, while the balance condition ensures nodal conservation and stabilizes the velocity field without losing accuracy. More intriguing is that our nodal solver can be regarded as a natural high-order extension of the HLLC and the HLLC-2D [41] solvers. The comparison between these solvers better demonstrates our innovative approach in addressing the difficulties encountered in constructing 2D high-order Lagrangian schemes. A variety of numerical experiments are carried out to illustrate the accuracy and robustness of the algorithm.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper presents a third-order cell-centered Lagrangian scheme (LMCV) for 2D compressible hydrodynamics that combines the multi-moment constrained finite volume (MCV) method—evolving both point values at cell vertices and volume-integrated averages—with a new 2D nodal Riemann solver derived from EUCCLHYD. The solver introduces novel jump and balance conditions to provide high-accurate pressure-to-velocity linking at nodes and to enforce nodal conservation while stabilizing the velocity field. The scheme is shown to reduce to the HLLC and HLLC-2D solvers in appropriate limits, satisfies geometric conservation by construction, and is supported by truncation-error analysis together with numerical convergence studies on smooth isentropic flows on moving meshes.

Significance. If the central claims hold, the work supplies a concrete high-order Lagrangian framework that preserves the accuracy of MCV while adding the conservation and robustness properties of a compatible nodal solver. The explicit construction of the jump and balance conditions, their reduction to established solvers, the truncation-error verification, and the demonstration of design-order convergence on deforming meshes constitute verifiable strengths that could support further development of robust high-order Lagrangian hydrodynamics codes.

minor comments (3)
  1. [Abstract] The abstract states that 'a variety of numerical experiments are carried out' but does not list the specific test problems (smooth isentropic vortex, shock-tube, etc.); a one-sentence enumeration in the abstract or introduction would improve readability.
  2. [Nodal solver derivation] In the section deriving the nodal solver, the limiting process that recovers the HLLC/HLLC-2D solvers is asserted but not written out equation-by-equation; inserting the explicit reduction steps would make the comparison more transparent.
  3. [Numerical results] Convergence tables report L2 errors for density but omit corresponding norms for velocity or total energy; adding these would strengthen the claim that the full scheme retains third-order accuracy.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading and positive evaluation of our work, including the accurate summary of the LMCV scheme, its connection to MCV and the EUCCLHYD-based nodal solver, the new jump and balance conditions, the reduction to HLLC/HLLC-2D, geometric conservation, truncation-error analysis, and the convergence results on deforming meshes. We appreciate the recommendation for minor revision.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper explicitly constructs the new jump and balance conditions for the 2D nodal solver from EUCCLHYD with high-order surface data, demonstrates reduction to HLLC/HLLC-2D, and verifies third-order accuracy via truncation-error analysis plus convergence tests on smooth flows. All load-bearing steps (nodal velocity update, flux compatibility, GCL satisfaction) are derived directly in the text or from independent external references without reducing to self-definition, fitted inputs renamed as predictions, or self-citation chains. The central claim of preserved accuracy plus added conservation is independently verifiable from the provided derivations and numerical evidence.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 2 invented entities

Scheme rests on standard hydrodynamics equations and prior MCV accuracy claims; new conditions are introduced without independent verification in the abstract.

axioms (2)
  • domain assumption Compressible Euler equations govern the flow.
    Implicit foundation for any hydrodynamics scheme.
  • domain assumption MCV method preserves third-order accuracy on moving meshes.
    Taken from cited references [16,51,52].
invented entities (2)
  • New jump condition no independent evidence
    purpose: High-accuracy link between cell surface pressure and nodal velocity.
    Introduced in this work as part of the nodal solver.
  • New balance condition no independent evidence
    purpose: Nodal conservation and velocity stabilization.
    Introduced to ensure compatibility with mesh motion.

pith-pipeline@v0.9.0 · 5597 in / 1295 out tokens · 63307 ms · 2026-05-08T15:57:50.311371+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Batten, N

    P . Batten, N. Clarke, C. Lambert and D. Causon, On the choi ce of wavespeeds for the HLLC Riemann solver, J. Sci. Comput. 18(6) (1997) 1553–1570

    work page 1997

  2. [2]

    Burton, N

    D. Burton, N. Morgan, T. Carney and M. Kenamond, Reductio n of dissipation in Lagrange cell-centered hydrodynamics CCH through corner gradient reconstruction CGR, J. Comput. Phys. 299 (2015) 229–280

    work page 2015

  3. [3]

    Charest, C

    M. Charest, C. Groth and P . Gauthier, A high-order centra l ENO finite-volume scheme for three- dimensional low-speed viscous flows on unstructured mesh, C ommun. Comput. Phys. 17(3) (2015) 615- 656. 40

    work page 2015

  4. [4]

    Cheng and C.W

    J. Cheng and C.W. Shu, A third order conservative Lagrang ian type scheme on curvilinear meshes for the compressible Euler equations, Commun. Comput. Phys. 4 (200 8) 1008–1024

    work page

  5. [5]

    Cockburn and C.W

    B. Cockburn and C.W. Shu, TVB Runge-Kutta local projecti on discontinuous Galerkin finite element method for conservation laws. II. General framework, Math. Comput. 52(186) (1989) 411–435

    work page 1989

  6. [6]

    Cockburn, S.Y

    B. Cockburn, S.Y . Lin and C.W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: 1D systems, J. Com put. Phys. 84(1) (1989) 90–113

    work page 1989

  7. [7]

    Cockburn, S.C

    B. Cockburn, S.C. Hou and C.W. Shu, The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV . The multidimensi onal case, Math. Comput. 54(190) (1990) 545–581

    work page 1990

  8. [8]

    Cockburn and C.W

    B. Cockburn and C.W. Shu, The Runge–Kutta discontinuous Galerkin method for conservation laws V: Multidimensional systems, J. Comput. Phys. 141(2) (1998) 1 99–224

    work page 1998

  9. [9]

    Després and C

    B. Després and C. Mazeran, Lagrangian Gas Dynamics in Two Dimensions and Lagrangian systems, Arch. Rational Mech. Anal. 178 (2005) 327–372

    work page 2005

  10. [10]

    Dobrev, T

    V . Dobrev, T. Kolev and R. Rieben, High-order curviline ar finite element methods for Lagrangian hydro- dynamics, SIAM J. Sci. Comput. 34 (2012) B606–B641

    work page 2012

  11. [11]

    Dukowicz and B

    J. Dukowicz and B. Meltz, V orticity errors in multidime nsional Lagrangian codes, J. Comput. Phys. 99 (1992) 115–134

    work page 1992

  12. [12]

    Godunov, A di fference method for numerical calculation of discontinuous e quations of hydrodynam- ics, Math

    S.K. Godunov, A di fference method for numerical calculation of discontinuous e quations of hydrodynam- ics, Math. Sb. 47 (1959) 217

    work page 1959

  13. [13]

    Harten, B

    A. Harten, B. Engquist, S. Osher and S.R. Chakravarthy, Uniformly high order accurate essentially non- oscillatory schemes, III, J. Comput. Phys. 131(1) (1997) 3– 47

    work page 1997

  14. [14]

    Harten, P .D

    A. Harten, P .D. Lax, B. van Leer, On upstream di fferencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev. 25 (1983) 35–61

    work page 1983

  15. [15]

    Hirt, B.D

    C.W. Hirt, B.D. Nichols, V olume of fluid (VOF) method for the dynamics of free boundaries, J. Comput. Phys. 39 (1) (1981) 201–225

    work page 1981

  16. [16]

    Ii and F

    S. Ii and F. Xiao, High order multi-moment constrained fi nite volume method. Part I: Basic formulation, J. Comput. Phys. 228(10) (2009) 3669–3707

    work page 2009

  17. [17]

    Ivan and C

    L. Ivan and C. Groth, High-order solution-adaptive cen tral essentially non-oscillatory (CENO) method for viscous flows, J. Comput. Phys. 257 (2014) 830-862

    work page 2014

  18. [18]

    Jiang and C.W

    G.S. Jiang and C.W. Shu, E fficient implementation of weighted ENO schemes, J. Comput. Ph ys. 126(1) (1996) 202–228. 41

    work page 1996

  19. [19]

    Law and Philip T

    Timothy R. Law and Philip T. Barton, A cell-centred Eule rian volume-of-fluid method for compressible multi-material flows, J. Comput. Phys. 497 (2024) 112592

    work page 2024

  20. [20]

    Y . Liu, M. Vinokur and Z.J. Wang, Spectral di fference method for unstructured grids I: Basic formulation, J. Comput. Phys. 216(2) (2006) 780–801

    work page 2006

  21. [21]

    X.D. Liu, S. Osher and T. Chan, Weighted essentially non -oscillatory schemes, J. Comput. Phys. 115(1) (1994) 200–212

    work page 1994

  22. [22]

    X.D. Liu, N. Morgan and D. Burton, A high-order Lagrangi an discontinuous Galerkin hydrodynamic method for quadratic Cells using a subcell mesh stabilizati on scheme, J. Comput. Phys. 386 (2019) 101–157

    work page 2019

  23. [23]

    Loubère, V alidation Test Case Suite for Compressibl e Hydrodynamics Computation, Technical report, Los Alamos National Laboratory, 2005

    R. Loubère, V alidation Test Case Suite for Compressibl e Hydrodynamics Computation, Technical report, Los Alamos National Laboratory, 2005

    work page 2005

  24. [24]

    Maire, R

    P-H. Maire, R. Abgrall, J. Breil and J. Ovadia, A cell-ce ntered Lagrangian scheme for 2D compressible flow problems, SIAM J. Sci. Comput. 29 (2007) 1781–1824

    work page 2007

  25. [25]

    Maire, A high-order cell-centered Lagrangian sch eme for 2D compressible fluid flows on unstruc- tured mesh, J

    P-H. Maire, A high-order cell-centered Lagrangian sch eme for 2D compressible fluid flows on unstruc- tured mesh, J. Comput. Phys. 228 (2009) 2391–2425

    work page 2009

  26. [26]

    Maire, R

    P-H. Maire, R. Loubère, P . V achal, Staggered Lagrangia n discretization based on cell-centered Riemann solver and associated hydrodynamics scheme, Commun. Compu t. Phys. 10 (2011) 940–978

    work page 2011

  27. [27]

    Margolin, A Centered Artificial Viscosity for Cells w ith Large Aspect Ratios, Technical Report UCRL- S3882, Los Alamos National Laboratory, 1988

    L. Margolin, A Centered Artificial Viscosity for Cells w ith Large Aspect Ratios, Technical Report UCRL- S3882, Los Alamos National Laboratory, 1988

    work page 1988

  28. [28]

    Michalak and C

    C. Michalak and C. Ollivier-Gooch, Accuracy preservin g limiter for the high-order accurate solution of the Euler equations, J. Comput. Phys. 228(23) (2009) 8693-8 711

    work page 2009

  29. [29]

    Morgan, J

    N. Morgan, J. Waltz, D. Burton, M. Charest, T. Canfield an d J. Wohlbier, A Godunov-like point-centered essentially Lagrangian hydrodynamic approach, J. Comput. Phys. 281 (2014) 614–652

    work page 2014

  30. [30]

    Morgan, X

    N. Morgan, X. Liu and D. Burton, A Lagrangian discontinu ous Galerkin hydrodynamic method for higher- order triangular elements, in: 2018 AIAA Aerospace Science s Meeting, Kissimmee, Florida, 2018, AIAA- 2018-1092

    work page 2018

  31. [31]

    Morgan, X

    N. Morgan, X. Liu and D. Burton, Reducing spurious mesh m otion in Lagrangian finite volume and dis- continuous Galerkin hydrodynamic methods, J. Comput. Phys . 372 (2018) 35–61

    work page 2018

  32. [32]

    Noh, Errors for calculations of strong shocks using a n artificial viscosity and an artificial heat flux, J

    W. Noh, Errors for calculations of strong shocks using a n artificial viscosity and an artificial heat flux, J. Appl. Phys. 72 (1987) 78–120

    work page 1987

  33. [33]

    Osher, J.A

    S. Osher, J.A. Sethian, Fronts propagating with curvat ure-dependent speed: algorithms based on Hamilton- Jacobi formulations, J. Comput. Phys. 79 (1) (1988) 12–49. 42

    work page 1988

  34. [34]

    Park, S.H

    J.S. Park, S.H. Y oon, and C. Kim, Multi-dimensional lim iting process for hyperbolic conservation laws on unstructured grids, J. Comput. Phys. 229(3) (2010) 788-812

    work page 2010

  35. [35]

    Pederson, B

    C. Pederson, B. Brown and N. Morgan, The Sedov blast wave as a radial piston verification problem, J. V erification V alid. Uncertain. Quantification 1 (2016) 1–9

    work page 2016

  36. [36]

    Pilliod, E.G

    J.E. Pilliod, E.G. Puckett, Second-order accurate vol ume-of-fluid algorithms for tracking material inter- faces, J. Comput. Phys. 199 (2) (2004) 465–502

    work page 2004

  37. [37]

    Qiu and C.W

    J.X. Qiu and C.W. Shu, Hermite WENO Schemes and Their App lication as Limiters for Runge–Kutta Discontinuous Galerkin Method: 1D Case, J. Comput. Phys. 19 3(1) (2004) 115–135

    work page 2004

  38. [38]

    Qiu and C.W

    J.X. Qiu and C.W. Shu, Hermite WENO schemes and their app lication as limiters for Runge–Kutta dis- continuous Galerkin method II: Two dimensional case, Compu t. Fluids 34(6) (2005) 642–663

    work page 2005

  39. [39]

    Roe, Approximate Riemann solvers, parameter vect or and di fference schemes, J

    P .L. Roe, Approximate Riemann solvers, parameter vect or and di fference schemes, J. Comput. Phys. 43 (1981) 357–372

    work page 1981

  40. [40]

    Sedov, Similarity and Dimensional Methods in Mechan ics, Academic Press, 1959

    L. Sedov, Similarity and Dimensional Methods in Mechan ics, Academic Press, 1959

    work page 1959

  41. [41]

    Z.J. Shen, W. Y an and G.W. Y uan, A robust and contact reso lving Riemann solver on unstructured mesh, Part II, ALE method, J. Comput. Phys. 268 (2014) 456–484

    work page 2014

  42. [42]

    Shu, Essentially non-oscillatory and weighted es sentially non-oscillatory schemes for hyperbolic conservation laws, ICASE Report No

    C.W. Shu, Essentially non-oscillatory and weighted es sentially non-oscillatory schemes for hyperbolic conservation laws, ICASE Report No. 97-65, 1997

    work page 1997

  43. [43]

    Shu, Total-variation-diminishing time discreti zations, SIAM J

    C.W. Shu, Total-variation-diminishing time discreti zations, SIAM J. Sci. Stat. Comput. 9 (1988) 1073–1084

    work page 1988

  44. [44]

    G. A. Sod, A survey of several finite di fference methods for systems of non-linear hyperbolic conser vation laws, J. Comput. Phys. 27 (1978) 1–31

    work page 1978

  45. [45]

    van Leer, Towards the ultimate conservative di fference scheme

    B. van Leer, Towards the ultimate conservative di fference scheme. V . A second-order sequel to Godunov /AP s method, J. Comput. Phys. 32(1) (1979) 101–136

    work page 1979

  46. [46]

    Vilar, Cell-centered discontinuous Galerkin discr etization for two-dimensional Lagrangian hydrody- namics, Comput

    F. Vilar, Cell-centered discontinuous Galerkin discr etization for two-dimensional Lagrangian hydrody- namics, Comput. Fluids 64 (2012) 64–73

    work page 2012

  47. [47]

    Vilar, P-H

    F. Vilar, P-H. Maire and R. Abgrall, A discontinuous Gal erkin discretization for solving the 2D gas dy- namics equations written under total Lagrangian formulati on on general unstructured grids, J. Comput. Phys. 276 (2014) 188–234

    work page 2014

  48. [48]

    Vilar, C.W

    F. Vilar, C.W. Shu and P-H. Maire, Positivity-preservi ng cell-centered Lagrangian schemes for multi- material compressible flows: From first-order to high-order s. Part II: The two-dimensional case, J. Comput. Phys. 312 (2016) 416-442. 43

    work page 2016

  49. [49]

    Z.J. Wang, Y . Liu, G. May and A. Jameson, Spectral difference method for unstructured grids II: Extension to the Euler equations, J. Sci. Comput. 32(1) (2007) 45–71

    work page 2007

  50. [50]

    Wendro ff, A two-dimensional HLLE Riemann solver and associated Godu nov-type difference scheme for gas dynamics, Comput

    B. Wendro ff, A two-dimensional HLLE Riemann solver and associated Godu nov-type difference scheme for gas dynamics, Comput. Math. 38 (1999) 175–185

    work page 1999

  51. [51]

    B. Xie, S. Ii, A. Ikebata and F. Xiao, A multi-moment finit e volume method for incompressible Navier–Stokes equations on unstructured grids: V olume-av erage/point-value formulation, J. Comput. Phys. 227 (2014) 138–162

    work page 2014

  52. [52]

    Xie and F

    B. Xie and F. Xiao, A multi-moment constrained finite vol ume method on arbitrary unstructured grids for incompressible flows, J. Comput. Phys. 327 (2016) 747-778

    work page 2016

  53. [53]

    Zhu and J.X

    J. Zhu and J.X. Qiu, Hermite WENO schemes and their appli cation as limiters for Runge-Kutta discontin- uous Galerkin method, III: Unstructured meshes, J. Sci. Com put. 39(2) (2009) 293–321. 44

    work page 2009