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,171 records · Page 65

Enabling Metric-Based Mesh Adaptation for Advanced Compressible Flow Simulations Using US3D

This work presents a metric-based mesh adaptation capability added to the flow solver US3D that is used to improve the accuracy of atmospheric entry simulations. Flow past an entry vehicle in free-flight is usually unsteady and characterized by anisotropic flow features like a strong bow shock and small-scale isotropic flow features like turbulent eddies in the wake. The results presented in this paper demonstrate that we are currently able to iteratively apply metric-based mesh adaptation to improve the solution accuracy by simultaneously targeting isotropic and anisotropic flow features. A simplified workflow of running CFD simulations for entry vehicles is presented where we start off with a coarse initial mesh to compute a first approximation of the flow. This coarse mesh is then iteratively refined using the newly developed metric to improve the resolution of the CFD computation. Furthermore, we demonstrate that we are able to adapt large mesh sizes in a timely fashion which allows us to start thinking about applying this adaptive framework to unsteady atmospheric entry CFD simulations.

Mesh adaptation↗

Flow Visualization for Plume-Surface Interaction Testing Within Large-Scale Vacuum Environments at Conditions Relevant to Lunar and Martian Landers

Interactions between rocket exhaust plumes and the landing surface during powered spacecraft descent on the Moon and Mars pose significant risks to the landing vehicle, landing site, and nearby infrastructure. Understanding the underlying plume-surface interaction phenomena through ground test data can provide critical insights on the sensitivities of parameters such as the spacecraft altitude and thrust. In the present work, we summarize a scaled ground test recently conducted within a 20-ft vacuum chamber environment located in the historical East Test Area at the NASA Marshall Space Flight Center. The ground test featured a Mach 5.3 inert gas plume impinging upon an instrumented flat plate. Planar laser-induced fluorescence, which is a 2D laser-based flow field measurement technique, was performed at this test area for the first time to visualize salient flow features such as the barrel shock, stagnation shock, and wall jet. Measurements were obtained at discrete lander altitudes, made dimensionless using the nozzle exit diameter, corresponding to h/D = 10, 8, 5, 4, 3, and 2. The stagnation pressures of the plume were varied from approximately 0.04 to 1 MPa, providing sensitivity to the engine thrust of a powered spacecraft. The unique test facility allowed for near-lunar conditions to be obtained at initial vacuum chamber pressures less than 0.1 Pa. Martian-relevant measurements were also performed at ambient pressures near 600 Pa. Furthermore, a third set of measurements were obtained at so called lunar-relevant conditions near 3 Pa chamber pressure

Plume Surface Interaction↗

Flow Visualization for Plume-Surface Interaction Testing Within Large-Scale Vacuum Environments at Conditions Relevant to Lunar and Martian Landers

Interactions between rocket exhaust plumes and the landing surface during powered spacecraft descent on the Moon and Mars pose significant risks to the landing vehicle, landing site, and nearby infrastructure. Understanding the underlying plume-surface interaction phenomena through ground test data can provide critical insights on the sensitivities of parameters such as the spacecraft altitude and thrust. In the present work, we summarize a scaled ground test recently conducted within a 20-ft vacuum chamber environment located in the historical East Test Area at the NASA Marshall Space Flight Center. The ground test featured a Mach 5.3 inert gas plume impinging upon an instrumented flat plate. Planar laser-induced fluorescence, which is a 2D laser-based flow field measurement technique, was performed at this test area for the first time to visualize salient flow features such as the barrel shock, stagnation shock, and wall jet. Measurements were obtained at discrete lander altitudes, made dimensionless using the nozzle exit diameter, corresponding to h/D = 10, 8, 5, 4, 3, and 2. The stagnation pressures of the plume were varied from approximately 0.04 to 1 MPa, providing sensitivity to the engine thrust of a powered spacecraft. The unique test facility allowed for near-lunar conditions to be obtained at initial vacuum chamber pressures less than 0.1 Pa. Martian-relevant measurements were also performed at ambient pressures near 600 Pa. Furthermore, a third set of measurements were obtained at so called lunar-relevant conditions near 3 Pa chamber pressure.

Plume Surface Interaction↗

The Evolution of the Viking Landing Gear

The primary function of the landing gear is to dissipate touchdown energy in a controlled fashion that minimizes the landing shock for onboard components while maximizing ground clearance and the probability of landing stability over the range of possible surface and touchdown parameter variations. Many other mission requirements and constraints were translated into the evolving Viking landing gear design as they arose. Hence, design considerations included such factors as prelaunch heat sterilization and noncontamination of the Martian landing site, gear stowage and deployment, terminal descent engine shutdown initiation, structural load attenuation, hard/soft landing surface capability, reliability, weight, and post-landed stability.

Pohlen, J. C.↗

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗

Impact-generated winds on Venus: Causes and effects

The pressure of the dense atmosphere of Venus significantly changes the appearance of ejecta deposits relative to craters on the Moon and Mercury. Conversely, specific styles and sequences of ejecta emplacement can be inferred to represent different intensities of atmospheric response winds acting over different timescales. Three characteristic timescales can be inferred from the geologic record: surface scouring and impactor-controlled (angle and direction) initiation of the long fluidized run-out flows; nonballistic emplacement of inner, radar-bright ejecta facies and radar-dark outer facies; and very late reworking of surface materials. These three timescales roughly correspond to processes observed in laboratory experiments that can be scaled to conditions on Venus (with appropriate assumptions): coupling between the atmosphere and earlytime vapor/melt (target and impactor) that produces an intense shock that subsequently evolves into blast/response winds; less energetic dynamic response of the atmosphere to the outward-moving ballistic ejecta curtain that generates nonthermal turbulent eddies; and late recovery of the atmosphere to impact-generated thermal and pressure gradients expressed as low-energy but long-lived winds. These different timescales and processes can be viewed as the atmosphere equivalent of shock melting, material motion, and far-field seismic response in the target. The three processes (early Processes, Atmospheric Processes, and Late Recovery Winds) are discussed at length.

Schultz, Pater H.↗

Radiative signals from impact of Shoemaker-Levy on Jupiter

The temperature and internal energy fields calculated by Takata et al. in the plume are used to calculate the greybody thermal radiation emitted versus wavelength to predict what might be observed by several spectral sensors operating from different platforms when fragments of Comet Shoemaker-Levy 9 (SL-9) impact Jupiter in July 1994. A SPH code was used by Takata et al. to calculate the full three dimensional flow and thermodynamic fields in the comet fragment and the atmosphere of Jupiter. We determined the fragment penetration depth, energy partitioning between the atmosphere and the impactor, and energy density deposited per unit length over the trajectory. Once the impactor had disintegrated and stopped, and the strong atmospheric shock decayed, the flow is driven by buoyancy effects. We then used our SPH code to calculate the flow and thermodynamic fields: pressure, article velocity, temperature, and internal energy distributions in the plume. The calculations for 2 and 10 km cometary fragments yield maximum deposition depths of approximately 175 and 525 km, respectively (1 bar = 0 km depth). We also calculated that 0.7 and 0.6 of the initial kinetic energy of the 10 and 2 km bolides, respectively, are deposited as internal energy in Jupiter's atmosphere.

Ahrens, Thomas J.↗

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗

Project DarkStar: Vision for LLNL in 2030

DarkStar was a Strategic Initiative (FY2021-FY2024) to investigate applications of Artificial Intelligence (AI) and Machine Learning (ML) to scientific problems of complex hydrodynamics, shockwave physics and energetic materials. The research focused on physics and engineering design as a process that can be tremendously accelerated through merging AI with advanced physics simulation on exascale-class platforms, and to experimentally validate this revolutionary new approach through dynamic materials campaigns. A central thread of scientific inquiry was in the application of AI to enable human understanding of how to control hydrodynamic instability (which has impacts to areas such as inertial confinement fusion) via engineering features and time-dependent sources. Motivated by an unfinished line of research started by Dr. Johnny von Neumann, AI-enabled simulation approaches were developed that allowed DarkStar researchers to uncover several ground-breaking discoveries regarding hydrodynamic instability, including how to completely suppress Richtmyer-Meshkov instability (RMI). These S&T discoveries, along with other advances, have shown the way for an entirely new approach to time-dependent problems known as inverse design – the idea that complex systems can be developed directly from a final state that is to be achieved and resolve the initial design via satisfying several constraints simultaneously via AI/ML. Through experimental campaigns conducted across a wide range of facilities in the NNSA complex (the High Explosive Application Facility at LLNL, the Dynamic Compression Sector/Advanced Photon Source at Argonne National Lab, and Special Technologies Laboratory at MSTS) the radical new AI/ML approach to engineering complex material dynamics was verified, establishing a new field of study within the realm of shock physics. As advanced manufacturing capabilities continue to develop, the great importance of inverse design as a means to apply that technology effectively for NNSA missions will feature prominently over this decade. DarkStar has positioned NNSA as a world-leader in this newly emerging cross-disciplinary area of AI methods for advanced physics simulation and pioneered multiple novel approaches that have enabled the broader scientific community. By allowing us to see past the horizon, to 2030 and beyond, DarkStar has illuminated the vast potential of AI/ML to impact a wide range of new national security missions and, consequently, multiple areas of further research have already emerged across the NNSA and DOD complex.

42 ENGINEERING↗

Study of exothermic processes in shock ignited gases by the use of laser shear interferometry.

The paper reports on the measurements of maximum exothermic power pulses attainable from a given chemical system. Experimental tests involved the use of a shock tube technique whereby the exothermic process of combustion was controlled by reflected shock so that it occurred under virtually inviscid flow conditions, the data having been obtained at a resolution commensurate with the actual rate of chemical reaction. Experimental observations were made by means of a novel method of laser shear interferometry. Chemical systems treated in this manner comprised of mixtures of stoichiometric hydrogen-oxygen with 80 and 90% argon, while the measurements covered the full permissible range of initial thermodynamic conditions.

Meyer, J. W.↗

A throat-bypass stability-bleed system using relief valves to increase the transient stability of a mixed-compression inlet

A stability-bleed system was installed in a YF-12 flight inlet that was subjected to internal and external airflow disturbances in the NASA Lewis 10 by 10 foot supersonic wind tunnel. The purpose of the system is to allow higher inlet performance while maintaining a substantial tolerance (without unstart) to internal and external disturbances. At Mach numbers of 2.47 and 2.76, the inlet tolerance to decreases in diffuser-exit corrected airflow was increased by approximately 10 percent of the operating-point airflow. The stability-bleed system complemented the terminal-shock-control system of the inlet and did not show interaction problems. For disturbances which caused a combined decrease in Mach number and increase in angle of attack, the system with valves operative kept the inlet started 4 to 28 times longer than with the valves inoperative. Hence, the stability system provides additional time for the inlet control system to react and prevent unstart. This was observed for initial Mach numbers of 2.55 and 2.68. For slow increase in angle of attack at Mach 2.47 and 2.76, the system kept the inlet started beyond the steady-state unstart angle. However, the maximum transient angles of attack without unstart could not be determined because wind-tunnel mechanical-stop limits for angle of attack were reached.

Neiner, G. H.↗

Time scales for formation and spreading of velocity shells of pickup ions in the solar wind

This paper discusses the process of assimilation (pickup) by the solar wind of newly ionized atoms and molecules. Generally, the pickup process is considered to evolve in three stages: (1) the initial interaction of newly created ions with the interplanetary magnetic field to form the ring-beam distribution; (2) pitch angle scattering of the ring beam to form a hollow shell; and (3) slower velocity diffusion to form a partially filled-in shell distribution. Using numerical simulations of turbulence such as would occur naturally in the solar wind and such as would be encountered near cometary bow shocks, the processes of shell formation and evolution are studied, and the results are used to estimate the time scales for shell formation and diffusion in several situations of recent observational interest, the interstellar He data obtained by AMPTE and cometary ion pickup distributions obtained by various spacecraft at comets Giacobini-Zinner and Halley.

Gaffey, J. D., Jr.↗

Lunar and Planetary Science XXXV: Effects of Impacts: Shock and Awe

This document discusses the following topics: Zircon as a Shock Indicator in Impactites of Drill Core Yaxcopoil-1, Chicxulub Impact Structure, Mexico; Experimental Investigation of Shock Effects in a Metapelitic Granulite; Experimental Reproduction of Shock Veins in Single-Crystal Minerals; Post-Shock Crystal-Plastic Processes in Quartz from Crystalline Target Rocks of the Charlevoix Impact Structure; Shock Reequilibration of Fluid Inclusions; How Does Tektite Glass Lose Its Water?; Assessing the Role of Anhydrite in the KT Mass Extinction: Hints from Shock-loading Experiments; A Mineralogical and Geochemical Study of the Nonmarine Permian/Triassic Boundary in the Southern Karoo Basin, South Africa; Extraterrestrial Chromium in the Permian-Triassic Boundary at Graphite Peak, Antarctica; Magnetic Fe,Si,Al-rich Impact Spherules from the P-T Boundary Layer at Graphite Peak, Antarctica; A Newly Recognized Late Archean Impact Spherule Layer in the Reivilo Formation, Griqualand West Basin, South Africa; Initial Cr-Isotopic and Iridium Measurements of Concentrates from Late Eocene Cpx-Spherule Deposits; An Ordinary Chondrite Impactor Composition for the Bosumtwi Impact Structure, Ghana, West Africa: Discussion of Siderophile Element Contents and Os and Cr Isotope Data.

Kyte, F. T.↗

The Adaptively-Refined MHD Solver (ARMS)

ARMS is a robust code for solving the single-fluid, fully compressible equations of magnetohydrodynamics. In this presentation I will describe how to access and run ARMS and how to use HelioSpace to analyze the results. The main features of the code are the following: 1. Equations are solved in either Cartesian or Spherical coordinates. For the latter, either a linear or a Log(r) coordinate can be used. 2. The latest flux-corrected-transport schemes developed at NRZ, are used to advance the convective equations. All waves and shocks are accurately captured. 3. A fully adaptive grid can be used. The code incorporates the PARAMESH technology for block-adaptation. Arbitrary adaptation criteria can be specified. 4. The code maintains a divergence-free magnetic field to machine accuracy throughout all refinement and de-refinement. 5. A uniform resistivity and thermal conduction can be included. Solution is either explicit or implicit. 6. Arbitrary boundary and initial conditions can be used in either Cartesian or spherical coordinates. 7. An accompanying visualization package, HelioSpace, is available for viewing and analyzing the native ARMS output, (without need for interpolating onto a structured grid).

Antichos, S. K.↗

Parameter Dependencies of Early-Stage Tangential Discontinuity-Driven Foreshock Bubbles in Local Hybrid Simulations

Foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere-ionosphere system. In the kinetic formation model, foreshock ions interact with the discontinuity by performing partial gyrations to generate currents that change the magnetic field topology around the discontinuity. However, how different foreshock ion properties affect the growth of the field variations is not well understood. Therefore, we use 2-D local hybrid simulations to study the effects of different foreshock ion distributions and properties on the growth of tangential discontinuity (TD)-driven FBs. We discover that for a gyrophase-bunched distribution with an initial phase where the guiding center is on the other side of the TD, the foreshock ions gyrate together across the TD, causing more foreshock ions to cross the TD and leading to a faster expansion of the structure than for a Maxwellian distribution. A ring distribution also yields higher expansion speeds because of the higher projected velocity into the new perpendicular direction. For Maxwellian distributions, there are positive and linear correlations of the FB expansion speeds with the initial foreshock ion densities, thermal speeds, parallel speeds, and sine of the TD magnetic shear angles. These parameter dependencies grow in strength as the structures evolve with time. The foreshock ion distributions and properties that lead to stronger currents produce more significant magnetic field variations and higher expansion speeds. Our study helps quantify the formation and expansion of FBs to forecast their space weather effects and contribution to shock acceleration.

Andrew Vu↗

Development of pebble-based extruded carbon rods for extreme plasma heat flux environments

This work presents first experiments toward the development of continuously renewable (extrudable) pebble-based carbon rods for use as plasma-facing components in extreme steady-state plasma flux environments. The primary envisioned application of this work is a first wall that can survive long-term in future magnetic fusion power reactors while also improving recovery of the reactor fuel (tritium and deuterium atoms). Bench tests applying extreme steady-state front-surface heat loads of up to 50 MW/m 2 are presented. Continuous pebble rod front-surface recession and intact pebble recovery are successfully demonstrated, at a rate of order 0.2 cm/s. Numerical simulations of the pebble rod front-surface recession are able to match observations reasonably well, indicating that the recession mechanism can be understood as occurring due to pebble thermal expansion and resulting shock and cracking of the inter-pebble binder. Tests of the pebble rod extrusion demonstrate that friction between the rods and the stainless-steel extrusion channel is tolerably low (<50 N for the expected channel length) over a wide range of temperatures. Front-surface outgassing rates below 1000 Torr L/s/m 2 are achieved, believed to be sufficiently low for use in magnetic fusion reactors. In conclusion, initial parametric scans over pebble rod size and binder fraction to vary front-surface recession rates are presented.

36 MATERIALS SCIENCE↗

High speed viscous flow calculations about complex configurations

Applications of the NASA Ames Parabolized Navier-Stokes (PNS) code to a variety of complex generic configurations is presented. The algorithm, boundary conditions, initial conditions, and grid generators are discussed as applied to these configurations. The PNS code was used as the mainline procedure to numerically simulate the viscous supersonic flow over these generic configurations. The turbulence model that was used in this study is the Baldwin-Lomax model. The boundary conditions are the usual viscous no slip at the wall, and a characteristic procedure is used to fit the bow shock wave which is the outermost boundary. An elliptic grid generator is employed to discretize the flow domain. In addition, an equilibrium air capability has been incorporated into the code. It uses the curve fits of Tannehill, et al. The flow regimes vary from a Mach number of 2 up to 25. Both laminar and turbulent flow are considered. Varying angles of attack have also been computed. Configurations vary from simple cone-type bodies to lifting winged bodies, such as the space shuttle or the generic supersonic cruise fighter.

Chaussee, D. S.↗

Re-forming supercritical quasi-parallel shocks. II - Mechanism for wave generation and front re-formation

This paper continues the study of Thomas et al. (1990) in which hybrid simulations of quasi-parallel shocks were performed in one and two spatial dimensions. To identify the wave generation processes, the electromagnetic structure of the shock is examined by performing a number of one-dimensional hybrid simulations of quasi-parallel shocks for various upstream conditions. In addition, numerical experiments were carried out in which the backstreaming ions were removed from calculations to show their fundamental importance in reformation process. The calculations show that the waves are excited before ions can propagate far enough upstream to generate resonant modes. At some later times, the waves are regenerated at the leading edge of the interface, with properties like those of their initial interactions.

Winske, D.↗