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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1967

Fluid physics

Three-dimensional viscous disturbances, shock structure of helium-argon mixtures, and magnetohydrodynamic experiments with large interaction parameter

MAGNETOHYDRODYNAMIC STABILITY↗
NASA NTRS1988

Accretion shock geometries in the magnetic variables

The first self consistent shock models for the AM Herculis-type systems successfully identified the dominant physical processes and their signatures. These homogenous shock models predict unpolarized, Rayleigh-Jeans optical spectra with sharp cutoffs and rising polarizations as the shocks become optically thin in the ultraviolet. However, the observed energy distributions are generally flat with intermediate polarizations over a broad optical band. These and other observational evidence support a non-homogenous accretion profile which may extend over a considerable fraction of the stellar surface. Both the fundamental assumptions underlying the canonical 1-D shock model and the extension of this model to inhomogenous accretion shocks were identified, for both radial and linear structures. The observational evidence was also examined for tall shocks and little evidence was found for relative shock heights in excess of h/R(1) greater than or equal to 0.1. For several systems, upper limits to the shock height can be obtained from either x ray or optical data. These lie in the region h/R(1) is approximately 0.01 and are in general agreement with the current physical picture for these systems. The quasi-periodic optical variations observed in several magnetic variables may eventually prove to be a major aid in further understanding their accretion shock geometries.

Stockman, H. S.↗
NASA NTRS1986

Configurations of corotating shocks in the outer heliosphere

Configurations of corotating shocks in the outer heliosphere are computed on the assumption that the shocks move at constant speeds between collisions. The basic physical process is the interaction of a forward shock with a reverse shock. The points at which this interaction takes place are determined algebraically in terms of the initial shock positions and speeds. Between the points at which the shocks interact, the shocks have the form of Archimedian spirals. Shock configurations are determined for the case of two corotating shock pairs originating at 2.5 AU and for a single shock pair at 2.5 AU. In both cases the heliosphere is divided into several distinct regions distinguished by the number of times the particles have passed through a shock. Since the shock strengths decrease when a forward shock interacts with a reverse shock, corotating shocks are expected to be weaker at larger distances from the sun.

Burlaga, L. F.↗
NASA NTRS1986

Why do interstellar grains exist?

There exists a discrepancy between calculated destruction rates of grains in the interstellar medium and postulated sources of new grains. This problem was examined by modelling the global life cycle of grains in the galaxy. The model includes: grain destruction due to supernovae shock waves; grain injection from cool stars, planetary nebulae, star formation, novae, and supernovae; grain growth by accretion in dark clouds; and a mixing scheme between phases of the interstellar medium. Grain growth in molecular clouds is considered as a mechanism or increasing the formation rate. To decrease the shock destruction rate, several new physical processes, such as partial vaporization effects in grain-grain collisions, breakdown of the small Larmor radius approximation for betatron acceleration, and relaxation of the steady-state shock assumption are included.

Seab, C. G.↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1949

Correlation of physical properties of ceramic materials with resistance to fracture by thermal shock

An analysis is made to determine which properties of materials affect their resistance to fracture by thermal stresses.From this analysis, a parameter is evaluated that is correlated with the resistance of ceramic materials to fracture by thermal shock as experimentally determined. This parameter may be used to predict qualitatively the resistance of a material to fracture by thermal shock. Resistance to fracture by thermal shock is shown to be dependent upon the following material properties: thermal conductivity, tensile strength, thermal expansion, and ductility modulus. For qualitative prediction of resistance of materials to fracture by thermal shock, the parameter may be expressed as the product of thermal conductivity and tensile strength divided by the product of linear coefficient of thermal expansion and ductility modulus of the specimen.

Lidman, W G↗
DOE OSTI2020

Electron Energy Partition across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Anomalous Shocks on the Measured Near-Field Pressure Signatures of Low-Boom Wind-Tunnel Models

Unexpected shocks on wind-tunnel-measured pressure signatures prompted questions about design methods, pressure signature measurement techniques, and the quality of measurements in the flow fields near lifting models. Some of these unexpected shocks were the result of component integration methods. Others were attributed to the three-dimension nature of the flow around a lifting model, to inaccuracies in the prediction of the area-ruled lift, or to wing-tip stall effects. This report discusses the low-boom model wind-tunnel data where these unexpected shocks were initially observed, the physics of the lifting wing/body model's flow field, the wind-tunnel data used to evaluate the applicability of methods for calculating equivalent areas due to lift, the performance of lift prediction codes, and tip stall effects so that the cause of these shocks could be determined.

Mack, Robert J.↗
DOE OSTI2021

A comparison of Raman and pyrometry dynamic temperature measurements of shocked cyclohexane

Optical pyrometry is a well-developed, broadly applicable method to determine temperature in dynamic compression experiments. However, it measures radiation from only the skin depth of the radiant material, is highly susceptible to spurious backgrounds, and requires an often-unmeasurable free parameter—the sample dynamic emissivity. Raman spectroscopy offers a method to determine temperature directly from fundamental theory, but it is practically limited in applications to a small subset of shock experiments. In this work, we describe an experiment allowing simultaneous application of Raman and pyrometric temperature techniques to benchmark a specific instance of pyrometry. We attempted both measurements on multiple experiments and compare the successful temperature results between the separate but highly repeatable shots. Here, we find that dynamic shock temperatures inferred from Raman spectroscopy and pyrometry generally agree within the uncertainties of each measurement, although those uncertainties are larger than typical for a specially designed pyrometry experiment on a shocked metal sample.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2014

Development and Transition of the Radiation, Interplanetary Shocks, and Coronal Sources (RISCS) Toolset

We outline a plan to develop and transition a physics based predictive toolset called The Radiation, Interplanetary Shocks, and Coronal Sources (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 Astronomical Units; 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. The described transition plan is based on a well established approach developed in the Earth Science discipline that ensures that the customer has a tool that meets their needs

Spann, James F.↗
NASA NTRS1994

Dense molecular shocks and accretion onto protostellar disks

We have examined the physical and chemical processes which accompany the propagation of a J-type shock through molecular astrophysical gas of density 10(exp 7.5) - 10(exp 12) H nuclei per cu cm. Our study is primarily relevant to the high-density accretion shocks which are associated with the supersonic infall of material during the collapse of a molecular cloud core to form a protostar. We have carried out a general parameter study to determine the characteristic infrared emission spectrum of dense molecular shocks and the extent of grain destruction within such shocks as a function of relevant shock parameters. Accretion shocks with mass accretion rates approximately greater than 10(exp -5) solar masses per y are plausible sources of the CO v = 2-0 vibrational band emission observed in many protostars, at least for sources of moderate CO v = 2-0 luminosity (approximately less than few x 10(exp 23) W). Rovibrational emissions from H2O and OH are also predicted to be luminous in dense molecular accretion shocks. Vaporization is the dominant mechanism of grain destruction within dense molecular shocks, leading to the efficient removal of those grains which are heated to their vaporization temperature by the combined effects of radiative heating, drag heating, and thermal heating by the hot postshock gas. The criteria for grain vaporization are obtained for grains composed of several different materials: metallic iron, silicate, troilite, refractory organic materials, volatile organic materials, and water ice.

Neufeld, David A.↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1993

Prependicular transport in 1- and 2-dimensional shock simulations

We consider the foundations of 1- and 2-dimensional shock simulations in which the physical quantities are independent of a coordinate which is not parallel to the magnetic field. We show analytically that in such simulations the ions are effectively tied to the convected magnetic lines of force because of the presence of an ignorable spatial coordinate. This conclusion has important consequences. In particular we conclude that the acceleration of energetic charged particles at quasi-perpendicular shocks cannot be properly studied in such simulations because the role of perpendicular diffusion cannot be properly evaluated.

Jokipii, J. R.↗
NASA NTRS1998

Ionospheric Instabilities and Shock Capturing Algorithms

We describe the basic physics of the gradient-drift instability and the collisional Rayleigh-Taylor instability, the two plasma instabilities which dominate the production of small-scale structure in the earth's ionosphere. We then focus our attention on the nonlinear evolution of these instabilities for two particular cases: 1) the recursive bifurcation and "freezing" of artificial ion clouds released in the ionosphere; and 2) equatorial spread F. Finally, we discuss the numerical challenges faced in attempting to numerically simulate these instabilities, and the role played by shock-capturing algorithms therein.

Zalesak, Steven T.↗
DOE OSTI2022

Subcell limiting strategies for discontinuous Galerkin spectral element methods

Here, we present a general family of subcell limiting strategies to construct robust high-order accurate nodal discontinuous Galerkin (DG) schemes. The main strategy is to construct compatible low order finite volume (FV) type discretizations that allow for convex blending with the high-order variant with the goal of guaranteeing additional properties, such as bounds on physical quantities and/or guaranteed entropy dissipation. For an implementation of this main strategy, four main ingredients are identified that may be combined in a flexible manner: (i) a nodal high-order DG method on Legendre–Gauss–Lobatto nodes, (ii) a compatible robust subcell FV scheme, (iii) a convex combination strategy for the two schemes, which can be element-wise or subcell-wise, and (iv) a strategy to compute the convex blending factors, which can be either based on heuristic troubled-cell indicators, or using ideas from flux-corrected transport methods. By carefully designing the metric terms of the subcell FV method, the resulting methods can be used on unstructured curvilinear meshes, are locally conservative, can handle strong shocks efficiently while directly guaranteeing physical bounds on quantities such as density, pressure or entropy. We further show that it is possible to choose the four ingredients to recover existing methods such as a provably entropy dissipative subcell shock-capturing approach or a sparse invariant domain preserving approach. We test the versatility of the presented strategies and mix and match the four ingredients to solve challenging simulation setups, such as the KPP problem (a hyperbolic conservation law with non-convex flux function), turbulent and hypersonic Euler simulations, and MHD problems featuring shocks and turbulence.

97 MATHEMATICS AND COMPUTING↗