Search

This work is concerned with the numerical simulation of ablation of geological materials using a millimetre wave source. To this end, a new mathematical model is developed for a thermal approach to the problem, allowing for large scale simulations, whilst being able to include the strong temperature dependence of material parameters to ensure accurate modelling of power input into the rock. The model presented is implemented within an adaptive meshing framework, such that resolution can be placed where needed, for example at the borehole wall, to further improve the computational efficiency of large scale simulations. This approach allows for both the heating of the rock, and the removal of evaporated material, allowing rate of penetration and the shape of the resulting borehole to be quantified. The model is validated against experimental results, which indicates that the approach can accurately predict temperatures, and temperature gradients within the rock. The validated model is then exercised to obtain initial results demonstrating its capabilities for simulating the millimetre wave drilling process. The effects of the conditions at the surface of the rock are investigated, highlighting the importance of understanding the physical processes which occur between the wave guide and the rock. Additionally, the absorptivity of the rock, and the impact this has on the evaporation behaviour is considered. Simulations are carried out both for isotropic rock, and also for a multi-strata configuration. It is found that strata between similar rock types, such as granite and basalt, absorptive properties pose little problem for uniform drilling. However, larger variations in material parameters are shown to have strong implications on the evaporation behaviour of the wellbore, and hence the resulting structure., Comment: 13 Pages, 12 figures

This work outlines a diffuse interface method for the study of fracture and fragmentation in ductile metals at high strain-rates in Eulerian finite volume simulations. The work is based on an existing diffuse interface method capable of simulating a broad range of different multi-physics applications, including multi-material interaction, damage and void opening. The work at hand extends this method with a technique to model realistic material inhomogeneities, and examines the performance of the method on a selection of challenging problems. Material inhomogeneities are included by evolving a scalar field that perturbs a material's plastic yield stress. This perturbation results in non-uniform fragments with a measurable statistical distribution, allowing for underlying defects in a material to be modelled. As the underlying numerical scheme is three dimensional, parallelisable and multi-physics-capable, the scheme can be tested on a range of strenuous problems. These problems especially include a three-dimensional explosively driven fracture study, with an explicitly resolved condensed phase explosive. The new scheme compares well with both experiment and previous numerical studies.

This work is concerned with the development of a numerical modelling approach for studying the time-accurate response of aerospace fasteners subjected to high electrical current loading from a simulated lightning strike. The electromagnetic, thermal and elastoplastic response of individual fastener components is captured by this method allowing a critical analysis of fastener design and material layering. Under high electrical current loading, ionisation of gas filled cavities in the fastener assembly can lead to viable current paths across internal voids. This ionisation can lead to localised pockets of high pressure plasma through the Joule heating effect. The multi-physics approach developed in this paper extends an existing methodology that allows a two-way dynamic non-linear coupling of the plasma arc, the titanium aerospace fastener components, the surrounding aircraft panels, the internal supporting structure and internal plasma-filled cavities. Results from this model are compared with experimental measurements of a titanium fastener holding together carbon composite panels separated by thin dielectric layers. The current distribution measurements are shown to be accurately reproduced. A parameter study is used to assess the internal cavity modelling strategy and to quantify the relation between the internal cavity plasma pressure, the electrical current distribution and changes in the internal cavity geometry., Comment: 16 pages, 17 figures

This work outlines a new multi-physics-compatible immersed rigid body method for Eulerian finite-volume simulations. To achieve this, rigid bodies are represented as a diffuse scalar field and an interface seeding method is employed to mediate the interface boundary conditions. The method is based on an existing multi-material diffuse interface method that is capable of handling an arbitrary mixture of fluids and elastoplastic solids. The underlying method is general and can be extended to a range of different applications including high-strain rate deformation in elastoplastic solids and reactive fluid mixtures. As such, the new method presented here is thoroughly tested through a variety of problems, including fluid-rigid body interaction, elastoplastic-rigid body interaction, and detonation-structure interaction. Comparison is drawn between both experimental work and previous numerical results, with excellent agreement in both cases. The new method is straightforward to implement, highly local, and parallelisable. This allows the method to be employed in three dimensions with multiple levels of adaptive mesh refinement using complex immersed geometries. The rigid body field can be static or dynamic, with the THINC interface reconstruction method being used to keep the interface sharp in the dynamic case.

In this article we present a new formulation and an associated algorithm for the simultaneous numerical simulation of multiple condensed phase explosives in direct contact with each other, which may also be confined by (or interacting with one or more) compliant inert materials. Examples include composite rate-stick problems and interaction of shock waves with chemically-active particles in condensed-phase explosives. There are several formulations which address the compliant or structural response of confiners and particles due to detonations, but the direct interaction of explosives remains a challenge for most formulations and algorithms. The proposed formulation addresses this problem by extending the conservation laws and mixture rules of an existing hybrid formulation to model the interaction of multiple explosive mixtures. An algorithm for the solution of the resulting system of partial differential equations is presented, which includes a new robust method for the retrieval of the densities of the constituents of each explosive mixture. The algorithm is implemented in a hierarchical adaptive mesh refinement framework and validated against results from problems with known solutions. It is evaluated for robustness for rate-stick and shock-induced flows in particle-laden explosives case-studies. It is shown that the method can simulate the interaction of detonation waves produced by military grade and commercial explosives in direct contact, each with its own distinct equation of state and reaction rate law. The ability of the new model to simulate reactive particles which are explicitly resolved in a heterogeneous explosive is demonstrated by a case-study of a shock wave interacting with a high explosive bead embedded in liquid nitromethane.

The motivation of this work is to produce an integrated formulation for material response due to detonation wave loading. Here, we focus on elastoplastic structural response. In particular, we are interested to capture miscible and immiscible behaviour within condensed-phase explosives arising from the co-existence of a reactive carrier mixture of miscible materials, and several material interfaces due to the presence of immiscible impurities such as particles or cavities. The dynamic and thermodynamic evolution of the explosive is communicated to one or more inert confiners through their shared interfaces, which may undergo severe topological change. We also wish to consider elastic and plastic structural response of the confiners, rather than make a hydrodynamic assumption for their behaviour. Previous work by these authors has met these requirements by means of the simultaneous solution of appropriate systems of equations for the behaviour of the condensed-phase explosive and the elastoplastic behaviour of the confiners. In the present work, we employ a single system of partial differential equations (PDEs) proposed by Peshkov and Romenski, which is able to account for different states of matter by means of generalising the concept of distortion tensors beyond solids. We amalgamate that formulation with a single system of PDEs which meets the requirement of co-existing miscible and immiscible explosive mixtures. We present the mathematical derivation and construct appropriate algorithms for its solution. The resulting model is validated against exact solutions for several use-cases, including mechanically- and thermally-induced, inviscid and viscous detonations. Results indicate that the model can accurately simulate a very broad range of problems involving the nonlinear interaction between reactive and inert materials within a single framework., Comment: 28 pages, 9 figures

GPUs are now used for a wide range of problems within HPC. However, making efficient use of the computational power available with multiple GPUs is challenging. The main challenges in achieving good performance are memory layout, affecting memory bandwidth, effective use of the memory spaces with a GPU, inter-GPU communication, and synchronization. We address these problems with the Ripple library, which provides a unified view of the computational space across multiple dimensions and multiple GPUs, allows polymorphic data layout, and provides a simple graph interface to describe an algorithm from which inter-GPU data transfers can be optimally scheduled. We describe the abstractions provided by Ripple to allow complex computations to be described simply, and to execute efficiently across many GPUs with minimal overhead. We show performance results for a number of examples, from particle motion to finite-volume methods and the eikonal equation, as well as showing good strong and weak scaling results across multiple GPUs., Comment: Preprint submitted to the Journal of Parallel and Distributed Computing

This work is concerned with the development of a novel, accurate equation of state for describing partially ionised air plasma in local thermodynamic equilibrium. One key application for this new equation of state is the simulation of lightning strike on aircraft. Due to the complexities of species ionisation and interaction, although phenomenological curve fitting of thermodynamic properties is possible, these curves are intractable for practical numerical simulation. The large difference in size of the parameters (many orders of magnitude) and complexity of the equations means they are not straightforward to invert for conversion between thermodynamic variables. The approach of this paper is to take an accurate 19-species phenomenological model, and use this to generate a tabulated data set. Coupled with a suitable interpolation procedure this offers an accurate and computationally efficient technique for simulating partially ionised air plasma. The equation of state is implemented within a multiphysics methodology which can solve for two-way coupling between a plasma arc and an elastoplastic material substrate. The implementation is validated against experimental results, both for a single material plasma, and an arc coupled to a substrate. It is demonstrated that accurate, oscillation-free thermodynamic profiles can be obtained, with good results even close to material surfaces., Comment: 12 pages, 12 figures

This work outlines a new three-dimensional diffuse interface finite volume method for the simulation of multiple solid and fluid components featuring large deformations, sliding and void opening. This is achieved by extending an existing reduced-equation diffuse interface method by means of a number of novel flux-modifiers and interface seeding routines that enable the application of different material boundary conditions. The method allows for slip boundary conditions across solid interfaces, material-void interaction, and interface separation. The method is designed to be straightforward to implement, inexpensive and highly parallelisable. This makes it suitable for use in large, multi-dimensional simulations that feature many complex materials and physical processes interacting over multiple levels of adaptive mesh refinement. Furthermore, the new method allows for the generation of new interfaces in a conservative fashion and therefore naturally facilitates the simulation of high-strain rate fracture. Hence, the governing model is augmented to include ductile damage to allow for validation of the method against demanding physical experiments. The method is shown to give excellent agreement with both experiment and existing Eulerian interface tracking algorithms that employ sharp interface methods., Comment: arXiv admin note: text overlap with arXiv:2009.14555

This article presents a multi-physics methodology for the numerical simulation of physical systems that involve the non-linear interaction of multi-phase reactive fluids and elastoplastic solids, inducing high strain-rates and high deformations. Each state of matter is governed by a single system of non-linear, inhomogeneous partial differential equations, which are solved simultaneously on the same computational grid, and do not require special treatment of immersed boundaries. To this end, the governing equations for solid and reactive multiphase fluid mechanics are written in the same mathematical form and are discretised on a regular Cartesian mesh. All phase and material boundaries are treated as diffuse interfaces. An interface-steepening technique is employed at material boundaries to keep interfaces sharp whilst maintaining the conservation properties of the system. These algorithms are implemented in a highly-parallelised hierarchical adaptive mesh refinement platform, and are verified and validated using numerical and experimental benchmarks. Results indicate very good agreement with experiment and an improvement of numerical performance compared to certain existing Eulerian methods, without loss of conservation.

The sensitizing effect of cavities in the form of microbubbles on the shock initiation of a homogeneous liquid explosive is studied computationally. While the presence of voids in an explosive has long been known to induce so-called hot spots that greatly accelerate the global reaction rate, the ability to computationally resolve the details of the interaction of the shock front with heterogeneities existing on the scale of the detonation reaction zone has only recently become feasible. In this study, the influence of the spatial distribution of air-filled cavities has been examined, enabled by the use of graphic processing unit (GPU) accelerated computations that can resolve shock initiation and detonation propagation through an explosive while fully resolving features at the mesoscale. Different spatial distributions of cavities are examined in two-dimensional simulations, including regular arrays of cavities, slightly perturbed arrays, random arrays (with varying minimum spacing being imposed on the cavities), and randomly distributed clusters of cavities. The presence of the cavities is able to reduce the time required to initiate detonation---for a given input shock strength---by greater than 50%, in agreement with previous experimental results. Randomly distributing the cavities results in a 15-20% decrease in detonation initiation time in comparison to a regular array of cavities. Clustering the cavities---as would occur in the case of agglomeration---results in an additional 10% decrease in detonation initiation time in comparison to random arrays. The effect of clustering is shown not to be a result of the clusters forming an effectively larger cavity, but rather due to interactions between clusters upon shock loading occurring on the microscale.

This work is concerned with the numerical simulation of elastoplastic, electromagnetic and thermal response of aerospace materials due to their interaction with a plasma arc under lightning strike conditions. Current approaches treat the interaction between these two states of matter either in a decoupled manner or through one-way coupled 'co-simulation'. In this paper a methodology for multiphysics simulations of two-way interaction between lightning and elastoplastic materials is presented, which can inherently capture the non-linear feedback between these two states of matter. This is achieved by simultaneously solving the magnetohydrodynamic and the elastoplastic systems of equations on the same computational mesh, evolving the magnetic and electric fields dynamically. The resulting model allows for the topological evolution and movement of the arc attachment point coupled to the structural response and Joule heating of the substrate. The dynamic communication between the elastoplastic material and the plasma is facilitated by means of Riemann problem-based ghost fluid methods. This two-way coupling, to the best of the authors' knowledge, has not been previously demonstrated. The proposed model is first validated against experimental laboratory studies, demonstrating that the growth of the plasma arc can be accurately reproduced, dependent on the electrical conductivity of the substrate. It is then subsequently evaluated in a setting where the dynamically-evolved properties within the substrate feed back into the plasma arc attachment. Results are presented for multi-layered substrates of different materials, and for a substrate with temperature-dependent electrical conductivity. It is demonstrated that these conditions generate distinct behaviour due to the interaction between the plasma arc and the substrate., Comment: 20 pages, 16 figures

We present a methodology for simulating three-dimensional flow of incompressible viscoplastic fluids modelled by generalised Newtonian rheological equations. It is implemented in a highly efficient framework for massively parallelisable computations on block-structured grids. In this context, geometric features are handled by the embedded boundary approach, which requires specialised treatment only in cells intersecting or adjacent to the boundary. This constitutes the first published implementation of an embedded boundary algorithm for simulating flow of viscoplastic fluids on structured grids. The underlying algorithm employs a two-stage Runge-Kutta method for temporal discretisation, in which viscous terms are treated semi-implicitly and projection methods are utilised to enforce the incompressibility constraint. We augment the embedded boundary algorithm to deal with the variable apparent viscosity of the fluids. Since the viscosity depends strongly on the strain rate tensor, special care has been taken to approximate the components of the velocity gradients robustly near boundary cells, both for viscous wall fluxes in cut cells and for updates of apparent viscosity in cells adjacent to them. After performing convergence analysis and validating the code against standard test cases, we present the first ever fully three-dimensional simulations of creeping flow of Bingham plastics around translating objects. Our results shed new light on the flow fields around these objects.

We propose a numerical methodology for the simultaneous numerical simulation of four states of matter; gas, liquid, elastoplastic solids and plasma. The distinct, interacting physical processes are described by a combination of compressible, inert and reactive forms of the Euler equations, multiphase equations, elastoplastic equations and resistive MHD equations. Combinations of systems of equations are usually solved by coupling finite element for solid modelling and CFD models for fluid modelling or including material effects through boundary conditions rather than full material discretisation. Our simultaneous solution methodology lies on the recasting of all the equations in the same, hyperbolic form allowing their solution on the same grid with the same finite-volume numerical schemes. We use a combination of sharp and diffuse interface methods to track or capture material interfaces, depending on the application. The communication between the distinct systems of equations (i.e., materials separated by sharp interfaces) is facilitated by means of mixed-material Riemann solvers at the boundaries of the systems, which represent physical material boundaries. To this end we derive approximate mixed Riemann solvers for each pair of the above models based on characteristic equations. To demonstrate the applicability of the new methodology we consider a case study where we investigate the possibility of ignition of a combustible gas that lies over a liquid in a metal container that is struck by a plasma-arc akin to a lightning strike. We study the effect on the ignition of the metal container material and conductivity, of the presence of a dielectric coating, of insensitive combustible gases and sealed and pre-damaged metal surfaces., Comment: arXiv admin note: text overlap with arXiv:1710.01643

Two-dimensional, meso-resolved numerical simulations are performed to investigate the complete shock-to-detonation transition (SDT) process in a mixture of liquid nitromethane (NM) and air-filled, circular cavities. The shock-induced initiation behaviors resulting from the cases with neat NM, NM with an array of regularly spaced cavities, and NM with randomly distributed cavities are examined. For the case with randomly distributed cavities, hundreds of cavities are explicitly resolved in the simulations using a diffuse-interface approach to treat two immiscible fluids and GPU-enabled parallel computing. Without invoking any empirically calibrated, phenomenological models, the reaction rate in the simulations is governed by Arrhenius kinetics. For the cases with neat NM, the resulting SDT process features a superdetonation that evolves from a thermal explosion after a delay following the passage of the incident shock wave and eventually catches up with the leading shock front. For the cases wherein mesoscale heterogeneities are explicitly considered, a gradual SDT process is captured. These two distinct initiation behaviors for neat NM and heterogeneous NM mixtures agree with experimental findings. Via examining the global reaction rate of the mixture, a unique time scale characterizing the SDT process, i.e., the overtake time, is measured for each simulation. For an input shock pressure less than approximately $9.4~\mathrm{GPa}$, the overtake time resulting from a heterogeneous mixture is shorter than that for neat NM. This sensitizing effect is more pronounced for lower input shock pressures. A random distribution of cavities is found to be more effective in enhancing the SDT process than a regular array of cavities. Statistical analysis on the meso-resolved simulation data provides more insights into the mechanism of energy release underlying the SDT process., Comment: 23 pages, 17 figures

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR code (Incompressible Adaptive Mesh Refinement), which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve on adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

Direct numerical simulations of closely interacting infinite circular cylinders in a Bingham fluid are presented, and results compared to asymptotic solutions based on lubrication theory in the gap. Unlike for a Newtonian fluid, the macroscopic flow outside of the gap between the cylinders is shown to have a large effect on the pressure profile within the gap and the resulting lubrication force on the cylinders. The presented results indicate that the asymptotic lubrication solution can be used to predict the lubrication pressure only if the surrounding viscoplastic matrix is yielded by a macroscopic flow. This has implications for the use of sub-grid-scale lubrication models in simulations of non-colloidal particulate suspensions in viscoplastic fluids.

We present a numerical study of detonation propagation in unconfined explosive charges shaped as an annular arc (rib). Steady detonation in a straight charge propagates at constant speed but when it enters an annular section, it goes through a transition phase and eventually reaches a new steady state of constant angular velocity. This study examines the speed of the detonation wave along the annular charge during the transition phase and at steady state, as well as its dependence on the dimensions of the annulus. The system is modeled using a recently proposed diffuse-interface formulation which allows for the representation of a two-phase explosive and of an additional inert material. The explosive considered is the polymer-bonded TATB-based LX-17 and is modeled using two JWL equations of state and the Ignition and Growth reaction rate law. Results show that steady state speeds are in good agreement with experiment. In the transition phase, the evolution of outer detonation speed deviates from the exponential bounded growth function suggested by previous studies. We propose a new description of the transition phase which consists of two regimes. The first is caused by local effects at the outer edge of the annulus and leads to a dependence of outer detonation speed on angular position along the arc. The second regime is induced by effects originating from the inner edge of the annular charge and leads to the deceleration of the outer detonation until steady state is reached. The study concludes with a parametric study where the dependence of the steady state and the transition phase on the dimensions of the annulus is investigated.

The objective of this work is to improve the robustness and accuracy of numerical simulations of both ideal and non-ideal explosives by introducing temperature dependence in mechanical equations of state for reactants and products. To this end, we modify existing mechanical equations of state to appropriately approximate the temperature in the reaction zone. Mechanical equations of state of Mie-Gr\"{u}neisen form are developed with extensions, which allow the temperature to be evaluated appropriately, and the temperature equilibrium condition to be applied robustly. Furthermore the snow plow model is used to capture the effect of porosity on the reactants equation of state. We apply the methodology to predict the velocity of compliantly confined detonation waves. Once reaction rates are calibrated for unconfined detonation velocities, simulations of confined rate sticks and slabs are performed, and the experimental detonation velocities are matched without further parameter alteration, demonstrating the predictive capability of our simulations. We apply the same methodology to both ideal (PBX9502, a high explosive with principal ingredient TATB) and non-ideal (EM120D, an ANE or ammonium nitrate based emulsion) explosives.

We study effect of cavity collapse in non-ideal explosives as a means of controlling their sensitivity. The main aim is to understand the origin of localised temperature peaks (hot spots) that play a leading order role at early ignition stages. Thus, we perform 2D and 3D numerical simulations of shock induced single gas-cavity collapse in nitromethane. Ignition is the result of a complex interplay between fluid dynamics and exothermic chemical reaction. In part I of this work we focused on the hydrodynamic effects in the collapse process by switching off the reaction terms in the mathematical model. Here, we reinstate the reactive terms and study the collapse of the cavity in the presence of chemical reactions. We use a multi-phase formulation which overcomes current challenges of cavity collapse modelling in reactive media to obtain oscillation-free temperature fields across material interfaces to allow the use of a temperature-based reaction rate law. The mathematical and physical models are validated against experimental and analytic data. We identify which of the previously-determined (in part I of this work) high-temperature regions lead to ignition and comment on their reactive strength and reaction growth rate. We quantify the sensitisation of nitromethane by the collapse of the cavity by comparing ignition times of neat and single-cavity material; the ignition occurs in less than half the ignition time of the neat material. We compare 2D and 3D simulations to examine the change in topology, temperature and reactive strength of the hot spots by the third dimension. It is apparent that belated ignition times can be avoided by the use of 3D simulations. The effect of the chemical reactions on the topology and strength of the hot spots in the timescales considered is studied by comparing inert and reactive simulations and examine maximum temperature fields and their growth rates., Comment: Due to the limitation "The abstract field cannot be longer than 1,920 characters", the abstract appearing here is slightly shorter than the one in the PDF file

We study the effect of cavity collapse in non-ideal explosives as a means of controlling their sensitivity. The aim is to understand the origin of localised temperature peaks (hot spots) which play a key role at the early stages of ignition. Thus we perform 2D and 3D numerical simulations of shock induced gas-cavity collapse in nitromethane. Ignition is the result of a complex interplay between fluid dynamics and exothermic chemical reaction. To understand the relative contribution between these two processes we consider in this first part of the work the evolution of the physical system in the absence of chemical reactions. We employ a multi-phase mathematical formulation which accounts for the large density difference across the gas-liquid interface without generating spurious temperature peaks. The mathematical and physical models are validated against experimental, analytic and numerical data. Previous studies identified the impact of the upwind side of the cavity wall to the downwind one as the main reason for the generation of a hot-spot outside of the cavity; this is also observed in this work. However, it is apparent that the topology of the temperature field is more complex than previously thought and additional hot spots locations exist, arising from the generation of Mach stems rather than jet impact. To explain the generation mechanisms and topology of the hot spots we follow the complex wave patterns generated and identify the temperature elevation or reduction generated by each wave. This allows to track each hot spot back to its origins. We show that the highest hot spot temperatures can be more than twice the post-incident shock temperature of the neat material and can thus lead to ignition. By comparing the maximum temperature observed in the domain in 2D and 3D simulations we show that 3D calculations are necessary to avoid belated ignition times in reactive scenarios., Comment: Due to the limitation "The abstract field cannot be longer than 1,920 characters", the abstract appearing here is slightly shorter than the one in the PDF file

We propose a numerical methodology for the numerical simulation of distinct, interacting physical processes described by a combination of compressible, inert and reactive forms of the Euler equations, multiphase equations and elastoplastic equations. These systems of equations are usually solved by coupling finite element and CFD models. Here we solve them simultaneously, by recasting all the equations in the same, hyperbolic form and solving them on the same grid with the same finite-volume numerical schemes. The proposed compressible, multiphase, hydrodynamic formulation can employ a hierarchy of five reactive and non-reactive flow models, which allows simple to more involved applications to be directly described by the appropriate selection. The communication between the hydrodynamic and elastoplastic systems is facilitated by means of mixed-material Riemann solvers at the boundaries of the systems, which represent physical material boundaries. To this end we derive approximate mixed Riemann solvers for each pair of the above models based on characteristic equations. The components for reactive flow and elastoplastic solid modelling are validated separately before presenting validation for the full, coupled systems. Multi-dimensional use cases demonstrate the suitability of the reactive flow-solid interaction methodology in the context of impact-driven initiation of reactive flow and structural response due to violent reaction in automotive (e.g. car crash) or defence (e.g. explosive reactive armour) applications. Several types of explosives (C4, Deetasheet, nitromethane, gaseous fuel) in gaseous, liquid and solid state are considered.

Detonation propagation in a compressible medium wherein the energy release has been made spatially inhomogeneous is examined via numerical simulation. The inhomogeneity is introduced via step functions in the reaction progress variable, with the local value of energy release correspondingly increased so as to maintain the same average energy density in the medium, and thus a constant Chapman Jouguet (CJ) detonation velocity. A one-step Arrhenius rate governs the rate of energy release in the reactive zones. The resulting dynamics of a detonation propagating in such systems with one-dimensional layers and two-dimensional squares are simulated using a Godunov-type finite-volume scheme. The resulting wave dynamics are analyzed by computing the average wave velocity and one-dimensional averaged wave structure. In the case of sufficiently inhomogeneous media wherein the spacing between reactive zones is greater than the inherent reaction zone length, average wave speeds significantly greater than the corresponding CJ speed of the homogenized medium are obtained. If the shock transit time between reactive zones is less than the reaction time scale, then the classical CJ detonation velocity is recovered. The spatio-temporal averaged structure of the waves in these systems is analyzed via a Favre averaging technique, with terms associated with the thermal and mechanical fluctuations being explicitly computed. The analysis of the averaged wave structure identifies the super-CJ detonations as weak detonations owing to the existence of mechanical non-equilibrium at the effective sonic point embedded within the wave structure. The correspondence of the super-CJ behavior identified in this study with real detonation phenomena that may be observed in experiments is discussed.