We extend a higher-order finite volume cell-centered hydrodynamic (CCH) formulation to include an interface-aware subscale closure model and a multi-material remap for simulating 3D compressible hydrodynamic problems within an arbitrary Lagrangian-Eulerian (ALE) framework. This CCH formulation involves a multidirectional approximate Riemann solution using quadratic polynomial reconstructions of the stress tensor and the velocity. At the subscale level, we determine pair-wise material interactions by solving a distinct approximate Riemann problem at the common interface, using the volume of fluids (VOF) method to find the interface. Material interactions are constrained to ensure smooth pressure equilibration among materials. The accuracy and robustness of the ALE method is demonstrated by simulating a suite of 3D Cartesian multi-material problems covering both gas and solid dynamics, where each test case has two or more materials.

We present a multi-material cell model (closure model) for demanding arbitrary Lagrangian-Eulerian (ALE) simulations of fluids and solids. It is based on the interface-aware sub-scale dynamics (IASSD) approach which utilizes the exact material interface geometry within the computational cell to calculate internal material interactions. Our formulation of the closure model also aims to improve the accuracy in low-speed impact events. Voids are used to represent ambient vacuum and internal free boundaries of the distinct materials. Void regions can close and open at contact surfaces, allowing a transition from contact physics to free motion in vacuum. The coupling of void closure and opening with a new formulation of the IASSD model for solids is tested on several one- and two-dimensional numerical examples, ranging from gas expansion in vacuum to planar and round object impacts at various speeds. (C) 2020 Published by Elsevier Ltd.

We present the new discrete optimization-based interface-aware sub-scale dynamics (IA-SSD) closure model for multimaterial cells for Lagrangian cell-centered hydrodynamics. For the multimaterial cell, the kinematic and thermodynamic properties ( e.g., velocity, density, pressure and internal energy) will typically vary between the materials. The discrete closure model is responsible for an accurate update of the thermodynamic states of the individual material components in the multimaterial cell, and for determining the nodal forces that move the vertices of the cell. The IA-SSD closure model consists of two stages — a bulk stage followed by a sub-scale stage. During the bulk stage, the total change in the volume of the cell, total force applied to the cell, and total work done on the cell are distributed between the materials to update their volume, velocity and total energy. This distribution is performed using volume fractions of the materials. During the second stage, sub-scale interactions of the materials inside the multimaterial cell are taken into account. At this stage, information about the topology of the materials inside the multimaterial cell is used, allowing the orientations of internal interfaces to be included in the model. Each material interacts in a pair-wise fashion with the materials with which it has a common boundary. The interactions are based on the solution of the acoustic Riemann problem between each pair of materials and are limited using physically justified constraints: positivity of volume, positivity of internal energy, and controlled rate of pressure relaxation. To determine the values of the limiter coefficients, a constrained-optimization framework is employed using a quadratic objective function with linear constraints. It is a first-of-its kind application of constrained optimization to develop discrete closure models in a more rigorous fashion. The pair-wise interaction between materials is essentially one dimensional in the direction that is normal to interface. For this reason, we demonstrate in this paper the performance of our new model on one dimensional numerical examples.

In hydrocodes, voids are used to represent vacuum and model free boundaries between vacuum and real materials. We give a systematic description of a new treatment of void closure in the framework of the multimaterial arbitrary Lagrangian–Eulerian (ALE) methods. This includes a new formulation of the interface-aware sub-scale-dynamics (IA-SSD) closure model for multimaterial cells with voids, which is used in the Lagrangian stage of our indirect ALE scheme. The results of the comprehensive testing of the new model are presented for one- and two-dimensional multimaterial calculations in the presence of voids. We also present a sneak peek of a realistic shaped charge calculation in the presence of voids and solids.

This paper reviews recent developments in Arbitrary Lagrangian Eulerian (ALE) methods for modeling high speed compressible multimaterial flows in complex geometry on general polygonal meshes. We only consider the indirect ALE approach which consists of three key stages: a Lagrangian stage, in which the solution and the computational mesh are updated; a rezoning stage, in which the nodes of the computational mesh are moved to improve grid quality; and a remapping stage, in which the Lagrangian solution is transferred to the rezoned mesh.

Arbitrary Lagrangian Eulerian methods are well suited to modelling high speed compressible multimaterial flow problems in complex geometry. In this talk three stage ALE schemes will be considered, which are typically used for practical calculations of multimaterial flows. These three stages are: a Lagrangian stage, in which the solution and the computational mesh are updated; a rezoning stage, in which the nodes of the computational mesh are moved to more optimal positions; and a remapping stage, in which the Lagrangian solution is transferred to the rezoned mesh. This talk will give a brief overview of these methods [1] and then focuss on two topics in more detail; cell centred Lagrangian schemes and closure models for multimaterial cells.

A systematic description of the new interface-aware sub-scale-dynamics (IA-SSD) closure model for the Lagrangian stage of multimaterial arbitrary Lagrangian–Eulerian methods is presented. The IA-SSD closure model consists of two stages. During the first, bulk, stage, the well known equal compressibility model is used. During the second stage, sub-scale interactions of the materials inside the multimaterial cell are taken into account. At this stage, information about the topology of the materials inside the multimaterial cell is utilized, allowing the orientations of internal interfaces to be included in the model. Each material interacts in a pair-wise fashion with the materials with which it has a common boundary. The interactions are based on the solution of the acoustic Riemann problem between each pair of materials and is limited using physically justified constraints: positivity of volume, positivity of internal energy and controlled rate of pressure relaxation. To determine the values of the limiter coefficients, a constrained-optimization framework is employed using a quadratic objective function with linear constraints. The algorithm guarantees the positivity of the material volume and internal energy as well as the smooth relaxation of the pressure – this allows a significant increase in the robustness of the overall algorithm.The results of comprehensive testing of the new model have been presented for one- and two-dimensional multimaterial Lagrangian hydrodynamics along with representative results for 2D multimaterial arbitrary Lagrangian–Eulerian (ALE) calculations. The numerical tests have shown that in most cases the new IA-SSD closure model produces better results compared to the well known Tipton's closure model.

stage, the total change in the volume of the cell, total force applied to the cell, and total work done on the cell are distributed between the materials to update their volume, velocity and total energy. This distribution is performed using volume fractions of the materials. During the second stage, sub-scale interactions of the materials inside the multimaterial cell are taken into account. At this stage, information about the topology of the materials inside the multimaterial cell is used, allowing the orientations of internal interfaces to be included in the model. Each material interacts in a pair-wise fashion with the materials with which it has a common boundary. The interactions are based on the solution of the acoustic Riemann problem between each pair of materials and are limited using physically justified constraints: positivity of volume, positivity of internal energy, and controlled rate of pressure relaxation. To determine the values of the limiter coefficients, a constrained-optimization framework is employed using a quadratic objective function with linear constraints. It is a first-of-its kind application of constrained optimization to develop discrete closure models in a more rigorous fashion. The pair-wise interaction between materials is essentially one dimensional in the direction that is normal to interface. Finally, for this reason, we demonstrate in this paper the performance of our new model on one dimensional numerical examples.