High-order numerical simulation and modelling of the interaction of energetic and inert materials

We present an integrated algorithm on a Eulerian grid, for multimaterial simulations of energetic and inert materials modelled by non-ideal equations of state. We employ high-resolution shock capturing numerical algorithms for each material inside its domain and use an overlap domain method across the interface, maintained by a recently developed, hybrid, level-set algorithm. For applications to condensed explosives we implement a non-ideal, wide-ranging equation of state and reaction rate law. For inert materials, like plastic, metal, water, etc., we implement a (linear in the pressure) Mie–Grüneisen, (U _p ?U _s ), equation of state. We present a series of verifications of the integrated multimaterial code and show validations against experiment. We show examples of simulations of various experiments associated with real or planned experiments, some of which contain energetic materials (specifically the condensed explosives PBX-9502 and PBX-9501).     

Three-phase plane composites of minimal elastic stress energy: High-porosity structures

The paper establishes exact lower bound on the effective elastic energy of two-dimensional, three-material composite subjected to the homogeneous, anisotropic stress. It is assumed that the materials are mixed with given volume fractions and that one of the phases is degenerated to void, i.e., the effective composite is porous. Explicit formula for the energy bound is obtained using the translation method enhanced with additional inequality expressing certain property of stresses. Sufficient optimality conditions of the energy bound are used to set the requirements which have to be met by the stress fields in each phase of optimal effective material regardless of the complexity of its microstructural geometry. We show that these requirements are fulfilled in a special class of microgeometries, so-called laminates of a rank. Their optimality is elaborated in detail for structures with significant amount of void, also referred to as high-porosity structures. It is shown that geometrical parameters of optimal multi-rank, high-porosity laminates are different in various ranges of volume fractions and anisotropy level of external stress. Non-laminate, three-phase microstructures introduced by other authors and their optimality in high-porosity regions is also discussed by means of the sufficient conditions technique. Conjectures regarding low-porosity regions are presented, but full treatment of this issue is postponed to a separate publication. The corresponding “G-closure problem” of a three-phase isotropic composite is also addressed and exact bounds on effective isotropic properties are explicitly determined in these regions where the stress energy bound is optimal.     

Two-step hybrid conservative remapping for multimaterial arbitrary Lagrangian–Eulerian methods

We present a new hybrid conservative remapping algorithm for multimaterial Arbitrary Lagrangian–Eulerian (ALE) methods. The hybrid remapping is performed in two steps. In the first step, only nodes of the grid that lie inside subdomains occupied by single materials are moved. At this stage, computationally cheap swept-region remapping is used. In the second step, nodes that are vertices of mixed cells (cells containing several materials) and vertices of some cells in a buffer zone around mixed cells are moved. At this stage, intersection-based remapping is used. The hybrid algorithm results in computational expense that lies between swept-region and intersection-based remapping We demonstrate the performance of our new method for both structured and unstructured polygonal grids in two dimensions, as well as for cell-centered and staggered discretizations.