Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock capturing”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computational shock formation & development: An arbitrary Lagrangian-Eulerian characteristics approach [Slides]

We propose a new computational shock formation-development algorithm. We use a “good” geometry adapted to the evolving solution, along with a good set of variables defined in this geometry. We are able to: (1) Accurately capture the pre-shock; (2) Track distinguished characteristics and capture weak discontinuities; (3) Approximate solutions to classical Riemann problems; and (4) Accurately solve challenging problems for which standard methods fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2021

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth’s bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the “foot” region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock’s rest frame. That conversion allows for a unique study of the relative spatial scales of the shock’s various features, including the shock’s growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Evidence is also presented of nonlinear wave decay from electromagnetic whistler-mode “lion roars” to electrostatic solitary waves in the downstream plasma regime. Results highlight the critical cross-scale coupling between electron-kinetic-and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

D L Turner↗
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Time-split finite-volume method for three-dimensional blunt-body flow

An efficient numerical method for calculating plane, axisymmetric, and fully three-dimensional blunt-body flow is presented. It is a second-order-accurate, time-dependent finite-volume procedure that solves the Euler equations in integral conservation-law form. These equations are written with respect to a Cartesian coordinate system in which an embedded mesh adjusts in time to the motion of the bow shock that is automatically captured as part of the weak solution. With such an adjusting mesh, oscillations in flow properties near the shock are shown to be virtually eliminated. The scheme uses a time-splitting concept that accelerates the convergence appreciably. Comparisons are made between computed and experimental results.

Rizzi, A. W.↗
NASA NTRS2001

Adaptive Numerical Dissipation Controls for High Order Methods

A numerical scheme for direct numerical simulation of shock-turbulence interactions of high speed compressible flows would ideally not be significantly more expensive than the standard fourth or sixth-order compact or non-compact central differencing scheme. It should be possible to resolve all scales down to scales of order of the Kolmogorov scales of turbulence accurately and efficiently, while at the same time being able to capture steep gradients occurring at much smaller scales efficiently. The goal of this lecture is to review the progress and new development of the low dissipative high order shock-capturing schemes proposed by Yee et al. Comparison on the efficiency and accuracy of this class of schemes with spectral and the fifth-order WENO (weighted essentially nonoscillatory) scheme will be presented. A new approach to dynamically sense the appropriate amount of numerical dissipation to be added at each grid point using non-orthogonal wavelets will be discussed.

Yee, Helen C.↗
NASA NTRS1972

Numerical computation of multishocked, three-dimensional supersonic flow fields with real gas effects.

A computational procedure is presented which is capable of determining the supersonic flow field surrounding three-dimensional wing-body configurations such as a delta-wing space shuttle. The governing equations in conservation-law form are solved by a finite difference method using a second-order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically, and the intersection of these shocks with the bow shock posed no difficulty. Resulting flow fields about typical blunt nose shuttle-like configurations at angle of attack are presented. The differences between perfect and real gas effects for high Mach number flows are shown.

Kutler, P.↗
NASA NTRSDate not supplied

Zemax Modeling of a Self-aligned Focusing Schlieren System with Simulated Compressible Flow

A self-aligned focusing schlieren system is simulated in the ANSYS Zemax OpticStudio software environment. The model uses a custom scatter-defining dynamic link library on the retroreflective background in order to scatter rays in the direction of the incident ray rather than the specular ray. Simulated compressible flow refractive index fields are imported into Zemax by use of the grid gradient object and a workaround to the maximum point limit is introduced. The Ronchi ruling is built as a high resolution bitmap image in MATLAB and set in Zemax as a slide object allowing the bitmap image to either transmit or terminate rays at each line pair. The Zemax simulations are compared to computational numerical schlieren and experimental conventional, path-integrated schlieren results from a hypersonic stagnation point injection flow configuration and an experimental image of a shock wave from a laser-induced breakdown. In the stagnation point injection simulations, the Zemax schlieren images capture complex flow features such as the Mach disk, shear layers, contact surface, and bow shock well, but contains long, curved features (striations) between the model and the bow shock that are not present in the synthetic schlieren image from the CFD or the experimental image. The image from Zemax of the spark-induced shock wave captures the shape and curvature of the shock but appears blurry due to grid discretization in the computational data.

Nicholas A Mejia↗
NASA NTRS2009

Shock Layer Radiation Measurements and Analysis for Mars Entry

NASA's In-Space Propulsion program is supporting the development of shock radiation transport models for aerocapture missions to Mars. A comprehensive test series in the NASA Antes Electric Arc Shock Tube facility at a representative flight condition was recently completed. The facility optical instrumentation enabled spectral measurements of shocked gas radiation from the vacuum ultraviolet to the near infrared. The instrumentation captured the nonequilibrium post-shock excitation and relaxation dynamics of dispersed spectral features. A description of the shock tube facility, optical instrumentation, and examples of the test data are presented. Comparisons of measured spectra with model predictions are also made.

Bose, Deepak↗
NASA NTRS2009

CFD Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy, grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex. Although the one-dimensional isentropic nozzle plume results look reasonable, they fail to capture the sonic boom shock train in the highly underexpanded nozzle flow.

Bui, Trong T.↗
NASA NTRS1975

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Computational methods for shock waves in three-dimensional supersonic flow

A central artificial-viscosity and an upwind-biased difference method are contrived to solve the Euler equations for flowfields over typical spacecrafts. The spatial discretization is based on either nodal or cell-vertex formulation in the domain extending from free stream to the end of the vehicle. The outer boundary is treated as a bow shock in the first method but is placed in the free stream in the second, which captures both bow and internal shocks using an approximate Riemann solver based on high-order extrapolation to the cell face. These methods were tested for the Shuttle and Hermes orbiters at wind-tunnel conditions and angles of attack ranging from 0 to 60 deg. The artificial-viscosity method incorporated with a shock-fitting procedure shows smeared crossflow and wing-shock positions and required 15 percent more CPU per node than the upwind method. Greater flexibility and robustness is demonstrated by the latter on a fixed grid for all cases considered.

Li, Chien-Peng↗
DOE OSTI2025

Rippled shock propagation in a laser-driven target at multimegabar pressures

The evolution of non-uniform shocks produced by modulated laser irradiation or surface perturbations is relevant to studies of inertial confinement fusion and material properties at high-energy-density conditions. We present results from an experiment conducted at the OMEGA EP laser facility, where a 300 GPa shock was driven into a fused silica sample with pre-fabricated single-mode surface modulations. Using time-resolved optical velocimetry, we captured the continuous evolution of rippled shock motion, enabling a comprehensive mapping of the spatial amplitude history from formation to phase reversal in a single experiment. Initially, the ablation-driven shock inherits a fraction of the surface modulation amplitude from the sample, which subsequently grows before decaying, ultimately leading to the flattening of the rippled shock and a phase reversal. We find that two-dimensional inviscid hydrodynamic simulation of the experiment is able to qualitatively capture many aspects of the rippled shock evolution but over-predicts the initial amplitude growth. This experimental platform, capable of accommodating varying ripple wavelengths, lays the groundwork for a potential viscometry method at extreme pressures, where viscous effects manifest as differences in shock flattening times between rippled shocks of two distinct wavelengths propagating through the sample.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Towards Understanding the Fluid Dynamic Phenomenon of Interest to Rocket Base Heating: A Review

The significance of the base heating problem for rockets during ascent is due to the complex interaction between the rocket nozzle plumes and the external-flow which can change the flow characteristics in the base region dramatically. At lower altitudes the external-flow merges with the plume-flow, without the formation of a large separated flow region, and the cooler external-flow promotes convective cooling of the base wall. Under these conditions the majority of the base heating is due to radiative heating from the shock heated plume gases. At higher altitudes, however, the process of base heating is not so straightforward. The plume and the base flow expands dramatically and separated flow regions occur in the base area. Hot exhaust gases from the rocket nozzle will be entrained into the separated flow regions and produce a convective component to the base wall heating. Further, if the rocket exhaust-gas contains soot, the soot can increase the emission from the gas and dramatically increase the wall absorption coefficient for radiative heating if it is deposited on the walls . In addition, if the rocket exhaust gas is fuel rich, the fuel can bum in the separated flow regions and further increase the base heating. The base burning phenomenon, and the increased base heating caused by it at higher altitudes, have been observed for the Space Shuttle and Saturn Rocket. Under these conditions, the total heating is significantly higher than the heating without separated flow in the base region, and the increase in heating is directly attributable to the fluid dynamic complexity of the base region. Realistic simulation of the base heating requires that the calculated flow environment reproduce the fluid dynamic flow features accurately. Thus, it will be necessary to introduce into the CFD codes the capability for the flow to respond to the complex vehicle geometry, the effect of turbulence, the ability to accurately reproduce the plume shock/shear layer structures and capture such phenomenon as shock induced base separation and base-burning phenomenon. A survey of experimental, theoretical and computational work that details the fluid dynamics of the base flow environment will be presented in the proposed paper. CFD simulations of rocket base flows using standard CFD codes such as OVERFLOW or GASP will be explored to capture these phenomenon accurately. Merits and limitations of these codes for base flow environment predictions will be explored.

Venkatapathy, E.↗
DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

High-Speed Edge-Detecting Line Scan Smart Camera

A high-speed edge-detecting line scan smart camera was developed. The camera is designed to operate as a component in a NASA Glenn Research Center developed inlet shock detection system. The inlet shock is detected by projecting a laser sheet through the airflow. The shock within the airflow is the densest part and refracts the laser sheet the most in its vicinity, leaving a dark spot or shadowgraph. These spots show up as a dip or negative peak within the pixel intensity profile of an image of the projected laser sheet. The smart camera acquires and processes in real-time the linear image containing the shock shadowgraph and outputting the shock location. Previously a high-speed camera and personal computer would perform the image capture and processing to determine the shock location. This innovation consists of a linear image sensor, analog signal processing circuit, and a digital circuit that provides a numerical digital output of the shock or negative edge location. The smart camera is capable of capturing and processing linear images at over 1,000 frames per second. The edges are identified as numeric pixel values within the linear array of pixels, and the edge location information can be sent out from the circuit in a variety of ways, such as by using a microcontroller and onboard or external digital interface to include serial data such as RS-232/485, USB, Ethernet, or CAN BUS; parallel digital data; or an analog signal. The smart camera system can be integrated into a small package with a relatively small number of parts, reducing size and increasing reliability over the previous imaging system..

Prokop, Norman F.↗
NASA NTRS1982

Computation of viscous hypersonic flow over control surfaces

A parabolized Navier-Stokes (PNS) computer code has been used to compute hypersonic flow over control surfaces located on the aft end of generic re-entry configurations. The Beam-Warming implicit numerical algorithm as modified by Schiff and Steger is used to march in the spatial direction. The outer bow shock is treated as a sharp discontinuity; within the shock layer discontinuities are captured. Several different turbulent flow solutions were generated for flap angles varying from -7 deg to 15 deg and angles of attack varying from 0 deg to 20 deg. The present results are qualitatively correct when compared with previous compression corner studies.

Chaussee, D. S.↗