Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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,009 records · Page 56

CFD application to supersonic/hypersonic inlet airframe integration

Supersonic external compression inlets are introduced, and the computational fluid dynamics (CFD) codes and tests needed to study flow associated with these inlets are outlined. Normal shock wave turbulent boundary layer interaction is discussed. Boundary layer control is considered. Glancing sidewall shock interaction is treated. The CFD validation of hypersonic inlet configurations is explained. Scramjet inlet modules are shown.

Benson, Thomas J.↗

Aerodynamic tailoring of the Learjet Model 60 wing

The wing of the Learjet Model 60 was tailored for improved aerodynamic characteristics using the TRANAIR transonic full-potential computational fluid dynamics (CFD) code. A root leading edge glove and wing tip fairing were shaped to reduce shock strength, improve cruise drag and extend the buffet limit. The aerodynamic design was validated by wind tunnel test and flight test data.

Chandrasekharan, Reuben M.↗

Overview of Results from the Cassini Magnetospheric Imaging Instrument (MIMI) During the First Year of Operations

The MIMI investigation comprises three sensors covering the indicated energy ranges: the Ion and Neutral Camera (INCA) -- 7 keV/nuc 3 Mev) and protons (1.6 < E < 160 Mev) from the back end of the dual field-of-view telescope. The Saturn observation sequences began in January, 2004 and culminated in Saturn Orbit Insertion on July 1, 2004. The MIMI sensors observed substantial activity in interplanetary space for several months prior to SOI, including several interplanetary shocks associated with corotating interaction regions, numerous increases most likely originating from particle streams in the vicinity of the Saturnian bow shock and energetic neutral atoms (ENA) emanating from Saturn s magnetosphere. Results following SOI revealed: a dynamical magnetosphere with a day-night asymmetry and an 11-hour periodicity; several water-product ions (O+, OH+, H2O+), but little N+; inferred quantities of neutral gas sufficient to cause major losses in the trapped ions and electrons in the middle and inner magnetosphere; a Titan exosphere that is a copious source of ENA; INCA imaging through ENA has also revealed a previously unknown radiation belt residing inward of the D-ring that is most likely the result of double charge-exchange between the main radiation belt and the upper layers of Saturn s exosphere. Finally, there is ample evidence for the presence of substorm-like injections of plasma that subsequently corotates for a number of days before dissipating on the night-side magnetotail. The observations will be presented and discussed in the context of current theoretical models.

Krimigis, S. M.↗

Application of Micro-ramp Flow Control Devices to an Oblique Shock Interaction

Tests are planned in the 15cm x 15cm supersonic wind tunnel at NASA Glenn to demonstrate the applicability of micro-ramp flow control to the management of shock wave boundary layer interactions. These tests will be used as a database for computational fluid dynamics (CFD) validation and Design of Experiments (DoE) design information. Micro-ramps show potential for mechanically simple and fail-safe boundary layer control.

Hirt, Stefanie↗

When a Standard Candle Flickers

The Crab is the only bright steady source in the X-ray sky. The Crab consists of a pulsar wind nebula, a synchrotron nebula, and a cloud of expanding ejecta. On small scales, the Crab is extremely complex and turbulent. X-ray astronomers have often used the Crab as a standard candle to calibrate instruments, assuming its spectrum and overall flux remains constant over time. Four instruments (Fermi/GBM, RXTE/PCA, Swift/BAT, INTEGRAL/ISGRI) show a approx.5% (50 m Crab) decline in the Crab from 2008-2010. This decline appears to be larger with increasing energy and is not present in the pulsed flux, implying changes in the shock acceleration, electron population or magnetic field in the nebula. The Crab is known to be dynamic on small scales, so it is not too surprising that its total flux varies as well. Caution should be taken when using the Crab for in-orbit calibrations.

Wilson-Hodge, Colleen A.↗

Type II Radio Bursts and Their Association With Coronal Mass Ejections in Solar Cycles 23 and 24

Context. Meter-wavelength type II solar radio bursts are thought to be the signatures of shock-accelerated electrons in the corona. Studying these bursts can give information about the initial kinematics, dynamics, and energetics of coronal mass ejections (CMEs) in the absence of white-light observations. Aims. We investigate the occurrence of type II bursts in solar cycles 23 and 24 and their association with CMEs. We also explore whether type II bursts might occur in the absence of a CME. Methods. We performed a statistical analysis of type II bursts that occurred between 200 and 25 MHz in solar cycles 23 and 24 and determined the temporal association of these radio bursts with CMEs. We categorized the CMEs based on their linear speed and angular width and studied the distribution of type II bursts with fast (≥500 km s−1 ), slow (<500 km s−1), wide (≥60◦), and narrow (<60◦) CMEs. We explored the dependence of type II bursts occurrence on the phases of the solar cycle. Results. Our analysis shows that during solar cycles 23 and 24, 768 and 435 type II bursts occurred, respectively. Of these, 79% were associated with CMEs in solar cycle 23, and 95% were associated with CMEs in solar cycle 24. However, only 4% and 3% of the total number of CMEs were accompanied by type II bursts in solar cycle 23 and 24, respectively. Most of the type II bursts in both cycles were related to fast and wide CMEs (48%). We also determined the typical drift rate and duration for type II bursts, which is 0.06 MHz s−1 and 9 min. Our results suggest that type II bursts dominate at heights ≈1.7−2.3±0.3 R . A clear majority have an onset height around 1.7 ± 0.3 R assuming the four-fold Newkirk model. Conclusions. The results indicate that most of the type II bursts had a white-light CME counterpart, but a few type II bursts lacked a clear CME association. There were more CMEs in cycle 24 than in cycle 23. However, cycle 24 contained fewer type II radio bursts than cycle 23. The onset heights of type II bursts and their association with wide CMEs reported in this study indicate that the early lateral expansion of CMEs may play a key role in the generation of these radio bursts.

Sun: activity↗

Shock capturing schemes for multidimensional flow

Progress made in the development of a genuinely multidimensional finite volume algorithm for problems in inviscid gas dynamics is presented. The approach entails: (1) the reconstruction of flowfield data using a planar wave pattern in which the strengths and orientations of the component waves are derived independently of the mesh geometry; and (2) the development of a flux formula which provides a numerical approximation to the flux at a finite volume cell face during the passage of waves which are in general oblique to the face. Several algorithms are outlined, and the most recent developments are included. The results of several numerical cases are also included.

Parpia, Ijaz H.↗

The Three-Dimensional Bow Shock of Mars as Observed by MAVEN

The Martian magnetosphere is a product of the interaction of Mars with the interplanetary magnetic field and the supersonic solar wind. The location of the bow shock has been previously modeled as conic sections using data from spacecraft such as Phobos 2, Mars Global Surveyor, and Mars Express. The Mars Atmosphere and Volatile EvolutioN (MAVEN) mission spacecraft arrived in orbit about Mars in November 2014 resulting in thousands of crossings to date. We identify over 1,000 bow shock crossings. We model the bow shock as a three-dimensional surface accommodating asymmetry caused by crustal magnetic fields. By separating MAVEN's bow shock encounters based on solar condition, we also investigate the variability of the surface. We find that the shock surface varies in shape and location in response to changes in the solar radiation, the solar wind Mach number, dynamic pressure of the solar wind, and the relative local time location of the strong crustal magnetic fields (i.e., whether they are on the dayside or on the nightside).

Supersonic↗

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Refractive index of lithium fluoride at high dynamic stresses

Alkali halides are prototypical ionic solids whose refractive index is a fundamental property related to their lattice structure and ionic polarizability. In particular, lithium fluoride (LiF) has the largest band gap of any known transparent material and maintains transparency at high pressures, making it well suited for use as an optical window in dynamic compression experiments. While empirical fits to the measured density dependence of the refractive index have been provided, a model that is valid over the entire pressure range of experiments and based on a polarization-based theoretical description is lacking. We present a refractive index model for dynamically compressed LiF based on the Lorentz-Lorenz equation, where the molecular polarizability is determined using a single oscillator model with a strain polarizability parameter Λ=0.73. Here, we show that our modeling approach provides an excellent match to the LiF refractive index data for both shock and ramp compression experiments to 900 GPa (density compression greater than threefold). Additionally, updated fits are provided to determine the refractive index correction for LiF windows used in dynamic compression experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

Second coefficient of viscosity in air

Acoustic attenuation measurements in air were analyzed in order to estimate the second coefficient of viscosity. Data over a temperature range of 11 C to 50 C and at relative humidities between 6 percent and 91 percent were used. This analysis showed that the second coefficient of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of viscosity over the temperature and humidity range of the data. In addition, the data showed that the molecular relaxation effects, which are responsible for the magnitude of the second coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic equations of state. Compressible flows containing large streamwise velocity gradients, like shock waves, which cause significant changes in particle properties to occur during time intervals shorter than hundredths of seconds, must be modeled using dynamic equations of state. The dynamic model approach is described briefly.

Ash, Robert L.↗

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s -1 convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗

Temporal and Spatial Dynamics of Nitric Oxide Production at High Latitudes During an ICME-driven Storm on Dec. 14, 2006

During geomagnetic storms driven by interplanetary coronal mass ejections (ICMEs), Earth’s upper atmosphere neutral mass density is directly driven by the geomagnetic storm cycle from intense driving and subsequent recovery. The neutral density heating and cooling times are shorter for stronger storms and longer for weaker storms. The influx of energy (Joule Heating and precipitation) into Earth’s upper atmosphere during the storm main phase allows for the enhanced production of nitric oxide (NO) at high latitudes, which in turn cools the thermosphere by radiating away excess energy. As a result, greater NO production results in quicker thermospheric cooling. Yet while the production of NO on a global scale has been examined in detail, the space and time evolution of NO is not well-known. In this study, we investigate the effects of a particular ICME-driven storm, which occurred on Dec. 14th, 2006, on the production of NO at high latitudes. We compare NO measurements from the Thermosphere, Ionosphere, Mesosphere Dynamics (TIMED) spacecraft to neutral mass density measurements from the Challenging Minisatellite Payload (CHAMP) spacecraft and find that the impact of the shock prior to the storm, not just the onset of the storm itself, is responsible for an increase in NO production. We also observe a dawn-dusk asymmetry in high-latitude NO production and discuss how the geometry of the solar wind driver may play a role in creating this asymmetry. Finally, we discuss how the improved hemispheric and local time coverage provided by the GDC and DYNAMIC missions will enable more robust studies of the spatiotemporal evolution of NO during storms. Since NO is an important factor in regulating the neutral density of the thermosphere, such studies will be relevant for more accurate orbital drag calculations.

Kevin Delano↗

Electron Dynamics in a Subproton-Gyroscale Magnetic Hole

Magnetic holes are ubiquitous in space plasmas, occurring in the solar wind, downstream of planetary bow shocks, and inside the magnetosphere. Recently, kinetic-scale magnetic holes have been observed near Earth's central plasma sheet. The Fast Plasma Investigation on NASA's Magnetospheric Multiscale (MMS) mission enables measurement of both ions and electrons with 2 orders of magnitude increased temporal resolution over previous magnetospheric instruments. Here we present data from MMS taken in Earth's nightside plasma sheet and use high-resolution particle and magnetometer data to characterize the structure of a subproton-scale magnetic hole. Electrons with gyroradii above the thermal gyroradius but below the current layer thickness carry a current sufficient to account for a 10-20 depression in magnetic field magnitude. These observations suggest that the size and magnetic depth of kinetic-scale magnetic holes is strongly dependent on the background plasma conditions.

Gershman, Daniel J.↗

The bow shock structure of IRS 7 - Wind-wind collision near the Galactic center

New structural details of ionized gas associated with IRS 7, a mass-losing supergiant near the Galactic center, are reported, and radio continuum and theoretical maps are compared. High-resolution multiconfiguration VLA observations of IRS 7 reveal ionized gas with a bow shock structure at a wavelength of 2 cm. The apex of the bow shock is facing more toward the cluster of blue stellar objects, known as IRS 16, than Sgr A*, the compact nonthermal source near the Galaxy's dynamical center. It is inferred that the shape of the mass-losing envelope of IRS 7 is influenced by Galactic winds from IRS 16. It is concluded from a comparison of the morphology of the ionized envelope of IRS 7 with simple theoretical modeling of the bow shock structure that the shape of the shock results from the collision between the Galactic center and stellar winds rather than from the motion of IRS 7 through the interstellar medium.

Yusef-Zadeh, Farhad↗

A survey of grid generation techniques in computational fluid dynamics

The use of numerically-generated body-conforming coordinate systems has become well-established in computational fluid dynamics, allowing general codes to be written which can treat boundaries of arbitrary shape. There is now a rapidly increasing interest in the dynamic coupling of the grid with the physical solution so that the grid lines continually move to achieve concentration in regions of strong variation of the physical solution and alignment with shocks, flame fronts, etc. Such adaptive coordinate systems have been developed using a variety of techniques for sensing and reacting to the developing regions of concentration and/or alignment. This paper reviews the techniques and applications of these dynamically adaptive systems and other systems.

Thompson, J. F.↗