Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Wind and Rotation Enhanced Escape from the Early Terrestrial Atmospheres

The earliest atmospheres of the terrestrial planets are thought to have been hotter, have stronger winds and rotate faster than atmospheres of today. Since these primitive atmospheres were weakly bound, they evolved rapidly because atmospheric escape was very strong, often referred to as "blowoff." Such escape has been treated as hydrodynamic, transonic flow; similar to solar wind flow dynamics. However, in many cases, although the outward flow is hydrodynamic at low altitudes, it becomes collisionless at higher altitudes, before sonic speeds are ever attained. Recent models dealing with the transition from fluid to kinetic flow have applied the Jeans escape flux at the exobase. This approach leads to escape rates that are too low, because thermospheric winds and planetary rotation are known to increase the escape flux above the corresponding Jeans flux. In particular, for a given density and temperature at the exobase, the escape flux increases as the wind speed and/or the rotation rate increase. Also, for a given wind speed and rotation rate, the escape flux enhancement over the Jeans flux increases as the mass of an escaping constituent increases, an important factor in isotope fractionation, especially the enrichment of deuterium on Mars. Accounting for a range of possible temperatures, thermospheric wind speeds and planetary rotation rates in the primitive atmospheres of the terrestrial planets, estimates are made of light constituent escape flux increases over the corresponding Jeans fluxes.

Hartle, Richard E.↗
DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Stochastic reacceleration of cosmic rays in the interstellar medium

The effects of reacceleration on cosmic rays have been studied over a wide charge and energy range using a model of reacceleration by the interstellar turbulence. We take into account only inevitable stochastic reacceleration of cosmic rays by the random hydrodynamic waves, which supposedly exist in the interstellar medium and provide a means for cosmic-ray scattering and spatial diffusion in the Galaxy. Our calculations reproduce not only the B/C ratio but also the H and He data over the entire energy range where the measurements are available. However, the sub-Fe to Fe ratio is not fitted as well as the B/C ratio, and the reacceleration effect does not seem to remove the need for truncation of short path lengths, which is apparently required by the standard leaky box model. This work demonstrates that the cosmic-ray data can be represented at least as well by a reacceleration model with a simple rigidity power-law escape length, which agrees with the Kolmogorov-type spectrum of hydromagnetic turbulence, as they can by the standard leaky box model with its ad hoc escape lengths.

Seo, E. S.↗
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 NTRS2004

Underexpanded Screeching Jets From Circular, Rectangular, and Elliptic Nozzles

The screech frequency and amplitude, the shock spacing, the hydrodynamic-acoustic standing wave spacing, and the convective velocity of large organized structures are measured in the nominal Mach number range of 1.1 less than or = Mj less that or = l0.9 for supersonic, underexpanded jets exhausting from a circular, a rectangular and an elliptic nozzle. This provides a carefully measured data set useful in comparing the importance of various physical parameters in the screech generation process. The hydrodynamic-acoustic standing wave is formed between the potential pressure field of large turbulent structures and the acoustic pressure field of the screech sound. It has been demonstrated earlier that in the currently available screech frequency prediction models replacement of the shock spacing by the standing wave spacing provides an exact expression. In view of this newly found evidence, a comparison is made between the average standing wavelength and the average shock spacing. It is found that there exists a small, yet important, difference, which is dependent on the azimuthal screech mode. For example, in the flapping modes of circular, rectangular, and elliptic jets, the standing wavelength is slightly longer than the shock spacing, while for the helical screech mode in a circular jet the opposite is true. This difference accounts for the departure of the existing models from predicting the exact screech frequency. Another important parameter, necessary in screech prediction, is the convective velocity of the large organized structures. It is demonstrated that the presence of the hydrodynamic-acoustic standing wave, even inside the jet shear layer, becomes a significant source of error in the convective velocity data obtained using the conventional methods. However, a new relationship, using the standing wavelength and screech frequency is shown to provide more accurate results.

Panda, J.↗
NASA NTRS1995

Cooling Panel Optimization for the Active Cooling System of a Hypersonic Aircraft

Optimization of cooling panels for an active cooling system of a hypersonic aircraft is explored. The flow passages are of rectangular cross section with one wall heated. An analytical fin-type model for incompressible flow in smooth-wall rectangular ducts with coupled wall conduction is proposed. Based on this model, the a flow rate of coolant to each design minimum mass flow rate or coolant for a single cooling panel is obtained by satisfying hydrodynamic, thermal, and Mach number constraints. Also, the sensitivity of the optimal mass flow rate of coolant to each design variable is investigated. In addition, numerical solutions for constant property flow in rectangular ducts, with one side rib-roughened and coupled wall conduction, are obtained using a k-epsilon and wall function turbulence model, these results are compared with predictions of the analytical model.

Youn, B.↗
DOE OSTI2021

The relaxation behavior of supercooled and glassy imidacloprid

Employing dielectric spectroscopy, oscillatory shear rheology, and calorimetry, the present work explores the molecular dynamics of the widely used insecticide imidacloprid above and below its glass transition temperature. In its supercooled liquid regime, the applied techniques yield good agreement regarding the characteristic structural (alpha) relaxation times of this material. In addition, the generalized Gemant–DiMarzio–Bishop model provides a good conversion between the frequency-dependent dielectric and shear mechanical responses in its viscous state, allowing for an assessment of imidacloprid’s molecular hydrodynamic radius. In order to characterize the molecular dynamics in its glassy regime, we employ several approaches. These include the application of frequency–temperature superposition (FTS) to its isostructural dielectric and rheological responses as well as use of dielectric and calorimetric physical aging and the Adam–Gibbs–Vogel model. While the latter approach and dielectric FTS provide relaxation times that are close to each other, the other methods predict notably longer times that are closer to those reflecting a complete recovery of ergodicity. In conclusion, this seemingly conflicting dissimilarity demonstrates that the molecular dynamics of glassy imidacloprid strongly depends on its thermal history, with high relevance for the use of this insecticide as an active ingredient in technological applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

Parameterizations of Chromospheric Condensations in dG and dMe Model Flare Atmospheres

The origin of the near-ultraviolet and optical continuum radiation in flares is critical for understanding particle acceleration and impulsive heating in stellar atmospheres. Radiative-hydrodynamic (RHD) simulations in 1D have shown that high energy deposition rates from electron beams produce two flaring layers at T approximately 10 (exp 4) K that develop in the chromosphere: a cooling condensation (downflowing compression) and heated non-moving (stationary) flare layers just below the condensation. These atmospheres reproduce several observed phenomena in flare spectra, such as the red-wing asymmetry of the emission lines in solar flares and a small Balmer jump ratio in M dwarf flares. The high beam flux simulations are computationally expensive in 1D, and the (human) timescales for completing NLTE models with adaptive grids in 3D will likely be unwieldy for some time to come. We have developed a prescription for predicting the approximate evolved states, continuum optical depth, and emergent continuum flux spectra of RHD model flare atmospheres. These approximate prescriptions are based on an important atmospheric parameter: the column mass (m(sub ref)) at which hydrogen becomes nearly completely ionized at the depths that are approximately in steady state with the electron beam heating. Using this new modeling approach, we find that high energy flux density (>F11) electron beams are needed to reproduce the brightest observed continuum intensity in IRIS data of the 2014 March 29 X1 solar flare, and that variation in m(sub ref) from 0.001 to 0.02 g cm (exp -2) reproduces most of the observed range of the optical continuum flux ratios at the peak of M dwarf flares.

Kowalski, Adam F.↗
NASA NTRS2001

Hydrodynamic Stability of Liquid-Propellant Combustion: Landau's Problem Revisited

Hydrodynamic, or Landau, instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. As its name suggests, it stems from hydrodynamic effects connected with thermal expansion across the reaction region. In the context of liquid-propellant combustion, the classical models that originally predicted this phenomenon have been extended to include the important effects that arise from a dynamic dependence of the burning rate on the local pressure and temperature fields. Thus, the onset of Landau instability has now been shown to occur for sufficiently small negative values of the pressure sensitivity of the burning rate, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. It has also been shown that the onset of instability occurs for decreasing values of the disturbance wave number as the gravitational-acceleration parameter decreases. Consequently, in an appropriate weak-gravity limit, Landau instability becomes a long-wave phenomena associated with the formation of large cells on the liquid-propellant surface. Additionally, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. This instability occurs for sufficiently large negative values of the pressure sensitivity, and is enhanced by increasing values of the burning-rate temperature sensitivity. It is further shown that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wave number, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗
DOE OSTI2022

Hot Circumsingle Disks Drive Binary Black Hole Mergers in Active Galactic Nucleus Disks

Abstract Binary black hole (BBH) mergers, particularly those with component masses in the pair-instability gap, may be produced by hierarchical mergers in the disks surrounding Active Galactic Nuclei (AGNs). While the interaction of an embedded BBH with an AGN disk is typically assumed to facilitate a merger, recent high-resolution hydrodynamical simulations challenge this assumption. However, these simulations often have simplified treatments for gas thermodynamics. In this work, we model the possible consequence of various feedback from an embedded BBH with a simple model that maintains an enhanced temperature profile around each binary component. We show that when the minidisks around each BH become hotter than the background by a factor of three, the BBH orbital evolution switches from expansion to contraction. By analyzing the gravitational torque profile, we find that this change in direction is driven by a weakening of the minidisk spirals and their positive torque on the binary. Our results highlight the important role of thermodynamics around BBHs and its effect on their orbital evolution, suggesting that AGN disks could be efficient factories for BBH mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗
NASA NTRS1982

Aspects of Coulomb damping in rotors supported on hydrodynamic bearings

The paper is concerned with the effect of friction in drive couplings on the non-sychronous whirling of a shaft. A simplified model is used to demonstrate the effect of large coupling misalignments on the stability of the system. It is concluded that provided these misalignments are large enough, the system becomes totally stable provided the shaft is supported on bearings exhibiting a viscous damping capacity.

Morton, P. G.↗
NASA NTRS1984

Theory and modeling of atmospheric turbulence, part 1

The cascade transfer which is the only function to describe the mode coupling as the result of the nonlinear hydrodynamic state of turbulence is discussed. A kinetic theory combined with a scaling procedure was developed. The transfer function governs the non-linear mode coupling in strong turbulence. The master equation is consistent with the hydrodynamical system that describes the microdynamic state of turbulence and has the advantages to be homogeneous and have fewer nonlinear terms. The modes are scaled into groups to decipher the governing transport processes and statistical characteristics. An equation of vorticity transport describes the microdynamic state of two dimensional, isotropic and homogeneous, geostrophic turbulence. The equation of evolution of the macrovorticity is derived from group scaling in the form of the Fokker-Planck equation with memory. The microdynamic state of turbulence is transformed into the Liouville equation to derive the kinetic equation of the singlet distribution in turbulence. The collision integral contains a memory, which is analyzed with pair collision and the multiple collision. Two other kinetic equations are developed in parallel for the propagator and the transition probability for the interaction among the groups.

Source record↗
DOE OSTI2025

Cosmological Hydrodynamics at Exascale: A Trillion-Particle Leap in Capability

Resolving the most fundamental questions in cosmology requires simulations that match the scale, fidelity, and physical complexity demanded by next-generation sky surveys. To achieve the realism needed for this critical scientific partnership, detailed gas dynamics must be treated self-consistently with gravity for end-to-end modeling of structure formation. Exascale computing enables simulations that span survey-scale volumes while incorporating key astrophysical processes that shape complex cosmic structures. We present results from CRK-HACC, a cosmological hydrodynamics code built for extreme scalability. Using separation-of-scale techniques, GPU-resident tree solvers, in situ analysis pipelines, and multi-tiered I/O, CRK-HACCexecuted Frontier-E: a four trillion particle full-sky simulation, over an order of magnitude larger than previous efforts. The run achieved 513.1 PFLOPs peak performance, processing 46.6 billion particles per second and writing more than 100 PB of data in just over one week of runtime. Frontier-E marks a significant advance in predictive modeling for next-generation cosmological science.

Frontiere, Nicholas [Argonne National Laboratory (↗
DOE OSTI2021

Improved modeling of the solid-to-plasma transition of polystyrene ablator for laser direct-drive inertial confinement fusion hydrocodes

We report that the target performance of laser direct-drive inertial confinement fusion (ICF) can be limited by the development of hydrodynamic instabilities resulting from the nonhomegeneous laser absorption at the target surface, i.e., the laser imprint on the ablator. To understand and describe the formation of these instabilities, the early ablator evolution during the laser irradiation should be considered. In this work, an improved modeling of the solid-to-plasma transition of a polystyrene ablator for laser direct-drive ICF is proposed. This model is devoted to be implemented in hydrocodes dedicated to ICF which generally assume an initial plasma state. The present approach consists of the two-temperature model coupled to the electron, ion and neutral dynamics including the chemical fragmentation of polystyrene. The solid-to-plasma transition is shown to significantly influence the temporal evolution of both free electron density and temperatures, which can lead to different shock formation and propagation compared with an initial plasma state. The influence of the solid-to-plasma transition on the shock dynamics is evidenced by considering the scaling law of the pressure with respect to the laser intensity. The ablator transition is shown to modify the scaling law exponent compared with an initial plasma state.

36 MATERIALS SCIENCE↗
NASA NTRS2023

Optimizing the Resolution of Hydrodynamic Simulations for MCRaT Radiative Transfer Calculations

Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydro dynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affect the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

Jose Arita-Escalante↗
DOE OSTI2023

Collective dynamics of polarized spin-half fermions in relativistic heavy-ion collisions

Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2023

The F-CHROMA Grid of 1D RADYN Flare Models

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

hydrodynamics↗