Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Particle Hydrodynamics with Material Strength for Multi-Layer Orbital Debris Shield Design

Three dimensional simulation of oblique hypervelocity impact on orbital debris shielding places extreme demands on computer resources. Research to date has shown that particle models provide the most accurate and efficient means for computer simulation of shield design problems. In order to employ a particle based modeling approach to the wall plate impact portion of the shield design problem, it is essential that particle codes be augmented to represent strength effects. This report describes augmentation of a Lagrangian particle hydrodynamics code developed by the principal investigator, to include strength effects, allowing for the entire shield impact problem to be represented using a single computer code.

Fahrenthold, Eric P.↗
DOE OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗
NASA NTRS1972

Planetary geology, stellar evolution and galactic cosmology

Field studies of selected basalt flows in the Snake River Plain, Idaho, were made for comparative lunar and Mars geological investigations. Studies of basalt lava tubes were also initiated in Washington, Oregon, Hawaii, and northern California. The main effort in the stellar evolution research is toward the development of a computer code to calculate hydrodynamic flow coupled with radiative energy transport. Estimates of the rotation effects on a collapsing cloud indicate that the total angular momentum is the critical parameter. The study of Paschen and Balmer alpha lines of positronium atoms in the center of a galaxy is mentioned.

Source record↗
NASA NTRS1981

On the mass of the Cepheid X Cygni

Using a new non-Lagrangian hydrodynamic radiation transport code, models of the Cepheid X Cygni have been calculated and compared with new observations. Adopting an effective temperature of 5300 K, the luminosity was varied for stellar masses of 6, 8, and 9 solar masses in order to obtain a period of 16.4d. The results indicate that in order to find reasonable agreement between the light curves as calculated and observed, the mass used in the models must be nearly the evolutionary mass, not some 60% less as found in earlier work for Cepheids with periods from 7d to 10d. In addition, the new observations support the physical reality of the so-called 'artificial viscosity dip'.

Davis, C. G.↗
NASA NTRS2003

Studies of Solar Wind Interaction and Ionospheric Processes at Venus and Mars

This is the final report summarizing the work done during the last three years under NASA Grant NAG5-8946. Our efforts centered on a systematic development of a new generation of three dimensional magneto-hydrodynamic (MHD) numerical code, which models the interaction processes of the solar wind or fast flowing magnetospheric plasma with 'non-magnetic' solar system bodies (e.g. Venus, Mars, Europa, Titan). We have also worked on a number of different, more specific and discrete studies, as various opportunities arose. In the next few pages we briefly summarize these efforts.

Bogan, Denis↗
NASA NTRS2005

Computer Simulation of the VASIMR Engine

The goal of this project is to develop a magneto-hydrodynamic (MHD) computer code for simulation of the VASIMR engine. This code is designed be easy to modify and use. We achieve this using the Cactus framework, a system originally developed for research in numerical relativity. Since its release, Cactus has become an extremely powerful and flexible open source framework. The development of the code will be done in stages, starting with a basic fluid dynamic simulation and working towards a more complex MHD code. Once developed, this code can be used by students and researchers in order to further test and improve the VASIMR engine.

Garrison, David↗
NASA NTRS1995

The Production of Jets From Magnetic Accretion Disks: Simulation of the Blandford-Payne Mechanism

We have performed magnetohydrodynamic (MRD) simulations of the production of jets from magnetized accretion disks with a factor of 5 greater extent in space and time, and with more models, than any study published so far. We find that jets are produced by such disks in a broad range of parameter space, and by at least two different mechanisms. We also are able to follow the propagation of the jet well beyond the accretion disk into the region of hydrodynamic collimation. The code used is our MHD simulation code FLOW (K. Lind, D. Payne, D. Meier, and R. Blandford, 1989), converted to run on Caltech's massively parallel Intel Touchstone Delta supercomputer. Some of these models may be directly applicable to observed radio sources.

Blandford-Payne↗
DOE OSTI2025

The AGORA High-Resolution Galaxy Simulations Comparison Project

Context. Satellite galaxies experience multiple physical processes when interacting with their host halos, often leading to the quenching of star formation. In the Local Group, satellite quenching has been shown to be highly efficient, affecting nearly all satellites except the most massive ones. While recent surveys study Milky Way-analogs to assess how representative our Local Group is, the dominant physical mechanisms behind satellite quenching in Milky Way-mass halos remain under debate. Aims. We analyze satellite quenching within the same Milky Way-mass halo simulated using various widely used astrophysical codes, each using different hydrodynamic methods and implementing different supernovae feedback recipes. The goal is to determine whether quenched fractions, quenching timescales, and the dominant quenching mechanisms are consistent across codes or if they show sensitivity to the specific hydrodynamic method and supernovae feedback physics employed. Methods. We used a subset of high-resolution cosmological zoom-in simulations of a Milky Way-mass halo from the multiple-code AGORA CosmoRun suite. Our analysis focuses on comparing satellite quenching across the different models and against observational data. We also analyzed the dominant mechanisms driving satellite quenching in each model. Results. We find that the quenched fraction is consistent with the latest SAGA Survey results within its 1σ host-to-host scatter across all the models. Regarding quenching timescales, all the models reproduce the trend observed in the ELVES survey, Local Group observations, and previous simulations: The less massive the satellite, the shorter its quenching timescale. All of our models converge on the dominant quenching mechanisms: Strangulation halts cold gas accretion in all satellites, while ram pressure stripping is the predominant mechanism for gas removal, and it is particularly effective in satellites with M * <10 8 M ⊙ . Nevertheless, the efficiency of the stripping mechanisms differs among the codes, showing a strong sensitivity to the different supernovae feedback implementations and/or hydrodynamic methods employed.

Local Group↗
NASA NTRS2001

Formation and Early Evolution of Solar and Extra-Solar Giant Planets

This project investigates the origin of giant planets, both in the Solar System and around other stars. It is assumed that the planets form by the core accretion process: small solid particles in a disk surrounding a young star gradually coagulate into objects of a few kilometers in size, known as planetesimals, which then accumulate into solid protoplanetary cores. Once the cores have become large enough, they are able to attract gas from the surrounding disk to form the deep gaseous envelope of the giant planet. Our code simulates giant planet growth in a spherical approximation, and it has been quite successful in addressing a number of basic planetary properties. Further improvements to the code have been made to achieve a more realistic understanding of planetary formation. The computations of the models were based on an earlier version of our code and were stopped at the onset of runaway gas accretion. Now, improved boundary conditions have been incorporated into the code to allow for hydrodynamic inflow of gas and to handle the late stages of evolution when the planet evolves at constant mass. These changes were made to the version of the code that uses a constant accretion rate and to the version that uses a self-consistent method for calculating both the solid and gas accretion rates. The equation of state has been updated to incorporate the detailed tables of Saumon, Chabrier, and Van Horn. The opacities were updated to include the results of Alexander and Ferguson. The outer boundary conditions were modified. During the accretion phase when the planet's radius is between the accretion radius and the tidal radius, we set the outer boundary at a 'modified' accretion radius, which is the point where thermal energy is enough to bring gas to the edge of the Hill sphere.

Bodenheimer, P. H.↗
NASA NTRS1991

Three-dimensional compressible hydrodynamic convection in the sun and stars

A pseudospectral code is used to obtain numerical solutions of 3D compressible hydrodynamic convection in a stratified medium where the boundaries are open and where radiative losses are strongest at the top. Results are presented for a series of three simulations in progressively deeper boxes. The code does not make the anelastic approximation, but follows sound waves explicitly. Spatial organization is found in the form of vertical 'stacking' of smaller structures on top of larger ones, with each structure being roughly Hp to a vertical extent. Compressibility effects are apparent in the density: the snapshots indicate that, on a horizontal plane, the density at certain points may be of an order of 10 percent larger or smaller than elsewhere on that plane.

Hossain, Murshed↗
NASA NTRS1994

Two-dimensional hydrodynamical simulations of wind-compressed disks around rapidly rotating B stars

We use a two-dimensional piecewise parabolic method (PPM) code to simulate numerically the hydrodynamics of a radiation-driven stellar wind from a rapidly rotating Be star. The results generally confirm predictions of the semianalytic 'wind-compressed disk' model recently proposed by Bjorkman and Cassinelli to explain the circumstellar disks inferred observationally to exist around such rapidly rotating stars. However, this numerical simulation is able to incorporate several important effects not accounted for in the simple model, including a dynamical treatment of the outward radiative driving and gas pressure, as well as a rotationally distorted, oblate stellar surface. This enables us to model quantitatively the compressed wind and shock that forms the equatorial disk. The simulation results thus do differ in several important details from the simple method, showing, for example, an inner disk inflow not possible in the heuristic approach of assuming a fixed outward velocity law. There is also no evidence for the predicted detachment of the disk that arises in the fixed outflow picture. The peak equatorward velocity in the dynamical models is furthermore about a factor of 2 smaller than the lytically predicted value of approximately 50% of the stellar equatorial rotation speed. As a result, the dynamical disks are somewhat weaker than predicted, with a wider opening angle, lower disk/pole density ratio, and smaller shock velocity jump. The principal cause of these latter differences appears to be an artificially strong equatorward drift of the subsonic outflow in the original analytic model. Much better agreement with the dynamical results can be obtained, however, from a slightly modified, analytic wind-compression model with a more detailed specification of the fixed wind outflow and a lower boundary set to the sonic radius along a rotationally oblate stellar surface. Hence, despite these detailed differences, the general predicted effect of disk formation by wind compression toward the equator is substantially confirmed.

Owocki, Stanley P.↗
DOE OSTI2025

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

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↗