Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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 1,081 records · Page 60

Differentiable lagrangian shock hydrodynamics with application to stable shock acceleration of density interfaces

We develop a gradient based optimization approach for the equations of compressible, Lagrangian hydrodynamics and demonstrate how it can be employed to automatically uncover strategies to control hydrodynamic instabilities arising from shock acceleration of density interfaces. Strategies for controlling the Richtmyer-Meshkov instability (RMI) are of great benefit for inertial confinement fusion (ICF) where shock interactions with many small imperfections in the density interface lead to instabilities which rapidly grow over time. These instabilities lead to mixing which, in the case of laser driven ICF, quenches the runaway fusion process ruining the potential for positive energy return. Here, we demonstrate that control of these instabilities can be achieved by optimization of initial conditions with ( > 100) parameters. Optimizing over a large parameter space like this is not possible with gradient-free optimization strategies. This requires computation of the gradient of the outputs of a numerical solution to the equations of Lagrangian hydrodynamics with respect to the inputs. We show that the efficient computation of these gradients is made possible via a judicious application of (i) adjoint methods, the exact formal representation of sensitivities involving partial differential equations, and (ii) automatic differentiation (AD), the algorithmic calculation of derivatives of functions. Careful regularization of multiple operators including artificial viscosity and timestep control is required. We perform design optimization of > 100 parameter energy field driving the Richtmyer Meshkov instability showing significant suppression while simultaneously enhancing the acceleration of the interface relative to a nominal baseline case.

Hydrophysics↗

Deep Koopman Neural Network for Analyzing High-Energy-Density Simulations of Electrical Wire Explosions

Megaampere-scale electrical wire experiments (EWEs) provide a platform for studying magnetohydrodynamic (MHD) instability growth in magneto-inertial fusion (MIF) devices. Even when nonlinear simulations of these experiments can digitally reproduce much of the experimentally observed instability growth, interpreting the results and understanding mode growth and evolution can be non-trivial. As a first step toward providing better interpretation of these simulation features, this work investigates the use of a deep neural network that uses Koopman operator theory to analyze the dynamics of pulsed-power-driven explosions of EWEs. This deep neural network is trained on 1-D resistive MHD simulations of EWEs. This neural network learns to transform the nonlinear data into a lower-dimensional representation where the time dynamics are linear. Layers of this neural network are shown to learn features of the simulations, including the locations of shock waves and different physical regimes of the simulation. Using the learned features, the network can compress a time state of the simulation consisting of 5120 data point into a 36-parameter lower-dimensional latent space embedding. Furthermore, these embeddings are shown to be clustered in the latent space by initial radius and time state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis of Si-Fe Chondrule-like Dust Analogues in RF Discharge Plasmas

Chondrules are tiny particles that occur in stony meteorites and are considered as the building blocks of early asteroids and planets. It is believed that they were formed by the fast heating of the dust in the solar nebula. To date, there is no lab-scale experimental study of the formation of chondrules from the initial gas phase precursors following fast heating and crystallisation. The motivation of this work is a pre-trial study of the formation of chnodrule-like particles. The formation of meteorites in the space environment is associated with the aggregation of small particles or molecular clouds under the influence of shock waves or high-energy gas discharges in the solar nebula. In this work, the properties of product formation at the nanoscale-level were investigated using different feedstock materials which are the dominant elements in the meteorite. The structural and morphological properties of the synthesised Si-Fe nanomaterials were analysed by scanning/transmission electron microscopy (SEM/TEM), and chemical composition was analysed by X-ray energy-dispersive spectroscopy (EDS). The identification of crystalline phases was carried out by X-ray diffraction (XRD), whereas the presence of an Fe-Si system in the synthesised particles was demonstrated by Mössbauer spectroscopy. The obtained materials were exposed to the relatively high-energy pulsed plasma beam on the substrate with the aim to emulate the possible fast heating and melting of the formed nanoparticles. The formation steps of growing synthetic (engineered) chondro-like particles and nanostructures in laboratory conditions is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A viscous starter solution for shuttle flow field computations

A time dependent numerical procedure that calculates the viscous flow in the shock layer around the shuttle nose has been developed. It has been applied to the free stream conditions that may be encountered in the shuttle flight trajectory and found to be efficient and in satisfactory agreement with other steady and unsteady techniques. Two features associated with the present formulation are responsible for its successful applications. The first is that the flow field computation is made within the shock layer; the second is that the mapping of the computational region places proportionately more mesh points in the vicinity of the wall than in the remaining computational region. The numerical results may be used for low altitude flight to determine the edge conditions for a thin boundary layer, to provide initial boundary layer profiles for downstream boundary layer calculations, and to start supersonic computations. For high altitude flight, it may be used to determine the complete flow field in the nose region and also provide a starter solution for subsequent supersonic flow computations.

Li, C. P.↗

Shock-turbulence interaction and the generation of noise

Interaction of convected field of turbulence with shock wave is analyzed to yield modified turbulence, entropy spottiness, and noise generated downstream of the shock. Analysis is generalization of single-spectrum-wave treatment of NACA-TN-2864. Formulas for spectra and correlations are obtained. Numerical calculations yield curves of rms velocity components, temperature, pressure, and noise in db against Mach number for m = 1 to infinity; both isotropic and strongly axisymmetric (lateral/longitudinal = 36/1) initial turbulence are treated. In either case, turbulence of 0.1 percent longitudinal component generates about 120 dbs of noise.

Ribner, H S↗

Relaxation from Steady States Far from Equilibrium and the Persistence of Anomalous Shock Behavior in Weakly Ionized Gases

The decay of anomalous effects on shock waves in weakly ionized gases following plasma generator extinction has been measured in the anticipation that the decay time must correlate well with the relaxation time of the mechanism responsible for the anomalous effects. When the relaxation times cannot be measured directly, they are inferred theoretically, usually assuming that the initial state is nearly in thermal equilibrium. In this paper, it is demonstrated that relaxation from any steady state far from equilibrium, including the state of a weakly ionized gas, can proceed much more slowly than arguments based on relaxation from near equilibrium states might suggest. This result justifies a more careful analysis of the relaxation times in weakly ionized gases and suggests that although the experimental measurements of relaxation times did not lead to an unambiguous conclusion, this approach to understanding the anomalous effects may warrant further investigation.

Rubinstein, Robert↗

Application of Time-resolved PIV to Supersonic Hot Jets

This presentation lays out the ground-breaking work at bringing high-speed (25kHz) particle image velocimetry (PIV) to bear on measurements of noise-producing turbulence in hot jets. The work is still in progress in that the tremendous amount of data obtained are still be analyzed, but the method has been validated and initial results of interest to jet noise modeling have been obtained. After a brief demonstration of the validation process used on the data, results are shown for hot jets at different temperatures and Mach numbers. Comparisons of first order statistics show the relative indifference of the turbulence to the presence of shocks and independence to jet temperature. What does come out is that when the shock-containing jets are in a screech mode the turbulence is highly elevated, showing the importance of removing screech phenomena from model-scale jets before applying findings to full-scale aircraft which typically do not contain shocks.

Bridges, James↗

Progress on Transonic Flutter and Shock Buffet Computationsin Support of the Third Aeroelastic Prediction Workshop

This paper reports on the progress of the NASA Langley team contributions to the third Aeroelastic Pre-diction Workshop’s (AePW-3) High Angle Working Group (HAWG). The primary objectives of HAWG is to predict the fluter dynamic pressure of the NASA Benchmark Supercritical Wing (BSCW) configuration at Mach 0.8 and 5◦angle of attack. The secondary objective is to determine if a shock-buffet onset is present at or near that flow condition. The computational results are obtained using FUN3D, an unstructured grid Reynolds-averaged Navier-Stokes solver developed at the NASA Langley Research Center. The preliminary analysis results show a computationally-obtained flutter dynamic pressure of approximately 120 psf. Initial results describing unforced BSCW unsteady flow environment at flutter condition are also presented.

Pawel Chwalowski↗

Synchrotron radiography of Richtmyer–Meshkov instability driven by exploding wire arrays

We present a new technique for the investigation of shock-driven hydrodynamic phenomena in gases, liquids, and solids in arbitrary geometries. The technique consists of a pulsed power-driven resistive wire array explosion in combination with multi-MHz synchrotron radiography. Compared to commonly used techniques, it offers multiple advantages: (1) the shockwave geometry can be shaped to the requirements of the experiment, (2) the pressure (P > 300 MPa) generated by the exploding wires enables the use of liquid and solid hydrodynamic targets with well-characterized initial conditions (ICs), (3) the multi-MHz radiography enables data acquisition to occur within a single experiment, eliminating uncertainties regarding repeatability of the ICs and subsequent dynamics, and (4) the radiographic measurements enable estimation of compression ratios from the x-ray attenuation. In addition, the use of a synchrotron x-ray source allows the hydrodynamic samples to be volumetrically characterized at a high spatial resolution with synchrotron-based microtomography. This experimental technique is demonstrated by performing a planar Richtmyer–Meshkov instability (RMI) experiment on an aerogel–water interface characterized by Atwood number A0∼−0.8 and Mach number M∼1.5. The qualitative and quantitative features of the experiment are discussed, including the energy deposition into the exploding wires, shockwave generation, compression of the interface, startup phase of the instability, and asymptotic growth consistent with Richtmyer's impulsive theory. Additional effects unique to liquids and solids—such as cavitation bubbles caused by rarefaction flows or initial jetting due to small perturbations—are observed. It is also demonstrated that the technique is not shape dependent by driving a cylindrically convergent RMI experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plant heat stress: Concepts directing future research

Predicted increases in future global temperatures require us to better understand the dimensions of heat stress experienced by plants. In this review we highlight four key areas for improving our approach towards understanding plant heat stress responses. First, although the term "heat stress" is broadly used, that term encompasses heat shock, heat wave, and warming experiments, which vary in the duration and magnitude of temperature increase imposed. A greater integration of results and tools across these approaches is needed to better understand how heat stress associated with global warming will affect plants. Secondly, there is a growing need to associate plant responses to tissue temperatures. We review how plant energy budgets determine tissue temperature and discuss the implications of using leaf versus air temperature for heat stress studies. Third, we need to better understand how heat stress affects reproduction, particularly understudied stages such as floral meristem initiation and development. Fourth, we emphasize the need to integrate heat-stress recovery into breeding programs to complement recent progress in improving plant heat-stress tolerance. Taken together, we provide insights into key research gaps in plant heat stress and provide suggestions on addressing these gaps to enhance heat stress resilience in plants.

54 ENVIRONMENTAL SCIENCES↗

Initial observations of plasma electrons from the Pioneer 10 flyby of Jupiter

Initial results are presented from the electron measurements made by the Ames Research Center Plasma Analyzer during the inbound passage of Pioneer 10 in the vicinity of Jupiter. The observations indicate that as in the case of the earth's magnetosheath, there is an increase in electron temperature across the Jovian bow shock. During the second extended magnetosheath traversal (approximately 54 to 46.5 Jupiter radii) the electron temperatures were generally higher than those observed during the first magnetosheath traversal (109 to 96 Jupiter radii). These higher electron temperatures are in agreement with the measured higher magnitude of the magnetic field and ion density during this traversal. These observations are consistent with the contraction of the Jovian magnetosphere due to an increase in the solar wind dynamic pressure.

Intriligator, D. S.↗

Interstellar clouds - From a dynamical perspective on their chemistry

The possibility is examined that in the course of its dynamical evolution, a single mass of interstellar gas would exhibit properties of diffuse clouds, dense clouds and finally also of clouds perturbed by shocks or intense UV or X-ray radiation generated by a star of its own creation. This concept provides a common thread through the bewildering diversity of physical and chemical compositional properties shown by interstellar clouds. From this perspective, instead of being static objects, interstellar clouds are possibly incessantly evolving from initially diffuse to later dense state and then to star formation which ultimately restructures or disperses the remaining cloud material to begin the whole evolutionary process once again. Based on a simplified study of interstellar chemistry from a dynamical perspective, the ideas are presented as an heuristic: to encourage thought on the future direction of molecular astrophysics and the need to consider the chemical behavior of interstellar clouds in conjunction with, rather than in isolation from, their dynamical behavior. A physical basis must be sought for the semiempirical temperature formula which has been given a critical role in the collapse of diffuse clouds. Self-shielding effects in the chemistry of CO were neglected and this drawback should be removed; the ability of the model to explain the fractional abundances of more complex molecules, such as cyanopolyynes, should be examined.

Prasad, S. S.↗

Development and verification of NASA Standard Initiator-2 (NSI-2)

The results of verification testing under NASA Contract NAS9-17249 are summarized in this Final Report. An Appendix is published separately. The report contains a brief discussion of design elements, test results, and a discussion of the latter. The contract was awarded for development and demonstration of design modifications to the NASA Standard Initiator-1 (NSI-1), resulting in the development of NSI-2. The intent of the project was to develop and prove the following design features: (1) Protection from premature initiation due to electrostatic potential by means other than a spark gap; (2) Development of an integral insulated charge cup for use with a zirconium/potassium perchlorate type mixture; (3) Improved cost and reliability figures by the use of a deposited ridge; and (4) Elimination of dudding as a result of the pyrotechnic shock effect, the movement of the primer composition away from the bridge circuit.

Corbin-Roberts, Linda↗

Voyager and Pioneer Missions to the boundaries of the heliosphere

The interaction of the solar wind with the interstellar medium is expected to result in a complex, probably dynamic outer heliospheric boundary region. In increasing distance from the sun the boundary region includes the solar wind termination shock, the heliopause, and perhaps a heliobowshock. The continuing missions of Voyager 1 and 2 and Pioneer 10 and 11 provide a unique opportunity to make in situ particle and field observations of the boundaries of the heliosphere and associated phenomena. Observations already made by these spacecraft suggest that the termination shock may be located 60-100 AU from the sun. Voyager 1 will reach 60 AU in 1995 and 100 AU in 2006. NASA's Space Physics Division is planning with the Voyager and Pioneer Projects to ensure timely and maximum scientific return from spacecraft encounters with the solar wind termination region. Initial emphasis is on being able to reliably determine when an encounter is imminent, so that spacecraft tracking and spacecraft and instrument configurations can be changed to encounter modes. It is anticipated that multiple encounters with heliospheric boundaries will occur due to motion of the boundaries in response to large variations in the solar wind pressure.

Pesses, M. E.↗

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗

On the instability of hypersonic flow past a flat plate

The instability of hypersonic boundary-layer flows over flat plates is considered. The viscosity of the fluid is taken to be governed by Sutherland's law, which gives a much more accurate representation of the temperature dependence of fluid viscosity at hypersonic speeds than Chapman's approximate linear law; although at lower speeds the temperature variation of the mean state is less pronounced so that the Chapman law can be used with some confidence. Attention is focussed on the so-called (vorticity) mode of instability of the viscous hypersonic boundary layer. This is thought to be the fastest growing inviscid disturbance at hypersonic speeds; it is also believed to have an asymptotically larger growth rate than any viscous or centrifugal instability. As a starting point the instability of the hypersonic boundary layer which exists far downstream from the leading edge of the plate is investigated. In this regime the shock that is attached to the leading edge of the plate plays no role, so that the basic boundary layer is non-interactive. It is shown that the vorticity mode of instability of this flow operates on a significantly different lengthscale than that obtained if a Chapman viscosity law is assumed. In particular, it is found that the growth rate predicted by a linear viscosity law overestimates the size of the growth rate by O(M(exp 2). Next, the development of the vorticity mode as the wavenumber decreases is described, and it is shown that acoustic modes emerge when the wavenumber has decreased from it's O(1) initial value to O(M (exp -3/2). Finally, the inviscid instability of the boundary layer near the leading edge in the interaction zone is discussed and particular attention is focussed on the strong interaction region which occurs sufficiently close to the leading edge. It is found that the vorticity mode in this regime is again unstable, and that it is concentrated in the transition layer at the edge of the boundary layer where the temperature adjusts from its large, O(M(exp 2), value in the viscous boundary layer, to its O(1) free stream value. The existence of the shock indirectly, but significantly, influences the instability problem by modifying the basic flow structure in this layer.

Blackaby, Nicholas↗

Measuring the effectiveness of high-performance Co-Optima biofuels on suppressing soot formation at high temperature

Soot emissions in combustion are unwanted consequences of burning hydrocarbon fuels. The presence of soot during and following combustion processes is an indication of incomplete combustion and has several negative consequences including the emission of harmful particulates and increased operational costs. Efforts have been made to reduce soot production in combustion engines through utilizing oxygenated biofuels in lieu of traditional nonoxygenated feedstocks. The ongoing Co-Optimization of Fuels and Engines (Co-Optima) initiative from the US Department of Energy (DOE) is focused on accelerating the introduction of affordable, scalable, and sustainable biofuels and high-efficiency, low-emission vehicle engines. The Co-Optima program has identified a handful of biofuel compounds from a list of thousands of potential candidates. In this study, a shock tube was used to evaluate the performance of soot reduction of five high-performance biofuels downselected by the Co-Optima program. Current experiments were performed at test conditions between 1,700 and 2,100 K and 4 and 4.7 atm using shock tube and ultrafast, time-resolve laser absorption diagnostic techniques. The combination of shock heating and nonintrusive laser detection provides a state-of-the-art test platform for high-temperature soot formation under engine conditions. Soot reduction was found in ethanol, cyclopentanone, and methyl acetate; conversely, an α-diisobutylene and methyl furan produced more soot compared to the baseline over longer test times. For each biofuel, several reaction pathways that lead towards soot production were identified. The data collected in these experiments are valuable information for the future of renewable biofuel development and their applicability in engines.

09 BIOMASS FUELS↗

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Bui, Trong T.↗