Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

The interstellar medium on the Gamma Cassiopeiae line of sight

The interstellar medium on the Gamma Cas line of sight is studied through the observations of O I, H2, H I, D I, and Ar I absorption features with the Copernicus satellite. By using the velocity structure of the line of sight previously determined through atomic nitrogen, it is demonstrated that these neutrals are located in the same physical regions and that the O I and Ar I abundances could be solar, on the average, relative to N I in all components. The 1066.660 A Ar I line profile is contaminated by an unidentified line, already detected in the spectrum of Zeta Pup. Furthermore, a peculiar behavior of the Ar I abundance in the different components is reported, possibly associated with high-velocity interstellar gas present on the line of sight. It is concluded that the Gamma Cas line of sight could intercept a dense cloud which has been engulfed and disrupted by a shock.

Ferlet, R.↗
DOE OSTI2021

Investigation of UV filaments and their applications

UV laser pulses at 266 nm with 170 ps duration and up to 300 mJ energy create UV filaments. The effects of different beam preparation scenarios on UV filamentation are discussed. Beam preparation by focusing in vacuum and launching into the atmosphere through an aerodynamic window establishes the existence of filaments as self-guided non-diffracting entities. In addition, focusing through the aerodynamic window converts the initial super-Gaussian into a parabolic beam profile at the edge of the nonlinear medium (air). The shock-wave generated by the UV filament is investigated via shadowgraphy techniques. The shock-wave is cylindrically symmetric, leaving a low density channel behind, which is well-sustained for 1 ms, as it propagates outward. This low density channel is a waveguide that can be used for guiding purposes. Applications of UV filaments including guiding, remote sensing, and drag reduction for supersonic aircrafts are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Algorithm and code development for unsteady three-dimensional Navier-Stokes equations

In the last two decades, there have been extensive developments in computational aerodynamics, which constitutes a major part of the general area of computational fluid dynamics. Such developments are essential to advance the understanding of the physics of complex flows, to complement expensive wind-tunnel tests, and to reduce the overall design cost of an aircraft, particularly in the area of aeroelasticity. Aeroelasticity plays an important role in the design and development of aircraft, particularly modern aircraft, which tend to be more flexible. Several phenomena that can be dangerous and limit the performance of an aircraft occur because of the interaction of the flow with flexible components. For example, an aircraft with highly swept wings may experience vortex-induced aeroelastic oscillations. Also, undesirable aeroelastic phenomena due to the presence and movement of shock waves occur in the transonic range. Aeroelastically critical phenomena, such as a low transonic flutter speed, have been known to occur through limited wind-tunnel tests and flight tests. Aeroelastic tests require extensive cost and risk. An aeroelastic wind-tunnel experiment is an order of magnitude more expensive than a parallel experiment involving only aerodynamics. By complementing the wind-tunnel experiments with numerical simulations the overall cost of the development of aircraft can be considerably reduced. In order to accurately compute aeroelastic phenomenon it is necessary to solve the unsteady Euler/Navier-Stokes equations simultaneously with the structural equations of motion. These equations accurately describe the flow phenomena for aeroelastic applications. At Ames a code, ENSAERO, is being developed for computing the unsteady aerodynamics and aeroelasticity of aircraft and it solves the Euler/Navier-Stokes equations. The purpose of this contract is to continue the algorithm enhancements of ENSAERO and to apply the code to complicated geometries. During the last year, the geometric capability of the code was extended to simulate transonic flows, a wing with oscillating control surface. Single-grid and zonal approaches were tested. For the zonal approach, a new interpolation technique was introduced. The key development of the algorithm was an interface treatment between moving zones for a control surface using the virtual-zone concept. The work performed during the period, 1 Apr. 1992 through 31 Mar. 1993 is summarized. Additional details on the various aspects of the study are given in the Appendices.

Obayashi, Shigeru↗
DOE OSTI2026

Efficacy of Gadolinium Gallium Garnet (GGG) as a High-Impedance Optical Window for Shock Wave Experiments

We conducted a series of plate impact experiments to examine the efficacy of < 111>-oriented gadolinium gallium garnet (GGG) single crystals as high-impedance optical window for Photonic Doppler velocimetry (PDV) under shock and double-shock loading. At ~ 123 GPa, shocked GGG remains fully transparent to 1550 nm light for at least 250 ns without any signal degradation. Above 135 GPa, PDV data measured through GGG exhibit a gradual loss of fringe contrast following shock entrance, which eventually leads to transparency loss. The duration for which shocked GGG remains transparent decreases with increasing pressure, and at ~ 148 GPa, it becomes opaque within ~ 20–30 ns. This limits the use of GGG as an interferometry window between ~ 110–140 GPa under single shock loading. Within this pressure range, the refractive index of GGG increases linearly with density: n = 1.552 + 0.054ρ. In contrast to single shock loading, where GGG becomes opaque rapidly above ~ 140 GPa, double-shocked GGG remains optically transparent for over 100 ns when it is first shocked to ~ 123 GPa and then reshocked to significantly higher pressures (215–233 GPa). Our findings raise the exciting possibility of GGG being used as a high-impedance optical window in multi-shock and shock-ramp loading experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Why Are There so Few Reports of High-Energy Electron Drift Resonances? Role of Radial Phase Space Density Gradients

Models of monochromatic Pc5 (2–7 mHz) ultralow frequency (ULF) wave interactions with high energy (greater than ~1 MeV) electrons predict drift resonant interactions that can cause rapid radial transport and acceleration. There are few reports of electron drift resonance at energies greater than ~1 MeV, in contrast to lower energies; moreover, all previous reports occur in the aftermath of interplanetary shocks. These two facts are difficult to reconcile with theory and numerical simulations predicting that greater than ~1 MeV drift resonances should occur more often and in a wider variety of driving conditions. In this study, we show that a combination of observational sampling biases and nominal radial phase space density gradients is one explanation for this discrepancy between theory and observations. In particular, we examine electron dynamics in two case studies with very similar satellite coverage, solar wind conditions, and Pc5 wave properties, yet with different radial phase space density profiles. Using global wave and particle observations, we show that the events have vastly different particle responses despite having similar wave properties. Placing these results in context with past studies, we further show here that nominal radial PSD gradients near geostationary orbit can mask the expected drift resonance particle response and explain (1) the small number of past greater than ~1 MeV drift resonance reports and (2) the restriction of these reports to interplanetary shock events. We argue that future observational studies characterizing radial transport via drift resonance should examine global particle dynamics, including observations of the radial phase space density profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic simulations verifying reconnection rates measured in the laboratory, spanning the ion-coupled to near electron-only regimes

The rate of reconnection characterizes how quickly flux and mass can move into and out of the reconnection region. In the Terrestrial Reconnection EXperiment (TREX), the rate at which the antiparallel asymmetric reconnection occurs is modulated by the presence of a shock and a region of flux pileup in the high-density inflow. Simulations utilizing a generalized Harris-sheet geometry have tentatively shown agreement with TREX's measured reconnection rate scaling relative to system size, which is indicative of the transition from ion-coupled toward electron-only reconnection. Here, in this study, we present simulations tailored to reproduce the specific TREX geometry, which confirm both the reconnection rate scale and the shock jump conditions previously characterized experimentally in TREX. The simulations also establish an interplay between the reconnection layer and the Alfvénic expansions of the background plasma associated with the energization of the TREX drive coils; this interplay has not yet been experimentally observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
NASA NTRS2004

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolution, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗
NASA NTRS2004

DSMC Simulations of Hypersonic Flows With Shock Interactions and Validation With Experiments

The capabilities of a relatively new direct simulation Monte Carlo (DSMC) code are examined for the problem of hypersonic laminar shock/shock and shock/boundary layer interactions, where boundary layer separation is an important feature of the flow. Flow about two model configurations is considered, where both configurations (a biconic and a hollow cylinder-flare) have recent published experimental measurements. The computations are made by using the DS2V code of Bird, a general two-dimensional/axisymmetric time accurate code that incorporates many of the advances in DSMC over the past decade. The current focus is on flows produced in ground-based facilities at Mach 12 and 16 test conditions with nitrogen as the test gas and the test models at zero incidence. Results presented highlight the sensitivity of the calculations to grid resolutions, sensitivity to physical modeling parameters, and comparison with experimental measurements. Information is provided concerning the flow structure and surface results for the extent of separation, heating, pressure, and skin friction.

Moss, James N.↗
NASA NTRS1991

The effect of a type 3 and type 4 shock/shock interaction on heat transfer in the stagnation region

One of the major engineering challenges in designing the National Aerospace Plane, NASP, is to overcome augmented heating on the intake cowl lip from shock/shock interactions. The shock/shock interaction arises when the bow shock from the craft's nose interferes with the bow shock from the cowl lip. Considering only the region immediately around the cowl lip, the problem geometry may be simplified as that of an oblique shock impinging on a bow shock from a circular cylinder. Edney classified six different interference patterns resulting from an oblique-shock/curved bow-shock interaction. Of these six types, type 3 and 4 are most significant in that augmented surface heat transfer may be ten to thirty times greater than the case without the shock/shock interaction. The objective was to begin to develop a mathematical model which is capable of predicting the effect of a type 3 and 4 shock/shock interaction in the stagnation region of an arbitrary 2-D body. This model must be capable of predicting the maximum surface heat flux and the surface stagnation point pressure once the outer (effectively inviscid) flowfield is given. Therefore, it must capture the unsteady physics of the impinging shear layer.

Wilson, Dennis↗
NASA NTRS1996

Studies of Interstellar Pickup Ions in the Solar Wind

The work under this grant involves studies of the interaction of interstellar pickup ions with the solar wind, with the goal of a comprehensive model of the particle distributions and wave intensities to be expected throughout the heliosphere, as well as the interactions of those distributions with the solar wind termination shock. In the past year, we have completed a number of projects, including observations and modeling of the effects of a large scattering mean free path on the pickup He(+) seen at AMPTE, an analytical model of anisotropic pickup tons in a steady radial magnetic field, and a derivation of a reduced solar wind Mach number due to increased estimates on the inflowing hydrogen density allowing for a weak termination shock. In the next year, we plan to investigate in more detail the correspondence between our models of anisotropic pickup ions and the data on spectra, variations, and proton-He(+) correlation provided by AMPTE, Ulysses, and our instrument on SOHO. We will model the time-dependent pickup ion density resulting from finite periods of radial magnetic field. We will also incorporate the effects of a large mean free path into our analysis of the He(+) focusing cone, leading to more accurate parameter values for the interstellar helium gas. This progress report also includes a discussion of our Space Physics Educational Outreach activities in the past year and plans for the next year.

Isenberg, Philip A.↗
DOE OSTI2021

Shock to detonation transition of pentaerythritol tetranitrate (PETN) initially pressed to 1.65 g/cm 3

A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. In this work, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.

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

Multi-shock design tools in double shell models

The essence of the project was to implement a program in python that demonstrated certain models put forth by fusion physicists. Across the time spent at Los Alamos National Laboratory, much was spent programming. Though fundamentally programmatical in nature, the task contained considerable difficulty hid within mathematically deriving and manipulating the models while maintaining data discretization. The result of solving these complications gave invaluable experience in applied mathematics. The final product, a piece of lengthy code, ended up having the full range of desired capabilities, with plans for more. Given the laser power history from a pulse shot on NIF, Omega, or any other Inertial Confinement Fusion facility, the code maps the terrain of radiation temperature inside the hohlraum throughout the period of the shot. Additionally, shock wave convergence prediction was implemented as additional capability. As the name of the project suggests, this tool acts as a compass of direction for researchers designing laser power profiles. Further development of this tool will undoubtedly assist in the research conducted on single and double-shell indirect ICF, consequently expanding the capabilities of operating stewardship over the national nuclear stockpile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
DOE OSTI2020

Characterizing Dislocation Defects in Dark Field X-ray Microscopy Images of Bulk Aluminum

Material defects play a large role in material response under shock loading, yet the process of how these defects initiate, propagate, and annihilate is not well understood. Using a newly developed diagnostic, dark field X-ray microscopy (DFXM), we can now visualize the behavior of dislocation defects in materials at the mesoscale under varying conditions. Using data from the European Synchrotron Radiation Facility, we apply a variety of image processing techniques to capture relevant features to locate and characterize size and orientation of dislocation defects in bulk aluminum DFXM images. Beyond simple visualization, this analysis drives statistical characterization of the defects and their dynamic response to further improve the relevant physics models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Optimization Methodology of Polar Direct-Drive Illumination for the National Ignition Facility

Improved laser illumination uniformity drives shocks and implosions to create more extreme high energy density environments. Predominantly, the geometry of experiments that can be performed is dictated by the layout of beams at laser facilities, limiting inter-facility and multiscale investigations. This Letter presents the first automated, algorithmic approach for generating illumination configurations for high energy density experiments. The method is demonstrated in comparison to a polar direct drive solid target experiment at the National Ignition Facility. The new illumination configuration is simulated to create greater than ×3 higher peak pressure and almost ×2 higher density by maintaining better shock uniformity. Furthermore, the optimization process is performed with reduced computational expense and isotropic plasma profiles, while accounting for the impact of cross-beam energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗