Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The optimization of high resolution topographic data for 1D hydrodynamic models

The main focus of our research presented in this paper is to optimize and use high resolution topographical data (HRTD) for hydrological modelling. Optimization of HRTD is done by generating adaptive mesh by measuring distance of coarse mesh and the surface of the dataset and adapting the mesh from the perspective of keeping the geometry as close to initial resolution as possible. Technique described in this paper enables computation of very accurate 1-D hydrodynamic models. In the paper, we use HEC-RAS software as a solver. For comparison, we have chosen the amount of generated cells/grid elements (in whole discretization domain and selected cross sections) with respect to preservation of the accuracy of the computational domain. Generation of the mesh for hydrodynamic modelling is strongly reliant on domain size and domain resolution. Topographical dataset used in this paper was created using LiDAR method and it captures 5.9km long section of a catchment of the river Olše. We studied crucial changes in topography for generated mesh. Assessment was done by commonly used statistical and visualization methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Fast, High-fidelity Lyα Forests with Convolutional Neural Networks

Full-physics cosmological simulations are powerful tools for studying the formation and evolution of structure in the universe but require extreme computational resources. Here, we train a convolutional neural network to use a cheaper N-body-only simulation to reconstruct the baryon hydrodynamic variables (density, temperature, and velocity) on scales relevant to the Lyα forest, using data from Nyx simulations. We show that our method enables rapid estimation of these fields at a resolution of ~20 kpc, and captures the statistics of the Lyα forest with much greater accuracy than existing approximations. Because our model is fully convolutional, we can train on smaller simulation boxes and deploy on much larger ones, enabling substantial computational savings. Furthermore, as our method produces an approximation for the hydrodynamic fields instead of Lyα flux directly, it is not limited to a particular choice of ionizing background or mean transmitted flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Global Λ hyperon polarization in nuclear collisions

The extreme energy densities generated by ultra-relativistic collisions between heavy atomic nuclei produce a state of matter that behaves surprisingly like a fluid, with exceptionally high temperature and low viscosity. Non-central collisions have angular momenta of the order of 1,000ћ, and the resulting fluid may have a strong vortical structure that must be understood to describe the fluid properly. The vortical structure is also of particular interest because the restoration of fundamental symmetries of quantum chromodynamics is expected to produce novel physical effects in the presence of strong vorticity. But, no experimental indications of fluid vorticity in heavy ion collisions have yet been found. Since vorticity represents a local rotational structure of the fluid, spin–orbit coupling can lead to preferential orientation of particle spins along the direction of rotation. Here we present measurements of an alignment between the global angular momentum of a non-central collision and the spin of emitted particles (in this case the collision occurs between gold nuclei and produces Λ baryons), revealing that the fluid produced in heavy ion collisions is the most vortical system so far observed. (At high energies, this fluid is a quark–gluon plasma.) We find that Λ and hyperons show a positive polarization of the order of a few per cent, consistent with some hydrodynamic predictions. (A hyperon is a particle composed of three quarks, at least one of which is a strange quark; the remainder are up and down quarks, found in protons and neutrons.) A previous measurement that reported a null result, that is, zero polarization, at higher collision energies is seen to be consistent with the trend of our observations, though with larger statistical uncertainties. Furthermore, these data provide experimental access to the vortical structure of the nearly ideal liquid created in a heavy ion collision and should prove valuable in the development of hydrodynamic models that quantitatively connect observations to the theory of the strong force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Statistics of velocity fluctuations of Geldart A particles in a circulating fluidized bed riser

Here, the statistics of fluctuating velocity components are studied in the riser of a closed-loop circulating fluidized bed with fluid catalytic cracking catalyst particles. Our analysis shows distinct similarities as well as deviations compared to existing theories and bench-scale experiments. The study confirms anisotropic and non-Maxwellian distribution of fluctuating velocity components. The velocity distribution functions (VDFs) corresponding to transverse fluctuations exhibit symmetry, and follow a stretched-exponential behavior up to three standard deviations. The form of the transverse VDF is largely determined by interparticle interactions. The tails become more overpopulated with an increase in particle loading. The observed deviations from the Gaussian distribution are represented using the leading order term in the Sonine expansion, which is commonly used to approximate the VDFs in kinetic theory for granular flows. The vertical fluctuating VDFs are asymmetric and the skewness shifts as the wall is approached. In comparison to transverse fluctuations, the vertical VDF is determined by the local hydrodynamics. This is an observation of particle velocity fluctuations in a large-scale system and their quantitative comparison with the Maxwell-Boltzmann statistics.

42 ENGINEERING↗
DOE OSTI2014

LyMAS: Predicting large-scale Lyα forest statistics from the dark matter density field

We describe Lyα Mass Association Scheme (LyMAS), a method of predicting clustering statistics in the Lyα forest on large scales from moderate-resolution simulations of the dark matter (DM) distribution, with calibration from high-resolution hydrodynamic simulations of smaller volumes. We use the 'Horizon-MareNostrum' simulation, a 50 h {sup –1} Mpc comoving volume evolved with the adaptive mesh hydrodynamic code RAMSES, to compute the conditional probability distribution P(F{sub s} |δ {sub s}) of the transmitted flux F{sub s} , smoothed (one-dimensionally, 1D) over the spectral resolution scale, on the DM density contrast δ {sub s}, smoothed (three-dimensionally, 3D) over a similar scale. In this study we adopt the spectral resolution of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) at z = 2.5, and we find optimal results for a DM smoothing length σ = 0.3 h {sup –1} Mpc (comoving). In its simplest form, LyMAS draws randomly from the hydro-calibrated P(F{sub s} |δ {sub s}) to convert DM skewers into Lyα forest pseudo-spectra, which are then used to compute cross-sightline flux statistics. In extended form, LyMAS exactly reproduces both the 1D power spectrum and one-point flux distribution of the hydro simulation spectra. Applied to the MareNostrum DM field, LyMAS accurately predicts the two-point conditional flux distribution and flux correlation function of the full hydro simulation for transverse sightline separations as small as 1 h {sup –1} Mpc, including redshift-space distortion effects. It is substantially more accurate than a deterministic density-flux mapping ({sup F}luctuating Gunn-Peterson Approximation{sup )}, often used for large-volume simulations of the forest. With the MareNostrum calibration, we apply LyMAS to 1024{sup 3} N-body simulations of a 300 h {sup –1} Mpc and 1.0 h {sup –1} Gpc cube to produce large, publicly available catalogs of mock BOSS spectra that probe a large comoving volume. LyMAS will be a powerful tool for interpreting 3D Lyα forest data, thereby transforming measurements from BOSS and other massive quasar absorption surveys into constraints on dark energy, DM, space geometry, and intergalactic medium physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Effective field theory of dissipative fluids

We develop an effctive fi eld theory for dissipative fluids which governs the dynamics of long-lived gapless modes associated with conserved quantities. The resulting theory gives a path integral formulation of fluctuating hydrodynamics which systematically incorporates nonlinear interactions of noises. The dynamical variables are mappings between a "fluid spacetime" and the physical spacetime and an essential aspect of our formulation is to identify the appropriate symmetries in the fluid spacetime. The theory applies to nonlinear disturbances around a general density matrix. For a thermal density matrix, we require an additional Z2 symmetry, to which we refer as the local KMS condition. This leads to the standard constraints of hydrodynamics, as well as a nonlinear generalization of the Onsager relations. It also leads to an emergent supersymmetry in the classical statistical regime, and a higher derivative deformation of supersymmetry in the full quantum regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Scaling Invariance and Characteristics of the Cloud of Spherical Projectile Fragmentation Products upon High-Velocity Impact on a Thin Mesh Shield

In this paper, we consider the problem of the fragmentation of an aluminum projectile on a thin steel mesh shield at high-velocity impact in a three-dimensional (3D) setting. The numerical simulations are carried out by the smoothed particle hydrodynamics method applied to the equations of mechanics of deformable solids. Quantitative characteristics of the projectile fragmentation are obtained by studying statistics of the cloud of fragments. Considerable attention is given to scaling laws accompanying the fragmentation of the projectile. Scaling is carried out using the parameter K, which defines the number of the mesh cells in the projectile diameter. It is found that the dependence of the critical velocity V{sub c} of fragmentation on the parameter K consists of two branches that correspond to two modes of the projectile fragmentation associated with the “small” and “large” aperture of the mesh cell. We obtain the dependences of the critical velocity V{sub c} on the projectile diameter and the mesh parameters for both modes of the fragmentation. It is shown that the average cumulative mass distributions constructed at V{sub c} exhibit the property of scale invariance, splitting into two groups of distributions corresponding exactly to the modes of the projectile fragmentation. In each group, the average cumulative distributions show good coincidence in the entire mass region; moreover, in the intermediate mass region, each group of distributions has a power-law distribution with an exponent τ different from that in the other group. Conclusions about the dependence of the exponent of the power-law distribution τ on the fragmentation mode are made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2015

Characterizing the pressure smoothing scale of the intergalactic medium

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Lyα forest is degenerate with the 3D structure of the IGM on ∼100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T ≃ 10{sup 4} K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected filtering scale cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Lyα flux, F{sub real}, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The F{sub real} power spectrum is accurately described by a simple fitting function with cutoff at λ{sub F}, allowing us to rigorously quantify the pressure smoothing scale for the first time: we find λ{sub F} = 79 kpc (comoving) at z = 3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Lyα forest absorption in close quasar pairs, recently proposed as a method to measure the pressure smoothing scale. Our results enable one to quantify the pressure smoothing scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude T{sub 0} and slope γ of the temperature–density relation T=T{sub 0}(ρ/ ρ-bar ){sup γ−1} should be augmented with a third pressure smoothing scale parameter λ{sub F}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Elucidating size effects on the yield strength of single-crystal Cu via the Richtmyer–Meshkov instability

Capturing the dynamic response of a material under high strain-rate deformation often demands challenging and time consuming experimental effort. While shock hydrodynamic simulation methods can aid in this area, a priori characterizations of the material strength under shock loading and spall failure are needed in order to parameterize constitutive models needed for these computational tools. Moreover, parameterizations of strain-rate-dependent strength models are needed to capture the full suite of Richtmyer–Meshkov instability (RMI) behavior of shock compressed metals, creating an unrealistic demand for these training data solely on experiments. Herein, we sweep a large range of geometric, crystallographic, and shock conditions within molecular dynamics (MD) simulations and demonstrate the breadth of RMI in Cu that can be captured from the atomic scale. In this work, yield strength measurements from jetted and arrested material from a sinusoidal surface perturbation were quantified as Y RMI = 0.787 ± 0.374 GPa, higher than strain-rate-independent models used in experimentally matched hydrodynamic simulations. Defect-free, single-crystal Cu samples used in MD will overestimate Y RMI , but the drastic scale difference between experiment and MD is highlighted by high confidence neighborhood clustering predictions of RMI characterizations, yielding incorrect classifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2016

IMPROVED PERFORMANCES IN SUBSONIC FLOWS OF AN SPH SCHEME WITH GRADIENTS ESTIMATED USING AN INTEGRAL APPROACH

In this paper, we present results from a series of hydrodynamical tests aimed at validating the performance of a smoothed particle hydrodynamics (SPH) formulation in which gradients are derived from an integral approach. We specifically investigate the code behavior with subsonic flows, where it is well known that zeroth-order inconsistencies present in standard SPH make it particularly problematic to correctly model the fluid dynamics. In particular, we consider the Gresho–Chan vortex problem, the growth of Kelvin–Helmholtz instabilities, the statistics of driven subsonic turbulence and the cold Keplerian disk problem. We compare simulation results for the different tests with those obtained, for the same initial conditions, using standard SPH. We also compare the results with the corresponding ones obtained previously with other numerical methods, such as codes based on a moving-mesh scheme or Godunov-type Lagrangian meshless methods. We quantify code performances by introducing error norms and spectral properties of the particle distribution, in a way similar to what was done in other works. We find that the new SPH formulation exhibits strongly reduced gradient errors and outperforms standard SPH in all of the tests considered. In fact, in terms of accuracy, we find good agreement between the simulation results of the new scheme and those produced using other recently proposed numerical schemes. These findings suggest that the proposed method can be successfully applied for many astrophysical problems in which the presence of subsonic flows previously limited the use of SPH, with the new scheme now being competitive in these regimes with other numerical methods.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Pilot-wave hydrodynamics in a rotating frame: Exotic orbits

We present the results of a numerical investigation of droplets walking on a rotating vibrating fluid bath. The drop's trajectory is described by an integro-differential equation, which is simulated numerically in various parameter regimes. As the forcing acceleration is progressively increased, stable circular orbits give way to wobbling orbits, which are succeeded in turn by instabilities of the orbital center characterized by steady drifting then discrete leaping. In the limit of large vibrational forcing, the walker's trajectory becomes chaotic, but its statistical behavior reflects the influence of the unstable orbital solutions. The study results in a complete regime diagram that summarizes the dependence of the walker's behavior on the system parameters. Our predictions compare favorably to the experimental observations of Harris and Bush [“Droplets walking in a rotating frame: from quantized orbits to multimodal statistics,” J. Fluid Mech. 739, 444–464 (2014)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗