Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic codes”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Shocks in the Scrape-off Layer: Hydrodynamical Aspects (Summary) [Slides]

Hydrodynamical model provides a robust description of mass, momentum, and energy transfer. It can be used for benchmarking hydrodynamical modules of SOL codes. It clarifies the role of boundary conditions at the end-plates. It provides a “geometrical” model for partition of power flows between two divertor legs, both for the steady-state and pulsed flows. It indicates possible presence of shocks in the divertor legs. A framework for finding the shock location and strength has been developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗
NASA NTRS2003

The Formation of Asteroid Satellites in Catastrophic Impacts: Results from Numerical Simulations

We have performed new simulations of the formation of asteroid satellites by collisions, using a combination of hydrodynamical and gravitational dynamical codes. This initial work shows that both small satellites and ejected, co-orbiting pairs are produced most favorably by moderate-energy collisions at more direct, rather than oblique, impact angles. Simulations so far seem to be able to produce systems qualitatively similar to known binaries. Asteroid satellites provide vital clues that can help us understand the physics of hypervelocity impacts, the dominant geologic process affecting large main belt asteroids. Moreover, models of satellite formation may provide constraints on the internal structures of asteroids beyond those possible from observations of satellite orbital properties alone. It is probable that most observed main-belt asteroid satellites are by-products of cratering and/or catastrophic disruption events. Several possible formation mechanisms related to collisions have been identified: (i) mutual capture following catastrophic disruption, (ii) rotational fission due to glancing impact and spin-up, and (iii) re-accretion in orbit of ejecta from large, non-catastrophic impacts. Here we present results from a systematic investigation directed toward mapping out the parameter space of the first and third of these three collisional mechanisms.

Durda, D. D.↗
NASA NTRS2016

Patterns of Activity in A Global Model of A Solar Active Region

In this work we investigate the global activity patterns predicted from a model active region heated by distributions of nanoflares that have a range of frequencies. What differs is the average frequency of the distributions. The activity patterns are manifested in time lag maps of narrow-band instrument channel pairs. We combine hydrodynamic and forward modeling codes with a magnetic field extrapolation to create a model active region and apply the time lag method to synthetic observations. Our aim is not to reproduce a particular set of observations in detail, but to recover some typical properties and patterns observed in active regions. Our key findings are the following. (1) Cooling dominates the time lag signature and the time lags between the channel pairs are generally consistent with observed values. (2) Shorter coronal loops in the core cool more quickly than longer loops at the periphery. (3) All channel pairs show zero time lag when the line of sight passes through coronal loop footpoints. (4) There is strong evidence that plasma must be re-energized on a timescale comparable to the cooling timescale to reproduce the observed coronal activity, but it is likely that a relatively broad spectrum of heating frequencies are operating across active regions. (5) Due to their highly dynamic nature, we find nanoflare trains produce zero time lags along entire flux tubes in our model active region that are seen between the same channel pairs in observed active regions.

Bradshaw, S. J.↗
DOE OSTI2024

Hyperparameter Studies for Vision Transformers Trained on High-Fidelity Simulations

This library is a collection of python modules that define, train, and analyze vision-transformer (ViT) machine learning models. The code implements, with mild modifications, ViT models that have been made publicly available through publication and GitHub code. The training data for these models is hydrodynamic simulation output in the form of numpy arrays. This library contains code to train these ViT models on the hydrodynamic simulation output with a variety of hyperparameters, and to compare the results of such models. Furthermore, the library contains definitions of simple convolutional neural network (CNN) machine learning architectures which can be trained on the same hydrodynamic simulation output. These are included as a reference point to compare the ViT models to. Additionally, the library includes trained ViT and CNN models and example input data for demonstration purposes. The code is based on the PyTorch python library.

Callis, Skylar↗
NASA NTRS2014

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗
DOE OSTI2020

ALEGRA/Sceptre Code Coupling [Brief]

Researchers at Sandia have developed an advanced radiation hydrodynamic simulation capability by coupling the ALEGRA and Sceptre codes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1987

Analysis of planetary evolution with emphasis on differentiation and dynamics

In order to address the early stages of nebula evolution, a three-dimensional collapse code which includes not only hydrodynamics and radiative transfer, but also the effects of ionization and, possibly, magnetic fields is being addressed. As part of the examination of solar system evolution, an N-body code was developed which describes the latter stages of planet formation from the accretion of planetesimals. To test the code for accuracy and run-time efficiency, and to develop a stronger theoretical foundation, problems were studied in orbital dynamics. A regional analysis of the correlation in the gravity and topography fields of Venus was performed in order to determine the small and intermediate scale subsurface structure.

Kaula, William M.↗
DOE OSTI2020

Characterization of an imploding cylindrical plasma for electron transport studies using x-ray emission spectroscopy

We report here on the characterization of the conditions of an imploding cylindrical plasma by time-resolved x-ray emission spectroscopy. Knowledge about this implosion platform can be applied to studies of particle transport for inertial confinement fusion schemes or to astrophysical plasmas. A cylindrical Cl-doped CH foam within a tube of solid CH was irradiated by 36 beams (I total ~5×10 14 W/cm 2 , 1.5ns square pulse, and E total ~ 16.2kJ) of the OMEGA-60 laser to radially compress the CH toward the axis. The analysis of the time-resolved spectra showed that the compression can be described by four distinct phases, each presenting different plasma conditions. First the ablation of the cylinder is dominant; second, the foam is heated and induces a significant jump in emission intensities; third, the temperature and density of the foam reaches a maximum; and finally, the plasma expands. Ranges for the plasma temperature were inferred with the atomic physics code SCRAM (Spectroscopic Collisional-Radiative Atomic Model) and the experimental data have been compared to hydrodynamic simulations performed with the 2D code FLASH, which showed a similar implosion dynamic over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
NASA NTRS1989

Continuous infrared emission of proto- and young-planetary nebulae

Evolutionary sequences of planetary nebulae (PNe) were calculated beginning from the moment of shell ejection to its dissipation and IR spectra of the outgoing radiation were obtained for different stages of the evolution. To solve the coupled hydrodynamical and radiation transfer problem a computer code was used in which gas and dust are treated as two separate hydrodynamical components. The models, which are calculated for two grain materials (graphite and silicate) have a size distribution of particles based on that found for the interstellar medium. In the computations, the recent optical properties of graphite and astronomical silicate grains tabulated by Draine (1987) were employed. It was assumed that silicate of graphite grains exist below some critical temperature equal to 1300 or to 1700 K, respectively. Infrared model spectra are compared with the IR radiation emitted by some proto- and some young-planetary nebulae. The observed IR continua can be quite well matched with our models with grains having a reasonable size range.

Szczerba, Ryszard↗
DOE OSTI2025

Polarization of the 𝜙 meson in the hadronic phase with nucleon scatterings and a viscous hydrodynamic background

We extend our previous work on the spin alignment of the ϕ vector meson in the hadronic phase of the quark gluon plasma, by including effects of nucleon scatterings. The emission rates are calculated in a realistic hydrodynamic background simulated with the code FLUIDuM , for different beam energies. We find that all the effects taken into account cannot explain the out-of-plane spin alignment of the ϕ meson observed experimentally.

Polarization phenomena↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗
DOE OSTI2023

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. Furthermore, a terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Capture and migration of Jupiter and Saturn in mean motion resonance in a gaseous protoplanetary disc

Here, we study the dynamical evolution of Jupiter and Saturn embedded in a gaseous, solar nebula-type disc by means of hydrodynamics simulations with the FARGO2D1D code. We study the evolution for different initial separations of the planets’ orbits, Δa SJ , to investigate whether they become captured in mean motion resonance (MMR) and the direction of the subsequent migration of the planet (inwards or outwards). We also provide an assessment of the planet’s orbital dynamics at different epochs of Saturn’s growth. We find that the evolution of initially compact orbital configurations is dependent on the value of Δa SJ . This implies that an evolution as that proposed in the Grand Tack model depends on the precise initial orbits of Jupiter and Saturn and on the time-scales for their formation. Capture in the 1:2 MMR and inward or (nearly) stalled migration are highly favoured. Within its limits, our work suggests that the reversed migration, associated with the resonance capture of Jupiter and Saturn, may be a low-probability evolutionary scenario, so that other planetary systems with giant planets are not expected to have experienced a Grand Tack-like evolutionary path.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Studies of localized mix and asymmetry caused by the stalk in inertial fusion experiments at OMEGA

In an inertially confined fusion experiment, a capsule of fusion fuel is compressed to high temperatures and densities using lasers, creating a hot plasma in which thermonuclear fusion reactions take place. To achieve efficient conversion of laser energy to fuel heating, the implosion must be as spherically symmetric as possible. One common source of asymmetry that degrades performance is the stalk and glue spot that hold the capsule in place before the experiment. The exact nature and magnitude of this perturbation had not been previously quantified. In this work, inertially confined fusion experiments were executed with an additional stalk on the capsule to measure and quantify the effect of the stalk and glue spot on the symmetry and performance. It is found that the stalk drives low-mode asymmetry by shadowing the laser drive on the part of the shell near the glue spot, enhancing the hot-spot velocity toward the stalk by 40–80 km/s and elongating the hot-spot along the stalk axis. It is also found that localized mix due to the additional stalk has a negligible effect on implosion performance. Two-dimensional radiation-hydrodynamic simulations using the xRAGE code reproduce the direction of the experimentally observed change to hot-spot velocity and implosion performance metrics, but not the observed elongation, suggesting that this feature results from an interaction between stalk perturbations and other preexisting asymmetry seeds, such as laser drive asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Particle-in-cell simulations of expanding high energy density plasmas with laser ray tracing

The design and analysis of high energy density (HED) laser experiments typically rely on radiation hydrodynamics simulations. However, some laser–plasma interaction regimes are not collisional and cannot be adequately modeled with hydrodynamics. For example, strongly driven magnetic reconnection and magnetized collisionless shock experiments possess extended hydrodynamic or even kinetic properties, necessitating first-principles kinetic simulations. In this paper, we present the benchmarking and first results obtained with a laser-ray-tracing and inverse bremsstrahlung absorption module implemented in the particle-in-cell code PSC. The simulation results are compared to radiation hydrodynamic simulations using the FLASH code as well as analytical estimates. We successfully benchmark the energy deposition model and overall hydrodynamic evolution of the systems. We also consider possible kinetic effects that may be expected from laser-target ablation in the HED regime, including non-local transport and two-temperature effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗