Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

Abstract We present new Lagrangian discontinuous Galerkin (DG) hydrodynamic methods for compressible flows on unstructured meshes in axisymmetric coordinates. The physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with linear Taylor series polynomials. Two different approaches are used to discretize the evolution equations – the first one is the true volume approach and the second one is the area-weighted approach. For the true volume approach, the DG equations are derived using the true 3D volume that is consistent with the geometry conservation law (GCL). The Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem using surfaces areas for axisymmetric coordinates. This true volume approach conserves mass, momentum, and total energy and satisfies the GCL. However, it can not preserve spherical symmetry on an equal-angle polar grid with 1D radial flows. For the area-weighted approach, the DG equations are based on the 2D Cartesian geometry that is rotated about the axis of symmetry using a single, element average radius. With this approach, the Riemann velocity at the nodes on the surface of the element, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem in 2D Cartesian geometry. This area-weighted approach, in the limit of an infinitesimal mesh size, conserves physical momentum, and physical total energy. The area-weighted approach preserves spherical symmetry on an equal-angle polar grid for 1D radial flows, but it does not satisfy the GCL. A suite of test problems are calculated to demonstrate stable mesh motion, the expected second order accuracy of these methods, and that the new area-weighted DG method preserves spherical symmetry on 1D radial flow problems with equal-angle polar meshes.

[1]  Raphaël Loubère,et al.  Staggered Lagrangian Discretization Based on Cell-Centered Riemann Solver and Associated Hydrodynamics Scheme , 2011 .

[2]  Nathaniel R. Morgan,et al.  An approach for treating contact surfaces in Lagrangian cell-centered hydrodynamics , 2013, J. Comput. Phys..

[3]  Rémi Abgrall,et al.  A Cell-Centered Lagrangian Scheme for Two-Dimensional Compressible Flow Problems , 2007, SIAM J. Sci. Comput..

[4]  Hong Luo,et al.  A reconstructed discontinuous Galerkin method for the compressible Navier-Stokes equations on three-dimensional hybrid grids , 2017 .

[5]  Xiaodong Liu,et al.  Reducing spurious mesh motion in Lagrangian finite volume and discontinuous Galerkin hydrodynamic methods , 2018, J. Comput. Phys..

[6]  Xijun Yu,et al.  The cell-centered discontinuous Galerkin method for Lagrangian compressible Euler equations in two-dimensions , 2014 .

[7]  Donald E. Burton,et al.  Compatible, energy and symmetry preserving 2D lagrangian hydrodynamics in rz - cylindrical coordinates , 2010, ICCS.

[8]  R. Kidder,et al.  Laser-driven compression of hollow shells: power requirements and stability limitations , 1976 .

[9]  Juan Cheng,et al.  Positivity-preserving and symmetry-preserving Lagrangian schemes for compressible Euler equations in cylindrical coordinates , 2017 .

[10]  P. P. Whalen,et al.  Algebraic Limitations on Two-Dimensional Hydrodynamics Simulations , 1996 .

[11]  L. Sedov.  Similarity and Dimensional Methods in Mechanics , 1960 .

[12]  Nathaniel R. Morgan,et al.  A Lagrangian discontinuous Galerkin hydrodynamic method for 2D Cartesian and RZ axisymmetric coordinates , 2018 .

[13]  Shudao Zhang,et al.  A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions , 2011, J. Comput. Phys..

[14]  Chi-Wang Shu,et al.  Second order symmetry-preserving conservative Lagrangian scheme for compressible Euler equations in two-dimensional cylindrical coordinates , 2014, J. Comput. Phys..

[15]  Pavel Váchal,et al.  On preservation of symmetry in r-z staggered Lagrangian schemes , 2016, J. Comput. Phys..

[16]  Nathaniel R. Morgan,et al.  A cell-centered Lagrangian Godunov-like method for solid dynamics , 2013 .

[17]  Donald E. Burton,et al.  A cell-centered Lagrangian finite volume approach for computing elasto-plastic response of solids in cylindrical axisymmetric geometries , 2013, J. Comput. Phys..

[18]  Nathaniel R. Morgan,et al.  On the question of area weighting in cell-centered hydrodynamics , 2013 .

[19]  E. J. Caramana,et al.  Numerical Preservation of Symmetry Properties of Continuum Problems , 1998 .

[20]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[21]  Bernardo Cockburn,et al.  The Runge-Kutta local projection discontinous Galerkin finite element method for conservation laws , 1990 .

[22]  Pavel Váchal,et al.  A symmetry preserving dissipative artificial viscosity in an r-z staggered Lagrangian discretization , 2014, J. Comput. Phys..

[23]  Hong Luo,et al.  Development and Assessment of a Reconstructed Discontinuous Galerkin Method for the Compressible Turbulent Flows on Hybrid Grids , 2016 .

[24]  Bruno Després,et al.  Lagrangian Gas Dynamics in Two Dimensions and Lagrangian systems , 2005 .

[25]  Pavel Váchal,et al.  A symmetry preserving dissipative artificial viscosity in r‐z geometry , 2014 .

[26]  Chi-Wang Shu,et al.  A cell-centered Lagrangian scheme with the preservation of symmetry and conservation properties for compressible fluid flows in two-dimensional cylindrical geometry , 2010, J. Comput. Phys..

[27]  Pierre-Henri Maire,et al.  A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry , 2009, J. Comput. Phys..

[28]  Nathaniel R. Morgan,et al.  A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for elastic-plastic flows , 2019, Comput. Math. Appl..

[29]  Mark L. Wilkins,et al.  Use of artificial viscosity in multidimensional fluid dynamic calculations , 1980 .

[30]  Roger B. Lazarus,et al.  Self-Similar Solutions for Converging Shocks and Collapsing Cavities , 1981 .

[31]  François Vilar,et al.  Cell-centered discontinuous Galerkin discretization for two-dimensional Lagrangian hydrodynamics , 2012 .

[32]  W. F. Noh.  Errors for calculations of strong shocks using an artificial viscosity and artificial heat flux , 1985 .

[33]  M. Shashkov,et al.  The Construction of Compatible Hydrodynamics Algorithms Utilizing Conservation of Total Energy , 1998 .

[34]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[35]  Christian Rohde,et al.  An Introduction to Recent Developments in Theory and Numerics for Conservation Laws: Proceedings of the International School on Theory and Numerics for Conservation Laws, Freiburg/Littenweiler, Germany, October 20-24, 1997 , 1999, Theory and Numerics for Conservation Laws.

[36]  G. Yuan,et al.  A cell-centered lagrangian scheme in two-dimensional cylindrical geometry , 2008 .

[37]  Nathaniel R. Morgan,et al.  A Lagrangian discontinuous Galerkin hydrodynamic method , 2018 .

[38]  Nathaniel R. Morgan,et al.  Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme , 2018, J. Comput. Phys..

[39]  David C. Look,et al.  Defect Models in Electron-Irradiated N-Type GaAs , 1992 .

[40]  Len G. Margolin,et al.  Using a Curvilinear Grid to Construct Symmetry-Preserving Discretizations for Lagrangian Gas Dynamics , 1999 .

[41]  Rémi Abgrall,et al.  A discontinuous Galerkin discretization for solving the two-dimensional gas dynamics equations written under total Lagrangian formulation on general unstructured grids , 2014, J. Comput. Phys..

[42]  Nathaniel R. Morgan,et al.  A Lagrangian staggered grid Godunov-like approach for hydrodynamics , 2014, J. Comput. Phys..

[43]  M. Shashkov,et al.  A discrete operator calculus for finite difference approximations , 2000 .

[44]  C. L. Rousculp,et al.  A Compatible, Energy and Symmetry Preserving Lagrangian Hydrodynamics Algorithm in Three-Dimensional Cartesian Geometry , 2000 .

[45]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[46]  Pierre-Henri Maire,et al.  A high-order one-step sub-cell force-based discretization for cell-centered Lagrangian hydrodynamics on polygonal grids , 2011 .

[47]  Chuanjin Wang.  Reconstructed Discontinous Galerkin Method for the Compressible Navier-Stokes Equations in Arbitrary Langrangian and Eulerian Formulation , 2017 .

[48]  Pierre-Henri Maire,et al.  A high-order cell-centered Lagrangian scheme for two-dimensional compressible fluid flows on unstructured meshes , 2009, J. Comput. Phys..

[49]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[50]  Hong Luo,et al.  A Comparative Study of Rosenbrock-Type and Implicit Runge-Kutta Time Integration for Discontinuous Galerkin Method for Unsteady 3D Compressible Navier-Stokes equations , 2016 .

[51]  Donald E. Burton,et al.  Multidimensional discretization of conservation laws for unstructured polyhedral grids , 1994 .