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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Reacceleration of Galactic Cosmic Rays beyond the Knee at the Termination Shock of a Cosmic-Ray-driven Galactic Wind

The origin of cosmic rays (CRs) above the knee in the spectrum is an unsolved problem. We present a wind model in which interstellar gas flows along a nonrotating, expanding flux tube with a changing speed and cross-sectional area. CRs from Galactic sources, such as supernova remnants, which are coupled to the plasma via Alfvén waves, provide the main pressure source for driving this outflow. These CRs are then subject to diffusive shock reacceleration at the Galactic wind termination shock, which is located at a distance ~200 kpc. Some of the highest-energy reaccelerated particles propagate upstream against the wind and can contribute to the petaelectronvolt to exaelectronvolt range of the spectrum. We analyze the conditions under which efficient reacceleration can occur and find that rigidities ~10–40 PV can be obtained and that the termination shock may account for between 10% and 50% of the proton spectrum measured in IceCube/IceTop experiment. In our model, the termination shock is unable to fully explain the CR spectrum in the petaelectronvolt to exaelectronvolt range. The highest-energy particles that escape downstream from our termination shock, and similar shocks surrounding most galaxies, can be further accelerated by intergalactic shock fronts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Investigations of transonic trailing edge flows

An experimental study of several of the trailing edge and wake turbulence properties for a NACA 64A010 airfoil section was completed. The experiments were conducted at the Ohio State University Aeronautical and Astronautical Research Laboratory in the 6 inch x 22 inch transonic wind tunnel facility. The data were obtained at a free stream Mach number of 0.80 and a flow Reynolds number (based on chord length) of 5 million. The principal diagnostic tool was a dual-component laser Doppler velocimeter. The experimental data included surface static pressures, chordwise and vertical mean velocities, RMS turbulence intensities, local flow angles, and a determination of turbulence kinetic energy in the wake at chordwise locations between the transonic shock wave and the trailing edge, in the near wake just downstream of the trailing edge and in the far wake. At the two angles of attack tested (0 and 2 degrees), the turbulence intensities and turbulence kinetic energy were observed to decay in the streamwise direction. In the far wake, for the non-lifting case, the turbulence intensities were nearly isotropic. For the two degree case, the horizontal component of the turbulence intensity was observed to be substantially higher than the vertical component.

Petrie, S. L.↗
NASA NTRS1974

Observation of a driver gas-tangential discontinuity

A complete analysis of an interplanetary disturbance of Nov. 19, 1970 using the Apollo 12-SIDE (Suprathermal Ion Detector Experiment) is presented. The SIDE detectors were pointing at 26.3 degrees from the normal solar-wind direction during the observations. The data were least-squares fitted (using a parabolic hypersurface approximation) to a convected Maxwell-Boltzmann distribution function. The results of the fit combined with two other experiments showed a drastic change in the wind speed (from an steady 352 km/sec down to 219 km/sec), direction, and temperature. Except for a delta-function increase at the onset, the density remained constant. There was a considerable enhancement in the abundance of He and probably of heavier elements. The interplanetary magnetic field exhibited a jump of 21 gamma with a change in latitude from -56 to -76 degrees in solar ecliptic coordinates. It is concluded that the disturbance was due to the driver gas-tangential discontinuity of a solar flare-induced shock wave. The characteristic of the tangential discontinuity fit well with theoretical prediction.

Medrano, R. A.↗
NASA NTRS1986

An upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes (PNS) equations. This method does not require the addition of user specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming scheme in terms of accuracy, stability, computer time and storage, and programming effort. The new algorithm has been validated by applying it to three laminar test cases including flat plate boundary-layer flow, hypersonic flow past a 15 deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with the results obtained using the conventional Beam-Warming algorithm.

Lawrence, S. L.↗
NASA NTRS1944

The "limiting line" in mixed subsonic and supersonic flow of compressible fluids

It is well known that the vorticity for any fluid element is constant if the fluid is non-viscous and the change of state of the fluid is isentropic. When a solid body is placed in a uniform stream, the flow far ahead of the body is irrotational. Then if the flow is further assumed to be isentropic, the vorticity will be zero over the whole filed of flow. In other words, the flow is irrotational. For such flow over a solid body, it is shown by Theodorsen that the solid body experiences no resistance. If the fluid has a small viscosity, its effect will be limited in the boundary layer over the solid body and the body will have a drag due to the skin friction. This type of essentially isentropic irrotational flow is generally observed for a streamlined body placed in a uniform stream, if the velocity of the stream is kept below the so-called "critical speed." At the critical speed or rather at a certain value of the ratio of the velocity of the undisturbed flow and the corresponding velocity of sound, shock waves appear. This phenomenon is called the "compressibility bubble." Along a shock wave, the change of state of the fluid is no longer isentropic, although still adiabatic. This results in an increase in entropy of the fluid and generally introduces vorticity in an originally irrotational flow. The increase in entropy of the fluid is, of course, the consequence of changing part of the mechanical energy into heat energy. In other words, the part of fluid affected by the shock wave has a reduced mechanical energy. Therefore, with the appearance of shock waves, the wake of the streamline body is very much widened, and the drag increases drastically. Furthermore, the accompanying change in the pressure distribution over the body changes the aerodynamic moment acting on it and in the case of an airfoil decreases the lift force. All these consequences of the breakdown of isentropic irrotational flow are generally undesirable in applied aerodynamics. Its occurrence should be delayed as much as possible by modifying the shape or contour of the body. However, such endeavor will be very much facilitated if the cause or the criterion for the breakdown can be found first.

MACH NUMBER↗
NASA NTRS1995

Ulysses solar wind plasma observations from peak southerly latitude through perihelion and beyond

The Ulysses spacecraft trajectory includes a peak southern latitude of -80.2 deg, reached during September 1994, and perihelion in the ecliptic plane at 134 AU in March 1995. The near-perihelion mission phase features a rapid scan through solar latitude, with rates approaching one degree per day. We will present observations through mid-May 1995, when the spacecraft will be near 1.5 AU and +50 deg solar latitude. At the time of this writing, observations from the solar wind plasma experiment have been examined through -40 deg solar latitude. At that latitude Ulysses was still immersed in fast solar wind from the south polar coronal hole, with wind speeds of 700 to 800 km/s and with a variety of fine structure. Expectations for near-perihelion measurements include times of slow, dense wind characteristic of the near-equatorial heliomagnetic streamer belt. A non-zero tilt of the streamer belt would produce recurrent intervals of fast coronal hole wind and corotating interaction regions (CIRs) caused by fast wind overtaking slow wind. Forward and reverse shock waves bounding the CIRs, routinely observed at low and intermediate latitudes beyond approximately 2 AU, will likely be absent during the northward transit due to proximity to the Sun we will summarize solar wind plasma results concerning meridional gradients in fluid parameters, large-scale and fine structure, and transient events such as coronal mass ejections.

Phillips, J. L.↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Screech Noise Generation From Supersonic Underexpanded Jets Investigated

Many supersonic military aircraft and some of the modern civilian aircraft (such as the Boeing 777) produce shock-associated noise. This noise is generated from the jet engine plume when the engine nozzle is operated beyond the subsonic operation limit to gain additional thrust. At these underexpanded conditions, a series of shock waves appear in the plume. The turbulent vortices present in the jet interact with the shock waves and produce the additional shock-associated noise. Screech belongs to this noise category, where sound is generated in single or multiple pure tones. The high dynamic load associated with screech can damage the tailplane. One purpose of this study at the NASA Glenn Research Center at Lewis Field was to provide an accurate data base for validating various computational fluid dynamics (CFD) codes. These codes will be used to predict the frequency and amplitude of screech tones. A second purpose was to advance the fundamental physical understanding of how shock-turbulence interactions generate sound. Previously, experiments on shock-turbulence interaction were impossible to perform because no suitable technique was available. As one part of this program, an optical Rayleigh-scattering measurement technique was devised to overcome this difficulty.

Panda, Jayanta↗
DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗
NASA NTRS1985

Density and temperature structure over northern Europe

During the Energy Budget Campaign, a number of profiles of the density and temperature were obtained to study the structure and variability of the atmosphere. The measurements were made using rocketborne instrumentation launched from Esrange, Sweden, and Andoya Rocket Range, Norway, during November and December 1980. The techniques included meteorological temperature sondes, passive falling sphere, accelerometer instrumented falling spheres, density gauges, mass spectrometers and infrared emission experiments. The instruments provided data covering the altitude range from 20 to 150 km. The measurements were made during periods which have been grouped into three categories by level of geomagnetic activity. Analysis has been made to compare the results and to examine the wave features and variations in the vertical profiles for scales ranging between hundreds of meters and tens of kilometers. Most of the features observed fit qualitatively within the range expected for internal gravity waves. However, the features in the profiles during one of the measurement periods are unusual and may be due to aurorally generated shock waves. The geomagnetic storm conditions caused temperature increases in the lower thermosphere which maximized in the 120-140 km region.

Philbrick, C. R.↗
DOE OSTI2023

Non-Richtmyer–Meshkov instability ejecta production based on shallow bubble collapse

The study of shock-driven ejecta production has focused on Richtmyer–Meshkov instability (RMI) growth from geometric features of the material surface. Extensive study of this mechanism under both single- and multiple-shock conditions has found that the ejected mass tends to be closely associated with the shocked surface phase, and its temperature is not dramatically greater than the hydrodynamic shock temperature of the bulk. In this work, we propose and demonstrate a new ejecta production mechanism that can occur under multiple-shock conditions based on the collapse of bubbles near the free surface of the material. This mechanism produces ejected mass that is much greater in quantity than observed in the RMI case. The particles are much hotter than predicted by the shock Hugoniot state, and the ejected mass does not appear to be strongly dependent upon initial surface finish. The ejecta source extends into the material with no clear remaining free surface. We name this mechanism Shallow Bubble Collapse (SBC) and discuss the conditions under which it activates. We demonstrate resolved modeling methods that enable the calculation, design, and study of SBC as a mechanism and perform a series of experiments to compare with the models. Under some multiple-shock conditions, SBC ejection produces ten times more ejected mass than RMI growth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Understanding oblique impacts from experiments, observations, and modeling

Natural impacts in which the projectile strikes the target vertically are virtually nonexistent. Nevertheless, our inherent drive to simplify nature often causes us to suppose most impacts are nearly vertical. Recent theoretical, observational, and experimental work is improving this situation, but even with the current wealth of studies on impact cratering, the effect of impact angle on the final crater is not well understood. Although craters' rims may appear circular down to low impact angles, the distribution of ejecta around the crater is more sensitive to the angle of impact and currently serves as the best guide to obliquity of impacts. Experimental studies established that crater dimensions depend only on the vertical component of the impact velocity. The shock wave generated by the impact weakens with decreasing impact angle. As a result, melting and vaporization depend on impact angle; however, these processes do not seem to depend on the vertical component of the velocity alone. Finally, obliquity influences the fate of the projectile: in particular, the amount and velocity of ricochet are a strong function of impact angle.

NASA Discipline Exobiology↗
NASA NTRS1989

Upwind algorithm for the parabolized Navier-Stokes equations

A new upwind algorithm based on Roe's scheme has been developed to solve the two-dimensional parabolized Navier-Stokes equations. This method does not require the addition of user-specified smoothing terms for the capture of discontinuities such as shock waves. Thus, the method is easy to use and can be applied without modification to a wide variety of supersonic flowfields. The advantages and disadvantages of this adaptation are discussed in relation to those of the conventional Beam-Warming (1978) scheme in terms of accuracy, stability, computer time and storage requirements, and programming effort. The new algorithm has been validated by applying it to three laminar test cases, including flat-plate boundary-layer flow, hypersonic flow past a 15-deg compression corner, and hypersonic flow into a converging inlet. The computed results compare well with experiment and show a dramatic improvement in the resolution of flowfield details when compared with results obtained using the conventional Beam-Warming algorithm.

Lawrence, Scott L.↗
NASA NTRS2004

Modeling Nonlinear Acoustic Standing Waves in Resonators: Theory and Experiments

The overall goal of the cooperative research with NASA Glenn is to fundamentally understand, computationally model, and experimentally validate non-linear acoustic waves in enclosures with the ultimate goal of developing a non-contact acoustic seal. The longer term goal is to transition the Glenn acoustic seal innovation to a prototype sealing device. Lucas and coworkers are credited with pioneering work in Resonant Macrosonic Synthesis (RMS). Several Patents and publications have successfully illustrated the concept of Resonant Macrosonic Synthesis. To utilize this concept in practical application one needs to have an understanding of the details of the phenomenon and a predictive tool that can examine the waveforms produced within resonators of complex shapes. With appropriately shaped resonators one can produce un-shocked waveforms of high amplitude that would result in very high pressures in certain regions. Our goal is to control the waveforms and exploit the high pressures to produce an acoustic seal. Note that shock formation critically limits peak-to-peak pressure amplitudes and also causes excessive energy dissipation. Proper shaping of the resonator is thus critical to the use of this innovation.

Raman, Ganesh↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Scaling arguments for a plasma experiment relevant to an interplanetary coronal mass ejection

We present the design of a plasma experiment scaled to be relevant to an interplanetary coronal mass ejection (ICME) traveling from the Sun to the Earth. The experiment would launch a compact torus of plasma (the scaled analog of the coronal mass ejected by the Sun) into a magnetized background plasma (the scaled analog of the interplanetary medium, the diffuse magnetized plasma that fills the Solar System). We present the complete process of scaling the ICME system to a laboratory experiment: beginning with the ICME features that we aim to replicate with a laboratory experiment, explaining how these desired features are translated into specific design constraints, and presenting a set of experimental parameters that meet those constraints. We find that the experiment would be possible with plasma conditions typically achieved in large-scale facilities such as the Big Red Ball at the University of Wisconsin or the LArge Plasma Device at the University of California, Los Angeles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

On the quasi-conical flowfield structure of the swept shock wave-turbulent boundary layer interaction

The swept oblique shock-wave/turbulent-boundary-layer interaction generated by a 20-deg sharp fin at Mach 4 and Reynolds number 21,000 is investigated via a series of computations using both conical and three-dimensional Reynolds-averaged Navier-Stokes equations with turbulence incorporated through the algebraic turbulent eddy viscosity model of Baldwin-Lomax. Results are compared with known experimental data, and it is concluded that the computed three-dimensional flowfield is quasi-conical (in agreement with the experimental data), the computed three-dimensional and conical surface pressure and surface flow direction are in good agreement with the experiment, and the three-dimensional and conical flows significantly underpredict the peak experimental skin friction. It is pointed out that most of the features of the conical flowfield model in the experiment are observed in the conical computation which also describes the complete conical streamline pattern not included in the model of the experiment.

Knight, Doyle D.↗
NASA NTRS2018

Relative Sputtering Rates of FeS, MgS, and Mg Silicates: Implications for ISM Gas Phase Depletions of Rock-Forming Elements

Astronomical measurements of S abundances in the diffuse interstellar medium (ISM) indicate ionized S is a dominant species with little (< 5%) S residing in grains (e.g. Jenkins 2009). This is an enigmatic result, given that abundant Fe-sulfide grains are observed in dust around pre- and post-main sequence stars and are also observed as major components of primitive meteoritic and cometary samples. These disparate observations suggest that the lifetime of sulfide grains in the ISM is short because of destruction processes. Our previous work has shown that FeS and MgS retain their crystallinity and do not amorphize during radiation processing, whereas enstatite and forsterite are readily amorphized. We have extended this study to measure the relative sputtering rates of FeS and MgS compared to enstatite and forsterite. Irradiation of FeS with 4 keV He+ results in preferential sputtering of S and the formation of a thin 2-3 nm, compact Fe metal layer that armors the surface. The zone of S loss extends to a depth of ~8-10 nm below the exposed surface. Despite this S loss, the FeS retains its crystallinity and shows no sign of incipient amorphization. Irradiation of FeS with 5kV Ga+ in a focused ion beam (FIB) instrument resulted in preferential sputtering of S and the formation of a 5-8-nm thick surface layer of nanophase Fe metal. X-ray mapping shows that the zone of S sputtering extends to a depth of nearly 20 nm, but there is no evidence for FeS amorphization, consistent with our previous work. The irradiation experiments show that the relative sputtering rate of FeS and MgS are much higher than olivine or enstatite. Sputtering experiments utilizing 30 kV and 5 kV Ga ions in the FIB produced volume loss in troilite that was ~4X greater than in enstatite or forsterite. The sputter yield under these conditions is such that for every Si atom sputtered from enstatite, ~14 S atoms are sputtered from FeS. We have performed similar sputtering experiments on Fe-bearing niningerite (MgS) and co-existing enstatite from the ALH 84170 EH3 chondrite. MgS also sputters much more rapidly than enstatite with a relative Si:S sputter yield of 1:8. For MgS, sulfur is highly depleted at the surface and the S-depletion zone extends to a depth of ~15 nm (using 5 kV Ga+). There is a corresponding zone of Mg and especially Fe enrichment that extends from 5 to ~10 nm below the surface, respectively. The dominant grain destruction mechanism in the ISM is sputtering from passage of supernova-generated shock waves. This process also results in the amorphization of crystalline silicates in the ISM. Our results indicate that FeS and MgS grains produced in evolved stars and injected into the ISM will be destroyed more rapidly than crystalline silicates. This process may account for the lack of significant depletion of S from the gas phase in the ISM. However, rare nanophase FeS grains occur as inclusions in circumstellar amorphous silicate grains found in comet dust particles analyzed in the laboratory. These results show that a finite amount of S in the ISM is sequestered in solid grains.

SPUTTERING RATES↗