Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shock wave experiments”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42

Particle Seeding Method for Small-Scale, High-Pressure Nozzles

NASA tests new launch and reentry vehicle configurations in wind tunnels, where flow visualizations and quantitative flowfield measurements are often desired. Some of these vehicles have rocket motors for propulsion, retro-propulsion, or reaction control. High-pressure air is used to supply these rocket motor plumes. However, it is difficult to make off-body measurements in these regions, for several reasons. First, the plumes themselves lack seeding particles for flow diagnostics, and at the low pressures after expansion, Rayleigh scattering or other molecular techniques yield insufficient signal for flow velocity measurements. Second, the plumes displace particle-seeded tunnel air, preventing measurements in the vicinity of the plume. Third, the plumes force shock waves ahead of the vehicle, which melts the ice crystal fog commonly used for visualization and measurement techniques in certain facilities. In the current work, a novel method for seeding the flow in these small-scale, high-pressure nozzles has been devised and initially demonstrated, potentially enabling quantitative and qualitative measurements with particle-based instruments such as Doppler global velocimetry or particle image velocimetry. The method involves a Venturi contraction to draw the seed liquid out of a reservoir and into the nozzle channel, wherein shearing forces atomize the seed into particles. The concept was tested with a laser sheet visualization, which demonstrated that the flow rate of liquid spray was controllable; a valve could be adjusted to drop the flow rate by up to 65%. This relatively inexpensive and simple technique may prove useful in wind tunnel experiments involving particle-based laser diagnostics and small, high-pressure nozzles.

Aditya S. Acharya↗

Particle Seeding Method for Small-Scale, High-Pressure Nozzles

NASA tests new launch and reentry vehicle configurations in wind tunnels, where flow visualizations and quantitative flowfield measurements are often desired. Some of these vehicles have rocket motors for propulsion, retro-propulsion, or reaction control. High-pressure air is used to supply these rocket motor plumes. However, it is difficult to make off-body measurements in these regions, for several reasons. First, the plumes themselves lack seeding particles for flow diagnostics, and at the low pressures after expansion, Rayleigh scattering or other molecular techniques yield insufficient signal for flow velocity measurements. Second, the plumes displace particle-seeded tunnel air, preventing measurements in the vicinity of the plume. Third, the plumes force shock waves ahead of the vehicle, which melts the ice crystal fog commonly used for visualization and measurement techniques in certain facilities. In the current work, a novel method for seeding the flow in these small-scale, high-pressure nozzles has been devised and initially demonstrated, potentially enabling quantitative and qualitative measurements with particle-based instruments such as Doppler global velocimetry or particle image velocimetry. The method involves a Venturi contraction to draw the seed liquid out of a reservoir and into the nozzle channel, wherein shearing forces atomize the seed into particles. The concept was tested with a laser sheet visualization, which demonstrated that the flow rate of liquid spray was controllable; a valve could be adjusted to drop the flow rate by up to 65%. This relatively inexpensive and simple technique may prove useful in wind tunnel experiments involving particle-based laser diagnostics and small, high-pressure nozzles.

lasers↗

Numerical simulation of shock-induced combustion/detonation

A numerical study is conducted to simulate the shock-induced combustion in premixed H2-air mixtures. Two types of bodies, blunt (spherical projectile) and sharp (wedge), are considered in the study. A nine-species, 18-step finite-rate H2-air chemical reaction mechanism coupled with the Navier-Stokes equations is solved. The flow field over the blunt body is found to be unsteady, when the projectile velocity is same as the Chapman-Jouget velocity of the mixture. The unsteadiness is caused by the periodic instabilities originating in the stagnation zone. Numerical results show good qualitative agreement with the ballistic range shadowgraph. In addition, the frequency of oscillations, determined by using the Fourier power spectrum, is found to be in good agreement with the experiment. The flow field over the wedge is found to be stable for the conditions considered in this study. The oblique detonation wave structure is investigated and the important flow features are discussed.

Singh, D. J.↗

Skin friction measurement in complex flows using thin oil film techniques

The NASA Grant NAG2-261 was initiated to support a program of research to study complex flows that occur in flight and laboratory experiments by building, testing and optimizing an on-board technique for direct measurement of surface shear stress using thin oil film techniques. The program of research has proceeded under the supervision of the NASA Ames Research Center and with further cooperation from the NASA Ames-Dryden and NASA Langley Research Centers. In accordance with the original statement of work, the following research milestones were accomplished: (1) design and testing of an internally mounted one-directional skin friction meter to demonstrate the feasibility of the concept; (2) design and construction of a compact instrument capable of measuring skin friction in two directions; (3) study of transitional and fully turbulent boundary layers over a flat plate with and without longitudinal pressure gradients utilizing the compact two-directional skin friction meter; (4) study of the interaction between a turbulent boundary layer and a shock wave generated by a compression corner using the two-directional meter; and (5) flight qualification of the compact meter and accompanying electronic and pneumatic systems, preliminary installation into flight test fixture.

Source record↗

A Modular Approach to Model Oscillating Control Surfaces Using Navier Stokes Equations

The use of active controls for rotorcraft is becoming more important for modern aerospace configurations. Efforts to reduce the vibrations of helicopter blades with use of active-controls are in progress. Modeling oscillating control surfaces using the linear aerodynamics theory is well established. However, higher-fidelity methods are needed to account for nonlinear effects, such as those that occur in transonic flow. The aeroelastic responses of a wing with an oscillating control surface, computed using the transonic small perturbation (TSP) theory, have been shown to cause important transonic flow effects such as a reversal of control surface effectiveness that occurs as the shock wave crosses the hinge line. In order to account for flow complexities such as blade-vortex interactions of rotor blades higher-fidelity methods based on the Navier-Stokes equations are used. Reference 6 presents a procedure that uses the Navier-Stokes equations with moving-sheared grids and demonstrates up to 8 degrees of control-surface amplitude, using a single grid. Later, this procedure was extended to accommodate larger amplitudes, based on sliding grid zones. The sheared grid method implemented in EulerlNavier-Stokes-based aeroelastic code ENS AERO was successfully applied to active control design by industry. Recently there are several papers that present results for oscillating control surface using Reynolds Averaged Navier-Stokes (RANS) equations. References 9 and 10 report 2-D cases by filling gaps with overset grids. Reference 9 compares integrated forces with the experiment at low oscillating frequencies whereas Ref. 10 reports parametric studies but with no validation. Reference II reports results for a 3D case by modeling the gap region with a deformed grid and compares force results with the experiment only at the mid-span of flap. In Ref. II grid is deformed to match the control surface deflections at the section where the measurements are made. However, there is no indication in Ref. II that the gaps are explicitly modeled as in Ref. 6. Computations using overset grids are reported in Ref. 12 for a case by adding moving control surface to an existing blade but with no validation either with an experiment or another computation.

Navier↗

Use of computer vision for analysis of image datasets from high temperature plasma experiments

Great strides have been made in improving the quality of x-ray radiographs in high energy density plasma experiments, enabled in part by innovations in engineering and manufacturing of integrated circuits and materials. As a consequence, the radiographs of today are filled with a great deal of detail, but few of these features are extracted in a systematic way. Analysis techniques familiar to plasma physicists tend toward brittle 1D lineout or Fourier transform type analyses. The techniques applied to process our data have not kept pace with improvements in the quality of our data. Fortunately, the field of computer vision has a wealth of tools to offer, which have been widely used in industrial imaging and, more recently, adopted in biological imaging. We demonstrate the application of computer vision techniques to the analysis of x-ray radiographs from high energy density plasma experiments, as well as give a brief tutorial on the computer vision techniques themselves. These tools robustly extract 2D contours of shocks, boundaries of inhomogeneities, and secondary flows, thereby allowing for increased automation of analysis, as well as direct and quantitative comparisons with simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗

A novel, magnetically driven convergent Richtmyer–Meshkov platform

Here, we introduce a novel experimental platform for the study of the Richtmyer–Meshkov instability in a cylindrically converging geometry using a magnetically driven cylindrical piston. Magnetically driven solid liner implosions are used to launch a shock into a liquid deuterium working fluid and, ultimately, into an on-axis rod with a pre-imposed perturbation. The shock front trajectory is tracked through the working fluid and up to the point of impacting the rod through the use of on axis photonic Doppler velocimetry. This configuration allows for precise characterization of the shock state as it impacts the perturbed rod interface. Monochromatic x-ray radiography is used to measure the post-shock interface evolution and rod density profile. The ALEGRA MHD model is used to simulate the dynamics of the experiment in one dimension. We show that late in time the perturbation growth becomes non-linear as evidenced by the observation of high-order harmonics, up to n = 5. Two dimensional simulations performed using a combination of the GORGON MHD code and the xRAGE radiation hydrodynamics code suggest that the late time non-linear growth is modified by convergence effects as the bubbles and spikes experience differences in the pressure of the background flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards a Comprehensive Model of Jet Noise Using an Acoustic Analogy and Steady RANS Solutions

An acoustic analogy is developed to predict the noise from jet flows. It contains two source models that independently predict the noise from turbulence and shock wave shear layer interactions. The acoustic analogy is based on the Euler equations and separates the sources from propagation. Propagation effects are taken into account by calculating the vector Green's function of the linearized Euler equations. The sources are modeled following the work of Tam and Auriault, Morris and Boluriaan, and Morris and Miller. A statistical model of the two-point cross-correlation of the velocity fluctuations is used to describe the turbulence. The acoustic analogy attempts to take into account the correct scaling of the sources for a wide range of nozzle pressure and temperature ratios. It does not make assumptions regarding fine- or large-scale turbulent noise sources, self- or shear-noise, or convective amplification. The acoustic analogy is partially informed by three-dimensional steady Reynolds-Averaged Navier-Stokes solutions that include the nozzle geometry. The predictions are compared with experiments of jets operating subsonically through supersonically and at unheated and heated temperatures. Predictions generally capture the scaling of both mixing noise and BBSAN for the conditions examined, but some discrepancies remain that are due to the accuracy of the steady RANS turbulence model closure, the equivalent sources, and the use of a simplified vector Green's function solver of the linearized Euler equations.

Miller, Steven A. E.↗

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Melting anisotropy in crystalline solids

Despite the long and extensive history of melting studies, experiments to determine the dependence of melting on crystal orientation are lacking. Using longitudinal sound speed measurements in aluminum single crystals shock compressed along < 100 > and < 110 > to 168 GPa, we address this need and show that the melting transition (both onset and completion) is strongly anisotropic. Furthermore, the present findings demonstrate the need to consider the role of deformation induced lattice defects on the melting transition and provide a general approach to examine melting anisotropy in crystalline solids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ti-6Al-4V to over 1.2 TPa: Shock Hugoniot experiments, ab initio calculations, and a broad-range multiphase equation of state

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia’s Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia’s STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 10 5 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experiments with a self-correcting wind tunnel

The feasibility of controlling the flow actively through the walls of a transonic, porous wall wind tunnel in order to minimize wall interference effects on a test model is demonstrated. The method is based upon measuring the components of the disturbance velocity at discrete points along an imaginary surface in the flow field within the tunnel. A mathematical formulation of the flow field exterior to the surface including the boundary condition for unconfined flow, i.e., that all disturbance vanish at infinite, is used to determine if these measured velocity components are consistent with that boundary condition. If they are not, the theory provides a better approximation to the velocity component for unconfined flow, and the flow through the tunnel walls is readjusted iteratively until the measured quantities are consistent with unconfined flow. A brief review of theoretical methods is followed by a description of the Calspan self correcting wind tunnel design and operation, calibration with and without active wall control. Typical results obtained by approximating a conventional porous wall wind tunnel for an 0012 airfoil show that active wall control largely reproduces the correct shock wave position, eliminates wall interference of lift and drag, and reduces the interference effects on pitching moment to 10%.

Vidal, R. J.↗

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The hypersonic shock wave-turbulent boundary layer interaction generated by a sharp fin at Mach 8.2

A combined experimental and numerical study has been conducted on the hypersonic shock-wave turbulent-boundary layer interaction at Mach 8.2 generated by a single fin of angles alpha = 10 and 15 deg. Three models are considered: (1) the 3D compressible Reynolds-averaged Navier-Stokes (RANS) equations using the k-epsilon turbulence model, (2) the 3D RANS using the Rodi turbulence model, and (3) the conical RANS using the Baldwin-Lomax algebraic turbulence model. The computations are compared with various experimental data. The computations using models (1) and (2) show quantitatively very similar results and very good agreement with experimental data for surface pressure and skin friction. Comparison with boundary layer profiles of pitot pressure and yaw angle are also generally good, but the peak surface heat transfer is overestimated by up to 48 percent. The effect of the laminar boundary layer on the fin is restricted to the immediate vicinity of the fin surface. Conical calculations using model (3) show substantially poorer agreement with experiment.

Knight, Doyle D.↗

On the importance of three-dimensional modeling for high-energy-density physics experiments

Laser-driven cylindrical implosion experiments enable direct measurements of hydrodynamic instability growth in convergent geometries, providing a wealth of validation data in the high-energy-density regime. These experiments are designed to be nearly axially invariant, allowing for modeling with complementary two-dimensional slices of the cylinder. Two distinct hydrodynamics codes are employed to model a subset of these experiments, and the results are shown to be in very good agreement with each other and the available experimental data. While this 2D modeling approach adequately captures most of the physics of the implosion and ensuing instability growth, there are crucial aspects from the three-dimensional nature of the experiments that are missed in 2D. The first fully 3D simulations of these experiments are presented, and small but significant differences are found to arise from both the axial and azimuthal non-uniformity in the laser drive. Recent experimental results confirming the drive asymmetry are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dependence of 50-keV upstream ion events at IMP 7 and 8 upon magnetic field bow shock geometry

Results are presented for a statistical study of four years (1972-1976) of IMP 7 and 8 observations at approximately 40 R(E) of 50-200 keV upstream ion events measured with the Energetic Particle Experiment. A monotonic increase in the probability of observing upstream particle events is found with a decrease in the angle (theta-Bn) between the interplanetary magnetic field (IMF) and the local shock normal at the point where the IMF intersects the bow shock, independent of the length of time of bow shock connections exceeding approximately 10 min. Approximately equal probabilities are found for observing an event above a given flux from any portion of the bow shock with the same value of theta-Bn, a growth time of the 50-200 keV events of approximately 10 min, a maximum attainable flux of approximately 25,000/sq cm s sr, and a positive correlation between the probability of exceeding a given flux and the 3 hr Kp index. It is concluded that the local structure of the bow shock in the immediate vicinity of the field line connection is the dominant influence in the generation process of energetic upstream particle events and that wave-particle interactions produce a self-throttling mechanism that limits the maximum flux of ions escaping the upstream foreshock.

Mitchell, D. G.↗

Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501

We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. Overall, we conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.

42 ENGINEERING↗