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 541 records · Page 30

Simulation of the Flow Field Associated with a Rocket Thruster Having an Attached Panel

Two-dimensional inviscid and viscous numerical simulations are performed to predict the flow field induced by a H2-O2 rocket thruster and to provide insight into the heat load on the articles placed in the hot gas exhaust of the thruster under a variety of operating conditions, using the National Combustion Code (NCC). The simulations have captured physical details of the flow field, such as the plume formation and expansion, formation of the shock waves and their effects on the temperature and pressure distributions on the walls of the apparatus and the flat panel. Comparison between the computed results for 2-D and adiabatic walls and the related experimental measurements for 3-D and cooled walls shows that the results of the simulations are consistent with those obtained from the related rig tests.

Davoudzadeh, Farhad↗

Globular cluster origin of X-ray bursters

X-ray bursters and galactic bulge X-ray sources, or the most luminous X-ray sources in the Galaxy, are reasonably well constrained in their basic nature but not in their origin. It is suggested they may all have been produced by tidal capture in high density cores of globular clusters, which have now largely been disrupted by tidal stripping and shocking in the galactic plane. General arguments are presented for cluster disruption by the possible ring of giant molecular clouds in the Galaxy. Tests of the cluster disruption hypothesis are in progress and preliminary results are summarized here. The G-K star 'companions' previously noted for at least four bursters have spectra (in the two cases observed) consistent with metal-rich cluster giants. Several possibilities are discussed, including the formation of hierarchical triples in the dissolving cluster or in the galactic plane.

Grindlay, J. E.↗

The shocking development of lithium (and boron) in supernovae

It is shown that significant amounts of Li-7 and B-11 are produced in Type 2 supernovae. The synthesis of these rare elements occurs as the supernova shock traverses the base of the hydrogen envelope burning He-3 to masses 7 and 11 via alpha capture. The yields in this process are sufficient to account for the difference in lithium abundance observed between Pop 2 and Pop 1 stars. Since lithium (and boron) would, in this manner, be created in the same stars that produce the bulk of the heavy elements, the lithium abundance even in old Pop 1 stars would be high (as observed). The B-11 production may remedy the long-standing problem of the traditional spallation scenario to account for the observed isotopic ratio of boron. Observational consequences of this mechanism are discussed, including the evolution of lithium and boron isotope ratios in the Galaxy and the possible use of the boron yields to constrain the number of blue progenitor Type 2 supernovae.

Dearborn, David S. P.↗

The shocking development of lithium (and boron) in supernovae

It is shown that significant amounts of Li-7 and B-11 are produced in Type 2 supernovae. The synthesis of these rare elements occurs as the supernova shock traverses the base of the hydrogen envelope burning He-3 to masses 7 and 11 via alpha capture. The yields in this process are sufficient to account for the difference in lithium abundance observed between Pop 2 and Pop 1 stars. Since lithium (and boron) would, in this manner, be created in the same stars that produce the bulk of the heavy elements, the lithium abundance even in old Pop 1 stars would be high (as observed). The B-11 production may remedy the long-standing problem of the traditional spallation scenario to account for the observed isotopic ratio of boron. Observational consequences of this mechanism are discussed, including the evolution of lithium and boron isotope ratios in the Galaxy and the possible use of the boron yields to constrain the number of blue progenitor Type 2 supernovae.

Dearborn, David S. P.↗

Computational Analysis of a Low-Boom Supersonic Inlet

A low-boom supersonic inlet was designed for use on a conceptual small supersonic aircraft that would cruise with an over-wing Mach number of 1.7. The inlet was designed to minimize external overpressures, and used a novel bypass duct to divert the highest shock losses around the engine. The Wind-US CFD code was used to predict the effects of capture ratio, struts, bypass design, and angles of attack on inlet performance. The inlet was tested in the 8-ft by 6-ft Supersonic Wind Tunnel at NASA Glenn Research Center. Test results showed that the inlet had excellent performance, with capture ratios near one, a peak core total pressure recovery of 96 percent, and a stable operating range much larger than that of an engine. Predictions generally compared very well with the experimental data, and were used to help interpret some of the experimental results.

Chima, Rodrick V.↗

An Experimental Investigation of Incompressible Richtmyer-Meshkov Instability

Richtmyer-Meshkov (RM) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and eventually become a turbulent flow. It is closely related to Rayleigh-Taylor instability, which is the instability of a planar interface undergoing constant acceleration, such as caused by the suspension of a heavy fluid over a lighter one in the earth's gravitational field. Like the well-known Kelvin-Helmholtz instability, RM instability is a fundamental hydrodynamic instability which exhibits many of the nonlinear complexities that transform simple initial conditions into a complex turbulent flow. Furthermore, the simplicity of RM instability (in that it requires very few defining parameters), and the fact that it can be generated in a closed container, makes it an excellent test bed to study nonlinear stability theory as well as turbulent transport in a heterogeneous system. However, the fact that RM instability involves fluids of unequal densities which experience negligible gravitational force, except during the impulsive acceleration, requires RM instability experiments to be carried out under conditions of microgravity. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The impulsive acceleration is generated by bouncing a rectangular tank containing two different density liquids off a retractable vertical spring. The initial perturbation is produced prior to release by oscillating the tank in the horizontal direction to produce a standing wave. The instability evolves in microgravity as the tank travels up and then down the vertical rails of a drop tower until hitting a shock absorber at the bottom. Planar Laser Induced Fluorescence (PLIF) is employed to visualize the flow. PLIF images are captured by a video camera that travels with the tank. Figure 1 is as sequence of images showing the development of the instability from the initial sinusoidal disturbance far into the nonlinear regime which is characterized by the appearance of mushroom structures resulting from the coalescence of baroclinic vorticity produced by the impulsive acceleration. At later times in this sequence the vortex cores are observed to become unstable showing the beginnings of the transition to turbulence in this flow. The amplitude of the growing disturbance after the impulsive acceleration is measured and found to agree well with theoretical predictions. The effects of Reynolds number (based on circulation) on the development of the vortices and the transition to turbulence are also determined.

Jacobs, J. W.↗

Eucrite Impact Melt NWA 5218 - Evidence for a Large Crater on Vesta

Northwest Africa (NWA) 5218 is a 76 g achondrite that is classified as a eucrite [1]. However, an initial classification [2] describes it as a "eucrite shock-melt breccia...(in which) large, partially melted cumulate basalt clasts are set in a shock melt flow...". We explore the petrology of this clast-bearing impact melt rock (Fig. 1), which could be a characteristic lithology at large impact craters on asteroid Vesta [3]. Methods: Optical microscopy, scanning electronmicroscopy, and Raman spectroscopy were used on a thin section (Fig. 1) for petrographic characterization. The impact melt composition was determined by 20 m diameter defocused-beam analyses with a Cameca SX-100 electron microprobe. The data from 97 spots were corrected for mineral density effects [4]. Constituent mineral phases were analyzed with a focusedbeam. Bidirectonal visible and near-infrared (VNIR) and biconical FT-IR reflectance spectra were measured on the surface of a sample slab on its central melt area and on an eucrite clast, and from 125-500 m and <125 m powders of melt. Results: General petrography: The sample specimen is a coherent, medium dark-grey (N4), melt rock. The thin section captures a central, subophitic-textured melt that contains ~1 cm to tens of m-size subangular to rounded, variably-shocked eucrite clasts. Clasts >100 m are coarse-grained with equigranular ~1 mm size plagioclase, quartz, and clinopyroxene (Fig. 1). Single crystals of chromite, ilmenite, zircon, Ca-Mg phosphate, Fe-metal, and troilite are embedded in the melt. Polymineralic clasts are mostly compositionally similar to the above mentioned larger clasts but scarce granulitic fragments are observed as well.

Wittmann, Axel↗

Origin of the Moon

The similarities in siderophile abundances strongly suggest that the Moon was derived from the Earth's mantle after the Earth's core had formed. The energy required to remove material from the Earth's mantle and place it into geocentric orbit can be supplied most readily by impact processes during accretion of the Earth. Impacts of late-accreting, high-velocity planetismals would evaporate many times their masses of mantle material. These gases would be accompanied by a massive spray of shock-melted silicate droplets. It is suggested that the gases produced from such near-equatorial impacts were rapidly spun out into equatorial geocentric orbit. The evaporated material was selectively recondensed, and, accompanied by the shock melted, devolatilized silicate droplets, accreted to form a sediment ring of Earth-orbiting planetismals. This sediment ring also captured a significant proportion of Earth-bound planetismals. The Moon was formed by accretion from planetismals comprising the sediment ring.

Ringwood, A. E.↗

Uncertainty Assessments of 2D and Axisymmetric Hypersonic Shock Wave - Turbulent Boundary Layer Interaction Simulations at Compression Corners

This paper is one of a series of five papers in a special session organized by the NASA Fundamental Aeronautics Program that addresses uncertainty assessments for CFD simulations in hypersonic flow. Simulations of a shock emanating from a compression corner and interacting with a fully developed turbulent boundary layer are evaluated herein. Mission relevant conditions at Mach 7 and Mach 14 are defined for a pre-compression ramp of a scramjet powered vehicle. Three compression angles are defined, the smallest to avoid separation losses and the largest to force a separated flow engaging more complicated flow physics. The Baldwin-Lomax and the Cebeci-Smith algebraic models, the one-equation Spalart-Allmaras model with the Catrix-Aupoix compressibility modification and two-equation models including Menter SST, Wilcox k-omega 98, and Wilcox k-omega 06 turbulence models are evaluated. Each model is fully defined herein to preclude any ambiguity regarding model implementation. Comparisons are made to existing experimental data and Van Driest theory to provide preliminary assessment of model form uncertainty. A set of coarse grained uncertainty metrics are defined to capture essential differences among turbulence models. Except for the inability of algebraic models to converge for some separated flows there is no clearly superior model as judged by these metrics. A preliminary metric for the numerical component of uncertainty in shock-turbulent-boundary-layer interactions at compression corners sufficiently steep to cause separation is defined as 55%. This value is a median of differences with experimental data averaged for peak pressure and heating and for extent of separation captured in new, grid-converged solutions presented here. This value is consistent with existing results in a literature review of hypersonic shock-turbulent-boundary-layer interactions by Roy and Blottner and with more recent computations of MacLean.

Gnoffo, Peter A.↗

Simulation of Low-density Nozzle Plumes in Non-zero Ambient Pressures

The direct simulation Monte-Carlo (DSMC) method was applied to the analysis of low-density nitrogen plumes exhausting from a small converging-diverging nozzle into finite ambient pressures. Two cases were considered that simulated actual test conditions in a vacuum facility. The numerical simulations readily captured the complicated flow structure of the overexpanded plumes adjusting to the finite ambient pressures, including Mach disks and barrel shaped shocks. The numerical simulations compared well to experimental data of Rothe.

Chung, Chan-Hong↗

Axisymmetric Calculations of a Low-Boom Inlet in a Supersonic Wind Tunnel

This paper describes axisymmetric CFD predictions made of a supersonic low-boom inlet with a facility diffuser, cold pipe, and mass flow plug within wind tunnel walls, and compares the CFD calculations with the experimental data. The inlet was designed for use on a small supersonic aircraft that would cruise at Mach 1.6, with a Mach number over the wing of 1.7. The inlet was tested in the 8-ft by 6-ft Supersonic Wind Tunnel at NASA Glenn Research Center in the fall of 2010 to demonstrate the performance and stability of a practical flight design that included a novel bypass duct. The inlet design is discussed here briefly. Prior to the test, CFD calculations were made to predict the performance of the inlet and its associated wind tunnel hardware, and to estimate flow areas needed to throttle the inlet. The calculations were done with the Wind-US CFD code and are described in detail. After the test, comparisons were made between computed and measured shock patterns, total pressure recoveries, and centerline pressures. The results showed that the dual-stream inlet had excellent performance, with capture ratios near one, a peak core total pressure recovery of 96 percent, and a large stable operating range. Predicted core recovery agreed well with the experiment but predicted bypass recovery and maximum capture ratio were high. Calculations of offdesign performance of the inlet along a flight profile agreed well with measurements and previous calculations.

Chima, Rodrick V.↗

Three-dimensional supersonic vortex breakdown

Three-dimensional supersonic vortex-breakdown problems in bound and unbound domains are solved. The solutions are obtained using the time-accurate integration of the unsteady, compressible, full Navier-Stokes (NS) equations. The computational scheme is an implicit, upwind, flux-difference splitting, finite-volume scheme. Two vortex-breakdown applications are considered in the present paper. The first is for a supersonic swirling jet which is issued from a nozzle into a supersonic uniform flow at a lower Mach number than that of the swirling jet. The second is for a supersonic swirling flow in a configured circular duct. In the first application, an extensive study of the effects of grid fineness, shape and grid-point distribution on the vortex breakdown is presented. Four grids are used in this study and they show a substantial dependence of the breakdown bubble and shock wave on the grid used. In the second application, the bubble-type and helix-type vortex breakdown have been captured.

Kandil, Osama A.↗

Low Dissipative High Order Numerical Simulations of Supersonic Reactive Flows

The objective of this paper is to evaluate the performance of a newly developed low dissipative sixth-order spatial and fourth-order temporal scheme for viscous reactive flows interacting with shock waves that contain fine scale flow structures. The accuracy and efficiency of the scheme, and to what degree the scheme can capture the correct physical wave speeds of stiff reactive flows will be included.

Sjoegreen, B.↗

Representing Extremes in Agricultural Models

AgMIP and related projects are conducting several activities to understand and improve crop model response to extreme events. This involves crop model studies as well as the generation of climate datasets and scenarios more capable of capturing extremes. Models are typically less responsive to extreme events than we observe, and miss several forms of extreme events. Models also can capture interactive effects between climate change and climate extremes. Additional work is needed to understand response of markets and economic systems to food shocks. AgMIP is planning a Coordinated Global and Regional Assessment of Climate Change Impacts on Agricultural Production and Food Security with an aim to inform the IPCC Sixth Assessment Report.

Agricultural systems↗

Numerical Study of Interaction of a Vortical Density Inhomogeneity with Shock and Expansion Waves

We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.

Povitsky, A.↗

Measurement of Material Recession and Shock Standoff in Plasma Windtunnel using Neural Nets

Arcjets are plasma wind tunnels used to test the performance of heatshield materials for spacecraft atmospheric entry. These facilities present an extremely harsh flow environment with heat fluxes up to 109 W/m2 for up to 30 minutes. The plasma is low-temperature (∼1 eV) but high pressure (> 10 kPa) creating high-enthalpy supersonic flows similar to atmospheric entry conditions. Typically, material samples are measured before and after a test to characterize the total recession. However, this does not capture time-dependent effects such as material expansion and non-linear recession. This work will present new analysis of arcjet test videos which measure both the time-dependent 2D recession of the material samples and the shock standoff distance. The results show non-linear time-dependent effects are present for some conditions. The material and shock edges are extracted from the videos by training and applying a convolutional neural network. Due to the consistent camera settings, the machine learning model achieves high accuracy (± 2 px) relative to manually segmented images with only a small number of training frames (80).

Neural network↗

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

Two-dimensional planar and axisymmetric numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to develop a computational methodology to identify nozzle side load physics using simplified two-dimensional geometries, in order to come up with a computational strategy to eventually predict the three-dimensional side loads. The computational methodology is based on a multidimensional, finite-volume, viscous, chemically reacting, unstructured-grid, and pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system modeling. The side load physics captured in the low aspect-ratio, two-dimensional planar nozzle include the Coanda effect, afterburning wave, and the associated lip free-shock oscillation. Results of parametric studies indicate that equivalence ratio, combustion and ramp rate affect the side load physics. The side load physics inferred in the high aspect-ratio, axisymmetric nozzle study include the afterburning wave; transition from free-shock to restricted-shock separation, reverting back to free-shock separation, and transforming to restricted-shock separation again; and lip restricted-shock oscillation. The Mach disk loci and wall pressure history studies reconfirm that combustion and the associated thermodynamic properties affect the formation and duration of the asymmetric flow.

Wang, Ten-See↗

Measurement of material recession and shock standoff in plasma windtunnel using neural nets

Arcjets are plasma wind tunnels used to test the performance of heatshield materials for spacecraft atmospheric entry. These facilities present an extremely harsh flow environment with heat fluxes up to 10^9 W/m^2 for up to 30 minutes. The plasma is low-temperature (~1 eV) but high pressure (> 10 kPa) creating high-enthalpy supersonic flows similar to atmospheric entry conditions. Typically, material samples are measured before and after a test to characterize the total recession. However, this does not capture time-dependent effects such as material expansion and non-linear recession. This work will present new analysis of arcjet test videos which measure both the time-dependent 2D recession of the material samples and the shock standoff distance. New results showing non-linear material erosion rates will be highlighted. The material and shock edges are extracted from the videos by training and applying a convolutional neural network. Due to the consistent camera settings, the machine learning model achieves high accuracy (~99%) on new data with only a small number of training frames (~80). The new results will be discussed in the context of temperature dependent plasma-surface interaction.

Magnus A Haw↗