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 559 records · Page 31

Computational Relativistic Astrophysics Using the Flow Field-Dependent Variation Theory

We present our method for solving general relativistic nonideal hydrodynamics. Relativistic effects become pronounced in such cases as jet formation from black hole magnetized accretion disks which may lead to the study of gamma-ray bursts. Nonideal flows are present where radiation, magnetic forces, viscosities, and turbulence play an important role. Our concern in this paper is to reexamine existing numerical simulation tools as to the accuracy and efficiency of computations and introduce a new approach known as the flow field-dependent variation (FDV) method. The main feature of the FDV method consists of accommodating discontinuities of shock waves and high gradients of flow variables such as occur in turbulence and unstable motions. In this paper, the physics involved in the solution of relativistic hydrodynamics and solution strategies of the FDV theory are elaborated. The general relativistic astrophysical flow and shock solver (GRAFSS) is introduced, and some simple example problems for computational relativistic astrophysics (CRA) are demonstrated.

Richardson, G. A.↗

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerothermal Effects of Hypersonic Shock Impingement on a Highly Swept Leading Edge

A numerical study is conducted to investigate aerothermal implications of oblique shock impingement on a highly swept leading edge in hypersonic flow, with flow conditions and geometry relevant to non-propulsive leading edge heating on hypersonic cruise vehicles. In contrast to propulsive inlet leading edge shock impingement, which has been subject to extensive historical study, bow shock impingement on a highly swept wing leading edge has received comparatively little attention in the literature. Distinguishing characteristics of this type of flow include a potentially dominant surface heating contribution from crossflow-induced viscous dissipation, possible wall velocity slip and temperature jump effects, and various modes of shock layer thermal nonequilibrium that can interact with shock impingement flowfield structures. One overall goal of this work is to assess the impact of bow shock impingement and related nonequilibrium phenomena for wing leading edge design, materials selection and thermal management. Another goal is to improve understanding – and thereby reduce associated design margins – for the physical processes and mechanisms that may influence wing leading edge heating in the presence of an impinging shock. A combination of direct simulation Monte Carlo and Reynolds-averaged Navier-Stokes analysis is used to assess surface and flowfield characteristics for a series of hypersonic shock impingement problems, based on historical ground tests and related flight conditions. Favorable comparisons are performed with experimental data, and various thermal nonequilibrium effects are quantified.

direct simulation Monte Carlo↗

Shock-driven three-fluid mixing with various chevron interface configurations

When a shock wave crosses a density interface, the Richtmyer–Meshkov instability causes perturbations to grow. Richtmyer–Meshkov instabilities arise from the deposition of vorticity from the misaligned density and pressure gradients at the shock front. In many engineering applications, microscopic surface roughness will grow into multi-mode perturbations, inducing mixing between the fluid on either side of an initial interface. Applications often have multiple interfaces, some of which are close enough to interact in the later stages of instability growth. In this study, we numerically investigate the mixing of a three-layer system with periodic zigzag (or chevron) interfaces, calculating the dependence of the width and mass of mixed material on properties such as the shock timing, chevron amplitude, multi-mode perturbation spectrum, density ratio, and shock mach number. The multi-mode case is also compared with a single-mode perturbation. The Flash hydrodynamic code is used to solve the Euler equations in three dimensions with adaptive grid refinement. Key results include a significant increase in mixed mass when changing from a single-mode to a multi-mode perturbation on one of the interfaces. The mixed width is mainly sensitive to the density ratio and chevron amplitude, whereas the mixed mass also depends on the multi-mode spectrum. In conclusion, steeper initial perturbation spectra have lower mixed mass at early times but a greater mixed mass after the reflected shock transits back across the layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗

Introduction to Special Issue on the Early History of Nuclear Fusion

This introductory paper to the special issue of Fusion Science and Technology commemorates early research on fusion conducted at Los Alamos (the singular entity denoted Los Alamos Laboratory/Los Alamos Scientific Laboratory/Los Alamos National Laboratory at different times is designated “Los Alamos” in this paper) in support of the eventual H-bomb program. We survey the historical origins of the thermonuclear program, what was known of fusion reactions at the outbreak of the war, and the remarkable breakthroughs involving particularly the prospect of deuterium-tritium (DT) reactions conducted during the war, and we summarize the papers in this volume. Much of the nuclear fusion technical history presented herein has not been previously reported. Papers describe aspects of fusion science during these days, on shock hydrodynamics and on electron-radiation coupling, and on nuclear physics including the discoveries of resonances in both the DT cross section and in the lithium tritium-breeding cross section. Three papers follow our colleague Mark Paris’s finding Arthur Ruhlig’s 1938 paper on the first observation of DT fusion: one on how it influenced subsequent Manhattan Project research, another on a modern calculation of that historic experiment, and a third that has repeated the experiment using modern experimental capabilities. Other papers discuss how the first H-bomb test, Ivy Mike, led to the discovery of the new elements einsteinium and fermium and how the DT fusion processes played a key role in our universe’s development after the Big Bang. We also present a paper that analyzes the pioneering Cambridge University 1934 experiment by Marcus Oliphant, Paul Harteck, and Ernest Rutherford where deuterium-deuterium fusion was first observed and that describes how Ernest Lawrence missed identifying fusion in 1933. Finally, we present a summary of early concepts for controlled fusion energy that grew out of wartime discussions at Los Alamos. The papers show how J. Robert Oppenheimer played a leading technical role in the early developments of the H-bomb, before his later opposition—our first paper in this issue addresses the U.S. Department of Energy’s 2022 vacation of the earlier 1954 decision to revoke his security clearance.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Aerodynamic interaction of heat transfer with steady transonic flow

An experimental study investigated the interactive effect between the heat transfer and the steady transonic flow in a simple convergent-divergent duct. The results show that the position of the shock is considerably affected by the effect of wall cooling. The physical explanation for the shift is also offered. The results may have implications in the other internal flow problems.

Kurosaka, M.↗

Mach number effects on transonic aeroelastic forces and flutter characteristics

Transonic aeroelastic stability analysis and flutter calculations are presented for a generic transport-type wing based on the use of the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code. The CAP-TSD code was recently developed for transonic unsteady aerodynamic and aeroelastic analysis of complete aircraft configurations. A binary aeroelastic system consisting of simple bending and torsion modes was used to study aeroelastic behavior at transonic speeds. Generalized aerodynamic forces are presented for a wide range of Mach number and reduced frequency. Aeroelastic characteristics are presented for variations in freestream Mach number, mass ratio, and bending-torsion frequency ratio. Flutter boundaries are presented which have two transonic dips in flutter speed. The first dip is the usual transonic dip involving a bending-dominated flutter mode. The second dip is characterized by a single degree-of-freedom torsion oscillation. These aeroelastic results are physically interpreted and shown to be related to the steady state shock location and changes in generalized aerodynamic forces due to freestream Mach number.

Mohr, Ross W.↗

Mach number effects on transonic aeroelastic forces and flutter characteristics

Transonic aeroelastic stability analysis and flutter calculations are presented for a generic transport-type wing based on the use of the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code. The CAP-TSD code was recently developed for transonic unsteady aerodynamic and aeroelastic analysis of complete aircraft configurations. A binary aeroelastic system consisting of simple bending and torsion modes was used to study aeroelastic behavior at transonic speeds. Generalized aerodynamic forces are presented for a wide range of Mach number and reduced frequency. Aeroelastic characteristics are presented for variations in freestream Mach number, mass ratio, and bending-torsion frequency ratio. Flutter boundaries are presented which have two transonic dips in flutter speed. The first dip is the usual transonic dip involving a bending-dominated flutter mode. The second dip is characterized by a single degree-of-freedom torsion oscillation. These aeroelastic results are physically interpreted and shown to be related to the steady state shock location and changes in generalized aerodynamic forces due to freestream Mach number.

Mohr, Ross W.↗

On the 'flip-flop' instability of Bondi-Hoyle accretion flows

A simple physical interpretation is advanced by means of an analysis of the shock cone in the accretion flows past a compact object and with an examination of the accretion-line stability analyses. The stability of the conical shock is examined against small angular deflections with attention given to several simplifying assumptions. A line instability is identified in the Bondi-Hoyle accretion flows that leads to the formation of a large opening-angle shock. When the opening angle becomes large the instability becomes irregular oscillation. The analytical methodology is compared to previous numerical configurations that demonstrate different shock morphologies. The Bondi-Hoyle accretion onto a compact object is concluded to generate a range of nonlinear instabilities in both homogeneous and inhomogeneous cases with a quasiperiodic oscillation in the linear regime.

Livio, Mario↗

Automated Fluid Feature Extraction from Transient Simulations

In the past, feature extraction and identification were interesting concepts, but not required to understand the underlying physics of a steady flow field. This is because the results of the more traditional tools like iso-surfaces, cuts and streamlines were more interactive and easily abstracted so they could be represented to the investigator. These tools worked and properly conveyed the collected information at the expense of much interaction. For unsteady flow-fields, the investigator does not have the luxury of spending time scanning only one 'snap-shot' of the simulation. Automated assistance is required in pointing out areas of potential interest contained within the flow. This must not require a heavy compute burden (the visualization should not significantly slow down the solution procedure for co-processing environments like pV3). And methods must be developed to abstract the feature and display it in a manner that physically makes sense. The following is a list of the important physical phenomena found in transient (and steady-state) fluid flow: Shocks; Vortex ores; Regions of Recirculation; Boundary Layers; Wakes.

Haimes, Robert↗

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗

A Physics-Based Approach to Hypersonic Nonequilibrium Chemistry

NASA has undertaken a concerted effort to generate a complete model for the chemical reactions and other collisional processes that take place in the bow shock layer and wake of spacecraft entering the atmospheres of Earth and Mars. This new model is based on first-principles physics calculations: interaction potentials between colliding species obtained from accurate quantum mechanics calculations and reaction rate coefficients and cross sections based on classical mechanics simulations of heavy particle collisions.This presentation will give an introduction to the underlying theory and computational methods used in this work and describe the range of chemical and physical processes that take place in the close vicinity of spacecraft during atmospheric entry. I will also present some of our recent results for reactions involving nitrogen, carbon monoxide and CO2.

Jaffe, Richard↗

Molecular line emission models of Herbig-Haro objects. I - H2 emission

A comprehensive model for molecular hydrogen emssion in Herbig-Haro objects that are associated with the heads of radiative stellar jets is presented by using a simple representation of the jet head as a comprising a leading bow shock and a trailing jet shock, separated by a dense layer of cool shocked gas. Attention is given to collisional excitation in a nondissociative shock and formation pumping in the molecular reformation zone behind a dissociative shock, employing detailed shock and photodissociation-region emission models that incorporate most of the relevant atomic physics and chemistry. The conditions under which each of these excitation mechanisms may be expected to contribute to the observed emission are discussed, and a general diagnostic scheme for discriminating among them is constructed. Applying this scheme to the HH 1-2 system, strong evidence for excitation by the radiation field of a fast shock is found. It is inferred that FUV pumping contributes a significant fraction of the H2 line emission, and it is shown that this can occur only if the UV pump lines are not strongly self-shielded.

Wolfire, Mark G.↗

arcjetCV: automating recession extraction from video

Arc jet Computer Vision (arcjetCV)[1][2] is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗

ArcjetCV: Automating Recession Tracking

Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗

ArcjetCV: Automating Arc Jet Analysis

Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗