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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cooling flow regime of a plasma thermal quench

A large class of Laboratory, Space, and Astrophysical plasmas is nearly collisionless. When a localized energy or particle sink, for example, in the form of a radiative cooling spot or a black hole, is introduced into such a plasma, it can trigger a plasma thermal collapse, also known as a thermal quench in tokamak fusion. Here we show that the electron thermal conduction in such a nearly collisionless plasma follows the convective energy transport scaling in itself or in its spatial gradient, due to the constraint of ambipolar transport. As a result, a robust cooling flow aggregates mass toward the cooling spot and the thermal collapse of the surrounding plasma takes the form of four propagating fronts that originate from the radiative cooling spot, along the magnetic field line in a magnetized plasma. The slowest one, which is responsible for deep cooling, is a shock front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shocks in the Scrape-off Layer: Hydrodynamical Aspects (Summary) [Slides]

Hydrodynamical model provides a robust description of mass, momentum, and energy transfer. It can be used for benchmarking hydrodynamical modules of SOL codes. It clarifies the role of boundary conditions at the end-plates. It provides a “geometrical” model for partition of power flows between two divertor legs, both for the steady-state and pulsed flows. It indicates possible presence of shocks in the divertor legs. A framework for finding the shock location and strength has been developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

On the numerical solution of two-dimensional, laminar compressible flows with imbedded shock waves.

The complete, time-dependent Navier-Stokes equations are expressed in conservation form and solved by employing an explicit finite difference numerical technique which incorporates artificial viscosity terms of the form first suggested by Rusanov for numerical stability in the vicinity of shock waves. Surface boundary conditions are developed in a consistent and unique manner through the use of a physically oriented extrapolation procedure. From numerical experimentation an extended range for the explicit stability parameter is established. Also employed is an additional convergence parameter which relates incremental spatial steps. Convergence of the transient solution to a steady state flow was obtained after 400 to 500 time steps.

Goodrich, W. D.↗
NASA NTRS1992

Practical aspects of spatially high accurate methods

The computational qualities of high order spatially accurate methods for the finite volume solution of the Euler equations are presented. Two dimensional essentially non-oscillatory (ENO), k-exact, and 'dimension by dimension' ENO reconstruction operators are discussed and compared in terms of reconstruction and solution accuracy, computational cost and oscillatory behavior in supersonic flows with shocks. Inherent steady state convergence difficulties are demonstrated for adaptive stencil algorithms. An exact solution to the heat equation is used to determine reconstruction error, and the computational intensity is reflected in operation counts. Standard MUSCL differencing is included for comparison. Numerical experiments presented include the Ringleb flow for numerical accuracy and a shock reflection problem. A vortex-shock interaction demonstrates the ability of the ENO scheme to excel in simulating unsteady high-frequency flow physics.

Godfrey, Andrew G.↗
NASA NTRS1993

Unsteady Phenomena in Shock Wave/Boundary Layer Interaction

A brief review is given of the unsteadiness of shock wave/turbulent boundary layer interaction. The focus is on interactions generated by swept and unswept compression ramps, by flares, steps and incident shock waves, by cylinders and blunt fins, and by glancing shock waves. The effects of Mach number, Reynolds number, and separated flow scale are discussed as are the physical causes of the unsteadiness. The implications that the unsteadiness has for interpreting time-average surface and flowfield data, and for comparisons of such experimental data with computation, is also briefly discussed. Finally, some suggestions for future work are given. It is clear that there are large gaps in the data base and that many aspects of such phenomena are poorly understood. Much work remains to be done.

Dolling, D. S.↗
NASA NTRS2000

Solar Spicules: A Review of Recent Models and Targets for Future Observations

Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfven speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis on the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere, and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most- commonly -quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating TRACE satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.

Sterling, Alphonse C.↗
NASA NTRS2000

Solar Spicules: A Review of Recent Models and Targets for Future Observations

Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfven speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis an the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most commonly quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating Transition Region and Coronal Explorer (TRACE) satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.

Sterling, Alphonse C.↗
NASA NTRS2003

Effects of Wing Leading Edge Penetration with Venting and Exhaust Flow from Wheel Well at Mach 24 in Flight

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of: (1) holes through the leading edge into a vented cavity; and (2) a scarfed, conical nozzle directed toward the centerline of the vehicle from the forward, inboard corner of the landing gear door. The simulations were generated relatively quickly and early in the investigation because simplifications were made to the leading edge cavity geometry and an existing utility to merge scarfed nozzle grid domains with structured baseline external domains was implemented. These simplifications in the breach simulations enabled: (1) a very quick grid generation procedure; and (2) high fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a two-inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the six-inch diameter breach the boundary layer is fully ingested. The intent of externally directed jet simulations in the second scenario was to approximately model aerodynamic effects of a relatively large internal wing pressure, fueled by combusting aluminum, which deforms the corner of the landing gear door and directs a jet across the windside surface. These jet interactions, in and of themselves, were not sufficiently large to explain observed aerodynamic behavior.

Gnoffo, Peter A.↗
NASA NTRS2004

Simulation of Flow Through Breach in Leading Edge at Mach 24

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of holes through the leading edge into a vented cavity. The simulations were generated relatively quickly and early in the investigation by making simplifications to the leading edge cavity geometry. These simplifications in the breach simulations enabled: 1) A very quick grid generation procedure; 2) High fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a 2 inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the 6 inch diameter breach the boundary layer is fully ingested.

Gnoffo, Peter A.↗
NASA NTRS2012

Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-ray Observations of the Coma Cluster of Galaxies with VERITAS and FERMI

Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated and, on account of weak energy losses, can accumulate, suggesting that clusters may also be sources of very high energy (VHE; E greater than100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cerenkov telescopes, with complementing Fermi Large Area Telescope observations at GeV energies. No significant gamma-ray emission from the Coma Cluster was detected. Integral flux upper limits at the 99 confidence level were measured to be on the order of (2-5) x 10(sup -8) photons m(sup -2) s(sup -1) (VERITAS,greater than 220 GeV) and approximately 2 x 10(sup -6) photons m(sup -2) s(sup -1) (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain cosmic rays (CRs) and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be less than 16% from VERITAS data and less than 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be 50. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of approximately (2-5.5)microG, depending on the radial magnetic field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this renders the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, (sigma upsilon)

cosmic↗
DOE OSTI2021

Exceptionally high spallation strength for a high-entropy alloy demonstrated by experiments and simulations

High-entropy alloys are materials with an increasing number of technological applications. Amongst them, the Cantor alloy, FeMnCoCrNi, shows desirable mechanical properties at normal loading conditions. In this study we focus on the performance of the Cantor alloy at the ultra-high deformation rates of shock waves. We study shock-induced spallation using both experiments and atomistic simulations. Experimental loading is achieved using high power laser, with VISAR to obtain velocity profiles and spall strength, followed by transmission electron microscopy of the recovered samples. Molecular Dynamics (MD) simulations of shock-induced spallation are compared with experiments. Both experiments and simulations show a high spall strength which would be beneficial for certain applications, with experiments giving ~8 GPa at ~10 7 s -1 and MD giving almost ~30 GPa at ~10 9 s -1 . The difference between experiments and simulations can be explained by the difference in strain rate. Post-mortem analysis of the experimental samples shows nanotwins near the spall plane, while MD simulations show a highly disordered region giving rise to void nucleation and spall during loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The colliding planar shocks platform to study warm dense matter at the National Ignition Facility

We have developed an experimental platform at the National Ignition Facility that employs colliding planar shocks to produce warm dense matter with uniform conditions and enable high-precision equation of state measurements. The platform uses simultaneous x-ray Thomson scattering and x-ray radiography to measure the density, electron temperature, and ionization state in warm dense matter. The experimental platform is designed to create a large volume of uniform plasma (approximately 700×700×150μm 3 ) at pressures approaching 100 Mbar and minimize the distribution of plasma conditions in the x-ray scattering volume, significantly improving the precision of the measurements. Here, in this study, we present the experimental design of the platform and compare hydrodynamic simulations to x-ray radiography data from initial experiments studying hydrocarbons, producing uniform densities within ±25% of the average probed condition. We show that the platform creates a homogeneous plasma that can be characterized using x-ray Thomson scattering. Thus, the new platform enables accurate measurements of plasma conditions necessary to test models for the equation of state and ionization potential depression in the warm dense matter regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1968

Simple test for physical stability of cryogenic tank insulation

Qualitative test determines the ability of insulation liners used on liquid hydrogen tanks to withstand stresses produced by the thermal shocks imparted to the insulation during tank filling and drainage. Test specimens are bonded to metal plates with a low thermal expansion coefficient and are immersed in liquid hydrogen.

Rossello, D.↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
DOE OSTI2021

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗