Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Magnetized Plasma Expansion and its Interaction with a Plasma Stream

Expansion of magnetized plasma in the magnetic field of a solenoid is studied by means of simulations using a 3-dimensional hybrid code. The plasma expands against a high- density and slow plasma stream (PS). The expansion causes inflation of the magnetic field; near the solenoid the magnetic field variation with increasing distance (r) remains as B(alpha)r(sup -3), but at farther distances B(alpha)r(sup -p), where the exponent p is found in the range 0.5 approx. less than p approx. less than 1.2 forming a plateau in the magnetic field distribution B(r). At the start of injection of plasma from the ends of the solenoid, the PS interacts with the solenoid magnetic field and creates a bow shock at a distance where the Larmor radius (r(sub il)) of the PS ions in the solenoid magnetic field nearly equals the scale length (L) of B(r), that is, r(sub il) approx. L = (B(sup -1)dB L /dr)(sup -1). As the injected plasma accumulates in the solenoid magnetic field, it expands inflating the magnetic field and the bow shock moves outward. The speed of the expansion front and the shock progressively decreases and eventually a stand-off occurs when the PS dynamic pressure is eventually balanced by the magnetic and kinetic pressures of the expanding plasma and the inflating magnetic field. The inflating field shows wave-like behavior, with considerable structures in the field and current distributions.

Singh, Nagendra↗
NASA NTRS2004

Recent Enhancements to USM3D Unstructured Flow Solver for Unsteady Flows

The NASA USM3D unstructured flow solver is undergoing extensions to address dynamic flow problems in support of NASA and NAVAIR efforts to study the applicability of Computational Fluid Dynamics tools for the prediction of aircraft stability and control characteristics. The initial extensions reported herein include two second-order time stepping schemes, Detached-Eddy Simulation, and grid motion. This paper reports the initial code verification and validation assessment of the dynamic flow capabilities of USM3D. The cases considered are the classic inviscid shock-tube problem, low Reynolds number wake shedding from a NACA 0012 airfoil, high Reynolds number DES-based wake shedding from a 4-to-1 length-to-diameter cylinder, and forced pitch oscillation of a NACA 0012 airfoil with inviscid and turbulent flow.

Pandya, Mohagna J.↗
NASA NTRS2005

Large-Scale Experimental Planetary Science Meets Planetary Defense: Deorbiting an Asteroidal Satellite

Other than remote-sensing and spacecraft-derived data, the only information that exists regarding the physical and chemical properties of asteroids is that inferred through calculations, numerical simulations, extrapolation of experiments, and meteorite studies. Our understanding of the dynamics of accretion of planetesimals, collisional disruption of asteroids, and the macroscopic, shock-induced modification of the surfaces of such small objects is also, for the most part, founded on similar inferences. While considerable strides have been made in improving the state of asteroid science, too many unknowns remain to assert that we understand the parameters necessary for the more practical problem of deflecting an asteroid or asteroid pair on an Earth-intersecting trajectory. Many of these deficiencies could be reduced or eliminated by intentionally deorbiting an asteroidal satellite and monitoring the resulting collision between it and the primary asteroid, a capability that is well within the limitations of current technology.

Cintala, M. J.↗
NASA NTRS2003

Influence of Shock Wave on the Flutter Behavior of Fan Blades Investigated

Modern fan designs have blades with forward sweep; a lean, thin cross section; and a wide chord to improve performance and reduce noise. These geometric features coupled with the presence of a shock wave can lead to flutter instability. Flutter is a self-excited dynamic instability arising because of fluid-structure interaction, which causes the energy from the surrounding fluid to be extracted by the vibrating structure. An in-flight occurrence of flutter could be catastrophic and is a significant design issue for rotor blades in gas turbines. Understanding the flutter behavior and the influence of flow features on flutter will lead to a better and safer design. An aeroelastic analysis code, TURBO, has been developed and validated for flutter calculations at the NASA Glenn Research Center. The code has been used to understand the occurrence of flutter in a forward-swept fan design. The forward-swept fan, which consists of 22 inserted blades, encountered flutter during wind tunnel tests at part speed conditions.

Srivastava, Rakesh↗
NASA NTRS2009

Lunar Surface Potential Increases during Terrestrial Bow Shock Traversals

Since the Apollo era the electric potential of the Moon has been a subject of interest and debate. Deployed by three Apollo missions, Apollo 12, Apollo 14 and Apollo 15, the Suprathermal Ion Detector Experiment (SIDE) determined the sunlit lunar surface potential to be about +10 Volts using the energy spectra of lunar ionospheric thermal ions accelerated toward the Moon. We present an analysis of Apollo 14 SIDE "resonance" events that indicate the lunar surface potential increases when the Moon traverses the dawn bow shock. By analyzing Wind spacecraft crossings of the terrestrial bow shock at approximately this location and employing current balancing models of the lunar surface, we suggest causes for the increasing potential. Determining the origin of this phenomenon will improve our ability to predict the lunar surface potential in support of human exploration as well as provide models for the behavior of other airless bodies when they traverse similar features such as interplanetary shocks, both of which are goals of the NASA Lunar Science Institute's Dynamic Response of the Environment At the Moon (DREAM) team.

Collier, Michael R.↗
DOE OSTI2024

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

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

Time-resolved velocity and ion sound speed measurements from simultaneous bow shock imaging and inductive probe measurements

We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (“b-dot”) probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large (R M > 10), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach–Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the sound speed to infer the value of Z̄T e , where Z̄ is the average ionization and T e is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic (M S ~ 8), super-Alfvénic (M A ~ 2) aluminum plasma generated during the ablation stage of an exploding wire array on the Magpie generator (1.4 MA, 250 ns). The velocity and Z̄T e measurements agree well with the optical Thompson scattering measurements reported in the literature and with 3D resistive magnetohydrodynamic simulations in GORGON.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
DOE OSTI2021

Observation of Fundamental Mechanisms in Compression-Induced Phase Transformations Using Ultrafast X-ray Diffraction

As theoretically hypothesized for several decades in group IV transition metals, we have discovered a dynamically stabilized body-centered cubic (bcc) intermediate state in Zr under uniaxial loading at sub-nanosecond timescales. Under ultrafast shock wave compression, rather than the transformation from α-Zr to the more disordered hex-3 equilibrium ω-Zr phase, in its place we find the formation of a previously unobserved nonequilibrium bcc metastable intermediate. We probe the compression-induced phase transition pathway in zirconium using time-resolved sub-picosecond x-ray diffraction analysis at the Linac Coherent Light Source. We also present molecular dynamics simulations using a potential derived from first-principles methods which independently predict this intermediate phase under ultrafast shock conditions. In contrast with experiments on longer timescale (> 10 ns) where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time dependence in the kinetics of phase transformations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Direct observations of shock-induced melting in a porous solid using time-resolved x-ray diffraction

Here in this work we provide direct evidence of shock-induced melting and associated kinetics in a porous solid (aluminum powder) using time-resolved x-ray diffraction. Unambiguous evidence of melting in 50% porous aluminum (Al) powder samples, shocked to peak pressures between ~13 - 19 GPa, was provided by the broadening of the Debye-Scherrer ring corresponding to the (111) peak. Shocked Al powder did not melt completely in any of our experiments within the durations of measurement. Incomplete (partial) melting of the powder, even after several hundreds of nanoseconds of shock loading, provides insights into thermal transport with Al powder particles under high-pressure dynamic loading. Such insights are quite valuable for developing well-constrained melting models and thermodynamic equations of state for porous Al and other porous solids relevant to planetary and materials science.

36 MATERIALS SCIENCE↗
NASA NTRS1995

High-velocity, high-excitation neutral carbon in a cloud in the Vela supernova remnant

HD 72089 is situated behind the Vela supernova remnant, and the interstellar absorption lines in the spectrum of this star are remarkable for two reasons. First, there are six distinct velocity components that span the (heliocentric) velocity range -60 to +121 km/s in the lines of Na I and Ca II. Second, two of the components at high velocity, one at +85 km/s and another at +121.5 km/s, have densities that are large enough to produce observable lines from neutral carbon. The gas moving at +121.5 km/s has such a large pressure that the excited fine-structure levels of the ground electronic state of C I are collisionally populated nearly in proportion to their level degeneracies. This high-velocity gas exhibits unusually low column densities of Mg I and Na I, compared to that of C I. We propose that the +121.5 km/s component represents gas that has cooled and recombined in a zone that follows a shock driven into a cloud by the very recent passage of a supernova blast wave. A representative preshock density of n(sub H) approximately = 13/cc and velocity v(sub s) = 100 km/s is indicated by the strength of diffuse (O III) emission lines seen in directions very near HD 72089. The strong collisional population of excited C I and apparent absence of excited levels of O I give a most favorable fit to the conditions 1000 less than n(sub H) less than 2900/cc over a temperature range 300 less than T less than 1000 K. The fact that the compression is not substantially more than this indicates that the preshock gas may have had an embedded, transverse magnetic field with a strength B greater than or approximately = 1 micro-G. The large dynamical pressure of the supernova blast wave that would be needed to create the cloud shock that we describe implies that the energy of the supernova was 8 x 10(exp 51) ergs, if the Vela remnant is 500 pc away. We can bring this value much closer to typical supernova energies E less than or approximately = 10(exp 51) ergs if the distance to the remnant is revised downward by at least a factor of 2.

Jenkins, Edward B.↗
NASA NTRS1966

Dynamic Environment of the Ranger Spacecraft I Through IX (Final Report)

The dynamic environment of the Ranger spacecraft (I through IX)during the launch portion of flight is defined in this Report. Flight data from each of nine spacecraft launches have been reviewed and are included. The environments receiving emphasis herein include liftoff acoustics, liftoff and transonic vibration, and the transient vibrations of the various pyrotechnic and staging events in the launch sequence. The systems for data acquisition and analysis are briefly described. In addition, post-flight comparison of flight data and ground test specification levels is made, and the dynamic test program is briefly discussed.

SHOCK↗
NASA NTRS1981

On the complex state of the interplanetary medium of 28-29 July 1977

Observations of plasma and magnetic field variations in the near-Earth solar wind are discussed. Both a corotating stream and a driven shock are present. The driver gas seems to be enveloped in the rising speed phase of this stream; this appearance is attributed to a convoluted surface separating the two plasma domains. The magnetic field in the post shock flow (0030-1230 UT of July 29) has a large and geoeffective southward component at times; the energy coupling coefficient reaches approximately 5.4 x 10 to the 19th power ergs/s. In the driver gas (1230 UT of July 29 to 0110 of July 30) the magnetic field is dominantly northward. The density and dynamic pressure decrease by almost two orders of magnitude (100 to 2 cm/3) from just behind the interplanetary shock to approximately 3 hours into the driver gas flow. The dominant magnetic field variation in the driver gas is modeled by a cloud-like structure. Significant plasma parameter variations within the driver gas are attributed to structure in the parent solar mass ejection event and to interplanetary kinematics.

King, J. H.↗
NASA NTRS2008

Modeling of Fixed-Exit Porous Bleed Systems

A model has been developed to simulate a fixed-exit porous bleed system for supersonic inlets. The fixed-exit model allows the amount of bleed flow to vary according to local flow conditions and fixed-exit characteristics of the bleed system. This variation is important for the control of shock-wave/boundary-layer interactions within the inlet. The model computes the bleed plenum static pressure rather than requiring its specification. The model was implemented in the Wind-US computational fluid dynamics code. The model was then verified and validated against experimental data for bleed on a flat plate with and without an impinging oblique shock and for bleed in a Mach 3.0 axisymmetric, mixed-compression inlet. The model was able to accurately correlate the plenum pressures with bleed rates and simulate the effect of the bleed on the downstream boundary layer. Further, the model provided a realistic simulation of the initiation of inlet unstart. The results provide the most in-depth examination to date of bleed models for use in the simulation of supersonic inlets. The results also highlight the limitations of the models and aspects that require further research.

Slater, John W.↗
NASA NTRS1972

Digital-computer normal shock position and restart control of a Mach 2.5 axisymmetric mixed-compression inlet

Digital computer control of a mixed-compression inlet is discussed. The inlet was terminated with a choked orifice at the compressor face station to dynamically simulate a turbojet engine. Inlet diffuser exit airflow disturbances were used. A digital version of a previously tested analog control system was used for both normal shock and restart control. Digital computer algorithms were derived using z-transform and finite difference methods. Using a sample rate of 1000 samples per second, the digital normal shock and restart controls essentially duplicated the inlet analog computer control results. At a sample rate of 100 samples per second, the control system performed adequately but was less stable.

Neiner, G. H.↗
DOE OSTI2022

Atomistic deformation mechanism of silicon under laser-driven shock compression

Abstract Silicon (Si) is one of the most abundant elements on Earth, and it is the most widely used semiconductor. Despite extensive study, some properties of Si, such as its behaviour under dynamic compression, remain elusive. A detailed understanding of Si deformation is crucial for various fields, ranging from planetary science to materials design. Simulations suggest that in Si the shear stress generated during shock compression is released via a high-pressure phase transition, challenging the classical picture of relaxation via defect-mediated plasticity. However, direct evidence supporting either deformation mechanism remains elusive. Here, we use sub-picosecond, highly-monochromatic x-ray diffraction to study (100)-oriented single-crystal Si under laser-driven shock compression. We provide the first unambiguous, time-resolved picture of Si deformation at ultra-high strain rates, demonstrating the predicted shear release via phase transition. Our results resolve the longstanding controversy on silicon deformation and provide direct proof of strain rate-dependent deformation mechanisms in a non-metallic system.

36 MATERIALS SCIENCE↗