SEARCH · Search NASA

Results for “Statistical hydrodynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

NASA NTRS1987

Collapse of 9 solar mass stars

General relativistic hydrodynamical calculations of the collapse of O + Ne + Mg cores of a 9 solar mass star are reported. Collapse is induced by rapid electron captures as the O + Ne + Mg is burned to nuclear statistical equilibrium. The high entropy in the core after burning leads to a large abundance of free protons which readily capture electrons. This leads to large neutrino losses and a correspondingly small infalling homologous core. The hydrodynamic shock thus forms at a small mass point. The shock stalls before reaching the edge of the O + Ne + Mg core and thereby fails to produce a successful supernova explosion by the direct mechanism. No enhancement in the shock energy due to nuclear burning is found.

Baron, E.↗

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↗

NASA NTRS1994

Multivariate space - time analysis of PRE-STORM precipitation

This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Polyak, Ilya↗

NASA NTRS1997

Method for Reducing Pumping Damage to Blood

Methods are provided for minimizing damage to blood in a blood pump wherein the blood pump comprises a plurality of pump components that may affect blood damage such as clearance between pump blades and housing, number of impeller blades, rounded or flat blade edges, variations in entrance angles of blades, impeller length, and the like. The process comprises selecting a plurality of pump components believed to affect blood damage such as those listed herein before. Construction variations for each of the plurality of pump components are then selected. The pump components and variations are preferably listed in a matrix for easy visual comparison of test results. Blood is circulated through a pump configuration to test each variation of each pump component. After each test, total blood damage is determined for the blood pump. Preferably each pump component variation is tested at least three times to provide statistical results and check consistency of results. The least hemolytic variation for each pump component is preferably selected as an optimized component. If no statistical difference as to blood damage is produced for a variation of a pump component, then the variation that provides preferred hydrodynamic performance is selected. To compare the variation of pump components such as impeller and stator blade geometries, the preferred embodiment of the invention uses a stereolithography technique for realizing complex shapes within a short time period.

Bozeman, Richard J., Jr.↗

DOE OSTI2025

Hypertriton puzzle in relativistic heavy-ion collisions

The yields of hadrons and light nuclei in relativistic collisions of heavy-nuclei at a center of mass energy of $\sqrt{s_{NN}}$ = 2.6 TeV can be described remarkably well by a thermal distribution of an ideal gas of hadrons and light nuclei interacting only via the decay of resonances. Given the particularly small binding energy of hypertritons relative to the temperature describing the yields (about 156 MeV), one might naturally expect hypertrions to dissociate in medium, making the agreement of hypertriton yields with thermal predictions highly puzzling. The puzzle is compounded by the fact that small binding energy is associated with the large size of the hypertriton. This size is on a similar scale to the overall size of the fireball and much larger than the length scale over which temperatures in the fireball vary over phenomenologically relevant amounts. Here, this paper quantifies the tension this effect causes and shows that it is sufficiently large to render the thermal model inconsistent: its natural assumptions are in conflict with its outputs. The possibility that hypertritons are formed at freeze out as compact objects, quark droplets, that subsequently evolve into hypertritons is considered as a way to resolve the puzzle. It is noted that beyond making the assumption that compact quark droplets form, additional detailed dynamical assumptions which have not been justified are needed to make the thermal model work. The issue of why, despite these issues, the hypertriton is well described by a simple statistical description at freeze out is unresolved. Resolving the hypertriton puzzle is important as it may clarify whether the phenomenological success of the simple thermal model for yields accurately reflects the simple picture of the underlying physics on which it is based.

hydrodynamic models↗

DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗

DOE OSTI2024

Impact of baryonic feedback on HSC-Y1 weak lensing non-Gaussian statistics

Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude S 8 = σ 8 Ω m / 0.3 . We focus on the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) data which, with its high source number density, closely resembles those expected from the upcoming Euclid and Rubin Legacy Survey of Space and Time. We study four non-Gaussian statistics of convergence maps—peak counts, minimum counts, the probability distribution function, and the scattering transform coefficients—in addition to the usual power spectrum. We first estimate the biases in S 8 using mock observations built from the IllustrisTNG and BAHAMAS hydrodynamical simulations and theoretical models built from dark-matter-only simulations. We find up to 1 σ bias in S 8 when the smallest scales (2 arcmin) and the highest feedback level are considered. We then analyze the HSC-Y1 data and compare the S 8 obtained for each statistic with different smoothing scales or scale cuts. As we expect that baryons mostly affect the small scales, comparing the results obtained from including and excluding small scales can indicate the level of impact from baryons. With HSC data, we find only minor ( ≤ 0.5 σ ) differences in S 8 for all statistics, even when considering very small scales (2 arcmin). Our results suggest that the effect of baryons is insignificant at the level of HSC-Y1 down to 2 arcmin for all statistics examined here, or it is canceled by other scale-dependent systematics.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2008

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

NASA NTRS2012

Experimenting with the GMAO 4D Data Assimilation

The Global Modeling and Assimilation Office (GMAO) has been working to promote its prototype four-dimensional variational (4DVAR) system to a version that can be exercised at operationally desirable configurations. Beyond a general circulation model (GeM) and an analysis system, traditional 4DV AR requires availability of tangent linear (TL) and adjoint (AD) models of the corresponding GeM. The GMAO prototype 4DVAR uses the finite-volume-based GEOS GeM and the Grid-point Statistical Interpolation (GSI) system for the first two, and TL and AD models derived ITom an early version of the finite-volume hydrodynamics that is scientifically equivalent to the present GEOS nonlinear GeM but computationally rather outdated. Specifically, the TL and AD models hydrodynamics uses a simple (I-dimensional) latitudinal MPI domain decomposition, which has consequent low scalability and prevents the prototype 4DV AR ITom being used in realistic applications. In the near future, GMAO will be upgrading its operational GEOS GCM (and assimilation system) to use a cubed-sphere-based hydrodynamics. This versions of the dynamics scales to thousands of processes and has led to a decision to re-derive the TL and AD models for this more modern dynamics, thus taking advantage of a two-dimensional MPI decomposition and improved scalability properties. With the aid of the Transformation of Algorithms in FORTRAN (l'AF) automatic adjoint generation tool and some hand-coding, a version of the cubed-sphere-based TL and AD models, with a simplified vertical diffusion scheme, is now available, enabling multiple configurations of standard implementations of 4DV AR in GEOS. Concurrent to this development, collaboration with the National Centers for Environmental Prediction (NCEP) and the Earth System Research Laboratory (ESRL) has allowed GMAO to implement a hybrid-ensemble capability within the GEOS data assimilation system. Both 3Dand 4D-ensemble capabilities are presently available thus allowing GMAO to now evaluate the performance and benefit of various ensemble and variational assimilation strategies. This presentation will cover the most recent developments taking place at GMAO and show results from various comparisons from traditional techniques to more recent ensemble-based ones.

Todling, R.↗

NASA NTRS1997

Reactive Hydrodynamics in Rotating Spherical and Cylindrical Geometry

In turbulent spray combustion among many complex interactions between local flow structures called turbulent eddies and droplets are those associated with rotation of droplets. In general, for a complete statistical description of turbulent sprays, consideration of at least four degrees of freedom respectively associated with translational, rotational, vibrational (pulsational), and internal motions of the droplet are needed. Clearly the interactions between all degrees of freedom of the droplets and those for the gaseous background field will be exceedingly complex. For example, one type of interaction between the translational and the rotational velocity of droplets results in droplet helicity, H(d) = w(d).v(d), the significance of which in turbulent spray combustion is yet to be recognized. The role of droplet rotation in turbulent spray combustion modeling and its impact on the evaporation of liquid fuel droplets was recently investigated. Also, the impact of rotation on combustion of solid particles such as is encountered in pulverized coal combustion has been emphasized. The problem of viscous flow around a rotating sphere discussed above also occurs in other areas of physical sciences such as astrophysics and geophysics. Consequently, the subject has been addressed in many classical as well as more recent investigations. According to these investigations, the rotation of a rigid sphere in an otherwise quiescent, unconfined environment results in the motion of the fluid towards the poles. The polar flows from the northern and southern hemispheres move along helical trajectories towards the equatorial plane. Eventually, the polar flows collide at the equatorial plane, thus producing a sheet of rotating fluid that is radially ejected outward on this plane. Therefore, a droplet induces a strained flow field as a result of its rotation. Since the spatial extent of equatorial jets could easily exceed many droplet diameters, interactions between neighboring droplets are enhanced as a result of their rotation. Also, the equatorial jet substantially alters the spherical geometry of the diffusion flame surface that surrounds a rotating droplet. The objective of the research is to gain more knowledge about the hydrodynamics within and around rotating spherical and cylindrical body of fluid, and the behavior of diffusion or premixed flame surfaces that could surround such symmetric body of rotating fluid.

Sohrab, Siavash H.↗

DOE OSTI2026

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗

NASA NTRS1996

Chiral Symmetry Breaking in Crystal Growth: Is Hydrodynamic Convection Relevant?

The effects of mechanical stirring on nucleation and chiral symmetry breaking have been investigated for a simple inorganic molecule, sodium chlorate (NaClO3). In contrast to earlier findings, our experiment suggests that the symmetry breaking may have little to do with hydrodynamic convection. Rather the effect can be reasonably accounted for by mechanical damage to incipient crystals. The catastrophic events, creating numerous small 'secondary' crystals, produce statistical domination of one chiral species over the other. Our conclusion is supported by a number of observations using different mixing mechanisms.

Martin, B.↗

NASA NTRS1992

Statistics of the cosmic Mach number from numerical simulations of a cold dark matter universe

Results are presented of an analysis of the cosmic Mach number, M, the ratio of the streaming velocity, v, to the random velocity dispersion, sigma, of galaxies in a given patch of the universe, which was performed on the basis of hydrodynamical simulations of the cold dark matter scenario. Galaxy formation is modeled by application of detailed physical processes rather than by the ad hoc assumption of 'bias' between dark matter and galaxy fluctuations. The correlation between M and sigma is found to be very weak for both components. No evidence is found for a physical 'velocity bias' in the quantities which appear in the definition of M. Standard cold-dark-matter-dominated universes are in conflict, at a statistically significant level, with the available observation, in that they predict a Mach number considerably lower than is observed.

Suto, Yasushi↗

NASA NTRS2009

The GMAO ADVAR System: Preliminary Results

The Goddard Earth Observing System (GEOS-S) Data Assimilation System is now configurable to run strong constraint 4DVAR. The system combines the 4DVAR-capable Grid-point Statistical Interpolation analysis with the GEOS-S general circulation model and an early version of the tangent linear and adjoint models of the finite-volume hydrodynamics with simplified physics. The GMAO 4DVAR uses a Lanczos-based conjugate gradient algorithm in a nested resolution inner-loop setting. A range of tests and experiments will have been carried out by the time of the Symposium. Results on balance-related issues, resolution configuration, and choice of assimilation time window will be shown. A comparison with the Global Modeling and Assimilation Office (GMAO) 3DVAR and a discussion of future directions of development will also be presented.

Todling, Ricardo↗

DOE OSTI2025

Cosmology with second- and third-order shear statistics for the Dark Energy Survey: Methods and simulated analysis

We present a new pipeline designed for the robust inference of cosmological parameters using both second- and third-order shear statistics. We build a theoretical model for rapid evaluation of three-point correlations using our fastnc code and integrate it into the cosmosis framework. We measure the two-point functions 𝜉 ± and the full configuration-dependent three-point shear correlation functions across all auto- and cross-redshift bins. We compress the three-point functions into the mass aperture statistic ⟨ℳ$^{3}_{ap}$⟩ for a set of 796 simulated shear maps designed to model the Dark Energy Survey Year 3 data. We estimate from it the full covariance matrix and model the effects of intrinsic alignments, shear calibration biases and photometric redshift uncertainties. We apply scale cuts to minimize the contamination from the baryonic signal as modeled through hydrodynamical simulations. We find a significant improvement of 83% on the figure of merit in the Ω m − 𝑆 8 plane when we add the ⟨ℳ$^{3}_{ap}$⟩ data to 𝜉 ± . Here, we present our findings for all relevant cosmological and systematic uncertainty parameters and discuss the complementarity of third-order and second-order statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel’dovich effect (kSZ). We precisely measure this effect (𝑆/𝑁 ≈ 10) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over ≥ 4,000 deg 2 . We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in 𝑔, 𝑟, and 𝑧 bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed (PTE)⁢ = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris 𝑧 = 0.5, PTE⁢ = 0.37) over weaker-feedback cases (IllustrisTNG 𝑧 = 0.8, PTE⁢ = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast 𝑆/𝑁 ∼ 50 for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. Furthermore, these measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

The progenitors and products of supernovae

Much of the recent radio and X-ray observational work relevant to the questions of the progenitor populations and explosion products of supernovae is reviewed. Several illustrative examples are given of how X-ray data on shell-type supernova remnants, coupled with analytic and numerical hydrodynamic models, can constrain the mass and evolutionary state of the stellar precursor. The population of Crab-like remnants is examined, pointing out the restrictions imposed on explosion models by the statistics of collapsed remnant cores.

Helfand, David J.↗