Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock wave”

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 1,639 records · Page 91

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Castner, Raymond S.↗

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Caster, Raymond S.↗

Mixing and combustion of an underexpanded H2 jet in supersonic flow

Theoretical calculations are compared with experimental data in an investigation of the injection, mixing, and combustion of an underexpanded hydrogen jet in a supersonic test stream. The experimental tests in reacting and nonreacting flows were conducted by use of air and nitrogen as test media. Tests were conducted in a free jet and a ducted mode. Theoretical calculations were performed with a modified viscous characteristics analysis. The computer program employs finite rate chemistry and a turbulent viscosity model which varies radially and axially; it handles exit shock waves, expansion fans, and thermal compressions generated in the reacting flow fields. Comparison of the theory and experimental data indicates good agreement for the near flow field; this agreement is directly related to the explicit handling of the near-field shock and expansion waves. Less satisfactory agreement at downstream locations for both free jet and ducted flows is attributed to the lack of continued computation of the shocks and expansions in the far field.

Jenkins, R. V.↗

Magnetospheric Multiscale Observations of Earth’s Oblique Bow Shock Reformation by Foreshock Ultra-Low Frequency Waves

Collisionless shocks can be nonstationary with periodic reformation shown in many simulation results, but direct observations are still tenuous and difficult to conclusively interpret. In this study, using Magnetospheric Multiscale (MMS) observations, we report direct observational evidence of Earth’s oblique bow shock reformation driven by the foreshock ultra-low frequency (ULF)waves. When the four MMS spacecraft were in a string-of-pearls formation roughly along the bow shock normal, they observed that when each period of foreshock ULF waves encountered the bow shock, a new shock ramp formed. Meanwhile, in the magnetosheath, the old bow shock’s remnants were observed periodically convecting downstream. We propose that the reformation mechanism of the oblique bow shock is the variation of the upstream conditions by the periodic ULF waves as they encounter the bow shock. We also examine the nature of reflected ions during the reformation process.

Terry Z Liu↗

Effect of stratification and geometrical spreading on sonic boom rise time

The purpose of our investigation is to determine the effect of unsteadiness (not associated with turbulence) on rise time. The unsteadiness considered here is due to (1) geometrical spreading, (2) stratification, which includes variation in density, temperature, and relative humidity, and (3) N shaped waveform. A very general Burgers equation, which includes all these effects, is the propagation model for our study. The equation is solved by a new computational algorithm in which all the calculations are done in the time domain. The present paper is a progress report in which some of the factors contributing to unsteadiness are studied, namely geometrical spreading and variation in relative humidity. The work of Pierce and Kang, which motivated our study, is first reviewed. We proceed with a discussion of the Burgers equation model and the algorithm for solving the equation. Some comparison tests to establish the validity of the algorithm are presented. The algorithm is then used to determine the distance required for a steady-state shock, on encountering an abrupt change in relative humidity, to reach a new steady state based on the new humidity. It is found that the transition distance for plane shocks of amplitude 70 Pa is about 4 km when the change in relative humidity is 10 percent. Shocks of amplitude 140 Pa require less distance. The effect of spherical and cylindrical spreading is also considered. We demonstrate that a spreading shock wave never reaches steady state and that its rise time will be less than the equivalent steady state shock. Finally we show that an N wave has a slightly shorter rise time than a step shock of the same amplitude.

Cleveland, Robin O.↗

Particle Acceleration at the Sun and in the Heliosphere

Energetic particles are accelerated in rich profusion at sites throughout the heliosphere. They come from solar flares in the low corona, from shock waves driven outward by coronal mass ejections (CMEs), from planetary magnetospheres and bow shocks. They come from corotating interaction regions (CIRs) produced by high-speed streams in the solar wind, and from the heliospheric termination shock at the outer edge of the heliospheric cavity. We sample all these populations near Earth, but can distinguish them readily by their element and isotope abundances, ionization states, energy spectra, angular distributions and time behavior. Remote spacecraft have probed the spatial distributions of the particles and examined new sources in situ. Most acceleration sources can be "seen" only by direct observation of the particles; few photons are produced at these sites. Wave-particle interactions are an essential feature in acceleration sources and, for shock acceleration, new evidence of energetic-proton-generated waves has come from abundance variations and from local cross-field scattering. Element abundances often tell us the physics the source plasma itself, prior to acceleration. By comparing different populations, we learn more about the sources, and about the physics of acceleration and transport, than we can possibly learn from one source alone.

Reames, Donald V.↗

Numerical computation of multishocked, three-dimensional supersonic flow fields with real gas effects.

A computational procedure is presented which is capable of determining the supersonic flow field surrounding three-dimensional wing-body configurations such as a delta-wing space shuttle. The governing equations in conservation-law form are solved by a finite difference method using a second-order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically, and the intersection of these shocks with the bow shock posed no difficulty. Resulting flow fields about typical blunt nose shuttle-like configurations at angle of attack are presented. The differences between perfect and real gas effects for high Mach number flows are shown.

Kutler, P.↗

Experimental study of the shock-induced unsteady laminar boundary layer on a flat plate

The behavior of the unsteady laminar boundary layer induced by the incident shock wave passing over a flat plate mounted in a shock tube has been experimentally studied for shock speeds ranging from 2.35 to 7.34 km/sec by measuring unsteady heat transfer rates to the plate using thin-film heat-flux gages. Theoretical heat-transfer rates were predicted from analytical solutions in the literature which describe the unsteady flat plate boundary layer development for equilibrium real-gas flows. Experimental results obtained for both air and nitrogen were found to be in good agreement with the theoretical predictions.

Cook, W. J.↗

Turbulent viscous shock layer solutions for Jovian entry at small angles of attack

The equations governing the laminar and turbulent flows of reacting gas mixtures in chemical equilibrium over axially symmetric blunt bodies at small angles of attack are developed and presented in the unsteady conservative form. Solutions are obtained in the planes of symmetry of the flow field for the typical conditions encountered by a probe entering the Jovian atmosphere at small angles of attack. The eddy-viscosity is approximated by a two-layer model. The shock wave is treated as a discontinuity across which the shock relations are used to compute the flow conditions behind the shock. A time-asymptotic finite-difference method is used to solve the equations. The zero angle of attack results for a hyperboloid are compared with the existing results but neither experimental nor computational results are available for comparison with the present results at angle of attack.

Jumar, A.↗

Coronal element abundances derived from solar energetic particles

The large gradual solar-energetic-particle (SEP) events, where abundances are commonly measured, are produced when coronal mass ejections (CMEs) drive shock waves through the corona and the interplanetary medium. The shock accelerates particles from the highly-ionized, approximately 1.5 MK, plasma in a manner that depends only weakly upon the Q/A of the ion, except at very high energies. Averaging the approximately 1 MeV/amu abundances over many events compensates for the acceleration effects to produce abundances that appear to correspond directly to those in the coronal source for all observed elements, including H. The resulting abundances reflect the 4 x enhancement of ions with low values of first ionization potential (FIP) arising from ion-neutral fractionation that occurs as the atoms are transported up from the photosphere. A different pattern of fractionation is found for ions that are shock-accelerated from the high speed solar wind emerging from coronal holes.

Reames, Donald V.↗

Vortex distortion during vortex-surface interaction in a Mach 3 stream

An experimental study has been conducted in Mach 3 wind tunnel to investigate the behavior of supersonic vortices as they interact with a wedge surface placed in their passage. The experimental setup was arranged so that interactions resulted in a close encounter of the vortex core and the wedge leading edge. Spark shadow photographs of the flowfield along the pressure measurements on the wedge surface were used to study the interaction problem. In their most organized form, distortion of streamwise vortices upon interacting with the wedge was found to result in formation of symmetric detached shock fronts far upstream of the wedge leading edge followed by an apparent slip surface separating a subsonic region from a supersonic zone. Interaction experiments leading to substantial changes in the structure of vortices revealed that the supersonic vortex distortion has strong resemblances to the incompressible 'B-breakdown' reported in the literature. Experimental results also indicate that the interaction strongly depends on the vortex strength and vortex proximity to the wedge leading edge, and the generated flowfield was found to be highly unsteady. Interaction of concentrated streamwise vortices with the oblique shock formed over the wedge surface resulted in formation of a locally three-dimensional shock wave with a limited subsonic region downstream of the shock.

Kalkhoran, Iraj M.↗

Cones in Supersonic Flow

In the case of cones in axially symmetric flow of supersonic velocity, adiabatic compression takes place between shock wave and surface of the cone. Interpolation curves betwen shock polars and the surface are therefore necessary for the complete understanding of this type of flow. They are given in the present report by graphical-numerical integration of the differential equation for all cone angles and airspeeds.

Hantzsche, W.↗

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗

Standing waves at low Mach number laminar bow shocks

Explorer 43 data were used to study 34 bow shock crossings observed from 5 to 16 earth radii upstream of the average bow shock location. Waves with periods of 6 to 130 s having amplitudes up to delta-B/B = 1 were detected. Wave polarization for the low-frequency waves is right-handed in relation to the average field direction when the observer moves from the upstream to downstream direction but is left-handed when the observer moves in the opposite sense. This fact identified the waves as standing whistler waves in the coordinate system of the shock. The waves are in agreement with collisionless low Mach number laminar shock theory. When the measured parameters were used to calculate theoretical wavelengths, the observed wave frequencies could be used to calculate velocities for the shock-wave coordinate system past the spacecraft; such velocities are mostly between 10 and 30 km/s. It is suggested that the higher-frequency propagating whistler waves may evolve from the standing whistler waves through a decay instability.

Fairfield, D. H.↗

Flight Tests of a Supersonic Natural Laminar Flow Airfoil

IR thermography was used to characterize the transition front on a S-NLF test article at chord Reynolds numbers in excess of 30 million Changes in transition due to Mach number, Reynolds number, and surface roughness were investigated - Regions of laminar flow in excess of 80% chord at chord Reynolds numbers greater than 14 million IR thermography clearly showed the transition front and other flow features such as shock waves impinging upon the surface A series of parallel oblique shocks, of yet unknown origin, were found to cause premature transition at higher Reynolds numbers. NASA has a current goal to eliminate barriers to the development of practical supersonic transport aircraft Drag reduction through the use of supersonic natural laminar flow (S-NLF) is currently being explored as a means of increasing aerodynamic efficiency - Tradeoffs work best for business jet class at M<2 Conventional high-speed designs minimize inviscid drag at the expense of viscous drag - Existence of strong spanwise pressure gradient leads to crossflow (CF) while adverse chordwise pressure gradients amplifies and Tollmien-Schlichting (TS) instabilities Aerion Corporation has patented a S-NLF wing design (US Patent No. 5322242) - Low sweep to control CF - dp/dx < 0 on both wing surfaces to stabilize TS - Thin wing with sharp leading edge to minimize wave drag increase due to reduction in sweep NASA and Aerion have partnered to study S-NLF since 1999 Series of S-NLF experiments flown on the NASA F-15B research test bed airplane Infrared (IR) thermography used to characterize transition - Non-intrusive, global, good spatial resolution - Captures significant flow features well

laminar flow↗

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite Element Simulation of the Acoustic Pressure Inside a Beverage Container for Non-Thermal, Ultrasound-based Pasteurization

The purpose of this effort is to investigate whether large acoustic pressure waves can be transmitted inside beverage containers to enable pasteurization. Acoustic waves are known to induce large nonlinear compressive forces and shock waves in fluids, suggesting that compression waves may be capable of damaging bacteria inside beverage containers without appreciably increasingly the temperature or altering the freshness and flavor of the beverage contents. Although a combined process such as thermosonication (e.g., sonication with heating) is likely more efficient, it is instructive to compute the acoustic pressure field distribution inside the beverage container. The COMSOL simulations used two and three-dimensional models of beverage containers placed in a water bath to compute the acoustic pressure field. A limitation of these COMSOL models is that they cannot determine the bacterial lysis efficiency, rather the models provide an indirect metric of bacterial lysis based on the magnitude of the pressure field and its distribution.

42 ENGINEERING↗

Mach stem structure in exothermic systems.

Incident waves having grossly similar characteristics were obtained by passing a shock wave over a ramp mounted on the divergent part of a convergent-divergent nozzle. The range of triple point parameters in both the naturally occurring detonative Mach stems and those generated in the ramp experiments compared favorably, but some differences were observed in the triple point trajectories and the decay rates of the lead waves. Whereas the lead waves in detonations have been found to decay at rates between those of planar and cylindrical blast waves, the decay observed in the ramp experiments was more like that of planar blast waves. The dominant factor in determining a triple point trajectory was found to be the orientation of the lead wave at the triple point. Comparison of reactive and nonreactive cases indicated a faster rate of decay of the lead wave and a steeper triple point trajectory for the reactive cases than for the nonreactive ones. Also, a sensitivity of the trajectory to initial pressure was observed in the former, although none was present in the nonreactive cases.

Liaugminas, R.↗