Planar reflection of gaseous detonations

Pipes containing flammable gaseous mixtures may be subjected to internal detonation. When the detonation normally impinges on a closed end, a reflected shock wave is created to bring the flow back to rest. This study built on the work of Karnesky (2010) and examined deformation of thin-walled stainless steel tubes subjected to internal reflected gaseous detonations. A ripple pattern was observed in the tube wall for certain fill pressures, and a criterion was developed that predicted when the ripple pattern would form. A two-dimensional finite element analysis was performed using Johnson-Cook material properties; the pressure loading created by reflected gaseous detonations was accounted for with a previously developed pressure model. The residual plastic strain between experiments and computations was in good agreement.

During the examination of detonation-driven deformation, discrepancies were discovered in our understanding of reflected gaseous detonation behavior. Previous models did not accurately describe the nature of the reflected shock wave, which motivated further experiments in a detonation tube with optical access. Pressure sensors and schlieren images were used to examine reflected shock behavior, and it was determined that the discrepancies were related to the reaction zone thickness extant behind the detonation front. During these experiments reflected shock bifurcation did not appear to occur, but the unfocused visualization system made certainty impossible. This prompted construction of a focused schlieren system that investigated possible shock wave-boundary layer interaction, and heat-flux gauges analyzed the boundary layer behind the detonation front. Using these data with an analytical boundary layer solution, it was determined that the strong thermal boundary layer present behind the detonation front inhibits the development of reflected shock wave bifurcation.

I. Impact spallation experiments: fracture patterns and spall velocities. II. Craters in carbonate rocks: an electron paramagnetic resonance analysis of shock damage

This work is divided into two independent papers.

PAPER 1.

Spall velocities were measured for nine experimental impacts into San Marcos gabbro targets. Impact velocities ranged from 1 to 6.5 km/sec. Projectiles were iron, aluminum, lead, and basalt of varying sizes. The projectile masses ranged from a 4 g lead bullet to a 0.04 g aluminum sphere. The velocities of fragments were measured from high-speed films taken of the events. The maximum spall velocity observed was 30 m/sec, or 0.56 percent of the 5.4 km/sec impact velocity. The measured velocities were compared to the spall velocities predicted by the spallation model of Melosh (1984). The compatibility between the spallation model for large planetary impacts and the results of these small scale experiments are considered in detail.

The targets were also bisected to observe the pattern of internal fractures. A series of fractures were observed, whose location coincided with the boundary between rock subjected to the peak shock compression and a theoretical "near surface zone" predicted by the spallation model. Thus, between this boundary and the free surface, the target material should receive reduced levels of compressive stress as compared to the more highly shocked region below.

PAPER 2.

Carbonate samples from the nuclear explosion crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron para- magnetic resonance, EPR. The first series of samples for OAK Crater were obtained from six boreholes within the crater, and the second series were ejecta samples recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to spectra of Solenhofen limestone, which had been shocked to known pressures.

The results of the OAK borehole analysis have identified a thin zone of highly shocked carbonate material underneath the crater floor. This zone has a maximum depth of approximately 200 ft below sea floor at the ground zero borehole and decreases in depth towards the crater rim. A layer of highly shocked material is also found on the surface in the vicinity of the reference bolehole, located outside the crater. This material could represent a fallout layer. The ejecta samples have experienced a range of shock pressures.

It was also demonstrated that the EPR technique is feasible for the study of terrestrial impact craters formed in carbonate bedrock. The results for the Meteor Crater analysis suggest a slight degree of shock damage present in the β member of the Kaibab Formation exposed in the crater walls.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

Anomalous electric dipole internal conversion

Electric dipole internal conversion has been experimentally studied for several nuclei in the rare earth region. Anomalies in the conversion process have been interpreted in terms of nuclear structure effects. It was found that all the experimental results could be interpreted in terms of the j ∙ r type of penetration matrix element; the j ∙ ∇ type of penetration matrix element was not important. The ratio λ of the El j ∙ r penetration matrix element to the El gamma-ray matrix element was determined from the experiments to be:

Lu¹⁷⁵,396 keV, λ = - 1000 ± 100;

282 keV, λ = 500 ± 100;

144 keV, λ = 500 ± 250;

Hf¹⁷⁷, 321 keV λ = - 1400 ± 200;

208 keV λ = - 90 ± 40;

72 keV |λ| ≤ 650;

Gd¹⁵⁵, 86 keV λ = - 150 ± 100;

Tm¹⁶⁹, 63 keV λ = - 100 ± 100;

W¹⁸², 152 keV, λ = - 160 ±80;

67 keV, λ = - 100 ± 100.

Predictions for λ are made using the unified nuclear model.

Nuclear structure effects in internal conversion

Experimental studies of nuclear effects in internal conversion in Ta¹⁸¹ and Lu¹⁷⁵ have been performed. Nuclear structure effects (“penetration” effects), in internal conversion are described in general. Calculation of theoretical conversion coefficients are outlined. Comparisons with the theoretical conversion coefficient tables of Rose and Sliv and Band are made. Discrepancies between our results and those of Rose and Sliv are noted. The theoretical conversion coefficients of Sliv and Band are in substantially better agreement with our results than are those of Rose. The ratio of the M1 penetration matrix element to the M1 gamma-ray matrix element, called λ, is equal to + 175 ± 25 for the 482 keV transition in Ta¹⁸¹ . The results for the 343 keV transition in Lu¹⁷⁵ indicate that λ may be as large as – 8 ± 5. These transitions are discussed in terms of the unified collective model. Precision L subshell measurements in Tm¹⁶⁹ (130keV), W¹⁸² (100 keV), and Ta¹⁸¹ (133 keV) show definite systematic deviations from the theoretical conversion coefficients. The possibility of explaining these deviations by penetration effects is investigated and is shown to be excluded. Other explanations of these anomalies are discussed.

I. Nuclear spin-internal-rotation-coupling. II. Fluorine spin-rotation interaction and magnetic shielding in fluorobenzene

Part I.

The interaction of a nuclear magnetic moment situated on an internal top with the magnetic fields produced by the internal as well as overall molecular rotation has been derived following the method of Van Vleck for the spin-rotation interaction in rigid molecules. It is shown that the Hamiltonian for this problem may be written

H_SR = Ῑ · M · Ĵ + Ῑ · M” · Ĵ”

Where the first term is the ordinary spin-rotation interaction and the second term arises from the spin-internal-rotation coupling.

The F¹⁹ nuclear spin-lattice relaxation time (T₁) of benzotrifluoride and several chemically substituted benzotrifluorides, have been measured both neat and in solution, at room temperature by pulsed nuclear magnetic resonance. From these experimental results it is concluded that in benzotrifluoride the internal rotation is crucial to the spin relaxation of the fluorines and that the dominant relaxation mechanism is the fluctuating spin-internal-rotation interaction.

Part II.

The radiofrequency spectrum corresponding to the reorientation of the F¹⁹ nuclear moment in flurobenzene has been studied by the molecular beam magnetic resonance method. A molecular beam apparatus with an electron bombardment detector was used in the experiments. The F¹⁹ resonance is a composite spectrum with contributions from many rotational states and is not resolved. A detailed analysis of the resonance line shape and width by the method of moments led to the following diagonal components of the fluorine spin-rotational tensor in the principal inertial axis system of the molecule:

F/Caa = -1.0 ± 0.5 kHz

F/Cbb = -2.7 ± 0.2 kHz

F/Ccc = -1.9 ± 0.1 kHz

From these interaction constants, the paramagnetic contribution to the F¹⁹ nuclear shielding in C₆H₅F was determined to be -284 ± ppm. It was further concluded that the F¹⁹ nucleus in this molecule is more shielded when the applied magnetic field is directed along the C-F bond axis. The anisotropy of the magnetic shielding tensor, σ_” - σ_⊥, is +160 ± 30 ppm.