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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Drivers and impacts of sediment deposition in Amazonian floodplains

Abstract The Amazon River carries enormous amounts of sediment from the Andes mountains, much of which is deposited in its floodplains. However, accurate quantification of the sediment sink at fine spatiotemporal scales is still challenging. Here, we present a high-resolution hydrodynamic-sediment model to simulate sediment deposition in a representative Amazon/Solimões floodplain. The process is found to be jointly driven by inundation, suspended sediment concentration in the Amazon River, and floodplain hydrodynamics and only weakly correlated with inundation level. By upscaling the sediment deposition rate (1.33 ± 0.24 kg m −2 yr −1 ), we estimate the trapping of 77.3 ± 13.9 Mt (or 6.1 ± 1%) of the Amazon River sediment by the Amazon/Solimões floodplains every year. Widespread deforestation would reduce the trapping efficiency of the floodplains over time, exacerbating downstream river aggradation. Additionally, we show that the deposition of sediment-associated organic carbon plays a minor role in fueling carbon dioxide and methane emissions in the Amazon.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Impacts of Water Clarity Variability on Temperature and Biogeochemistry in the Chesapeake Bay

Estuarine water clarity depends on the concentrations of aquatic constituents, such as colored dissolved organic matter, phytoplankton, inorganic suspended solids, and detritus, which are influenced by variations in riverine inputs. These constituents directly affect temperature because when water is opaque, sunlight heats a shallower layer of the water compared to when it is clear. Despite the importance of accurately predicting temperature variability, many numerical modeling studies do not adequately account for this key process. In this study, we quantify the effect of water clarity on heating by comparing two simulations of a hydrodynamic-biogeochemical model of the Chesapeake Bay for the years 2001–2005, in which (1) water clarity is constant in space and time for the computation of solar heating, compared to (2) a simulation where water clarity varies with modeled concentrations of light-attenuating materials. In the variable water clarity simulation, the water is more opaque, particularly in the northern region of the Bay. This decrease in water clarity reduces the total heat, phytoplankton, and nitrate throughout the Bay. During the spring and summer months, surface temperatures in the northern Bay are warmer by 0.1 °C and bottom temperatures are colder by 0.2 °C in the variable light attenuation simulation. Warmer surface temperatures encourage phytoplankton growth and nutrient uptake near the head of the Bay, and fewer nutrients are transported downstream. These impacts are greater during higher river flow years, when differences in temperature, nutrients, phytoplankton, and zooplankton extend further seaward compared to other years. This study demonstrates the consequences of utilizing different light calculations for estuarine heating and biogeochemistry.

water clarity↗
NASA NTRS1994

A simulation of the intracluster medium with feedback from cluster galaxies

We detail method and report first results from a three-dimensional hydrodynamical and N-body simulation of the formation and evolution of a Coma-sized cluster of galaxies, with the intent of studying the history of the hot, X-ray emitting intracluster medium. Cluster gas, galaxies, and dark matter are included in the model. The galaxies and dark matter fell gravitational forces; the cluster gas also undergoes hydrodynamical effects such as shock heating and PdV work. For the first time in three dimensions, we include modeling of ejection of processed gas from the simulated galaxies by winds, including heating and heavy element enrichment. For comparison, we employ a `pure infall' simulation using the same initial conditions but with no galaxies or winds. We employ an extreme ejection history for galactic feedback in order to define the boundary of likely models. As expected, feedback raises the entropy of the intracluster gas, preventing it from collapsing to densities as high as those attained in the infall model. The effect is more pronounced in subclusters formed at high redshift. The cluster with feedback is always less X-ray luminous, but experiences more rapid luminosity evolution, than the pure infall cluster. Even employing an extreme ejection model, the final gas temperature is only approximately 15% larger than in the infall model. The radial temperature profile is very nearly isothermal within 1.5 Mpc. The cluster galaxies in the feedback model have a velocity dispersion approximately 15% lower than the dark matter. This results in the true ratio of specific energies in galaxies to gas being less than one, beta(sub spec) approximately 0.7. The infall model predicts beta(sub spec) approximately 1.2. Large excursions in these values occur over time, following the complex dynamical history of the cluster. The morphology of the X-ray emission is little affected by feedback. The emission profiles of both clusters are well described by the standard beta-model with beta(sub fit) approximately equal to 0.7 - 0.9. X-ray mass estimates based on the assumptions of hydrostatic equilibrium and the applicability of the beta-model are quite accurate in both cases. A strong, radial iron abundance gradient is present, which develops as a consequence of the steepening of the galaxy density profile over time. Spectroscopic observations using nonimaging detectors with wide (approximately 45 min) fields of view dramatically smear the gradient. Observations with arcminute resolution, made available with the ASCA satellite, would readily resolve the gradient.

Metzler, Christopher A.↗
DOE OSTI2025

Reaction-in-flight neutrons as a diagnostic for hydrodynamical mixing in double shell inertial confinement fusion capsules

We examine reaction-in-flight (RIF) neutrons as diagnostics for hydrodynamical mixing of high-Z shell material into the hotspot of double shell capsules that are designed for the National Ignition Facility. In particular, we consider the effects of different levels of mixing of tungsten shell material into the DT gas on RIF spectra. Using a set of 1D simulations from the radiation hydrodynamic code xRAGE to determine the temperature and density profiles of the mixed W-DT gas, we find that increasing the mass of mix systematically reduces the ratio of RIF neutrons to primary 14 MeV neutrons. The shape of the RIF spectrum also changes with mix, with the predicted spectrum softening in energy with increasing mix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

Context. Solar flares are the result of the sudden release of magnetic energy in the corona. Much of this energy goes into accelerating charged particles to high velocity. These particles travel along the magnetic field and the energy is dissipated when the density gets high enough, primarily in the solar chromosphere. Modelling this region is difficult because the radiation energy balance is dominated by strong, optically thick spectral lines. Aims. Our aim is to provide the community with realistic simulations of a flaring loop with an emphasis on the detailed treatment of the chromospheric energy balance. This will enable a detailed comparison of existing and upcoming observations with synthetic observables from the simulations, thereby elucidating the complex interactions in a flaring chromosphere. Methods. We used the 1D radiation hydrodynamics code RADYN to perform simulations of the effect of a beam of electrons injected at the apex of a solar coronal loop. A grid of models was produced, varying the total energy input, the steepness, and low-energy cutoff of the beam energy spectrum. Results. The full simulation results for a grid of models are made available online. Some general properties of the simulations are discussed.

hydrodynamics↗
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗
NASA NTRS1974

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗
DOE OSTI2024

Enhanced sensitivity to target offset when using cross-beam energy transfer mitigation techniques in direct-drive inertial confinement fusion implosions

In direct-drive inertial confinement fusion, target offset from the target chamber center (or center of beam convergence) may lead to significant implosion asymmetry and fusion yield degradation. In addition, cross-beam energy transfer (CBET) has been shown to be a significant source of laser energy scattering and leads to a reduction in implosion velocity and yield. To improve energy coupling and implosion performance, several techniques for CBET mitigation have been proposed. Recent simulations, however, have shown that CBET also substantially mitigates the effect of target offset on implosion asymmetry and yield [Anderson et al., Phys. Plasmas 27, 112713 (2020)]. Furthermore, the inclusion of CBET models in radiation-hydrodynamics codes was shown to greatly improve agreement between simulations and experiments involving substantial target offset distances. This paper explores the intensity dependence of this CBET–offset effect. In addition, it is shown that enhanced sensitivity to target offset can be expected when CBET-mitigation techniques are used in direct-drive implosions. This is shown through simulations of two such CBET-mitigation techniques on the OMEGA laser: (1) decreased beam-to-target radius, and (2) beam-to-beam frequency detuning. For the typical target offset distances (<15 μm) observed in experiments on OMEGA, however, overall yield is still anticipated to be substantially higher when CBET-mitigation techniques are employed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
DOE OSTI2026

Quark flavor equilibration of the quark-gluon plasma

The early stage of a heavy-ion collision is marked by rapid entropy production and the transition from a gluon saturated initial condition to a plasma of quarks and gluons that evolves hydrodynamically. However, during the early times of the hydrodynamic evolution, the chemical composition of the QCD medium is still largely unknown. We present a study of quark chemical equilibration in the (Q)GP using a novel model of viscous hydrodynamic evolution in partial chemical equilibrium. Motivated by the success of gluon saturated initial condition models, we initialize the QCD medium as a completely gluon dominated state. Local quark production during the hydrodynamic phase is then simulated through the evolution of time-dependent fugacities for each independent quark flavor, with the timescales set as free parameters to compare different rates of equilibration. We present the results of complete heavy-ion collision simulations using this partial chemical equilibrium model, and show the effects on hadronic and electromagnetic observables. In particular, we show that the development of flow is sensitive to the equilibration timescale, providing an empirical way to probe the chemical equilibration of the QCD medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2013

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗
DOE OSTI2025

Hydrodynamic characterization of the redox chemistry of crown-encapsulated uranyl complexes

The redox properties of actinide-containing species strongly influence their reactivity, speciation, and interfacial behavior, but the experimental quantification of the electrochemical characteristics of molecular actinide complexes in nonaqueous media has not received the attention it deserves. Here, results from hydrodynamic methods and electrochemical simulations of U(VI)/U(V) redox are reported, including quantification of heterogeneous electron-transfer kinetics and estimation of chemical reversibility of U(VI)/U(V) interconversion at electrodes in acetonitrile-based electrolyte. The complexes investigated are recently reported U(VI) and U(V) complexes in which the uranyl ion (UO 2 n+ ) is encapsulated in a macrocyclic 18-crown-6-like moiety templated by a Pt(II) center. These complexes feature the most positive value U VI /U V reduction potential yet reported and are thus particularly relevant to study of facile U(V) generation from U(VI) precursors as well as uranium electroanalysis. Rotating disk electrode (RDE) studies have been used to quantify the diffusion coefficients of the U(VI) and U(V) complexes, and standard heterogeneous electron transfer rate constants ( k 0 ) for the redox have been determined using a conventional Koutecký-Levich analysis. Rotating ring-disk electrode (RRDE) studies have been used to directly interrogate the chemical reversibility of U(VI)-U(V) interconversion, confirming that reduction of the U(VI) complex at an Au disk is associated with formation of the U(V) analogue that can be readily re-oxidized at a Pt ring under hydrodynamic (rotating) conditions. Because measurements of the type reported here are generally associated with current flows that are larger than those found in corresponding quiescent (unstirred) conditions, our findings suggest that hydrodynamic methods could be advantageous for design of electroanalytical approaches to detection of actinide species and study of their redox properties.

actinides↗
NASA NTRS1983

Modeling of transient disturbances in coronal-streamer configurations

Numerical simulations of the formation and propagation of mass ejection, loop transients in coronal streamers are discussed. The simulations are accomplished with numerical solutions of the single fluid, ideal MHD equations of motion in the meridional plane. The streamer is produced by simulating the relaxation of an initially radial hydrodynamic flow coupled with a dipole magnetic field. The simulated transient then results from an energy release at the base of the streamer. The legs of the loop transient produced remain essentially stationary while the loop expands mainly in the radial direction with velocities of 400 to 750 km s-1. Once the leading edge of the transient has passed out of the lower corona, the initial streamer configuration is restored after 15 to 24 hours. A second energy release 2 hours later than, and with an energy release identical to, the first does not produce a significant coronal disturbance.

Steinolfson, R. S.↗
DOE OSTI2024

Effects of stratification on overshooting and waves atop the convective core of M ⊙ main-sequence stars

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a $5\,\mathrm{ M}_{\odot }$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evidence of a toroidal magnetic field in the core of 3C 84

The spatial scales of relativistic radio jets, probed by relativistic magneto-hydrodynamic (RMHD) jet launching simulations and by most very long baseline interferometry (VLBI) observations differ by an order of magnitude. Bridging the gap between these RMHD simulations and VLBI observations requires selecting nearby active galactic nuclei (AGN), the parsec-scale region of which can be resolved. The radio source 3C 84 is a nearby bright AGN fulfilling the necessary requirements: it is launching a powerful, relativistic jet powered by a central supermassive black hole, while also being very bright. Using 22 GHz globe-spanning VLBI measurements of 3C 84 we studied its sub-parsec region in both total intensity and linear polarisation to explore the properties of this jet, with a linear resolution of ~0.1 parsec. We tested different simulation set-ups by altering the bulk Lorentz factor Γ of the jet, as well as the magnetic field configuration (toroidal, poloidal, helical). We confirm the persistence of a limb brightened structure, which reaches deep into the sub-parsec region. The corresponding electric vector position angles (EVPAs) follow the bulk jet flow inside but tend to be orthogonal to it near the edges. Our state-of-the-art RMHD simulations show that this geometry is consistent with a spine-sheath model, associated with a mildly relativistic flow and a toroidal magnetic field configuration.

3C 84 (NGC 1275)↗
DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗