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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Reconnection-driven Magnetohydrodynamic Turbulence in a Simulated Coronal-hole Jet

Extreme-ultraviolet and X-ray jets occur frequently in magnetically open coronal holes on the Sun, especially at high solar latitudes. Some of these jets are observed by white-light coronagraphs as they propagate through the outer corona toward the inner heliosphere, and it has been proposed that they give rise to microstreams and torsional Alfvén waves detected in situ in the solar wind. To predict and understand the signatures of coronal-hole jets, we have performed a detailed statistical analysis of such a jet simulated by an adaptively refined magnetohydrodynamics model. The results confirm the generation and persistence of three-dimensional, reconnection-driven magnetic turbulence in the simulation. We calculate the spatial correlations of magnetic fluctuations within the jet and find that they agree best with the Müller–Biskamp scaling model including intermittent current sheets of various sizes coupled via hydrodynamic turbulent cascade. The anisotropy of the magnetic fluctuations and the spatial orientation of the current sheets are consistent with an ensemble of nonlinear Alfvén waves. These properties also reflect the overall collimated jet structure imposed by the geometry of the reconnecting magnetic field. A comparison with Ulysses observations shows that turbulence in the jet wake is in quantitative agreement with that in the fast solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Electrostatic soliton and double layer structures in unmagnetized degenerate pair plasmas

The acoustic solitons and double layers are studied in unmagnetized quantum electron-positron plasmas in the presence of stationary ions. The quantum hydrodynamic model is employed and reductive perturbation method is used to derive the Korteweg-de Vries (KdV) and extended KdV equations for solitons and double layers, respectively. It is found that in the linear limit both slow acoustic and fast Langmuir waves can propagate in such type of quantum plasmas like in classical pair-ion or pair plasmas. The amplitude and width of the electrostatic solitons are found to be decreasing with the increase in concentration of positrons (or decrease in the concentration of ions) in degenerate electron-positron-ion plasmas. It is found that only rarefactive double layer can exist in such plasmas which depend on various parameters. The dependence of double layer structure on ion concentration and quantum diffraction effects of electrons and positrons are also discussed. The results are also elaborated graphically by considering dense plasma parameters in the outer layers of astrophysical objects such as white dwarfs and neutron stars.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Generalized hydrodynamic analysis of transport through a finite open nanopore for two-component single-file systems

Single-file diffusion (SFD) in finite open nanopores is characterized by nonzero spatially varying tracer diffusion coefficients within a generalized hydrodynamic description. This contrasts with infinite SFD systems where tracer diffusivity vanishes. In standard tracer counterpermeation (TCP) analysis, two reservoirs, each containing a different species, are connected to opposite ends of a finite pore. Here, we implement an extended TCP analysis to allow the two reservoirs to contain slightly different mixtures of the two species. Then, determination of diffusion fluxes through the pore allows extraction of diffusion coefficients for near-constant partial concentrations of the two species. This analysis is applied for a lattice-gas model describing two-component SFD through a finite linear pore represented by a one-dimensional array of cells. Two types of particles, A and B, can hop only to adjacent empty cells with generally different rates, $h_A$ and $h_B$. Particles are noninteracting other than exclusion of multiple cell occupancy. Results reveal generalized hydrodynamic tracer diffusion coefficients which adopt small values inversely proportional to pore length in the pore center, but which are strongly enhanced near pore openings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

On the lower mass limit for the carbon detonation scenario

The lower mass limit for carbon detonation (about 4 solar masses) is due to an argument by Paczynski (1970) that stars in the range between 1.4 and 3.5 solar masses will lose their envelopes prior to carbon detonation. This argument is based on two criteria for mass loss: (1) the envelope is unstable to adiabatic oscillations, and (2) the total energy of the envelope is positive. Recent hydrodynamic calculations of red-giant envelopes have shown that the envelopes are highly nonadiabatic and that the second criterion is not a sufficient condition for mass loss. Thus, a calculation which takes into account the nonadiabatic nature of the envelope is required in order to determine the lower mass limit for carbon detonation. This lower mass limit is directly related to the formation mechanisms and statistics of supernovae and pulsars.

Endal, A. S.↗
DOE OSTI2015

A deep search for faint galaxies associated with very low-redshift C iv absorbers. II. Program design, absorption-line measurements, and absorber statistics

To investigate the evolution of metal-enriched gas over recent cosmic epochs as well as to characterize the diffuse, ionized, metal-enriched circumgalactic medium, we have conducted a blind survey for C iv absorption systems in 89 QSO sightlines observed with the Hubble Space Telescope Cosmic Origins Spectrograph. We have identified 42 absorbers at z < 0.16, comprising the largest uniform blind sample size to date in this redshift range. Our measurements indicate an increasing C iv absorber number density per comoving path length (dN/dX= 7.5 ± 1.1) and modestly increasing mass density relative to the critical density of the universe (Ω{sub C iv} = 10.0 ± 1.5 × 10{sup −8}) from z ∼ 1.5 to the present epoch, consistent with predictions from cosmological hydrodynamical simulations. Furthermore, the data support a functional form for the column density distribution function that deviates from a single power law, also consistent with independent theoretical predictions. As the data also probe heavy element ions in addition to C iv at the same redshifts, we identify, measure, and search for correlations between column densities of these species where components appear to be aligned in velocity. Among these ion–ion correlations, we find evidence for tight correlations between C ii and Si ii, C ii and Si iii, and C iv and Si iv, suggesting that these pairs of species arise in similar ionization conditions. However, the evidence for correlations decreases as the difference in ionization potential increases. Finally, when controlling for observational bias, we find only marginal evidence for a correlation (86.8% likelihood) between the Doppler line width b(C iv) and column density N(C iv).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Transverse momentum distributions of strange hadrons produced in p–p collisions at √s{sub NN} = 200 GeV

The mid-rapidity transverse momentum spectra of strange hadrons (K{sup +},K{sup −},K{sup −}/K{sup +},Λ,Λ̄,Ξ{sup −},Ξ{sup +}, and Ω) produced in p–p collisions at the highest RHIC energy √s{sub NN} = 200 GeV have been studied using a statistical unified thermal freeze-out model. The calculated results are found to be in good agreement with the experimental data taken from STAR and BRAHMS experiments. The fits of the transverse momentum spectra to the model calculations provide the thermal freeze-out conditions in terms of the temperature and collective flow effect parameters for different particle species. The model incorporates a longitudinal and a transverse hydrodynamic flow. The rapidity distributions of kaons and their ratios are also reproduced successfully, which reveals the presence of partial nuclear transparency effects in p–p collisions at √s{sub NN} = 200 GeV. The contributions from heavier decay resonances are also taken into account.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The MFiX Particle-in-Cell Method (MFiX-PIC) Theory Guide

MFiX (Multiphase Flow with Interphase eXchanges) is an open-source multiphase flow solver developed at the National Energy Technology Laboratory. Within the code, users have access to a single phase or interpenetrating continua-based multiphase two-fluid model (TFM), a discrete element model (DEM), and a particle-in-cell model (PIC). TFM, DEM, and PIC can all be used to create multiphase simulations that include hydrodynamics, chemical reactions, and heat transfer.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2013

THE PROGENITOR DEPENDENCE OF THE PRE-EXPLOSION NEUTRINO EMISSION IN CORE-COLLAPSE SUPERNOVAE

We perform spherically symmetric general-relativistic simulations of core collapse and the postbounce pre-explosion phase in 32 presupernova stellar models of solar metallicity with zero-age main-sequence masses of 12-120 M {sub Sun }. Using energy-dependent three-species neutrino transport in the two-moment approximation with an analytic closure, we show that the emitted neutrino luminosities and spectra follow very systematic trends that are correlated with the compactness ({approx}M/R) of the progenitor star's inner regions via the accretion rate in the pre-explosion phase. We find that these qualitative trends depend only weakly on the nuclear equation of state (EOS), but quantitative observational statements will require independent constraints on the EOS and the rotation rate of the core as well as a more complete understanding of neutrino oscillations. We investigate the simulated response of water Cherenkov detectors to the electron antineutrino fluxes from our models and find that the large statistics of a galactic core collapse event may allow robust conclusions on the inner structure of the progenitor star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

The Poisson equation at second order in relativistic cosmology

We calculate the relativistic constraint equation which relates the curvature perturbation to the matter density contrast at second order in cosmological perturbation theory. This relativistic ''second order Poisson equation'' is presented in a gauge where the hydrodynamical inhomogeneities coincide with their Newtonian counterparts exactly for a perfect fluid with constant equation of state. We use this constraint to introduce primordial non-Gaussianity in the density contrast in the framework of General Relativity. We then derive expressions that can be used as the initial conditions of N-body codes for structure formation which probe the observable signature of primordial non-Gaussianity in the statistics of the evolved matter density field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Classical-Nova Contribution to the Milky Way’s Al 26 Abundance: Exit Channel of the Key Al 25 ( p , γ ) Si 26 Resonance

We present that classical novae are expected to contribute to the 1809-keV Galactic γ-ray emission by producing its precursor 26 Al, but the yield depends on the thermonuclear rate of the unmeasured 25 Al(p,γ) 26 Si reaction. Using the β decay of 26 P to populate the key J π =3 + resonance in this reaction, we report the first evidence for the observation of its exit channel via a 1741.6±0.6(stat)±0.3(syst) keV primary γ ray, where the uncertainties are statistical and systematic, respectively. By combining the measured γ-ray energy and intensity with other experimental data on 26 Si, we find the center-of-mass energy and strength of the resonance to be E r =414.9±0.6(stat)±0.3(syst)±0.6(lit.) keV and ωγ=23±6(stat)$+11\atop{-10}$(lit.) meV, respectively, where the last uncertainties are from adopted literature data. Finally, we use hydrodynamic nova simulations to model 26 Al production showing that these measurements effectively eliminate the dominant experimental nuclear-physics uncertainty and we estimate that novae may contribute up to 30% of the Galactic 26 Al.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2011

Modeling of a Flooding Induced Station Blackout for a Pressurized Water Reactor Using the RISMC Toolkit

In the Risk Informed Safety Margin Characterization (RISMC) approach we want to understand not just the frequency of an event like core damage, but how close we are (or are not) to key safety-related events and how might we increase our safety margins. The RISMC Pathway uses the probabilistic margin approach to quantify impacts to reliability and safety by coupling both probabilistic (via stochastic simulation) and mechanistic (via physics models) approaches. This coupling takes place through the interchange of physical parameters and operational or accident scenarios. In this paper we apply the RISMC approach to evaluate the impact of a power uprate on a pressurized water reactor (PWR) for a tsunami-induced flooding test case. This analysis is performed using the RISMC toolkit: RELAP-7 and RAVEN codes. RELAP-7 is the new generation of system analysis codes that is responsible for simulating the thermal-hydraulic dynamics of PWR and boiling water reactor systems. RAVEN has two capabilities: to act as a controller of the RELAP-7 simulation (e.g., system activation) and to perform statistical analyses (e.g., run multiple RELAP-7 simulations where sequencing/timing of events have been changed according to a set of stochastic distributions). By using the RISMC toolkit, we can evaluate how power uprate affects the system recovery measures needed to avoid core damage after the PWR lost all available AC power by a tsunami induced flooding. The simulation of the actual flooding is performed by using a smooth particle hydrodynamics code: NEUTRINO.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

RELAXATION PROCESSES IN SOLAR WIND TURBULENCE

Based on global conservation principles, magnetohydrodynamic (MHD) relaxation theory predicts the existence of several equilibria, such as the Taylor state or global dynamic alignment. These states are generally viewed as very long-time and large-scale equilibria, which emerge only after the termination of the turbulent cascade. As suggested by hydrodynamics and by recent MHD numerical simulations, relaxation processes can occur during the turbulent cascade that will manifest themselves as local patches of equilibrium-like configurations. Using multi-spacecraft analysis techniques in conjunction with Cluster data, we compute the current density and flow vorticity and for the first time demonstrate that these localized relaxation events are observed in the solar wind. Such events have important consequences for the statistics of plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Time-domain theory of transient heat conduction in the local limit

Ultrafast and nanoscale heat conduction demands a unified theoretical framework that rigorously bridges macroscopic transport equations with microscopic material properties derived from statistical physics. Existing empirical generalizations of Fourier's law often lack a solid microscopic foundation, failing to connect observed non-Fourier behavior with underlying atomic-scale mechanisms. In this work, we present a time-domain theory of transient heat conduction rooted in Zwanzig's statistical theory of irreversible processes. Central to this framework is the time-domain transport function $\overleftrightarrow{𝑍}$⁡(𝑡) defined through equilibrium time-correlation functions of heat fluxes. This function generalizes the conventional concept of steady-state thermal conductivity, governing the transition of conduction dynamics from onset second sound type wave propagation at finite speeds to diffusion-dominated behavior across broad temporal and spatial scales. Unlike phonon hydrodynamic models that rely on mesoscopic constructs such as phonon drift velocity, our approach provides a quantitative and microscopic description of intrinsic memory effects in transient heat fluxes and applies universally to bulk materials at any temperature or length scale. By integrating atomistic-scale first-principles calculations with continuum-level macroscopic equations, this framework offers a robust foundation for numerical simulations of transient temperature fields. Furthermore, it facilitates the interpretation and design of transient thermal grating experiments using nanometer-scale heat sources and ultrafast laser systems in the extreme ultraviolet and x-ray wavelength ranges, advancing our understanding of heat dissipation dynamics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Impact of late-time neutrino emission on the diffuse supernova neutrino background

In the absence of high-statistics supernova neutrino measurements, estimates of the diffuse supernova neutrino background (DSNB) hinge on the precision of simulations of core-collapse supernovae. Understanding the cooling phase of protoneutron star (PNS) evolution (≳ 1 s after core bounce) is crucial, since approximately 50% of the energy liberated by neutrinos is emitted during the cooling phase. We model the cooling phase with a hybrid method by combining the neutrino emission predicted by 3D hydrodynamic simulations with several cooling-phase estimates, including a novel two-parameter correlation depending on the final baryonic PNS mass and the time of shock revival. We find that the predicted DSNB event rate at Super-Kamiokande can vary by a factor of ~2–3 depending on the cooling-phase treatment. Here we also find that except for one cooling estimate, the range in predicted DSNB events is largely driven by the uncertainty in the neutrino mean energy. With a good understanding of the late-time neutrino emission, more precise DSNB estimates can be made for the next generation of DSNB searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Machine learning assisted bayesian inference of mix and hot-spot conditions in NIF implosions

Experiments on the National Ignition Facility (NIF) have provided clear evidence of ablator material mixing into the Hot-Spot, leading to degraded performance. However, inferring the amount of mix and Hot-Spot conditions from typical experimental observations (e.g. x-ray spectra and images) is highly challenging. Here, we have developed an analysis method that utilizes machine learning assisted Bayesian inference to find the probability distributions of the Hot-Spot and mix conditions. This approach uses a neural network, trained on an idealized 2-dimensional representation of the Hot-Spot and mix distribution, and Bayesian inference to find the statistical distributions of Hot-Spot conditions that provide a match with observations. We have tested this method with synthetic data from simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

SEVEN-YEAR WILKINSON MICROWAVE ANISOTROPY PROBE (WMAP ) OBSERVATIONS: COSMOLOGICAL INTERPRETATION

The combination of seven-year data from WMAP and improved astrophysical data rigorously tests the standard cosmological model and places new constraints on its basic parameters and extensions. By combining the WMAP data with the latest distance measurements from the baryon acoustic oscillations (BAO) in the distribution of galaxies and the Hubble constant (H{sub 0}) measurement, we determine the parameters of the simplest six-parameter {Lambda}CDM model. The power-law index of the primordial power spectrum is n{sub s} = 0.968 {+-} 0.012 (68% CL) for this data combination, a measurement that excludes the Harrison-Zel'dovich-Peebles spectrum by 99.5% CL. The other parameters, including those beyond the minimal set, are also consistent with, and improved from, the five-year results. We find no convincing deviations from the minimal model. The seven-year temperature power spectrum gives a better determination of the third acoustic peak, which results in a better determination of the redshift of the matter-radiation equality epoch. Notable examples of improved parameters are the total mass of neutrinos, {Sigma}M{sub {nu}} < 0.58 eV(95%CL), and the effective number of neutrino species, N{sub eff} = 4.34{sup +0.86}{sub -0.88} (68% CL), which benefit from better determinations of the third peak and H{sub 0}. The limit on a constant dark energy equation of state parameter from WMAP+BAO+H{sub 0}, without high-redshift Type Ia supernovae, is w = -1.10 {+-} 0.14 (68% CL). We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis by measuring Y{sub p} = 0.326 {+-} 0.075 (68% CL). We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z = 1090 and the dominance of adiabatic scalar fluctuations. The seven-year polarization data have significantly improved: we now detect the temperature-E-mode polarization cross power spectrum at 21{sigma}, compared with 13{sigma} from the five-year data. With the seven-year temperature-B-mode cross power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved by 38% to {Delta}{alpha} = -1.{sup o}1 {+-} 1.{sup o}4(statistical) {+-} 1.{sup o}5(systematic) (68% CL). We report significant detections of the Sunyaev-Zel'dovich (SZ) effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data on a cluster-by-cluster basis. However, it is a factor of 0.5-0.7 times the predictions from 'universal profile' of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically expected SZ power spectrum recently measured by the South Pole Telescope Collaboration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Torqued fireballs in relativistic heavy-ion collisions

We show that the fluctuations in the wounded-nucleon model of the initial stage of relativistic heavy-ion collisions, together with the natural assumption that the forward- (backward-) moving wounded nucleons emit particles preferably in the forward (backward) direction, lead to an event-by-event torqued fireball. The principal axes associated with the transverse shape are rotated in the forward region in the opposite direction than in the backward region. On the average, the standard deviation of the relative torque angle between the forward and backward rapidity regions is {approx}20 deg. for the central and 10 deg. for the midperipheral collisions. The hydrodynamic expansion of a torqued fireball leads to a torqued collective flow, yielding, in turn, torqued principal axes of the transverse-momentum distributions at different rapidities. We propose experimental measures, based on cumulants involving particles in different rapidity regions, which should allow for a quantitative determination of the effect from the data. To estimate the nonflow contributions from resonance decays we run Monte Carlo simulations with therminator, a thermal heavy-ion generator. If the event-by-event torque effect is found in the data, it will support the assumptions concerning the fluctuations in the early stage of the fireball formation, as well as the hypothesis of the asymmetric rapidity shape of the emission functions of the moving sources in the nucleus-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗