Intersection-distribution-based remapping between arbitrary meshes for staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics

JOURNAL OF COMPUTATIONAL PHYSICS

卷:429

Article

Kenamond, Mack¹  Kuzmin, Dmitri²  Shashkov, Mikhail³ 

[1] Alamos Natl Lab, X Computat Phys, XCP 1, Los Alamos, NM 87545 USA

[2] TU Dortmund Univ, Inst Appl Math LS III, Vogelpothsweg 87, D-44227 Dortmund, Germany

[3] Alamos Natl Lab, X Computat Phys, XCP 4, Los Alamos, NM 87545 USA

关键词: Staggered arbitrary-Lagrangian-Eulerian hydrodynamics;    Multi-material remapping;    Arbitrary polygonal meshes;   

DOI  :  10.1016/j.jcp.2020.110014

来源: Elsevier

【 摘 要 】

We present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution of all material masses in each new cell to the cell's corners and then collecting those corner masses at new nodes. This distribution is done using a local constrained optimization approach for each cell in the new mesh. In order to remap nodal momentum we first define cell-centered momentum for each cell in the old mesh, conservatively remap this to the new mesh and then conservatively distribute the new zonal momentum to each cell's bounding nodes, again using local constrained optimization. Our method also conserves total energy by applying a new nodal kinetic energy correction that relies on a process similar to that used for remapping nodal mass and momentum. Cell-centered kinetic energy is computed, conservatively remapped and then distributed to nodes. The discrepancy between this conservatively remapped and actual nodal kinetic energy is then conservatively distributed to the internal energies of the materials in the cells surrounding each node. Unlike conventional cell-based corrections of this type, this new nodal kinetic energy correction has not been observed to drive material internal energy negative in any of our testing. Unlike flux based remapping, our new intersection-distribution method can be applied to remapping between source and target meshes that are arbitrarily different, which provides superior flexibility in the rezoning strategy. Our method is accurate, essentially conservative and essentially bounds preserving. (C) 2020 Elsevier Inc. All rights reserved.

【 授权许可】

Free