Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple shocks”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

The Effect of Shock on the Amorphous Component in Altered Basalt

Investigation of the geochemical and mineralogical composition of the Martian surface provides insight into the geologic history of the predominantly basaltic crust. The Chemistry and Mineralogy (CheMin) instrument onboard the Curiosity rover has returned the first X-Ray diffraction data from the Martian surface. However, large proportions (27 +/- 14 with some estimates as high as 50 weight percentage) of an amorphous component have been reported. As a remedy to this problem, mass balance equations using geochemistry, volatile chemistry, and mineralogy have been employed to constrain the geochemistry of the amorphous component. However, "the nature and number of amorphous phases that constitute the amorphous component is not unequivocally known". Multiple hypotheses have been proposed to explain the origin of this amorphous component: Allophane (Al2O); Basaltic glass (Volcanic and impact); Palagonite (Altered basaltic glass); Hisingerite (Fe (sup 3 plus)-bearing phyllosilicate); S/Cl-rich component (sulfates and/or akaganeite); Nanophase ferric oxide component (npOx). Establishing a multi-phase amorphous component from a basaltic precursor that has undergone physical and chemical weathering within geochemical constraints is of paramount importance to better understand the composition of a large portion of the Martian surface (up to 50 weight percentage). Shocked basalts from Lonar Crater in India are valuable analogs for the Martian surface because it is a well-preserved impact crater in a basaltic target. Having undergone pre- and post-shock aqueous alteration, these rocks provide crucial data regarding the effect of shock on the amorphous component in altered basalt. By conducting mass balance equations similar to what has been performed for Gale crater materials, we attempt to calculate the geochemistry of the amorphous component in altered basalts ranging from unshocked to Class 5 (Table 1). This has the potential to reveal the nature and origin (i.e. primary igneous, shock metamorphic, and/or aqueous alteration occurring before or after the impact event) of the amorphous component in shocked basalt with the goal of unravelling the history of the Martian surface.

Eckley, S. A.↗
NASA NTRS1985

Numerical solution of the Euler equation for compressible inviscid fluids

An effort is made to develop a satisfactory numerical method for the calculation of steady solutions of the Euler equations for inviscid compressible gas flows. The intended application is the prediction of the aerodynamic properties of aircraft flying at transonic speeds. Particular consideration is given to the modification of the equations to improve convergence to a steady state; finite volume formulation; adaptive dissipation; schemes designed to improve the resolution of shock waves; hybrid multistage time stepping schemes; residual averaging; a multigrid scheme; and trials with Burgers' equation. Results clearly demonstrate that the convergence of a time-dependent hyperbolic system to a steady state can be substantially accelerated by the introduction of multiple grids.

Jameson, A.↗
NASA NTRS1990

Three-dimensional compressible convection at low Prandtl numbers

Numerical simulations are used to study fully compressible thermal convection at large Rayleigh numbers. Results are presented from a sequence of three-dimensional simulations that reveal a transition from gradually-evolving laminar convection to nearly turbulent convection as the Prandtl number is reduced from a value of unity to one-tenth. The convective flows form irregular cellular patterns near the upper surface, possesing a network of fast downflow at cell peripheries and gentler upflow at cell centers. At greater depths the curving sheets of downflow collapse into plumes which may twist and possess substantial vertical vorticity. For the lowest Prandtl number, the convection near the bottom of the layer appears to be turbulent, yet the rapidly varying small-scale flow structure there is accompanied by more ordered sites of wavering upflow, with the latter able to penetrate all the way to the upper boundary. Thus a significant component of the flow is able to extend over multiple density scale heights, in contrast to what is argued in formulating mixing-length models for stellar convection. Results are also shown from two-dimensional simulations carried out with very high spatial resolution, which reveal that supersonic convection with fluttering shock systems can be realized.

Toomre, Juri↗
NASA NTRS1992

Computational fluid dynamics analysis of Space Shuttle main engine multiple plume flows at high-altitude flight conditions

Computational fluid dynamics (CFD) analysis is providing verification of Space Shuttle flight performance details and is being applied to Space Shuttle Main Engine Multiple plume interaction flow field definition. Advancements in real-gas CFD methodology that are described have allowed definition of exhaust plume flow details at Mach 3.5 and 107,000 ft. The specific objective includes the estimate of flow properties at oblique shocks between plumes and plume recirculation into the Space Shuttle Orbiter base so that base heating and base pressure can be modeled accurately. The approach utilizes the Rockwell USA Real Gas 3-D Navier-Stokes (USARG3D) Code for the analysis. The code has multi-zonal capability to detail the geometry of the plumes based region and utilizes finite-rate chemistry to compute the plume expansion angle and relevant flow properties at altitude correctly. Through an improved definition of the base recirculation flow properties, heating, and aerodynamic design environments of the Space Shuttle Vehicle can be further updated.

Dougherty, N. S.↗
NASA NTRS2012

Space Shuttle Program (SSP) Shock Test and Specification Experience for Reusable Flight Hardware Equipment

As commercial companies are nearing a preliminary design review level of design maturity, several companies are identifying the process for qualifying their multi-use electrical and mechanical components for various shock environments, including pyrotechnic, mortar firing, and water impact. The experience in quantifying the environments consists primarily of recommendations from Military Standard-1540, Product Verification Requirement for Launch, Upper Stage, and Space Vehicles. Therefore, the NASA Engineering and Safety Center (NESC) formed a team of NASA shock experts to share the NASA experience with qualifying hardware for the Space Shuttle Program (SSP) and other applicable programs and projects. Several team teleconferences were held to discuss past experience and to share ideas of possible methods for qualifying components for multiple missions. This document contains the information compiled from the discussions

Larsen, Curtis E.↗
DOE OSTI2019

Graphite to diamond transformation under shock compression: Role of orientational order

To gain insight into the role of orientational order on the shock-induced graphite to diamond phase transformation, three pyrolytic graphite types having different orientational orders were shock-compressed along the average c -axis to peak stresses between 35 and 69 GPa. The materials studied were ZYB-grade highly oriented pyrolytic graphite (HOPG), ZYH-grade HOPG, and as-deposited pyrolytic graphite (PG) having mosaic spreads of 0.8° ± 0.2°, 3.5° ± 1.5°, and ~45°, respectively. Wave profiles, obtained using laser interferometry, show a multiple-wave structure with a distinct, rapid (<10 ns) rise to the high-pressure phase for each graphite type. Multiple-wave profiles, first observed in this study for the less ordered ZYH-grade HOPG and PG samples, show that somewhat poorly oriented pyrolytic graphites also undergo a well-defined phase transformation. Previously, rapid transformation was reported for ZYB-grade but not ZYH-grade HOPG. The measured wave profiles for both HOPG grades are very similar and both grades show a ~22 GPa transformation stress. In contrast, the PG wave profiles are quite different and show a ~46 GPa transformation stress. The continuum results (stress-density states) presented here cannot distinguish between the different high-pressure phases [hexagonal diamond (HD) or cubic diamond] reported in recent x-ray studies. Because ZYB-grade HOPG was recently shown to transform to HD and due to the similar peak states for both HOPG grades, it seems likely that ZYH-grade also transforms into HD. In conclusion, the very different shock responses of PG and HOPG suggest different transformation mechanisms for PG and HOPG, but the high-pressure PG phase remains unclear in the present work.

58 GEOSCIENCES↗
NASA NTRS2011

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗
NASA NTRS2010

Oxidation States of Grim Glasses in EET79001 Based on Vanadium Valence

Gas-rich impact-melt (GRIM) glasses in SNC meteorites are very rich in Martian atmospheric noble gases and sulfur suggesting a possible occurrence of regolith-derived secondary mineral assemblages in these samples. Previously, we have studied two GRIM glasses, 506 and 507, from EET79001 Lith A and Lith B, respectively, for elemental abundances and spatial distribution of sulfur using EMPA (WDS) and FE-SEM (EDS) techniques and for sulfur-speciation using K-edge XANES techniques. These elemental and FE-SEM micro-graph data at several locations in the GRIM glasses from Shergotty (DBS), Zagami 994 and EET79001, Lith B showed that FeO and SO3 are positively correlated (SO3 represents a mixture of sulfide and sulfate). FE-SEM (EDS) study revealed that the sulfur-rich pockets in these glasses contain numerous micron-sized iron-sulfide (Fe-S) globules sequestered throughout the volume. However, in some areas (though less frequently), we detected significant Fe-S-O signals suggesting the occurrence of iron sulfate. These GRIM glasses were studied by K-edge microXANES techniques for sulfur speciation in association with iron in sulfur-rich areas. In both samples, we found the sulfur speciation dominated by sulfide with minor oxidized sulfur mixed in with various proportions. The abundance of oxidized sulfur was greater in 506 than in 507. Based on these results, we hypothesize that sulfur initially existed as sulfate in the glass precursor materials and, on shock-impact melting of the precursor materials producing these glasses, the oxidized sulfur was reduced to predominately sulfide. In order to further test this hypothesis, we have used microXANES to measure the valence states of vanadium in GRIM glasses from Lith A and Lith B to complement and compare with previous analogous measurements on Lith C (note: 506 and 507 contain the largest amounts of martian atmospheric gases but the gas-contents in Lith C measured by are unknown). Vanadium is ideal for addressing this re-dox issue because it has multiple valence states and is a well-studied element. Ferrous-dominated iron valences determined by microXANES on the Lith A and Lith B glasses provide little redox sensitivity. Vanadium valence measurements for impact glass in Lith C at three different locations yielded valence values of 3.1, 3.2 and 3.4 with inferred fO2 values of IW-0.7, IW-0.1 and IW+0.7, respectively. This range of oxygen-fugacity values is understandable because the glasses are shock-molten impact glasses which are heterogeneous in nature. Oxygen fugacity values obtained from the analysis of Fe-Ti oxides and Eu partitioning in pyroxenes from EET79001 Lith A and Lith B (host lithologies) were in the range of IW+0.3 to IW+1.9 suggesting that V in the Lith C impact glass was reduced in the impact process. Here, we examine whether the 506 from Lith A and 507 from Lith B GRIM glasses yield similar or different fO2 values from those of Lith C using the vanadium K-edge microXANES technique.

Sutton, S. R.↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Interfacial Strength Evaluation of Oxide Films on Carbon Steel by Using the Laser Shock Adhesion Test

We investigate the adhesion strength (interfacial strength) of the oxide film formed on a carbon steel substrate by using a laser shock adhesion test (LaSAT). A rolled carbon steel (SS400) plate is first heated at high temperature (700 °C) under the atmosphere (heat treatment) to cause the formation of oxide films about 20-50 μm thick on the steel surface. The film consisted of multiple layers, including wustite (FeO), magnetite (Fe{sub 3}O{sub 4}), and hematite (Fe{sub 2}O{sub 3}). Next, the x-ray diffraction method was used to measure the residual stress in the film. We find that compression stress develops in the film because of a mismatch of thermal strain in between the film and substrate after the heat treatment. With LaSAT, strong elastic wave generates by a pulsed laser ablation of the substrate backside, and results in interfacial fractures in the oxide films. To monitor the normal (out-of-plane) displacement of the top surface film, a laser ultrasonic interferometer is simultaneously utilized, such that we detect film delamination and determine the critical laser energy for interfacial fracture to obtain the critical stress (interfacial strength). Furthermore, additional pulsed laser irradiation is conducted against a pre-delaminated specimen in order to examine the growth of delamination area and evaluate the interfacial fracture toughness. Finally, elastic wave propagation in the sample is computed using the finite element method (FEM) in order to evaluate the interfacial stress field, accurately. For the FEM, multiple layers in the oxide film are modeled, and the residual stress due to thermal strain mismatch is introduced. The FEM computations reveal the stress distribution around the interface and evaluate the critical stress and critical stress intensity factor for delamination. We thus quantitatively evaluate the interfacial strength and interfacial fracture toughness of the oxide film formed on carbon steel. The mechanism of interfacial fracture is also discussed on the basis of microstructures in a film with multiple layers. This result may be useful for understanding the adhesion quality of oxide film during the hot-rolling process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Testing, Characterization, and Modeling of the Resonant Plate Test Environment

The resonant plate shock test is a dynamic test of a mid-field pyroshock environment where a projectile is struck against a plate. The structure undergoing the simulated field shock is mounted to the plate. The plate resonates when struck and provides a two sided shock that is representative of the shock observed in the field. This test environment shock simulates a shock in a single coordinate direction for components looking to provide evidence that they will survive a similar or less shock when deployed in their operating environment. However, testing in one axis at a time provides many challenges. The true environment is a multi-axis environment. The test environment exhibits strong off-axis motion when only motion in one axis is desired. Multiple fixtures are needed for a single test series. It would be advantageous if a single test could be developed that tests the multi-axis environment simultaneously. In order to design such a test, a model must be developed and validated. The model can be iterated in design and configuration until the specified multi-axis environment is met. The test can then execute the model driven test design. This report discusses the resonant plate model needed to design future tests and the steps and methods used to obtain the model. This report also details aspects of the resonant plate test discovered during the process of model development that aids in our understanding of the test.

42 ENGINEERING↗
NASA NTRS1972

The simulation of high pressure hydrogen/oxygen rocket engines.

Initial experiments have established the feasibility of using a backward running detonation wave in a shock tube to produce the correct enthalpy and chemistry of the simulated combustion chamber gases. This technique produces plume Reynolds numbers higher than those achievable by any other known technique. The practical pressure limitations are being explored and detailed surveys of both single and multiple engine plumes and impingement data on a cylindrical body are being obtained. Facility details to allow safe handling of the high pressure hydrogen and oxygen up to a design chamber pressure of 20,000 psi are described.

Hopkins, H.↗
NASA NTRS1990

Efficient parametric analysis of performance measures for communication networks

Efficient techniques for estimating performance measures in communication networks in a steady-state or transient setting are developed. These techniques may be used in a simulation environment or in connection with real-time observations. For Markov chain models, the recently proposed standard clock approach is extended, and a class of real-time algorithms for simultaneously generating multiple sample paths under different parameter sets is presented. Attention is focused on the link crash time estimation problem, where a 'crash' is defined as the first time a buffer overflows, given some initial conditions. An algorithm for estimating crash times under various traffic shocks is derived, where all estimates are obtained in parallel to an actual network's normal operation. An algorithm for crash time estimation of a G/D/1 link model is also derived using a different (perturbation-analysis-based) approach. Finally, extensive simulation results are provided to validate the proposed algorithms and compare them to brute-force simulation.

Cassandras, Christos G.↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Multiple-cycle Simulation of a Pulse Detonation Engine Ejector

This paper presents the results of a study involving single and multiple-cycle numerical simulations of various PDE-ejector configurations utilizing hydrogen-oxygen mixtures. The objective was to investigate the thrust, impulse and mass flow rate characteristics of these devices. The results indicate that ejector systems can utilize the energy stored in the strong shock wave exiting the detonation tube to augment the impulse obtained from the detonation tube alone. Impulse augmentation ratios of up to 1.9 were achieved. The axial location of the converging-diverging ejectors relative to the end of the detonation tube were shown to affect the performance of the system.

Yungster, S.↗
DOE OSTI2016

Development of an FPGA-based multipoint laser pyroshock measurement system for explosive bolts

Pyroshock can cause failure to the objective of an aerospace structure by damaging its sensitive electronic equipment, which is responsible for performing decisive operations. A pyroshock is the high intensity shock wave that is generated when a pyrotechnic device is explosively triggered to separate, release, or activate structural subsystems of an aerospace architecture. Pyroshock measurement plays an important role in experimental simulations to understand the characteristics of pyroshock on the host structure. This paper presents a technology to measure a pyroshock wave at multiple points using laser Doppler vibrometers (LDVs). These LDVs detect the pyroshock wave generated due to an explosive-based pyrotechnical event. Field programmable gate array (FPGA) based data acquisition is used in the study to acquire pyroshock signals simultaneously from multiple channels. This paper describes the complete system design for multipoint pyroshock measurement. The firmware architecture for the implementation of multichannel data acquisition on an FPGA-based development board is also discussed. An experiment using explosive bolts was configured to test the reliability of the system. Pyroshock was generated using explosive excitation on a 22-mm-thick steel plate. Three LDVs were deployed to capture the pyroshock wave at different points. The pyroshocks captured were displayed as acceleration plots. The results showed that our system effectively captured the pyroshock wave with a peak-to-peak magnitude of 303 741 g. The contribution of this paper is a specialized architecture of firmware design programmed in FPGA for data acquisition of large amount of multichannel pyroshock data. The advantages of the developed system are the near-field, multipoint, non-contact, and remote measurement of a pyroshock wave, which is dangerous and expensive to produce in aerospace pyrotechnic tests.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗