A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme

Liu, Xiaodong; Morgan, Nathaniel R.; Burton, Donald E.

doi:10.1016/j.jcp.2019.02.008

Title: A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme

Full Record
Other Related Research

Abstract

We present here a Lagrangian discontinuous Galerkin (DG) hydrodynamic method that is up to third-order accurate using subcell mesh stabilization (SMS) for compressible flows on quadratic meshes in two-dimensional (2D) Cartesian coordinates. Similar to the second-order accurate Lagrangian DG method with linear meshes, the physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with Taylor series polynomials. The Riemann velocity at the vertices of a curvilinear cell, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem. Curvilinear cells (e.g., quadratic quadrilateral meshes in this work) have many deformational degrees of freedom, and with these cells, they can deform in unphysical ways. Likewise, the Riemann solution at an edge vertex differs from the one at the corner of a cell. With SMS, each quadratic quadrilateral cell is decomposed into four quadrilateral subcells, that move in a Lagrangian manner. The edge vertex is surrounded by four subcells so that it is similar to the vertex at the cell corner. SMS can detect inconsistent density fields between the cell and subcells. The difference between these two density fields is used to correct the stress (pressure) input to themore »« less

Authors:

Liu, Xiaodong ^[1];

^[1];

^[1]

Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Publication Date:: Wed Feb 27 00:00:00 EST 2019

Research Org.:: Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

Sponsoring Org.:: USDOE National Nuclear Security Administration (NNSA)

OSTI Identifier:: 1542835

Alternate Identifier(s):: OSTI ID: 1635997

Report Number(s):: LA-UR-18-26689
Journal ID: ISSN 0021-9991

Grant/Contract Number:: 89233218CNA000001

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 386; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 96 KNOWLEDGE MANAGEMENT AND PRESERVATION; Lagrangian shock hydrodynamics; discontinuous Galerkin; cell-centered hydrodynamics; third-order accurate; quadratic curvilinear cells

Citation Formats


                    Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme.  United States: N. p., 2019. 
Web.  doi:10.1016/j.jcp.2019.02.008.

Copy to clipboard


                    Liu, Xiaodong, Morgan, Nathaniel R., & Burton, Donald E. A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme.  United States.  https://doi.org/10.1016/j.jcp.2019.02.008

Copy to clipboard


                    Liu, Xiaodong, Morgan, Nathaniel R., and Burton, Donald E. Wed .  
"A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme".  United States.  https://doi.org/10.1016/j.jcp.2019.02.008.  https://www.osti.gov/servlets/purl/1542835.

Copy to clipboard


                    
@article{osti_1542835,

  title        = {A high-order Lagrangian discontinuous Galerkin hydrodynamic method for quadratic cells using a subcell mesh stabilization scheme},

  author       = {Liu, Xiaodong and Morgan, Nathaniel R. and Burton, Donald E.},

  abstractNote = {We present here a Lagrangian discontinuous Galerkin (DG) hydrodynamic method that is up to third-order accurate using subcell mesh stabilization (SMS) for compressible flows on quadratic meshes in two-dimensional (2D) Cartesian coordinates. Similar to the second-order accurate Lagrangian DG method with linear meshes, the physical evolution equations for the specific volume, velocity, and specific total energy are discretized using a modal DG method with Taylor series polynomials. The Riemann velocity at the vertices of a curvilinear cell, and the corresponding surface forces, are calculated by solving a multidirectional approximate Riemann problem. Curvilinear cells (e.g., quadratic quadrilateral meshes in this work) have many deformational degrees of freedom, and with these cells, they can deform in unphysical ways. Likewise, the Riemann solution at an edge vertex differs from the one at the corner of a cell. With SMS, each quadratic quadrilateral cell is decomposed into four quadrilateral subcells, that move in a Lagrangian manner. The edge vertex is surrounded by four subcells so that it is similar to the vertex at the cell corner. SMS can detect inconsistent density fields between the cell and subcells. The difference between these two density fields is used to correct the stress (pressure) input to the Riemann solver. This SMS scheme enables stable mesh motion and accurate solutions in the context of a Lagrangian high-order DG method that is up to third-order with quadratic cells. We also present effective limiting strategies that ensure monotonicity of the primitive variables with the high-order DG method. This Lagrangian DG hydrodynamic method with SMS conserves mass, momentum, and total energy. A suite of test problems are calculated to demonstrate the designed order of accuracy (up to third-order accurate) of this method, and that the Lagrangian DG method using SMS preserves cylindrical symmetry on 1D radial flow problems with an equal-angle polar quadratic mesh.},

  doi          = {10.1016/j.jcp.2019.02.008},

  journal      = {Journal of Computational Physics},

  number       = ,

  volume       = 386,

  place        = {United States},

  year         = {Wed Feb 27 00:00:00 EST 2019},

  month        = {Wed Feb 27 00:00:00 EST 2019}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2019.02.008

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 24 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

A fourth-order Lagrangian discontinuous Galerkin method using a hierarchical orthogonal basis on curvilinear grids

Journal Article Liu, Xiaodong ; Morgan, Nathaniel R. ; Lieberman, Evan J. ; ... - Journal of Computational and Applied Mathematics

The existing high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic methods are restricted to using quadratic meshes with quadratic polynomials (P2), which in turn, yield up to third-order accuracy. Additionally, these existing DG hydrodynamic schemes, when extended to work with cubic meshes and cubic polynomials (P3), can be unstable on strong-shock problems. Therefore, this paper presents a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method to simulate compressible material dynamics ( e.g., gasses, fluids, and solids) with strong-shocks using cubic meshes and cubic polynomials, and delivers up to fourth-order accuracy on smooth flows. The stability on shock problems is achieved using new hierarchicalmore »« less
https://doi.org/10.1016/j.cam.2021.113890

Full Text Available
Reducing spurious mesh motion in Lagrangian finite volume and discontinuous Galerkin hydrodynamic methods

Journal Article Morgan, Nathaniel Ray ; Liu, Xiaodong ; Burton, Donald E. - Journal of Computational Physics

The Lagrangian finite volume (FV) cell-centered hydrodynamic (CCH) method and the Lagrangian discontinuous Galerkin (DG) CCH method have been demonstrated to be quite stable and capable of producing very accurate solutions on many mesh topologies. However, some challenges can arise with higher-order elements and polygonal elements that have many deformational degrees of freedom. With these types of meshes, elements can deform in unphysical ways and the mesh can tangle. In this study, we present methods for obtaining more robust Lagrangian solutions on polygonal and higher-order elements. The robustness is achieved by (1) incorporating a new iterative method that modifies themore »« less
Cited by 12
https://doi.org/10.1016/j.jcp.2018.06.008

Full Text Available
Lagrangian discontinuous Galerkin hydrodynamic methods in axisymmetric coordinates

Journal Article Liu, Xiaodong ; Morgan, Nathaniel R. ; Burton, Donald E. - Journal of Computational Physics

In this paper, 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 velocitymore »« less
Cited by 14
https://doi.org/10.1016/j.jcp.2018.06.073

Full Text Available
A Lagrangian discontinuous Galerkin hydrodynamic method

Journal Article Liu, Xiaodong ; Morgan, Nathaniel Ray ; Burton, Donald E. - Computers and Fluids

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for solving the two-dimensional gas dynamic equations on unstructured hybrid meshes. The physical conservation laws for the momentum and total energy are discretized using a DG method based on linear Taylor expansions. Three different approaches are investigated for calculating the density variation over the element. The first approach evolves a Taylor expansion of the specific volume field. The second approach follows certain finite element methods and uses the strong mass conservation to calculate the density field at a location inside the element or on the element surface. The thirdmore »« less
Cited by 32
https://doi.org/10.1016/j.compfluid.2017.12.007

Full Text Available
A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

Journal Article Wang, Chuanjin ; Luo, Hong ; Shashkov, Mikhail Jurievich - Computers and Fluids

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. Amore »« less
Cited by 7
https://doi.org/10.1016/j.compfluid.2020.104522

Full Text Available

Similar Records