Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Diabatic dynamic initialization

A version of the dynamical adjustment procedure is generalized in order to make the procedure applicable to a diabatic model. In this new initialization procedure, backward adiabatic model integration is followed by forward diabatic model integration, with a high-frequency (low-pass) time filter in the form of the Euler backward time-differencing scheme applied throughout the integration. The new procedure was applied to the Goddard Laboratory for Atmospheres 4D data-assimilation system. It was found that, right from the beginning of the forecast integration, the forecast tendencies (and fields) were free of any noise due to imbalance in initial conditions, and the shocks related to an initial imbalance between model physics and dynamics, or the initial spinup effect, were practically removed.

Fox-Rabinovitz, Michael S.↗
NASA NTRS1994

Energetic ion distributions on both sides of the Earth's magnetopause

The Active Magnetospheric Particle Tracer Explorer/Charge Composition Explorer (AMPTE/CCE) spacecraft, with an apogee of approximately 8.8R(sub E) and an inclination of approximately 4.3 deg, sampled the outer dayside equatorial magnetosphere for extended time periods and often crossed into the magnetosheath whenever the solar wind pressure was sufficiently high to compress the magnetopause to less than 8.8R(sub E). We have analyzed ion distributions on both sides of the magnetopause in order to investigate any local time and energy dependence, giving information about physical processes at the magnetopause and the bow shock. Particle measurements are from the Charge Energy Mass (CHEM)(1.5 to 300 keV/e) and Medium-Energy Particle Analyzer (MEPA)(10 keV to 2 MeV) instruments. The wide total energy range (1.5 keV to 2 MeV) covered describes the magnetospheric distributions quite well, and for the purpose of this study, adequately describes the high-energy part of the shocked solar wind. Thus both solar wind and magnetospheric components can be recognized in a mixed particle distribution. Case studies of representative magnetopause crossings at dawn, noon, and dusk, as well as a survey of several other crossings, indicated: (1) a local time and energy dependence of magnetosheath spectra at energies greater than or equal to 50 KeV; spectra were harder at the duskside than at the dawnside and also correlated with magnetospheric activity, (2) constantly much higher intensities in the magnetosphere than in the magnetosheath at energies greater than 10 KeV and an earthward gradient in the subsolar magnetosheath. In addition to the steady state magnetosheath population there exists a burst-type component indicative of a magnetospheric source, and most of the time this is recognized as a flux transfer event. Overall, the results about the origin of the greater than or equal to 50 KeV magnetosheath ions are consistent with the continuous leakage of magnetospheric particles across a tangential discontinuity magnetopause, locally distributed according to magnetospheric drift paths. Magnetic reconnection, although present, should not be a dominant source on average, because it is not continuous in time. Fermi acceleration should not be dominant because it predicts the opposite local time asymmetry, and shock drift acceleration should be a minor contributor at E greater than or equal to 50 keV because of upper-energy cutoff limitations. Our observations also indicate a significant magnetospheric contribution to energies as low as approximately 10 KeV, where the magnetosphere-magnetosheath intensity gradient reverses. However, in orderto examine the relative strength and local time distribution of all possible sources at these energies, a detailed analysis is required.

Paschalidis, N. P↗
NASA NTRS2008

Swift Observations of Supernovae during and after Shock Breakout

Over the past few years, space-based observatories have allowed detailed studies of energetic supernova explosions in previously unexplored time domains and wavelength regimes. NASA's Swift observatory is playing an important role in probing the early emission of SNe during and after shock breakout due to its fast response, flexible scheduling capabilities, and large wavelength band coverage, ranging from the optical, W, and X-ray to the Gamma-ray bands. By studying the outgoing SN shocks with material in its surroundings, the explosion physics and nature of progenitor stars can be studied. Furthermore, monitoring the X-ray emission of SNe with space-based X-ray observatories is being used to map the density structure in SN environments out to large radii from the sites of the explosions (>E20 cm), the transition of a SN into an old supernova remnant can be studied, and the mass-loss rates of the progenitor stars are being probed over significant timescales (>E4 years) in the stellar wind history. In combination, these observations give unprecedented insights into the nature of energetic explosions and their environments. During this talk, I will present highlights from recent observations, among them the first observation of a SN DURING the actual explosion with Swift, and I will discuss the "naked eye" burst at a redshift of -1, which was the most distant object humans could ever see with their own eyes.

Immler, Stefan↗
NASA NTRS1993

Chondrule formation by clumpy accretion onto the solar nebula

Chondrule textures and compositions appear to require rapid heating of precursor grain aggregates to temperatures in the range 1500 K to 2100 K, cooling times on the order of hours, and episodic and variable intensity events in order to produce chondrule rims and chemically distinct groups. Nebula shock waves have been proposed by Hood and Horanyi as a physical mechanism that may be capable of meeting the meteoritical constraints. Motivated by astronomical observations of the close environments of young stars, we suggest that the source of the nebula shock waves may be clumpy accretion onto the solar nebula - that is, episodic impacts onto the nebula by discrete cloud clumps with masses of at least 10(exp 22) g. If the cloud clumps are massive enough (10(exp 26) g), the resulting shockwave may be able to propagate to the midplane and process precursor aggregates residing in a dust sub-disk.

Boss, A. P.↗
DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Three-Dimensional Structure and Modeling of a Normal Bifurcated Shock Train from Experimental Measurements

An asymmetric, bifurcated normal shock train in a Mach 2 constant area, rectangular duct is investigated in the University of Michigan Direct Connect Isolator facility. High-speed schlieren imaging, wall static pressure measurements, surface oil flow visualization, and particle image velocimetry of a shock train are synthesized into a three-dimensional representation of shock train structure. This visualization is then used to inform the underlying flow physics of the distributed fluid dynamical processes along the structure. The detailed, three-dimensional morphology of the flow profile entering the shock train is shown to have a significant impact on the separated flow morphology of the shock train. This results in skewed supersonic core flow with alternating separation bubbles in the surrounding boundary layers. Finally, the efficacy of the prevailing pseudo-shock models in the literature is analyzed with the support of the available flow measurements. None are found to simultaneously model both the pressure and streamwise Mach number profile of the UMDCI pseudo-shock.

Hypersonics↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Automated Extraction of Flow Features

Computational Fluid Dynamics (CFD) simulations are routinely performed as part of the design process of most fluid handling devices. In order to efficiently and effectively use the results of a CFD simulation, visualization tools are often used. These tools are used in all stages of the CFD simulation including pre-processing, interim-processing, and post-processing, to interpret the results. Each of these stages requires visualization tools that allow one to examine the geometry of the device, as well as the partial or final results of the simulation. An engineer will typically generate a series of contour and vector plots to better understand the physics of how the fluid is interacting with the physical device. Of particular interest are detecting features such as shocks, recirculation zones, and vortices (which will highlight areas of stress and loss). As the demand for CFD analyses continues to increase the need for automated feature extraction capabilities has become vital. In the past, feature extraction and identification were interesting concepts, but not required in understanding the physics of a steady flow field. This is because the results of the more traditional tools like; iso-surface, cuts and streamlines, were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of a great deal of interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snapshot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for (co-processing environments). Methods must be developed to abstract the feature of interest and display it in a manner that physically makes sense.

Dorney, Suzanne↗
NASA NTRS2005

Automated Extraction of Flow Features

Computational Fluid Dynamics (CFD) simulations are routinely performed as part of the design process of most fluid handling devices. In order to efficiently and effectively use the results of a CFD simulation, visualization tools are often used. These tools are used in all stages of the CFD simulation including pre-processing, interim-processing, and post-processing, to interpret the results. Each of these stages requires visualization tools that allow one to examine the geometry of the device, as well as the partial or final results of the simulation. An engineer will typically generate a series of contour and vector plots to better understand the physics of how the fluid is interacting with the physical device. Of particular interest are detecting features such as shocks, re-circulation zones, and vortices (which will highlight areas of stress and loss). As the demand for CFD analyses continues to increase the need for automated feature extraction capabilities has become vital. In the past, feature extraction and identification were interesting concepts, but not required in understanding the physics of a steady flow field. This is because the results of the more traditional tools like; isc-surface, cuts and streamlines, were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of a great deal of interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one "snapshot" of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments). Methods must be developed to abstract the feature of interest and display it in a manner that physically makes sense.

Dorney, Suzanne↗
NASA NTRS1977

Dynamical response of the solar corona. III - Numerical simulation of the 1973 June 10 coronal transient

A spherically symmetric adiabatic single-fluid model is outlined for simulating the nonlinear time-dependent response of the corona to solar events that are simulated by perturbations in the appropriate physical variables from their steady-state values at the coronal base. Several observed features of the coronal transient that occurred on June 10, 1973, are simulated by using a particular steady-state solar wind, a specific combination of density and temperature perturbations, and a particular time dependence of the perturbations. A different steady-state solar wind, a perturbation of shorter duration, and other perturbation combinations are also employed to simulate the same transient so that the effect of each quantity can be determined. It is found that the model cannot adequately simulate all the observational results for the investigated transient, that the steady-state solar wind is relatively unimportant in such numerical simulations, and that studies which attempt only to reproduce observed shock trajectories may lead to erroneous conclusions regarding the physics of the solar event that produced the transient.

Steinolfson, R. S.↗
NASA NTRS1981

Flares and dynamic aspects

Theory of the dynamic nature of solar flares, eruptive prominences, and coronal transients is surveyed. Two types of flare are considered: compact eruptions (small flares) which are probably loops, and large eruptions of double filament flares. The physical sense of this classification is justified, using solar observations. Even leaving the question unanswered of whether a coronal transient is a loop of ejected matter or a shock wave, the important, perhaps fundamental, role of mass motion in flare physics is evidenced.

Hanssen, E. T.↗
NASA NTRS2016

Blunt-Body Aerothermodynamic Database from High-Enthalpy CO2 Testing in an Expansion Tunnel

An extensive database of heating, pressure, and flow field measurements on a 70-deg sphere-cone blunt body geometry in high-enthalpy, CO2 flow has been generated through testing in an expansion tunnel. This database is intended to support development and validation of computational tools and methods to be employed in the design of future Mars missions. The test was conducted in an expansion tunnel in order to avoid uncertainties in the definition of free stream conditions noted in previous studies performed in reflected shock tunnels. Data were obtained across a wide range of test velocity/density conditions that produced various physical phenomena of interest, including laminar and transitional/turbulent boundary layers, non-reacting to completely dissociated post-shock gas composition and shock-layer radiation. Flow field computations were performed at the test conditions and comparisons were made with the experimental data. Based on these comparisons, it is recommended that computational uncertainties on surface heating and pressure, for laminar, reacting-gas environments can be reduced to +/-10% and +/-5%, respectively. However, for flows with turbulence and shock-layer radiation, there were not sufficient validation-quality data obtained in this study to make any conclusions with respect to uncertainties, which highlights the need for further research in these areas.

Hollis, Brian R.↗
NASA NTRS1977

Numerical investigation of turbulence models for shock separated boundary-layer flows

Numerical solutions of the Navier-Stokes equations for shock separated turbulent boundary-layer flows are presented. Several turbulence models are investigated and assessed by their ability to predict the physical phenomena associated with two extensively documented experiments. The experimental flows consist of shock-wave boundary-layer interactions in axisymmetric internal and external geometries at Mach numbers of 1.5 and 7, respectively. Algebraic and one-equation eddy viscosity models are used to describe the Reynolds shear stress. Calculated values of skin friction, wall pressure distribution, kinetic energy of turbulence, and heat transfer are compared with measurements.

Viegas, J. R.↗
NASA NTRS2019

Using Computational Flow Imaging to Optimize Filtered Rayleigh Scattering Measurements of an Isolator Shock Train

Filtered Rayleigh scattering (FRS) is a laser diagnostic where the intensity of elastically scattered light is measured after it passes through a molecular absorption filter. The filter removes background interference that overwhelms the relatively weak scattering from the gas molecules. However, with a filter, the measured light intensity depends on many of the scattering gas properties, including pressure, density, temperature, and velocity. In this work, CFD simulations of an isolator shock train flow field are input into a physics-based model to predict the values of FRS intensity measurements in a proposed experiment. The goal is to evaluate if a simplistic FRS setup (that utilizes one camera, laser, and absorption filter) can be used to accurately quantify number density despite the fact that the scattered light intensity also depends on other gas properties. It is found that the experimental setup can be optimized such that a linear relationship describes the number density with an average prediction error of 2%. The vast majority of the flow exhibits prediction errors of less than 3%, but small regions of the flow reach up to 11% error. A sensitivity analysis shows that the prediction error increases with the central wavelength of the incident laser light and decreases with the angle between the camera and laser propagation directions. The optimal experimental parameters are chosen based on a compromise between the prediction error, spatial resolution, and the amount of unfiltered light. In the future, the proposed FRS setup will be implemented to acquire new and valuable information on an isolator shock train, a flow field that has been traditionally studied using wall static pressure measurements and path-integrated visualization techniques, such as schlieren and shadowgraphy.

Hunt, Robin L.↗
NASA NTRS2014

The Radiation, Interplanetary Shocks, and Coronal Sources (RISCS) Toolset

We outline a plan to develop a physics based predictive toolset RISCS to describe the interplanetary energetic particle and radiation environment throughout the inner heliosphere, including at the Earth. To forecast and "nowcast" the radiation environment requires the fusing of three components: 1) the ability to provide probabilities for incipient solar activity; 2) the use of these probabilities and daily coronal and solar wind observations to model the 3D spatial and temporal heliosphere, including magnetic field structure and transients, within 10 AU; and 3) the ability to model the acceleration and transport of energetic particles based on current and anticipated coronal and heliospheric conditions. We describe how to address 1) - 3) based on our existing, well developed, and validated codes and models. The goal of RISCS toolset is to provide an operational forecast and "nowcast" capability that will a) predict solar energetic particle (SEP) intensities; b) spectra for protons and heavy ions; c) predict maximum energies and their duration; d) SEP composition; e) cosmic ray intensities, and f) plasma parameters, including shock arrival times, strength and obliquity at any given heliospheric location and time. The toolset would have a 72 hour predicative capability, with associated probabilistic bounds, that would be updated hourly thereafter to improve the predicted event(s) and reduce the associated probability bounds. The RISCS toolset would be highly adaptable and portable, capable of running on a variety of platforms to accommodate various operational needs and requirements.

Zank, G. P.↗
DOE OSTI2024

Electron Heating in the Transrelativistic Perpendicular Shocks of Tilted Accretion Flows

Abstract General relativistic magnetohydrodynamic (GRMHD) simulations of black hole tilted disks—where the angular momentum of the accretion flow at large distances is misaligned with respect to the black hole spin—commonly display standing shocks within a few to tens of gravitational radii from the black hole. In GRMHD simulations of geometrically thick, optically thin accretion flows, applicable to low-luminosity sources like Sgr A* and M87*, the shocks have transrelativistic speed, moderate plasma beta (the ratio of ion thermal pressure to magnetic pressure is β pi1 ∼ 1–8), and low sonic Mach number (the ratio of shock speed to sound speed is M s ∼ 1–6). We study such shocks with 2D particle-in-cell simulations, and we quantify the efficiency and mechanisms of electron heating for the special case of preshock magnetic fields perpendicular to the shock direction of propagation. We find that the postshock electron temperature T e2 exceeds the adiabatic expectation T e2,ad by an amount T e 2 / T e 2 , ad − 1 ≃ 0.0016 M s 3.6 , nearly independent of the plasma beta and of the preshock electron-to-ion temperature ratio T e1 / T i1 , which we vary from 0.1 to unity. We investigate the heating physics for M s ∼ 5–6 and find that electron superadiabatic heating is governed by magnetic pumping at T e1 / T i1 = 1, whereas heating by B -parallel electric fields (i.e., parallel to the local magnetic field) dominates at T e1 / T i1 = 0.1. Our results provide physically motivated subgrid prescriptions for electron heating at the collisionless shocks seen in GRMHD simulations of black hole accretion flows.

Astronomy & Astrophysics↗
NASA NTRS1993

Higher-order viscous shock-layer solutions for high altitude flows

A higher-order viscous shock-layer (HVSL) method has been developed and is used to obtain physically consistent results under varying degrees of low-density conditions for perfect-gas and nonequilibrium flows past long slender bodies. This method is implemented with higher-order body and shock slip conditions. Detailed comparisons with the direct simulation Monte Carlo (DSMC) method and Navier-Stokes (NS) calculations establish HVSL as a computationally efficient code for the low-density calculations. Further, it is shown that the deficiencies in the standard viscous shock-layer (VSL) method when predicting low-density flows lie in the governing equations and not with the slip boundary conditions as considered by earlier researchers.

Gupta, Roop N.↗