Simulation-Driven Optimization of High-Order Meshes in ALE Hydrodynamics

Abstract

In this paper we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP) of Knupp. The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

arXiv comment26 pages, 10 figures. Submitted to Computers & Fluids

References in corpus (4)

• MFEM: a modular finite element methods library
• Direct Arbitrary-Lagrangian-Eulerian finite volume schemes on moving nonconforming unstructured meshes
• Nonconforming Schwarz-Spectral Element Methods For Incompressible Flow
• $rp$-adaptation for compressible flows

Cited by in corpus (9)

• Efficient Exascale Discretizations: High-Order Finite Element Methods
• High-performance finite elements with MFEM
• hr-adaptivity for nonconforming high-order meshes with the target matrix optimization paradigm
• Accelerating High-Order Mesh Optimization Using Finite Element Partial Assembly on GPUs
• High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries
• Combining high-order metric interpolation and geometry implicitization for curved $r$-adaption
• High order well-balanced Arbitrary-Lagrangian-Eulerian ADER discontinuous Galerkin schemes on general polygonal moving meshes
• A globalized and preconditioned Newton-CG solver for metric-aware curved high-order mesh optimization
• Adaptive Surface Fitting and Tangential Relaxation for High-Order Mesh Optimization