Simulation-driven optimization of high-order meshes in ALE hydrodynamics

•Mesh adaptivity is used to control numerical dissipation in ALE hydrodynamics.•Algebraic problem formulation applicable to curved meshes in any dimension.•Adaptivity to discrete simulation features is done through node movement.•General target construction methods can adapt all geometric properties...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Computers & fluids 2020-08, Vol.208, p.104602, Article 104602
Hauptverfasser: Dobrev, Veselin, Knupp, Patrick, Kolev, Tzanio, Mittal, Ketan, Rieben, Robert, Tomov, Vladimir
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:•Mesh adaptivity is used to control numerical dissipation in ALE hydrodynamics.•Algebraic problem formulation applicable to curved meshes in any dimension.•Adaptivity to discrete simulation features is done through node movement.•General target construction methods can adapt all geometric properties.•Adaptive ALE triggers can improve both accuracy and simulation time. 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 [1]. 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 [2].
ISSN:0045-7930
1879-0747
DOI:10.1016/j.compfluid.2020.104602