Contact-surface tailoring in real shock tubes
Contact surface tailoring of reflected shock wave interaction with boundary layer in real shock tubes
SEARCH · Engineering Papers
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.
Contact surface tailoring of reflected shock wave interaction with boundary layer in real shock tubes
A version of the two step Lax-Wendroff difference method with second order accuracy is used to seek solutions of the unsteady magnetohydrodynamic (MHD) equations for the study of an explosion in a non-homogeneous medium with and without magnetic fields. The explosion is initiated by introducing a finite amount of energy in a small volume of gas which leads to an instantaneous increase of gas temperature. Two types of magnetic field configurations (i.e. open and closed) are considered to illustrate the dependence and differences of magnetic and a non-magnetic flow motion. Numerical results show the strong dependence of induced MHD flow field on the magnetic field configuration and strength. The development and decay of the fast and slow MHD shock waves, as well as the ordinary gasdynamic shock waves, are also represented well by the present computer simulations. In conclusion, we have demonstrated that the numerical scheme we have used is a reliable one for studying 2-dimensional MHD problems in non-homogeneous medium.
SR-71 sonic boom signatures were measured to validate sonic boom propagation prediction codes. An SR-71 aircraft generated sonic booms from Mach 1.25 to Mach 1.6, at altitudes of 31,000 to 48,000 ft, and at various gross weights. An F-16XL aircraft measured the SR-71 near-field shock waves from close to the aircraft to more than 8,000 ft below, gathering 105 signatures. A YO-3A aircraft measured the SR-71 sonic booms from 21,000 to 38,000 feet below, recording 17 passes. The sonic booms at ground level and atmospheric data were recorded for each flight. Data analysis is underway. Preliminary results show that shock wave patterns and coalescence vary with SR-71 gross weight, Mach number, and altitude. For example, noncoalesced shock wave signatures were measured by the YO-3A at 21,000 ft below the SR-71 aircraft while at a low gross weight, Mach 1.25, and 31,000-ft altitude. This paper describes the design and execution of the flight research experiment. Instrumentation and flight maneuvers of the SR-71, F-16XL, and YO-3A aircraft and sample sonic boom signatures are included.
A conceptual view for a GAIM energy/driver system to maximize shock-tube performance through efficient interfacing of the energy source with the gas dynamics of the arc driver is presented. Electrical and arc-chamber requirements are evaluated utilizing two new computer codes. One code calculates the shock wave generated for a selected time rate and magnitude of arc-energy input; the other computes the values of external circuit elements required to produce the selected energy input, with the driver represented as the load element of the electrical discharge circuit. Results indicate that the energy-storage capability and the driver arrangement needed to produce the highest shock Mach number can be achieved by means of driver gas addition and by impedance matching (GAIM). Design criteria are presented for arc energy requirements necessary to produce given shock-wave speeds. Shock velocities as high as the 70 km/sec required for simulating Jovian entry now seem possible in shock-tube operation. Practical implementation of a GAIM system is discussed.
Anharmonic models are used to describe the high temperatures and the electronic effects characteristic of shock wave phenomena and planetary interiors. It is found that a generalized Grueneisen function dependent on temperature and volume is necessary for the self-consistent application of the Mie-Grueneisen equation at the physical conditions encountered in shock wave phenomena or in the interior of planets. It is noted that at temperatures of a few thousand degrees Kelvin, by deriving the function in terms of anharmonic lattice dynamics, the difference between gamma and gamma(th) is of the same order of a magnitude as gamma. This leads to the conclusions that (1) gamma is not a function of temperature and (2) gamma?gamma(th) does not explain either the equation of state of planetary interiors or high-temperature shock wave phenomena.
Alfven waves effect on MHD fast shock from numerical solution of extended Rankine-Hugoniot equations
The generation of slow shock waves as the result of the interaction of Alfven discontinuities with solar wind contact surfaces is considered. The latter are taken to be the boundaries of proton and alpha-particle concentration inhomogeneities. It is found that the intensity of the Alfven discontinuity may be increased as the result of its interaction with the more dense plasma. The converse (i.e., decrease of the Alfven discontinuity's intensity following interaction with a less dense plasma) is also indicated. Also discussed is the generation of a magnetic cloud as the result of the interaction of a quasi-parallel Alfven discontinuity with a dense plasma contact surface. It is shown that the (solar-generated) Alfven discontinuity may then be transformed into non-flare fast and slow shock waves as the result of this interaction. Thus, it is indicated that some fast shock waves in the solar wind may have a nonsolar origin.
The paper reports the observation of a large coronal transient that can only be interpreted as a three-dimensional structure. Its form is one which has not been observed before: a gradually expanding, sun-centered disk of excess brightness, whose projected radius increased from 4 to 8 solar radii during 0832-0958 UT on November 27, 1979. This earth-directed transient was the source of an interplanetary shock wave that reached ISEE 3 at 0649 UT, November 30, and earth at 0738 UT, November 30. Fitting the shock speed at ISEE 3 and the average transit speed from the sun to ISEE 3 to a power law of the form V = (V sub 0)(r exp -n), it is found that V sub 0 = 1980 km/s and n = 0.294, in good agreement with shock wave models. The shock speed predicted by the power law at 10 solar radii is 1000 km/s, in good agreement with the estimated frontal speed of the transient.
The relaminarization of the boundary layer over a flat plate in the shock tube was investigated by using the partially reflected shock wave technique. The flow Mach number was approximately 0.14, which corresponds to the inleft flow Mach number for the first row of vanes in a gas turbine. The thin film platinum heat gauges were used to measure the heat transfer rate and the Stanton number was calculated from the oscilloscope voltage traces. The Reynolds number was varied by changing the operation pressure of the shock tube and the values varied from 2.3 x 10 to the 4th to 5.3 x 10 to the 5th. For a Reynolds number range of 7 x 10 to the 4th to 3.5 x 10 to the 5th, the relaminarization of the boundary layer was observed. This phenomenon is due to the decay of the turbulence level in the flow as the reflected shock wave moves upstream from the flat plate. As the Reynolds number increased, the relaminarization was delayed and the delay was related to the turbulence generated by the reflected shock wave.
Inertial confinement fusion (ICF) holds increasing promise as a potential source of abundant, clean energy, but has been impeded by defects such as micro-voids in the ablator layer of the fuel capsules. It is critical to understand how these micro-voids interact with the laser-driven shock waves that compress the fuel pellet. At the Matter in Extreme Conditions (MEC) instrument at the Linac Coherent Light Source (LCLS), we utilized an x-ray pulse train with ns separation, an x-ray microscope, and an ultrafast x-ray imaging (UXI) detector to image shock wave interactions with micro-voids. To minimize the high- and low-frequency variations of the captured images, we incorporated principal component analysis (PCA) and image alignment for flat-field correction. After applying these techniques we generated phase and attenuation maps from a 2D hydrodynamic radiation code (xRAGE), which were used to simulate XPCI images that we qualitatively compare with experimental images, providing a one-to-one comparison for benchmarking material performance. Moreover, we implement a transport-of-intensity (TIE) based method to obtain the average projected mass density (areal density) of our experimental images, yielding insight into how defect-bearing ablator materials alter microstructural feature evolution, material compression, and shock wave propagation on ICF-relevant time scales.
Observations of the plasma waves associated with collisionless shocks are reviewed, and the understanding of their generation mechanisms and their importance to shock physics are summarized. The emphasis is on waves generated directly at the shock, especially ion acoustic and lower-hybrid-like modes. The observations are discussed in the context of shock structure, with attention given to the distinctions between waves generated in the shock foot and ramp. The behavior of resistive, dispersive, and supercritical quasi-perpendicular shocks is contrasted. Evidence for the operation of various generation mechanisms, including interactions with cross-field currents, gyrating reflected ions, and field-aligned electron beams, are summarized. The various forms of plasma heating which are actually observed are outlined, and the role of the various wave modes in this heating is discussed. Conclusion, it is argued that, while plasma wave turbulence may play a vital role in plasma heating for some special shocks, it is of second-order importance in most cases.
This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.
For many internal transonic flows of practical interest, some of the relevant nondimensional parameters typically are small enough that a perturbation scheme can be expected to give a useful level of numerical accuracy. A variety of steady and unsteady transonic channel and cascade flows is studied with the help of systematic perturbation methods which take advantage of this fact. Asymptotic representations are constructed for small changes in channel cross-section area, small flow deflection angles, small differences between the flow velocity and the sound speed, small amplitudes of imposed oscillations, and small reduced frequencies. Inside a channel the flow is nearly one-dimensional except in thin regions immediately downstream of a shock wave, at the channel entrance and exit, and near the channel throat. A study of two-dimensional cascade flow is extended to include a description of three-dimensional compressor-rotor flow which leads to analytical results except in thin edge regions which require numerical solution. For unsteady flow the qualitative nature of the shock-wave motion in a channel depends strongly on the orders of magnitude of the frequency and amplitude of impressed wall oscillations or fluctuations in back pressure. One example of supersonic flow is considered, for a channel with length large compared to its width, including the effect of separation bubbles and the possibility of self-sustained oscillations. The effect of viscosity on a weak shock wave in a channel is discussed.
Apatite is a phosphate mineral relevant to shock metamorphism in planetary materials. Here, we report on the response of natural fluorapatite from Durango, Mexico under shock wave loading between 14.5 and 119.5 GPa. Wave profile measurements were obtained in plate-impact experiments conducted on [0001]-oriented fluorapatite single crystals. To GPa peak stresses, we observed a two-wave structure indicating an elastic-inelastic response with elastic wave amplitudes of 10.5–13.1 GPa. Between 39.1 and 62.1 GPa, a complex wave structure was observed involving the propagation of three waves. At and above 73.7 GPa, only a single shock wave was observed. The data above 73.7 GPa provided the following linear shock velocity – particle velocity relationship: U s = 6.5(2) + 0.78(6) u p , (mm/μs). Above 80 GPa, the densities in the shocked state exceed both the extrapolated 300-K density of fluorapatite and the predicted 300-K density for a mixture of the high-pressure assemblage, tuite and CaF 2 . Furthermore, this result indicates that fluorapatite undergoes a transition to a denser structure under shock loading at these conditions. The shock response of fluorapatite is observed to be similar to enstatite but stiffer than quartz and albite at the stresses examined in this work.
Steady-state gasdynamical studies, previously limited to tightly wound normal spiral galaxies, are extended to models of barred spirals with a 5% to 10% perturbing potential. The models show that a strong wave manifestation is an important constituent of the bar structure in many barred spirals and that a density-wave shock wave can form a bar structure as pronounced as the narrow bars often evident in optical photographs of barred spirals. The dark narrow dust lanes often observed along the leading edges of bar structures are identified as tracers of shocks, and it is found that strong shocks along a bar structure during even a small part of a galaxy's lifetime might easily deplete a large enough proportion of the gas to cause a lack of gas in the inner annuli encompassing the bar by the time of the present epoch. It is emphasized that even moderate-amplitude barlike perturbations in the disk can drive large noncircular gas motions, typically 50 to 150 km/s.
An investigation of the transonic flow over a circular arc airfoil was conducted to obtain basic information for turbulence modeling of shock-induced separated flows and to verify numerical computer codes which are being developed to simulate such flows. The investigation included the employment of a laser velocimeter to obtain data concerning the mean velocity, the shear stress, and the turbulent kinetic energy profiles in the flowfield downstream of the airfoil midchord where the flow features are more complex. Depending on the free-stream Mach number, the flowfield developed was either steady with shock-wave-induced separation extending from the foot of the shock wave to beyond the trailing edge or unsteady with a periodic motion also undergoing shock-induced separation. The experimental data were compared with the results of numerical simulations in which a computer code was employed that solved the time-dependent Reynolds' averaged compressible Navier-Stokes equations.
A numerical investigation of transonic turbulent flows separated by curvature and shock wave - boundary layer interaction is presented. The free stream Mach numbers considered are 0.4, 0.5, 0.6, 0.7, 0.8, 0.825, 0.85, 0.875, 0.90, and 0.925. In the numerical method, the conservation of mass equation is replaced by a pressure correction equation for compressible flows and thus incremental pressure is solved for instead of density. The turbulence is described by a multiple-time-scale turbulence model supplemented with a near-wall turbulence model. The present numerical results show that there exists a reversed flow region at all free stream Mach numbers considered whereas various k-epsilon turbulence models fail to predict such a reversed flow region at low free stream Mach numbers. The numerical results also show that the size of the reversed flow region grows extensively due to the shock wave - turbulent boundary layer interaction as the free stream Mach number is increased. These numerical results show that the turbulence model can resolve the turbulence field subjected to extra strains caused by the curvature and the shock wave - turbulent boundary layer interaction and that the numerical method yields a significantly accurate solution for the complex compressible turbulent flow.
High-entropy alloys are materials with an increasing number of technological applications. Amongst them, the Cantor alloy, FeMnCoCrNi, shows desirable mechanical properties at normal loading conditions. In this study we focus on the performance of the Cantor alloy at the ultra-high deformation rates of shock waves. We study shock-induced spallation using both experiments and atomistic simulations. Experimental loading is achieved using high power laser, with VISAR to obtain velocity profiles and spall strength, followed by transmission electron microscopy of the recovered samples. Molecular Dynamics (MD) simulations of shock-induced spallation are compared with experiments. Both experiments and simulations show a high spall strength which would be beneficial for certain applications, with experiments giving ~8 GPa at ~10 7 s -1 and MD giving almost ~30 GPa at ~10 9 s -1 . The difference between experiments and simulations can be explained by the difference in strain rate. Post-mortem analysis of the experimental samples shows nanotwins near the spall plane, while MD simulations show a highly disordered region giving rise to void nucleation and spall during loading.