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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2023

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Time-dependent evolution of cosmic-ray-modified shock structure: Transition to steady state

Steady state solutions to the two-fluid equations of cosmic-ray-modified shock structure were investigated first by Drury and Volk (1981). Their analysis revealed, among other properties, that there exist regions of upstream parameter space where the equations possess three different downstream solutions for a given upstream state. In this paper we investigate whether or not all these solutions can occur as time-asymptotic states in a physically realistic evolution. To do this, we investigate the time-dependent evolution of the two-fluid cosmic-ray equations in going from a specified initial condition to a steady state. Our results indicate that the time-asymptotic solution is strictly single-valued, and it undergoes a transition from weakly to strongly cosmic-ray-modified at a critical value of the upstream cosmic ray energy density. The expansion of supernova remnant shocks is considered as an example, and it is shown that the strong to weak transition is in fact more likely. The third intermediate solution is shown to influence the time-dependent evolution of the shock, but it is not found to be a stable time-asymptotic state. Timescales for convergence to these states and their implications for the efficiency of shock acceleration are considered. We also investigate the effects of a recently introduced model for the injection of seed particles into the shock accelerated cosmic-ray population. The injection is found to result in a more strongly cosmic-ray-dominated shock, which supports our conclusion that for most classes of intermediate and strong cosmic-ray-modified shocks, the downstream cosmic-ray pressure component is at least as large as the thermal gas pressure, independent of the upstream state. As a result, cosmic rays almost always play a significant role in determining the shock structure and dissipation and they cannot be regarded as test particles.

Donohue, D. J.↗

DOE OSTI2024

Impacts of aluminum-cerium alloy deployment on different cerium commodities

An imbalance problem exists in the rare earth element (REE) market with cerium and lanthanum being more abundant and less valuable than REEs used in permanent magnets such as neodymium, praseodymium, and dysprosium. Research has been conducted on new applications for abundant REEs to improve mining economics. However, potential market impacts of these applications on individual REE commodities are unknown. A system dynamics model was built to investigate the impacts of a new aluminum-cerium (Al-Ce) alloy that can theoretically utilize three separate cerium forms: cerium carbonate, cerium oxide, or cerium metal. Cerium carbonate currently possesses the lowest price of the three available forms making it initially the most favorable choice. However, model results showed that cerium oxide could provide the best economics for future Al-Ce alloy deployment because its larger market size buffers major shocks to price.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Effect on compression of lowering the design adiabat in the SQ-n campaign

A set of experiments in the SQ-n campaign attempts to increase the compression in layered implosions by reducing the design adiabat from 3.0 down to ∼ 2.6. The experiments utilize nanocrystalline high density carbon ablators driven by a ramped laser pulse drive. Contrary to expectations, the lower design adiabat layered implosion reaches the same level of compression, within the measurement uncertainty, as the previous set of SQ-n layered implosions at the higher design adiabat. Here, we examine potential explanations for this result, including incomplete melting or refreezing of the ablator after the first shock passes into the fuel layer, reproducibility, and dopant areal density. We also discuss additional potential paths to reaching lower adiabat including increased dopant areal density, ramped dopant profiles, and increased initial capsule size which are also based on experimental campaign results from the National Ignition Facility.

Physics - Plasma physics↗

DOE OSTI2025

Interplays In The Rare Earth Element Value Chain

Severson, Michael [Idaho National Laboratory (INL)↗

NASA NTRS1988

Impact-induced devolatilization and hydrogen isotopic fractionation of serpentine: Implications for planetary accretion

Impact-induced devolatilization of porous serpentine was investigated using two independent experimental methods, the gas recovery and the solid recovery method, each yielding nearly identical results. For shock pressures near incipient devolatilization, the hydrogen isotopic composition of the evolved H2O is very close to that of the starting material. For shock pressures at which up to 12 percent impact-induced devolatilization occurs, the bulk evolved gas is significantly lower in deuterium than the starting material. There is also significant reduction of H2O to H2 in gases recovered at these higher shock pressures, probably caused by reaction of evolved H2O with the metal gas recovery fixture. Gaseous H2O-H2 isotopic fractionation suggests high temperature isotopic equilibrium between the gaseous species, indicating initiation of devolatilization at sites of greater than average energy deposition. Bulk gas-residual solid isotopic fractionations indicate nonequilibrium, kinetic control of gas-solid isotopic ratios. Impact-induced hydrogen isotopic fractionation of hydrous silicates during accretion can strongly affect the long-term planetary isotopic ratios of planetary bodies, leaving the interiors enriched in deuterium. Depending on the model used for extrapolation of the isotopic fractionation to devolatilization fractions greater than those investigated experimentally can result from this process.

Tyburczy, James A.↗

NASA NTRS1993

When clusters collide - A numerical Hydro/N-body simulation of merging galaxy clusters

A 3D numerical simulation of two merging clusters of galaxies, using a hybrid Hydro/N-body code, is presented. The hydrodynamics of the code is solved by an Eulerian finite difference method. Initial results disclose that the X-ray emission of the dominant cluster becomes elongated and broadened; heating occurs at the core of the dominant cluster as a result of multiple shocks, and high velocity gas motions within the intracluster medium. It is predicted that clusters which have undergone recent mergers and do not have cooling flows will have high peculiar gas velocities and that the shocks and turbulence generated during the merger may power cluster-wide radio halos. Prolonged high-velocity gas motions through the dominant cluster core possibly play a major role in the formation and shaping of wide-angle tailed radio sources associated with central dominant galaxies. The N-body component of the simulation reveals the subcluster to be dispersed as it passes through the dominant cluster.

Roettiger, Kurt↗

NASA NTRS1995

The solar/interplanetary event of 14 April 1994 observed by Yohkoh/SXT

The polar crown event of April 14 1994 is one of the largest scale eruptive events observed by the Yohkoh/SXT. Associated with the formation of an arcade of soft X-ray loops at the Sun was the detection of an interplanetary forward/reverse shock event by the Ulysses spacecraft some 4-7 days later. The relationship between the coronal and interplanetary signatures of these events is important if we are to address fully the initialization and consequent acceleration of interplanetary phenomena, such as CMEs and counter-streaming electrons, originating at the Sun. From detailed analysis of the energetics of the arcade formed during the eruption of April 14 1994, we find peak temperatures and emission measures of approximately 5MK and approximately 10(exp 48)cm(exp -3) respectively. The total thermal content of the arcade loop structure observed in soft X-rays is calculated to be some 5 x 10(exp 29) ergs. The development of these parameters as the event proceeds and their relationship to the dynamics of the eruption are investigated. Although spanning a longitudinal range of some 150 degrees the April 14 event displayed the typical helmet streamer structure normally associated with coronal mass ejections These helmet streamers are thought to be related to the global solar magnetic field through the heliospheric current sheet (HCS). The arcade formation, together with the eruption of material into interplanetary space, signifies a large-scale reconfiguration of the coronal magnetic field. We examine the effects of the formation of such a coronal arcade structure on the HCS and discuss the dynamics involved with the passage of a large scale disturbance through the interplanetary magnetic field.

Alexander, D.↗

NASA NTRS2004

Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Pierazzo, E.↗

NASA NTRS2003

Unsteady Flowfield in a High-Pressure Turbine Modeled by TURBO

Forced response, or resonant vibrations, in turbomachinery components can cause blades to crack or fail because of the large vibratory blade stresses and subsequent high-cycle fatigue. Forced-response vibrations occur when turbomachinery blades are subjected to periodic excitation at a frequency close to their natural frequency. Rotor blades in a turbine are constantly subjected to periodic excitations when they pass through the spatially nonuniform flowfield created by upstream vanes. Accurate numerical prediction of the unsteady aerodynamics phenomena that cause forced-response vibrations can lead to an improved understanding of the problem and offer potential approaches to reduce or eliminate specific forced-response problems. The objective of the current work was to validate an unsteady aerodynamics code (named TURBO) for the modeling of the unsteady blade row interactions that can cause forced response vibrations. The three-dimensional, unsteady, multi-blade-row, Reynolds-averaged Navier-Stokes turbomachinery code named TURBO was used to model a high-pressure turbine stage for which benchmark data were recently acquired under a NASA contract by researchers at the Ohio State University. The test article was an initial design for a high-pressure turbine stage that experienced forced-response vibrations which were eliminated by increasing the axial gap. The data, acquired in a short duration or shock tunnel test facility, included unsteady blade surface pressures and vibratory strains.

Bakhle, Milind A.↗

NASA NTRS2008

High-Order Energy Stable WENO Schemes

A new third-order Energy Stable Weighted Essentially NonOscillatory (ESWENO) finite difference scheme for scalar and vector linear hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. A rigorous truncation error analysis is presented showing that the accuracy of the 3rd-order ESWENO scheme is drastically improved if the tuning parameters of the weight functions satisfy certain criteria. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO finite difference scheme of Jiang and Shu in terms of accuracy, while providing essentially nonoscillatory solutions near strong discontinuities.

Yamaleev, Nail K.↗

NASA NTRS2000

Experiments on the Richtmyer-Meshkov Instability of Incompressible Fluids

Richtmyer-Meshkov (R-M) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and possibly become turbulent given the proper initial conditions. R-M instability is similar to the Rayleigh-Taylor (R-T) instability, which is generated when the two fluids undergo a constant acceleration. R-M instability is a fundamental fluid instability that is important to fields ranging from astrophysics to high-speed combustion. For example, R-M instability is currently the limiting factor in achieving a net positive yield with inertial confinement fusion. The experiments described here utilize a novel technique that circumvents many of the experimental difficulties previously limiting the study of the R-M instability. A Plexiglas tank contains two unequal density liquids and is gently oscillated horizontally to produce a controlled initial fluid interface shape. The tank is mounted to a sled on a high speed, low friction linear rail system, constraining the main motion to the vertical direction. The sled is released from an initial height and falls vertically until it bounces off of a movable spring, imparting an impulsive acceleration in the upward direction. As the sled travels up and down the rails, the spring retracts out of the way, allowing the instability to evolve in free-fall until impacting a shock absorber at the end of the rails. The impulsive acceleration provided to the system is measured by a piezoelectric accelerometer mounted on the tank, and a capacitive accelerometer measures the low-level drag of the bearings. Planar Laser-Induced Fluorescence is used for flow visualization, which uses an Argon ion laser to illuminate the flow and a CCD camera, mounted to the sled, to capture images of the interface. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The amplitude of the disturbance during the experiment is measured and compared to theory. The results show good agreement (within 10%) with linear stability theory up to nondimensional amplitude ka = 0.7 (wavenumber x amplitude). These results hold true for an initial ka (before acceleration) of -0.7 less than ka less than -0.06, while the linear theory was developed for absolute value of ka much less than 1. In addition, a third order weakly nonlinear perturbation theory is shown to be accurate for amplitudes as large as ka = 1.3, even though the interface becomes double-valued at ka = 1.1. As time progresses, the vorticity on the interface concentrates, and the interface spirals around the alternating sign vortex centers to form a mushroom pattern. At higher Reynolds Number (based on circulation), an instability of the vortex cores has been observed. While time limitations of the apparatus prevent determination of a critical Reynolds Number, the lowest Reynolds Number this vortex instability has been observed at is 5000.

Jacobs, J.↗

DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗

DOE OSTI2020

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗

NASA NTRS1981

Plasma waves near Saturn: Initial results from Voyager 1

The Voyager 1 encounter with Saturn provided the first opportunity to investigate plasma wave interactions in the magnetosphere of Saturn. An overview of the principal results from the Voyager 1 plasma wave instrument is presented starting with the initial detection of Saturn and ending about four weeks after closest approach. A survey plot of the electric field intensities detected during the Saturn encounter is shown starting shortly before the inbound shock crossing and ending shortly after the outbound magnetopause crossing. Many intense waves were observed in the vicinity of Saturn. To provide a framework for presenting the observations, the results are discussed more or less according to the sequence in which the data were obtained.

Gurnett, D. A.↗

NASA NTRS1981

Surface photometry of spiral galaxies. I - Theoretical color variation and surface brightness across spiral arms

A method for the construction of color variations and surface brightness across a galactic spiral arm, using the concepts and results of the density-wave theory, is developed. The color variation measured by the index Q varies with the ratio of the light of young stars to that of the old stars, the initial mass function, and the pattern speed. The fact that the variation of brightness and color across a spiral arm is relatively smooth and symmetrical and not sharp as in the galactic shock picture is due to the smearing-out effects of star formation inside a cloud in the regions behind the galactic shock, and the tendency of new stars to fall back to small galactocentric distances after their formation because of their lower angular velocity at birth. Both effects arise naturally in the context of the density-wave theory.

Yuan, C.↗

NASA NTRS1992

Geomechanical models of impact cratering: Puchezh-Katunki structure

Impact cratering is a complex natural phenomenon that involves various physical and mechanical processes. Simulating these processes may be improved using the data obtained during the deep drilling at the central mound of the Puchezh-Katunki impact structure. A research deep drillhole (named Vorotilovskaya) has been drilled in the Puchezh-Katunki impact structure (European Russia, 57 deg 06 min N, 43 deg 35 min E). The age of the structure is estimated at about 180 to 200 m.y. The initial rim crater diameter is estimated at about 40 km. The central uplift is composed of large blocks of crystalline basement rocks. Preliminary study of the core shows that crystalline rocks are shock metamorphosed by shock pressure from 45 GPa near the surface to 15-20 GPa at a depth of about 5 km. The drill core allows the possibility of investigating many previously poorly studied cratering processes in the central part of the impact structure. As a first step one can use the estimates of energy for the homogeneous rock target. The diameter of the crater rim may be estimated as 40 km. The models elaborated earlier show that such a crater may be formed after collapse of a transient cavity with a radius of 10 km. The most probable range of impact velocities from 11.2 to 30 km/s may be inferred for the asteroidal impactor. For the density of a projectile of 2 g/cu cm the energy of the impact is estimated as 1E28 to 3E28 erg. In the case of vertical impact, the diameter of an asteroidal projectile is from 1.5 to 3 km for the velocity range from 11 to 30 km/s. For the most probable impact angle of 45 deg, the estimated diameter of an asteroid is slightly larger: from 2 to 4 km. Numerical simulation of the transient crater collapse has been done using several models of rock rheology during collapse. Results show that the column at the final position beneath the central mound is about 5 km in length. This value is close to the shock-pressure decay observed along the drill core. Further improvement of the model needs to take into account the blocky structure of target rocks revealed by drilling.

Ivanov, B. A.↗

NASA NTRS1959

Pressure Drag of Axisymmetric Cowls Having Large Initial Lip Angles at Mach Numbers from 1.90 to 4.90

The results of experimental and theoretical data on nine cowls are presented to determine the effect of initial lip angle and projected frontal area on the cowl pressure drag coefficient at Mach numbers from 1.90 to 4.90. The experimental drag coefficients were approximated well with two-dimensional shock-expansion theory at the lower cowl-projected areas, but the difference between theory and experiment increased as the cowl area ratio was increased or as shock detachment at the cowl lips was approached. An empirical chart is presented, which can be used to estimate the cowl pressure drag coefficient of cowls approaching an elliptic contour.

Samanich, Nick E.↗