3 Matching Results

Search Results

Advanced search parameters have been applied.

Moving Overlapping Grids with Adaptive Mesh Refinement for High-Speed Reactive and Non-reactive Flow

Description: We consider the solution of the reactive and non-reactive Euler equations on two-dimensional domains that evolve in time. The domains are discretized using moving overlapping grids. In a typical grid construction, boundary-fitted grids are used to represent moving boundaries, and these grids overlap with stationary background Cartesian grids. Block-structured adaptive mesh refinement (AMR) is used to resolve fine-scale features in the flow such as shocks and detonations. Refinement grids are added to base-level grids according to an estimate of the error, and these refinement grids move with their corresponding base-level grids. The numerical approximation of the governing equations takes place in the parameter space of each component grid which is defined by a mapping from (fixed) parameter space to (moving) physical space. The mapped equations are solved numerically using a second-order extension of Godunov's method. The stiff source term in the reactive case is handled using a Runge-Kutta error-control scheme. We consider cases when the boundaries move according to a prescribed function of time and when the boundaries of embedded bodies move according to the surface stress exerted by the fluid. In the latter case, the Newton-Euler equations describe the motion of the center of mass of the each body and the rotation about it, and these equations are integrated numerically using a second-order predictor-corrector scheme. Numerical boundary conditions at slip walls are described, and numerical results are presented for both reactive and non-reactive flows in order to demonstrate the use and accuracy of the numerical approach.
Date: August 30, 2005
Creator: Henshaw, W D & Schwendeman, D W
Partner: UNT Libraries Government Documents Department

A Study of Detonation Diffraction in the Ignition-and-Growth Model

Description: Heterogeneous high-energy explosives are morphologically, mechanically and chemically complex. As such, their ab-initio modeling, in which well-characterized phenomena at the scale of the microstructure lead to a rationally homogenized description at the scale of observation, is a subject of active research but not yet a reality. An alternative approach is to construct phenomenological models, in which forms of constitutive behavior are postulated with an eye on the perceived picture of the micro-scale phenomena, and which are strongly linked to experimental calibration. Most prominent among these is the ignition-and-growth model conceived by Lee and Tarver. The model treats the explosive as a homogeneous mixture of two distinct constituents, the unreacted explosive and the products of reaction. To each constituent is assigned an equation of state, and a single reaction-rate law is prescribed for the conversion of the explosive to products. It is assumed that the two constituents are always in pressure and temperature equilibrium. The purpose of this paper is to investigate in detail the behavior of the model in situations where a detonation turns a corner and undergoes diffraction. A set of parameters appropriate for the explosive LX-17 is selected. The model is first examined analytically for steady, planar, 1-D solutions and the reaction-zone structure of Chapman-Jouguet detonations is determined. A computational study of two classes of problems is then undertaken. The first class corresponds to planar, 1-D initiation by an impact, and the second to corner turning and diffraction in planar and axisymmetric geometries. The 1-D initiation, although interesting in its own right, is utilized here as a means for interpretation of the 2-D results. It is found that there are two generic ways in which 1-D detonations are initiated in the model, and that these scenarios play a part in the post-diffraction evolution as well. For the parameter ...
Date: April 14, 2006
Creator: Kapila, A K; Schwendeman, D W; Bdzil, J B & Henshaw, W D
Partner: UNT Libraries Government Documents Department

A High-Resolution Godunov Method for Compressible Multi-Material Flow on Overlapping Grids

Description: A numerical method is described for inviscid, compressible, multi-material flow in two space dimensions. The flow is governed by the multi-material Euler equations with a general mixture equation of state. Composite overlapping grids are used to handle complex flow geometry and block-structured adaptive mesh refinement (AMR) is used to locally increase grid resolution near shocks and material interfaces. The discretization of the governing equations is based on a high-resolution Godunov method, but includes an energy correction designed to suppress numerical errors that develop near a material interface for standard, conservative shock-capturing schemes. The energy correction is constructed based on a uniform pressure-velocity flow and is significant only near the captured interface. A variety of two-material flows are presented to verify the accuracy of the numerical approach and to illustrate its use. These flows assume an equation of state for the mixture based on Jones-Wilkins-Lee (JWL) forms for the components. This equation of state includes a mixture of ideal gases as a special case. Flow problems considered include unsteady one-dimensional shock-interface collision, steady interaction of an planar interface and an oblique shock, planar shock interaction with a collection of gas-filled cylindrical inhomogeneities, and the impulsive motion of the two-component mixture in a rigid cylindrical vessel.
Date: February 13, 2006
Creator: Banks, J W; Schwendeman, D W; Kapila, A K & Henshaw, W D
Partner: UNT Libraries Government Documents Department