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

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↗

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↗

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.↗

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↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗

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↗

Screech Noise Generation From Supersonic Underexpanded Jets Investigated

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

Panda, Jayanta↗

New Transient Finite Energy Shock Prediction Methodology

Shock prediction is one of the top loads and dynamics discipline technical challenges identified within NASA and industry programs and projects. The physics-based Transient Finite Energy (TFE) shock prediction methodology has been developed and compared favorably against test results. TFE can provide another approach to develop predictions of shock response spectra (SRS) for use in the analysis of structural margins.

NESC Technical Bulletin No. 23-03↗

Electron acceleration at nearly perpendicular collisionless shocks. I - One-dimensional simulations without electron scale fluctuations

Under certain conditions electrons can be reflected and effectively energized at quasi-perpendicular shocks. This process is most prominent close to the point where the upstream magnetic field is tangent to the curved shock. A theoretical explanation of the underlying physical mechanism has been proposed which assumes conservation of magnetic moment and a static, simplified shock profile are performed. Test particle calculations of the electron reflection process in order to examine the results of the theoretical analysis without imposing these restrictive conditions. A one-dimensional hybrid simulation code generates the characteristic field variations across the shock. Special emphasis is placed on the spatial and temporal length scales involved in the mirroring process. The simulation results agree generally well with the predictions from adiabatic theory. The effects of the cross-shock potential and unsteadiness are quantified, and the influence of field fluctuations on the reflection process is discussed.

Krauss-Varban, D.↗

Particle Acceleration, X-Rays, and Gamma-Rays From Winds

The instability of the line-driven winds of hot stars leads to the formation of strong shocks. These shocks not only emit thermal X-rays, but also accelerate a small fraction of the thermal electrons and ions to relativistic energies. Synchrotron radiation from these energetic particles can account for the non-thermal radio emission observed from some hot stars, and can also explain the hard X-rays detected in the Einstein X-ray spectra. Our calculations indicate that the gamma-ray emission from non-thermal particles should be detectable by GRO. The detection (or non-detection) of these emissions over a wide energy range, from the radio to gamma-rays, should provide a great deal of information on the structure of the unstable winds and the physics of particle acceleration by shocks.

White, Richard L.↗

Particle acceleration, X-rays, and gamma-rays from winds

The instability of the line-driven winds of hot stars leads to the formation of strong shocks. These shocks not only emit thermal X-rays, but also accelerate a small fraction of the thermal electrons and ions to relativistic energies. Synchrotron radiation from these energetic particles can account for the non-thermal radio emission observed from some hot stars, and can also explain the hard X-rays detected in the Einstein X-ray spectra. Our calculations indicate that the gamma-ray emission from non-thermal particles should be detectable by Gamma Ray Observatory (GRO). The detection (or non-detection) of these emissions over a wide energy range, from the radio to gamma-rays, should provide a great deal of information on the structure of the unstable winds and the physics of particle acceleration by shocks.

White, Richard L.↗

A viscous flow study of shock-boundary layer interaction, radial transport, and wake development in a transonic compressor

A numerical study based on the 3D Reynolds-averaged Navier-Stokes equation has been conducted to investigate the detailed flow physics inside a transonic compressor. 3D shock structure, shock-boundary layer interaction, flow separation, radial mixing, and wake development are all investigated at design and off-design conditions. Experimental data based on laser anemometer measurements are used to assess the overall quality of the numerical solution. An additional experimental study to investigate end-wall flow with a hot-film was conducted, and these results are compared with the numerical results. Detailed comparison with experimental data indicates that the overall features of the 3D shock structure, the shock-boundary layer interaction, and the wake development are all calculated very well in the numerical solution. The numerical results are further analyzed to examine the radial mixing phenomena in the transonic compressor. A thin sheet of particles is injected in the numerical solution upstream of the compressor. The movement of particles is traced with a 3D plotting package. This numerical survey of tracer concentration reveals the fundamental mechanisms of radial transport in this transonic compressor.

Hah, Chunill↗