A Moving Discontinuous Galerkin Finite Element Method with Interface Conservation Enforcement for Compressible Flows

A moving discontinuous Galerkin finite element method with interface conservation enforcement (MDG+ICE) is developed for solving the compressible Euler equations. The MDG+ICE method is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG+ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for compressible flows are conducted to assess the accuracy and robustness of the MDG+ICE method. Numerical results obtained indicate that the MDG+ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

[1]  Wanai Li,et al.  The multi-dimensional limiters for solving hyperbolic conservation laws on unstructured grids , 2011, J. Comput. Phys..

[2]  Harold L. Atkins,et al.  QUADRATURE-FREE IMPLEMENTATION OF DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS , 1996 .

[3]  Jean-Yves Trépanier,et al.  A Conservative Shock Fitting Method on Unstructured Grids , 1996 .

[4]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[5]  W. M. Keck,et al.  Adaptive grid generation by minimizing residuals , 2002 .

[6]  Michael Dumbser,et al.  A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes , 2008, J. Comput. Phys..

[7]  Renato Paciorri,et al.  Unsteady shock‐fitting for unstructured grids , 2016 .

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

[9]  Rainald Löhner,et al.  A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids , 2007, J. Comput. Phys..

[10]  Marco Luciano Savini,et al.  Discontinuous Galerkin solution of the Reynolds-averaged Navier–Stokes and k–ω turbulence model equations , 2005 .

[11]  Andrew T. Corrigan,et al.  The Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Unsteady Three-Dimensional Flows , 2019, AIAA Scitech 2019 Forum.

[12]  Jun Liu,et al.  A shock-fitting technique for cell-centered finite volume methods on unstructured dynamic meshes , 2017, J. Comput. Phys..

[13]  Andrew Corrigan,et al.  A moving discontinuous Galerkin finite element method for flows with interfaces , 2018, International Journal for Numerical Methods in Fluids.

[14]  David A. Kessler,et al.  Convergence of the Moving Discontinuous Galerkin Method with Interface Condition Enforcement in the Presence of an Attached Curved Shock , 2019, AIAA Aviation 2019 Forum.

[15]  Leif Kobbelt,et al.  A Robust Procedure to Eliminate Degenerate Faces from Triangle Meshes , 2001, VMV.

[16]  D. Kessler,et al.  The Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Viscous Flows , 2020 .

[17]  David L. Darmofal,et al.  p-Multigrid solution of high-order discontinuous Galerkin discretizations of the compressible Navier-Stokes equations , 2005 .

[18]  M. D. Salas,et al.  A Brief History of Shock-Fitting , 2011 .

[19]  Hong Luo,et al.  A Hermite WENO reconstruction-based discontinuous Galerkin method for the Euler equations on tetrahedral grids , 2012, J. Comput. Phys..

[20]  Rainald Löhner,et al.  THE NUMERICAL SIMULATION OF STRONGLY UNSTEADY FLOW WITH HUNDREDS OF MOVING BODIES , 1999 .

[21]  P. Roe,et al.  Space-Time Methods for Hyperbolic Conservation Laws , 1998 .

[22]  Dmitri Kuzmin,et al.  Slope limiting for discontinuous Galerkin approximations with a possibly non‐orthogonal Taylor basis , 2013 .

[23]  M. Liou.  The Root Cause of the Overheating Problem , 2017 .

[24]  Bernardo Cockburn,et al.  A class of embedded discontinuous Galerkin methods for computational fluid dynamics , 2015, J. Comput. Phys..

[25]  S. Rebay,et al.  High-Order Accurate Discontinuous Finite Element Solution of the 2D Euler Equations , 1997 .

[26]  Li Wang,et al.  Implicit solution of the unsteady Euler equations for high-order accurate discontinuous Galerkin discretizations , 2006, J. Comput. Phys..

[27]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[28]  O. C. Zienkiewicz,et al.  Automatic Mesh Generation , 2005 .

[29]  J. J. Moré,et al.  Levenberg--Marquardt algorithm: implementation and theory , 1977 .

[30]  Raytcho D. Lazarov,et al.  Unified Hybridization of Discontinuous Galerkin, Mixed, and Continuous Galerkin Methods for Second Order Elliptic Problems , 2009, SIAM J. Numer. Anal..

[31]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[32]  Christoph Lehrenfeld,et al.  The Nitsche XFEM-DG Space-Time Method and its Implementation in Three Space Dimensions , 2014, SIAM J. Sci. Comput..