Engineering Papers⌕ Search

Engineering topics

Quirk, James J.

Publications and source records attributed to Quirk, James J..

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Detonation propagation for shock-driven, subsonic and supersonic confiner flow

In this work, we study the compressible flow dynamics of two-dimensional, steady detonation wave propagation in a high explosive (HE) confined by aluminium (Al) or stainless steel (SS), outside of which is an air layer. We examine how the thickness of the confinement affects the subsonic detonation driving zone structure (DDZ) and the detonation speed (D 0 ), demonstrating a strong dependence on whether the oblique shock-driven flow in the confiner is supersonic, as for SS, or subsonic, as for Al. A characteristic path analysis is used to examine the information flow from the material boundaries through the supersonic flow regions in both the HE and confiner that can impact the sonic surfaces bounding the subsonic flow regions. It is shown that the nature of gas-dynamic wave reflection off the SS–air or Al–air boundary can significantly influence the DDZ and D 0 .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Towards an Analysis of the GhostFluid Method

The GhostFluid method is built on the premise that a material interface, such as that shown in Figure 1 (a), can be discretized as two distinct materials that communicate, indirectly, via ghostcells. The situation is shown pictorially in Figure 1 (b). The attraction of the approach is that the two halves of the problem can be integrated using a regular, single material scheme, with none of the complex, book-keeping normally associated with a multi-material solver. On the down side, the reliability of the approach rests with the filling of the ghost states. Specifically, the states must be chosen such that the resultant evolution is consistent with the full problem, and for non-linear systems, there is no guarantee that satisfactory states actually exist.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Propagation of a stable gaseous detonation in a circular arc configuration

Here, we examine numerically the effect that detonation wavefront curvature has on gaseous detonation cellular structure in a circular arc geometry with rigid walls. For the model considered, a planar detonation is weakly unstable. The detonation curvature is induced by diffraction and flow divergence as the wave sweeps around the arc. We examine the effect of arc thickness, inner arc radius and activation energy on the cellular pattern. A combination of pseudo-Schlieren, pseudo-soot foil records and detonation shock and sonic flow loci are used to explore the observed detonation structures. We find that for sufficiently wide arcs, the detonation wavefront curvature renders the detonation hydrodynamically stable, i.e. laminar flow with no cells. The stable configuration consists of a broadly curved wave originating from the inner arc surface, which is typically connected to the outer arc boundary by a stable Mach stem. The stable detonation is driven by a small region of subsonic flow attached to the inner arc surface, outside of which the flow is supersonic. The initiation dynamics of a stable propagating wave case are explored. The evolution dynamics between an unstable (cellular) and stable arc thickness are also examined. Finally, smaller inner arc radii and lower activation energies increase the range of stable solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Analysis of steady compaction waves in polyurea aerogel

In this work, steady compaction waves in an inert porous material are investigated using a p - α model. In a steady traveling wave reference frame, the one-dimensional Euler equations are reduced to a set of ordinary differential equations. A Mie-Grüneisen equation of state (EOS) is used with parameters calibrated for polyurea aerogel (PUA). Analytic solutions for non-equilibrium compaction are developed which compliment numerical simulations and are able to predict the complete wave structure, including the compaction wave speed, zone length, and final compacted solid volume fraction. The dynamic compaction of PUA is studied for a range of piston velocities. Three regions of behavior are identified: supersonic, subsonic-complete, and subsonic-partial compaction. Below a critical piston velocity, a subsonic compaction wave is produced without a leading shock. Finally, at even lower piston velocities, there is partial compaction and a greater dependence on the dynamic compaction relation. Some features and limitations of the current model are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Detonation diffraction in a circular arc geometry of the insensitive high explosive PBX 9502

We describe the details of an unconfined insensitive high explosive (PBX 9502) circular arc section experiment, in which, after a transient period, a detonation sweeps around the arc with constant angular speed. The arc section is sufficiently wide that the flow along the centerline of the arc section remains two-dimensional. Data includes time-of-arrival diagnostics of the detonation along the centerline inner and outer arc surfaces, which is used to obtain the angular speed of the steadily rotating detonation. We also obtain the lead shock shape of the detonation as it sweeps around the arc. Reactive burn model simulations of the PBX 9502 arc experiment are then conducted to establish the structure of the detonation driving zone, i.e. the region enclosed between the detonation shock and flow sonic locus (in the frame of the steady rotating detonation). It is only the energy released in this zone which determines the speed at which the steady detonation sweeps around the arc. We show that the sonic flow locus of the detonation driving zone largely lies at the end of, or within, the fast reaction stage of the PBX 9502 detonation, with the largest section of the detonation driving zone lying close to the inner arc surface. In conclusion, we also demonstrate that the reactive burn model provides a good prediction of both the angular speed of the detonation wave and the curved detonation front shape.

42 ENGINEERING↗
DOE OSTI2018

High Explosive Detonation–Confiner Interactions

The primary purpose of a detonation in a high explosive (HE) is to provide the energy to drive a surrounding confiner, typically for mining or munitions applications. The details of the interaction between an HE detonation and its confinement are essential to achieving the objectives of the explosive device. For the high pressures induced by detonation loading, both the solid HE and confiner materials will flow. The structure and speed of a propagating detonation, and ultimately the pressures generated in the reaction zone to drive the confiner, depend on the induced flow both within the confiner and along the HE–confiner material interface. The detonation–confiner interactions are heavily influenced by the material properties and, in some cases, the thickness of the confiner. This review discusses the use of oblique shock polar analysis as a means of characterizing the possible range of detonation–confiner interactions. Computations that reveal the fluid mechanics of HE detonation–confiner interactions for finite reaction-zone length detonations are discussed and compared with the polar analysis. Furthermore this includes cases of supersonic confiner flow; subsonic, shock-driven confiner flow; subsonic, but shockless confiner flow; and sonic flow at the intersection of the detonation shock and confiner material interface. We also summarize recent developments, including the effects of geometry and porous material confinement, on detonation–confiner interactions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1995

Some observations on mesh refinement schemes applied to shock wave phenomena

This workshop's double-wedge test problem is taken from one of a sequence of experiments which were performed in order to classify the various canonical interactions between a planar shock wave and a double wedge. Therefore to build up a reasonably broad picture of the performance of our mesh refinement algorithm we have simulated three of these experiments and not just the workshop case. Here, using the results from these simulations together with their experimental counterparts, we make some general observations concerning the development of mesh refinement schemes for shock wave phenomena.

Quirk, James J.↗
NASA NTRS1994

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗
NASA NTRS1994

A Cartesian grid approach with hierarchical refinement for compressible flows

Many numerical studies of flows that involve complex geometries are limited by the difficulties in generating suitable grids. We present a Cartesian boundary scheme for two-dimensional, compressible flows that is unfettered by the need to generate a computational grid and so it may be used, routinely, even for the most awkward of geometries. In essence, an arbitrary-shaped body is allowed to blank out some region of a background Cartesian mesh and the resultant cut-cells are singled out for special treatment. This is done within a finite-volume framework and so, in principle, any explicit flux-based integration scheme can take advantage of this method for enforcing solid boundary conditions. For best effect, the present Cartesian boundary scheme has been combined with a sophisticated, local mesh refinement scheme, and a number of examples are shown in order to demonstrate the efficacy of the combined algorithm for simulations of shock interaction phenomena.

Quirk, James J.↗
NASA NTRS1993

A parallel adaptive mesh refinement algorithm

Over recent years, Adaptive Mesh Refinement (AMR) algorithms which dynamically match the local resolution of the computational grid to the numerical solution being sought have emerged as powerful tools for solving problems that contain disparate length and time scales. In particular, several workers have demonstrated the effectiveness of employing an adaptive, block-structured hierarchical grid system for simulations of complex shock wave phenomena. Unfortunately, from the parallel algorithm developer's viewpoint, this class of scheme is quite involved; these schemes cannot be distilled down to a small kernel upon which various parallelizing strategies may be tested. However, because of their block-structured nature such schemes are inherently parallel, so all is not lost. In this paper we describe the method by which Quirk's AMR algorithm has been parallelized. This method is built upon just a few simple message passing routines and so it may be implemented across a broad class of MIMD machines. Moreover, the method of parallelization is such that the original serial code is left virtually intact, and so we are left with just a single product to support. The importance of this fact should not be underestimated given the size and complexity of the original algorithm.

Quirk, James J.↗
NASA NTRS1992

A contribution to the great Riemann solver debate

The aims of this paper are threefold: to increase the level of awareness within the shock capturing community to the fact that many Godunov-type methods contain subtle flaws that can cause spurious solutions to be computed; to identify one mechanism that might thwart attempts to produce very high resolution simulations; and to proffer a simple strategy for overcoming the specific failings of individual Riemann solvers.

Quirk, James J.↗
NASA NTRS1992

An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies

In this paper we describe an approach for dealing with arbitrary complex, two dimensional geometries, the so-called cartesian boundary method. Conceptually, the cartesian boundary method is quite simple. Solid bodies blank out areas of a background, cartesian mesh, and the resultant cut cells are singled out for special attention. However, there are several obstacles that must be overcome in order to achieve a practical scheme. We present a general strategy that overcomes these obstacles, together with some details of our successful conversion of an adaptive mesh algorithm from a body-fitted code to a cartesian boundary code.

Quirk, James J.↗