Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic simulation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Investigating the impact of intermediate-mode perturbations on diagnosing plasma conditions in DT cryogenic implosions via synthetic x-ray Thomson scattering

The pursuit of inertial confinement fusion ignition target designs requires precise experimental validation of the conditions within imploding capsules, in particular the density and temperature of the compressed shell. Previous work has identified x-ray Thomson scattering (XRTS) as a viable diagnostic tool for inferring the in-flight compressed deuterium-tritium shell conditions during capsule implosions (Poole et al 2022 Phys. Plasmas 29 072703). However, this study focused on one-dimensional simulations, which do not account for the growth of hydrodynamic instabilities. In this work, two-dimensional DRACO simulations incorporating intermediate-mode perturbations up to Legendre mode l = 50 were used to generate synthetic XRTS spectra with the SPECT3D code. The analysis employed Markov-Chain Monte Carlo techniques to infer plasma conditions from these spectra. The results demonstrate that the XRTS diagnostic platform can effectively discern the in-flight compressed shell conditions for targets with varying adiabats, even in the presence of intermediate-mode perturbations. This work underscores the potential of XRTS for realistic inertial confinement fusion experiments, providing a robust method for probing the complex dynamics of fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

An Axisymmetric Hydrodynamical Model for the Torus Wind in AGN: X-ray Excited Funnel Flow - 2

We have calculated a series of models of outflows from the obscuring torus in active galactic nuclei (AGN). Our modeling assumes that the inner face of a rotationally supported torus is illuminated and heated by the intense X-rays from the inner accretion disk and black hole. As a result of such heating a strong biconical outflow is observed in our simulations. We calculate 3-dimensional hydrodynamical models, assuming axial symmetry, and including the effects of X-ray heating, ionization, and radiation pressure. We discuss the behavior of a large family of these models, their velocity fields, mass fluxes and temperature, as functions of the torus properties and X-ray flux. Synthetic warm absorber spectra are calculated, assuming pure absorption, for sample models at various inclination angles and observing times. We show that these models have mass fluxes and flow speeds which are comparable to those which have been inferred from observations of Seyfert 1 warm absorbers, and that they can produce rich absorption line spectra.

Dorodnitsyn, A.↗
NASA NTRS1997

Modeling of Water Injection into a Vacuum

A loosely coupled two-phase vacuum water plume model has been developed. This model consists of a source flow model to describe the expansion of water vapor, and the Lagrangian equations of motion for particle trajectories. Gas/Particle interaction is modeled through the drag force induced by the relative velocities. Particles are assumed traveling along streamlines. The equations of motion are integrated to obtain particle velocity along the streamline. This model has been used to predict the mass flux in a 5 meter radius hemispherical domain resulting from the burst of a water jet of 1.5 mm in diameter, mass flow rate of 24.2 g/s, and stagnation pressure of 21.0 psia, which is the nominal Orbiter water dump condition. The result is compared with an empirical water plume model deduced from a video image of the STS-29 water dump. To further improve the model, work has begun to numerically simulate the bubble formation and bursting present in a liquid stream injected into a vacuum. The technique of smoothed particle hydrodynamics was used to formulate this simulation. A status and results of the on-going effort are presented and compared to results from the literature.

Alred, John W.↗
DOE OSTI2024

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves↗
NASA NTRS1992

Giant impacts on a primitive Uranus

Simulations of collisions are conducted between a model of the primitive Uranus and 1-3 earth-mass impactors, using smooth-particle hydrodynamics. A series of collisions was simulated for each impactor while varying the total angular momentum of the system. Most of the simulation runs left ices in orbit; a subset of the runs also left rock or iron (from the impactor). It is concluded on the basis of these results that there is a wide range of giant impacts which could have produced the current period and inclination of the spin axis relative to the plane of the ecliptic. A subset of these could have deposited the material in orbit from which the regular satellites of Uranus were assembled.

Slattery, Wayne L.↗
NASA NTRS1993

Numerical simulation of plasma processes driven by transverse ion heating

The plasma processes driven by transverse ion heating in a diverging flux tube are investigated with numerical simulation. The heating is found to drive a host of plasma processes, in addition to the well-known phenomenon of ion conics. The downward electric field near the reverse shock generates a doublestreaming situation consisting of two upflowing ion populations with different average flow velocities. The electric field in the reverse shock region is modulated by the ion-ion instability driven by the multistreaming ions. The oscillating fields in this region have the possibility of heating electrons. These results from the simulations are compared with results from a previous study based on a hydrodynamical model. Effects of spatial resolutions provided by simulations on the evolution of the plasma are discussed.

Singh, Nagendra↗
NASA NTRS1988

Numerical simulation of nonaxisymmetric adiabatic accretion flow

The hydrodynamics of gas flow past a finite-sized gravitating central object is studied in the adiabatic approximation in two spatial dimensions. The flow morphology is shown to depend on the ratio (epsilon) of the axisymmetric accretion radius to the length scale of the density variation. For small to intermediate epsilon values the flows are highly time dependent, and for large epsilon values a quasi-steady state is obtained.

Fryxell, B. A.↗
NASA NTRS1993

Galaxy tracers in N-body simulations

Using the method of smoothed particle hydrodynamics, we have modeled the formation of a compact group of galaxies with sufficient resolution to trace galaxies. Radiative cooling allows the baryons to dissipate their thermal energy and collapse to overdensities characteristic of real galaxies. With their cross section greatly reduced, these galaxy tracers remain distinct during cluster formation while their dark matter halos merge. In addition, the number density, the mass distribution function, and even the morphology of these objects are similar to those of observed galaxies. A viable population of galaxy tracers can be unambiguously defined.

Summers, F. J.↗
NASA NTRS2012

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Dubey, A.↗
NASA NTRS1974

Wind flow over simulated structures for Martian and terrestrial atmospheres

A two-dimensional nonlinear hydrodynamic calculation is presented of a 50-m/sec wind flowing past a simulated Martian crater, and an attempt is made to reproduce a wind flow observed in a California valley. The computer program and boundary conditions used in the calculation are described, and streamlines of the flow over the simulated crater are plotted. The attempt to calculate the wind flow over Owens Valley, California, is made to test the validity of the results for the Martian wind. It is shown that the Martian calculation yields a quasi-steady flow similar to that produced in a wind tunnel simulation of a similar problem. The California wind could not be reproduced with a coarsely resolved calculation, leaving the validity of the Martian calculation and the wind tunnel simulation in doubt.

Jones, E. M.↗
DOE OSTI2026

The Effects of Realistic Impact-processed Asteroid Topology on Atmospheric Airbursts

When small asteroids (∼10–200 m) enter Earth’s atmosphere they can break up before hitting the ground, releasing a considerable amount of energy in what is known as an atmospheric airburst. To date, numerous studies have examined how an asteroid’s strength, size, shape, and entry characteristics affect the breakup process; but few have assessed the effects of heterogenous internal topologies. We present a two-step high-fidelity simulation approach, which we use to model the atmospheric breakup of asteroids with realistic impact-processed internal topologies and their associated ground effects. In the first step, we model the asteroid as it is impacted to generate an “impact-processed” topology. We then map the resulting fractures (i.e., damage field) from that simulation to the initial state of an atmospheric breakup simulation. For both phases, we use the smoothed particle hydrodynamics code Spheral++. As a demonstrative example, we apply our new simulation pipeline to a Tunguska-scale impactor (50 m diameter, 20 km s −1 entry velocity, and 45° angle) in a case study comparing impact-processed and fully intact asteroids in both planar 2D and full 3D. For sufficiently strong asteroids, we find that impact processing generally causes the asteroid to break up at higher altitudes. We also observe that impact-processed asteroids break up into smaller fragments than objects with equivalent properties that have not been impact-processed.

Astronomy and AstroPhysics↗
DOE OSTI2026

Capturing Sub‐Kilometer Flood Inundation Dynamics During the California Rain‐on‐Snow Events of 2017

The severe impacts of rain‐on‐snow (ROS) extreme events have been widely recognized and studied. However, unlike hydrological processes, flood inundation dynamics and the relative contribution of rainfall and snowmelt during ROS events remain under‐investigated. We diagnosed and documented sub‐kilometer spatio‐temporal dynamics of flood inundation during the 2017 California ROS events, simulated by a 2‐dimensional hydrodynamic model, River Dynamical Core (RDycore). RDycore shows good performance in capturing fine‐scale flood inundation dynamics against gauge measurements (with a median correlation coefficient of 0.81) and satellite observations. On top of rainfall, snowmelt not only increases the mean maximum inundation depth (12.0%–25.1%), but also expands the total flooded area (19.9%–31.9%) and prolongs the mean flood duration (3.4%–7.1%) across the events. Our findings offer an explicit and accurate picture of when and where ROS flooding could occur and how snowmelt increases flood hazard, valuable for risk assessment and infrastructure planning.

Flood inundation↗
DOE OSTI2025

Mediation of Colloidal Encounter Dynamics by Surface Roughness

Rigorous understanding on the self-assembly of colloidal nanocrystals is crucial to develop tailored nanostructured materials for energy storage, sensing, and optical fields. Despite extensive studies on the self-assembly, a mechanistic understanding of self-assembly under an external field still remains an ongoing challenge. Here, in this work, we used optical tweezers that impose an external attractive force field, resulting in the self-assembly of alpha-phase sodium yttrium fluoride nanocrystals. The dynamic force that is strongly dependent on the surface roughness of the nanocrystals is shown to be a decisive factor for a direct contact between the nanocrystals, manifested by the roughness-dependent hydrodynamic resistivity and Langevin dynamic simulations. Our study provides direct evidence that the role of dynamics is equally important in understanding the self-assembly, which has been rarely observed in the self-assembly in contrast with many studies on the equilibrium forces. These results may have further impact in other fields, such as having explanatory power for the probability of nonclassical crystal growth or the structure of viral spike proteins that affect the probability of viruses binding with cells.

Felsted, Robert G. [Univ. of Washington, Seattle, ↗
NASA NTRS1984

The interplanetary electric field, cleft currents and plasma convection in the polar caps

The relationship between the pattern of plasma convection in the polar cleft and the dynamics of the interplanetary electric field (IEF) is examined theoretically. It is shown that owing to the geometrical properties of the magnetosphere, the East-West component of the IEF will drive field-aligned currents which connect to the ionosphere at points lying on either side of noon, while currents associated with the North-South component of the IEF will connect the two polar caps as sheet currents, also centered at 12 MLT. In order to describe the consequences of the Interplanetary Magnetic Field (IMF) effects upon high-latitude electric fields and convection patterns, a series of numerical simulations was carried out. The simulations were based on a solution to the steady-state equation of current continuity in a height-integrated ionospheric current. The simulations demonstrate that a simple hydrodynamical model can account for the narrow 'throats' of strong dayside antisunward convection observed during periods of southward interplanetary IMF drift, as well as the sunward convection observed during periods of strongly northward IMF drift.

Banks, P. M.↗