Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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

Thermonuclear 28 P(p, γ ) 29 S reaction rate and astrophysical implication in ONe nova explosion

An accurate 28 P(p, γ) 29 S reaction rate is crucial to defining the nucleosynthesis products of explosive hydrogen burning in ONe novae. Using the recently released nuclear mass of 29 S, together with a shell model and a direct capture calculation, we reanalyzed the 28 P(p, γ) 29 S thermonuclear reaction rate and its astrophysical implication. We focus on improving the astrophysical rate for 28 P(p, γ) 29 S based on the newest nuclear mass data. Our goal is to explore the impact of the new rate and associated uncertainties on the nova nucleosynthesis. We evaluated this reaction rate via the sum of the isolated resonance contribution instead of the previously used Hauser-Feshbach statistical model. The corresponding rate uncertainty at different energies was derived using a Monte Carlo method. Nova nucleosynthesis is computed with the 1D hydrodynamic code SHIVA. The contribution from the capture on the first excited state at 105.64 keV in 28 P is taken into account for the first time. We find that the capture rate on the first excited state in 28 P is up to more than 12 times larger than the ground-state capture rate in the temperature region of 2.5 × 10 7 K to 4 × 10 8 K, resulting in the total 28 P(p, γ) 29 S reaction rate being enhanced by a factor of up to 1.4 at ~1 × 10 9 K. In addition, the rate uncertainty has been quantified for the first time. It is found that the new rate is smaller than the previous statistical model rates, but it still agrees with them within uncertainties for nova temperatures. The statistical model appears to be roughly valid for the rate estimation of this reaction in the nova nucleosynthesis scenario. Using the 1D hydrodynamic code SHIVA, we performed the nucleosynthesis calculations in a nova explosion to investigate the impact of the new rates of 28 P(p, γ) 29 S. Our calculations show that the nova abundance pattern is only marginally affected if we use our new rates with respect to the same simulations but statistical model rates. Finally, the isotopes whose abundance is most influenced by the present 28 P(p, γ) 29 S uncertainty are 28 Si, 33,34 S, 35,37 Cl, and 36 Ar, with relative abundance changes at the level of only 3% to 4%.

Astronomy & Astrophysics↗
DOE OSTI2023

Sucrose transport inside the phloem: Bridging hydrodynamics and geometric characteristics

In plants, the delivery of the products of photosynthesis is achieved through a hydraulic system labeled as phloem. This semi-permeable plant tissue consists of living cells that contract and expand in response to fluid pressure and flow velocity fluctuations. The Münch pressure flow theory, which is based on osmosis providing the necessary pressure gradient to drive the mass flow of carbohydrates, is currently the most accepted model for such sucrose transport. When this hypothesis is combined with the conservation of fluid mass and momentum as well as sucrose mass, many simplifications must be invoked to mathematically close the problem and to resolve the flow. Here this study revisits such osmotically driven flows by developing a new two-dimensional numerical model in cylindrical coordinates for an elastic membrane and a concentration-dependent viscosity. It is demonstrated that the interaction between the hydrodynamic and externally supplied geometrical characteristic of the phloem has a significant effect on the front speed of sucrose transport. These results offer a novel perspective about the evolutionary adaptation of plant hydraulic traits to optimize phloem soluble compounds transport efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Simultaneous inference of equation of state parameters and unknown data errors with uncertainty quantification via hierarchical Bayesian posterior maximization

Equations of state (EOSs) are a key component in running hydrodynamic simulations as they relate the thermodynamic states for the material. The Davis reactants EOS is commonly used for modeling high explosives (HEs), and the EOS model parameters are calibrated using material specific data. The calibrations are often performed with uncertainty quantification via Bayesian inference to account for uncertainty in the data and generate ensembles of likely parameters. However, there are relatively few HE data sets to use for calibration and many are historical and lack error information. In this work, we simultaneously calibrate the Davis reactants EOS model parameters and unknown data error terms for the high explosive PBX 9501. To quantify the uncertainty in the models and the data, we use a Bayesian framework for the calibration and compute the hierarchical Bayesian posterior distribution with both a posteriori maximization approach and Markov Chain Monte Carlo. In general, we find that, given our assumptions, the two approaches result in similar calibrated parameters, posterior covariance matrices, and insights about the parameters but that the posterior maximization requires far less computational resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Galaxy mergers in eagle do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

ABSTRACT We investigate the connection between galaxy–galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical eagle simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems, and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of ≈3 depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass (M* ∼ 1010 M⊙) gas rich (fgas > 0.2) central galaxies with lower mass BHs (MBH ∼ 107 M⊙) at lower redshifts (z < 1). We find no enhancement of AGN triggered via mergers in more massive galaxies (M* ≳ 1011 M⊙). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early remnants of merging systems (typically lagging ≈300 Myr post-coalescence of the two galaxies at z = 0.5). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \ge \frac{1}{4}$) nor minor mergers ($\frac{1}{10} \lt M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \lt \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2–3 times that of isolated galaxies, the majority of mass that is accreted on to BHs occurs outside the periods of a merger. We compute that on average no more than 15 per cent of a BHs final day mass comes from the enhanced accretion rates triggered via a merger.

McAlpine, Stuart↗
DOE OSTI2021

The protoplanetary disc around HD 169142: circumstellar or circumbinary?

ABSTRACT Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas distributions. In this work, we suggest that the structures shown in the central region of the protoplanetary disc HD 169142 are produced by the presence of an inner stellar binary and a circumbinary (P-type) planet. We find that a companion with a mass ratio of 0.1, semimajor axis of 9.9 au, eccentricity of 0.2, and inclination of 90°, together with a 2 MJ coplanar planet on a circular orbit at 45 au reproduce the structures at the innermost ring observed at 1.3 mm and the shape of spiral features in scattered light observations. The model predicts changes in the disc’s dust structure, and star’s astrometric parameters, which would allow testing its veracity by monitoring this system over the next 20 yr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 NTRS1980

Protostellar formation in rotating interstellar clouds. I - Numerical methods and tests

Attention is given to numerical methods and tests of a series of gravitational hydrodynamics computer codes constructed in order to numerically follow the dynamic collapse of spherically symmetric (1-D), axisymmetric (2-D), and non-axisymmetric (3-D) isothermal interstellar clouds. A spherical harmonic expansion is used to solve the Poisson equation for the gravitational potential. It is shown that the use of explicit donor-cell hydrodynamics on a moving spherical coordinate grid ensures mass and momentum conservation and allows the grid to follow the collapse of the fluid. Finally, the performance of the codes is examined.

Boss, A. P.↗
NASA NTRS1983

Filtering flight data prior to aerodynamic system identification

An algorithm for processing flight-test data to provide state estimates and instrument calibrations for aerodynamic or hydrodynamic system identification by the equation error estimation method is developed and demonstrated on synthesized data. The extended-Kalman-filter algorithm employs a locally level, north-pointing frame of reference, accounts for rotating ellipsoidal earth effects, and estimates sensor bias, scale factors, wind components, and process noise levels by maximum-likelihood parameters. The method is found to be most effective with navigation quality inertial input data. The algorithm is applied to data from a six-degree-of-freedom F-4 aircraft simulation and shown to produce state estimates in good agreement with the simulation values.

Trankle, T. L.↗
NASA NTRS1995

The Benard problem: A comparison of finite difference and spectral collocation eigen value solutions

The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Skarda, J. Raymond Lee↗
NASA NTRS1999

Photochemistry of Pluto's Atmosphere

This work include studies of two problems: (1) Modeling thermal balance, structure. and escape processes in Pluto's upper atmosphere. This study has been completed in full. A new method, of analytic solution for the equation of hydrodynamic flow from in atmosphere been developed. It was found that the ultraviolet absorption by methane which was previously ignored is even more important in Pluto's thermal balance than the extreme ultraviolet absorption by nitrogen. Two basic models of the lower atmosphere have been suggested, with a tropopause and a planetary surface at the bottom of the stellar occultation lightcurve, respectively, Vertical profiles, of temperature, density, gas velocity, and the CH4 mixing ratio have been calculated for these two models at low, mean, and high solar activity (six models). We prove that Pluto' " s atmosphere is restricted to 3060-4500 km, which makes possible a close flyby of future spacecraft. Implication for Pluto's evolution have also been discussed. and (2) Modeling of Pluto's photochemistry. Based on the results of (1), we have made some changes in the basic continuity equation and in the boundary conditions which reflect a unique can of hydrodynamic escape and therefore have not been used in modeling of other planetary atmospheres. We model photochemistry of 44 neutral and 23 ion species. This work required solution of a set of 67 second-order nonlinear ordinary differential equations. Two models have been developed. Each model consists of the vertical profiles for 67 species, their escape and precipitation rates. These models predict the chemical structure and basic chemical processes in the current atmosphere and possible implication of these processes for evolution. This study has also been completed in full.

Krasnopolsky, Vladimir A.↗
DOE OSTI2021

Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules

Advances in machine learning provide the ability to leverage data from expensive simulations of high-energy-density experiments to significantly cut down on computational time and costs associated with the search for optimal target designs. This study presents an application of cutting-edge Bayesian optimization methods to the one-dimensional (1D) design optimization of double shell graded layer targets for inertial confinement fusion experiments. This investigation attempts to reduce hydrodynamic instabilities while retaining high yields for future NIF experiments. Machine learning methods can use predictive physics simulations to identify graded layer designs from within the vast design space that demonstrate high predicted performance, including novel designs with high uncertainty in performance that may hold unexpected promise. By applying machine learning tools to the simulation design, we map the trade-off between 1D yield and instability, specifically isolating parameter ranges, which maintain high performance while showing significantly improved Rayleigh–Taylor stability over the point design. Furthermore, the groundwork laid in this study will be a useful design tool for future NIF experiments with graded layer targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗
DOE OSTI2023

Analytic insights into nonlocal energy transport: Steady state Fokker Planck theory in arbitrary Z plasmas

The generation of energetic electrons in laser fusion in an important issue. The electrons may either arise from a laser plasma instability, or from the uncoupled high temperature tail of a Maxwellian distribution. To study these in a laser fusion context, it is important to find a method accurate enough to be useful, and simple enough to be incorporated into a radiation hydrodynamics numerical simulation, the main workhorse for studying the laser fusion target. That is why analytic insights become important, they allow one to simplify the Fokker Planck theory so that a solution of it can be incorporated into a radiation hydrodynamic simulation. This work develops and analyzes a steady state Fokker Planck theory for plasmas of arbitrary Z. It develops a method of solving the simplified Fokker Planck method with a technique called sparse eigenfunction analysis. As a result, this method appears to work reasonably well when compared to the experimental results from the Rochester/NIF on plastic spherical targets with and without a silicon layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Film cooling on a convex wall: Heat transfer and hydrodynamic measurements for full and partial coverage

Turbine-blade cooling is an important issue for high-efficiency turbine engines, and discrete-hole injection is widely used as a cooling method. In the present study, detailed measurements were made of the heat transfer and hydrodynamics of a film-cooled flow on a convex wall, both for full and partial coverage. Two important parameters were altered: the blowing ratio, m, and the number of rows of injection holes. Three values of m were tested: m = 0.2, 0.4, and 0.6. In the blown region, m = 0.4 results in the lowest Stanton numbers of the three blowing ratios tested. This indicates that the value of m = 0.4 is near optimum on the convex wall from the point of view of cooling effect by injection. In the recovery region, Stanton numbers gradually approach the no injection values. Although the heat-transfer behavior during recovery from injection looks relatively complicated, the behavior of Stanton number can be explained in terms of two mechanisms: recovery from the thermal effect of injection and recovery from the turbulence augmentation. This interpretation of the data is supported by the hydrodynamic and temperture-profile measurements. For partial blowing cases, the data follow the full-coverage values inside the blown region. In the unblown region, both in the curved and in the flat plate, the effect of the number of blown rows is clearly seen. Hydrodynamic boundary-layer profiles were measured with the aid of a triple hot-water probe. Three mean-velocity components and six turbulence quantities were simultaneously measured, and inside the blown region strong three-dimensionality was observed.

Furuhama, K.↗