Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

JOURNAL OF COMPUTATIONAL PHYSICS

卷:355

Article

Burton, D. E.¹  Morgan, N. R.¹  Charest, M. R. J.¹  Kenamond, M. A.¹  Fung, J.¹ 

[1] Los Alamos Natl Lab, X Computat Phys Div, POB 1663, Los Alamos, NM 87545 USA

关键词: ALE;    Remap;    Reconstruction;    Bounds-preserving;    Energy-conserving;   

DOI  :  10.1016/j.jcp.2017.11.017

来源: Elsevier

【 摘 要 】

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian-Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. In particular, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. The paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order. (c) 2017 Elsevier Inc. All rights reserved.

【 授权许可】

Free