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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The geocentric particulate distribution: Cometary, asteroidal, or space debris?

Definition of the Low Earth Orbit (LEO) particulate environment has been refined considerably with the analysis of data from NASA's Long Duration Exposure Facility (LDEF). Measurements of the impact rates from particulates ranging from sub-micron to millimetres in dimension and, especially, information on their directionality has permitted new scrunity of the sources of the particulates. Modelling of the dynamics of both bound (Earth orbital) and unbound (hyperbolic interplanetary) particulates intercepting LDEF's faces leads to the conclusion that the source is dominantly interplanetary for particle dimensions of greater than some 5 microns diameter; however the anisotropy below this dimension demands lower velocities and is compatible with an orbital component. Characteristics of the LDEF interplanetary component are compatible with familiar meteoroid sources and deep space measurements. Understanding of the orbital component which exceeds the interplanetary flux by a factor of 4 is less clear; although the very small particulates in orbit have been associated with space debris (Lawrance and Brownlee, 1986) this data conflicts with other measurements (McDonnell, Carey and Dixon, 1984) at the same epoch. By analysis of trajectories approaching the Earth and its atmosphere, we have shown that a significant contribution could be captured by aerocapture, i.e., atmospheric drag, from either asteroidal or cometary sources; such enhancement is unlikely however to provide the temporal and spatial fluctuations observed by the LDEF Interplanetary Dust Experiment (Mullholland et al. 1992). A further new mechanism is also examined, that of aerofragmentation capture, where an atmospheric grazing trajectory, which would not normally lead to capture, leads to fragmentation by thermal or mechanical shock; the microparticulates thus created can be injected in large numbers, but only into short-lifetime orbits. The concentration in one particular orbit plane, could explain the temporal fluctuations seen on LDEF; space debris could also explain the phenomenon.

Mcdonnell, J. A. M.↗
NASA NTRS1993

A new procedure for dynamic adaption of three-dimensional unstructured grids

A new procedure is presented for the simultaneous coarsening and refinement of three-dimensional unstructured tetrahedral meshes. This algorithm allows for localized grid adaption that is used to capture aerodynamic flow features such as vortices and shock waves in helicopter flowfield simulations. The mesh-adaption algorithm is implemented in the C programming language and uses a data structure consisting of a series of dynamically-allocated linked lists. These lists allow the mesh connectivity to be rapidly reconstructed when individual mesh points are added and/or deleted. The algorithm allows the mesh to change in an anisotropic manner in order to efficiently resolve directional flow features. The procedure has been successfully implemented on a single processor of a Cray Y-MP computer. Two sample cases are presented involving three-dimensional transonic flow. Computed results show good agreement with conventional structured-grid solutions for the Euler equations.

Biswas, Rupak↗
NASA NTRS1993

Numerical analysis of base flowfield at high altitude for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a pressure based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution and turbulence model are two primary factors that influence the accuracy of the base flowfield prediction.

Wang, Ten-See↗
NASA NTRS1993

Numerical study of base pressure characteristic curve for a four-engine clustered nozzle configuration

The objective of this study is to benchmark a four-engine clustered nozzle base flowfield with a computational fluid dynamics (CFD) model. The CFD model is a three-dimensional pressure-based, viscous flow formulation. An adaptive upwind scheme is employed for the spatial discretization. The upwind scheme is based on second and fourth order central differencing with adaptive artificial dissipation. Qualitative base flow features such as the reverse jet, wall jet, recompression shock, and plume-plume impingement have been captured. The computed quantitative flow properties such as the radial base pressure distribution, model centerline Mach number and static pressure variation, and base pressure characteristic curve agreed reasonably well with those of the measurement. Parametric study on the effect of grid resolution, turbulence model, inlet boundary condition and difference scheme on convective terms has been performed. The results showed that grid resolution had a strong influence on the accuracy of the base flowfield prediction.

Wang, Ten-See↗
NASA NTRS2004

Transient Three-Dimensional Analysis of Side Load in Liquid Rocket Engine Nozzles

Three-dimensional numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, and pressure-based computational fluid dynamics formulation, and a simulated inlet condition based on a system calculation. Finite-rate chemistry was used throughout the study so that combustion effect is always included, and the effect of wall cooling on side load physics is studied. The side load physics captured include the afterburning wave, transition from free- shock to restricted-shock separation, and lip Lambda shock oscillation. With the adiabatic nozzle, free-shock separation reappears after the transition from free-shock separation to restricted-shock separation, and the subsequent flow pattern of the simultaneous free-shock and restricted-shock separations creates a very asymmetric Mach disk flow. With the cooled nozzle, the more symmetric restricted-shock separation persisted throughout the start-up transient after the transition, leading to an overall lower side load than that of the adiabatic nozzle. The tepee structures corresponding to the maximum side load were addressed.

Wang, Ten-See↗
NASA NTRS1972

Studies of engine-airframe integrated hypersonic vehicles

A parametric study of an integrated airframe and engine is presented for a hypersonic transport at an altitude of 70,000 feet and a free stream Mach number of 6. The engine considered is a subsonic combustion ramjet using conventional hydrocarbon fuels. The lift-to-drag ratio of the aircraft for two configurations, one with full capture and accelerated flight and the other allowing spillage of the leading shock and in unaccelerated flight, is studied. The parameters varied are the engine efficiencies, the angle of attack, the combustion rates, as well as the captured mass flow. Lift-to-drag ratios on the order of 6.5 are obtained.

Saland, H.↗
NASA NTRS1993

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗
DOE OSTI2021

In Situ X-ray Diffraction of High Density Polyethylene During Dynamic Compression: Polymer Chain Compression and Phase Transformations

Polymers are widely implemented materials often used for safety and security in our daily lives. However, we do not have a thorough understanding of the phase transformations of these materials when exposed to extreme condition, high strain rate, events. To study polymers in these severe environments, we performed plate impact experiments that impart a one-dimensional shock wave, at high pressures and temperatures, to observe the materials response. Generally, velocimetry is the experimental diagnostic used to capture the movement of high velocity interfaces to then calculate the mechanical variables associated with shock compression. However, when a material undergoes a phase transition, such as between two crystalline phases or melting to decomposition during shock compression, velocimetry does not explicitly indicate what structural change occurred. At the Dynamic Compression Sector at the Advanced Photon Source gas gun driven compression has been paired with x-ray diffraction, providing the ability to study both the shock parameters and structural changes concurrently. We explored high density polyethylene (HDPE) both on its reactants Hugoniot (<10.12 GPa) and within its transition region (24.3 GPa) where HDPE is hypothesized to decompose from reactants to products. With x-ray diffraction we observed both polymer chain compression and crystal structure reorientation from orthorhombic to monoclinic above 6 GPa, as well as either a possible phase change or chemical decomposition at 24.3 GPa. With these experiments we were able to determine the crystalline and amorphous Hugoniot loci at pressures below the polymer transition region.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Effects of Fin Leading Edge Sweep on Shock-Shock Interaction at Mach 6

The effects of fin leading edge sweep on peak heating rates due to shock-shock interaction have been experimentally examined in the Langley 20-Inch Mach 6 Tunnel. The shock interaction was produced by the intersection of a planar incident shock (16.8 deg shock angle relative to the freestream, generated by a 9 deg wedge) with the bow shock formed around a O.5-inch diameter cylindrical leading edge fin. Heating distributions along the leading edge stagnation line have been obtained using densely spaced thin film resistive-type sensors. Schlieren images were obtained to illustrate the very complex shock-shock interactions. The fin leading edge sweep angle was varied from 15-degrees swept back to 45-degrees swept forward for a freestream unit Reynolds number of 2 x 10(exp 6)/ft. Two models were utilized during the study, one with 0.025-inch spacing between gage centers, and the other 0.015-inch spacing. Gage spatial resolution on the order of 0.015-in appeared to accurately capture the narrow spike in heating. Peak heating due to shock interaction was maximized when the fin was swept forward 15 deg and 25 deg, both promoting augmentations about 7 times the baseline value. The schlieren images for these cases revealed Type 4 and Type 3 interactions, respectively.

Berry, Scott A.↗
NASA NTRS1992

Internal shock interactions in propulsion/airframe integrated three-dimensional sidewall compression scramjet inlets

The advantages and design requirements of propulsion/airframe integration for high Mach number flight are studied in terms of the 3D sidewall compression scramjet inlet. The present work addresses in a parametric fashion the inviscid effects of leading edge sweep, sidewall compression, and inflow Mach number on the internal shock structure in terms of inlet compression and mass capture. The source of the Mach number invariance with leading edge sweep for a constant sidewall compression class of inlet is identified, and a previously undocumented spillage phenomenon in a constant effective wedge angle class of inlets is discussed.

Holland, Scott D.↗
NASA NTRS1992

Internal Shock Interactions in Propulsion/Airframe Integrated Three-Dimensional Sidewall Compression Scramjet Inlets

The advantages and design requirements of propulsion/airframe integration for high Mach number flight have led to extensive study of the three-dimensional sidewall compression scramjet inlet in recent years. Recent research publications have indicated testing over a broad range of Mach number (2 to 18) in a variety of test gases, such as air, helium, and tetrafluoromethane. Multiple experimental techniques have been employed to obtain detailed internal shock interaction data, performance data, and inlet starting limits. Computational fluid dynamics has been effectively used for preliminary parametric studies as well as in parallel with experiments to aid in the explanation of unusual or unexpected flow phenomena. Inlets of this genre afford a relatively simple, generic geometry while producing a highly complex, three-dimensional flow field dominated by shock/shock and shock/boundary layer interactions. While the importance of the viscous effects in high speed inlet interactions is recognized, the present work addresses in a parametric fashion the inviscid effects of leading edge sweep, sidewall compression, and inflow Mach number on the internal shock structure in terms of inlet compression and mass capture. In the process, the source of the of the Mach number invariance with leading edge sweep for a constant sidewall compression class of inlet is identified, and a previously undocumented spillage phenomenon in a constant effective wedge angle class of inlets is discussed.

Holland, Scott D.↗
NASA NTRS1996

Interplanetary and Interstellar Dust Particles: Reentry Heating and Capture in Aerogel

Studies of interplanetary and interstellar dust can provide significant information on the evolution of the solar system or stars, respectively. However, for reliable analysis it is crucial to know how the particles have been modified during reentry (in the case of interplanetary dust particles, or [DP's) and impact into the capture medium. In the case of stratospheric capture, particles will be heated by atmospheric drag. Subsequent capture of the particles will result in heating, ablation, accretion of the capture medium and possible fragmentation. Modeling of these processes is a useful way of assessing their effects on the interpretation of the compositional data for these particles. Previous work on reentry heating has shown that heat diffusion alone cannot adequately account for temperature gradients observed in IDP's. In fact, for any reasonable thermal parameters, calculations show the particles to be nearly isothermal. Here we extend those studies to include phase transitions. These preliminary results are promising and show a significant temperature gradient for a 100 micron (diameter) particle. The actual capture of the particles in silica aerogel is being modeled using a comprehensive shock hydrodynamic code (called CTH). Various options of this code were investigated to attempt to make the most appropriate choice of methods of impact, equations of state, and processes of energy transfer from capture material to particle. The initial calculations with the code used only 'reasonable' estimates for the physical parameters of silica aerogel. Through the literature searches and personal contacts with the knowledgeable scientists, the best possible mechanical and thermal data have been made available for these computations.

Szydlik, Paul P.↗
NASA NTRS2015

Computational and Experimental Study of Supersonic Nozzle Flow and Shock Interactions

This study focused on the capability of NASA Tetrahedral Unstructured Software System's CFD code USM3D capability to predict the interaction between a shock and supersonic plume flow. Previous studies, published in 2004, 2009 and 2013, investigated USM3D's supersonic plume flow results versus historical experimental data. This current study builds on that research by utilizing the best practices from the early papers for properly capturing the plume flow and then adding a wedge acting as a shock generator. This computational study is in conjunction with experimental tests conducted at the Glenn Research Center 1'x1' Supersonic Wind Tunnel. The comparison of the computational and experimental data shows good agreement for location and strength of the shocks although there are vertical shifts between the data sets that may be do to the measurement technique.

Carter, Melissa B.↗
NASA NTRS2012

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2022

Spatial and Temporal Interpolation Analysis Process of Shock Loading

The following document serves to describe the current process that is utilized by W-13 analysts to approximate the pressure-time boundary condition seen by a test object from shock tube loading. The intent of this report is to capture the current capability and to encourage continued growth and development. Information on the required input variables, mapping process, and generated output information is covered in the following sections. Throughout the document, gaps and weaknesses of the mapping process are noted to encourage future development efforts.

42 ENGINEERING↗