A discontinuous Galerkin discretization for solving the two-dimensional gas dynamics equations written under total Lagrangian formulation on general unstructured grids

JOURNAL OF COMPUTATIONAL PHYSICS

卷:276

Article

Vilar, Francois¹  Maire, Pierre-Henri²  Abgrall, Remi³ 

[1] Brown Univ, Providence, RI 02912 USA

[2] CEA, CESTA, F-33116 Le Barp, France

[3] Univ Zurich, Inst Math, CH-8057 Zurich, Switzerland

关键词: Discontinuous Galerkin discretization;    Total Lagrangian formulation;    Updated Lagrangian formulation;    Cell-centered scheme;    Godunov-type method;    Unstructured moving grid;    Curvilinear grid;    Gas dynamics;   

DOI  :  10.1016/j.jcp.2014.07.030

来源: Elsevier

【 摘 要 】

Based on the total Lagrangian kinematical description, adiscontinuous Galerkin (DG) discretization of the gas dynamics equations is developed for two-dimensional fluid flows on general unstructured grids. Contrary to the updated Lagrangian formulation, which refers to the current moving configuration of the flow, the total Lagrangian formulation refers to the fixed reference configuration, which is usually the initial one. In this framework, the Lagrangian and Eulerian descriptions of the kinematical and the physical variables are related by means of the Piola transformation. Here, we describe a cell-centered high-order DG discretization of the physical conservation laws. The geometrical conservation law, which governs the time evolution of the deformation gradient, is solved by means of a finite element discretization. This approach allows to satisfy exactly the Piola compatibility condition. Regarding the DG approach, it relies on the use of a polynomial space approximation which is spanned by a Taylor basis. The main advantage in using this type of basis relies on its adaptability regardless the shape of the cell. The numerical fluxes at the cell interfaces are computed employing a node-based solver which can be viewed as an approximate Riemann solver. We present numerical results to illustrate the robustness and the accuracy up to third-order of our DG method. First, we show its ability to accurately capture geometrical features of a flow region employing curvilinear grids. Second, we demonstrate the dramatic improvement in symmetry preservation for radial flows. (C) 2014 Elsevier Inc. All rights reserved.

【 授权许可】

Free