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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗
DOE OSTI2021

Modeling soot growth in Comp-B using a detonation history

A chemical kinetics procedure for soot formation calculation is summarized. By importing Lagrangian hydrodynamic calculations, species history is consistently accounted for. Detonation product composition is assumed in region with the burnt high explosive (HE). Soot formation is calculated using the density and energy obtained in the hydrodynamics calculation of HE detonation. Ideal gas equation-of-state (EOS) with temperature-dependent specific heat is used in the HE product region. Sample calculations show reasonably-expected temperature and soot formation.

42 ENGINEERING↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes

Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Towards a quantitative understanding of bonding in supersonic single particle impacts: A three-dimensional FIB-SEM exploration

Particle bonding is crucial to coating quality in cold spray, but it has been a challenge to accurately quantify bonding even in single particle impacts. This paper uses FIB-SEM to explicitly map the particle-substrate interface for Cu-on-Cu single microparticle impacts in a full 3D rendering that spans a wide range of impact velocities. This approach permits a detailed quantification of the total bonding area and all of its associated components. In addition to revealing why prior 2D characterization efforts have missed important details about impact bonding, these data quantitatively reveal the evolution of bonding from its onset at the “critical velocity” V cr (where bonding is generally poor, ∼6 %) to its peak at around 1.3‧V cr (where almost 90 % of the particle bonds). Further increase in the velocity to 1.5‧V cr and beyond finds the onset of hydrodynamic penetration and a decrease in bonding. These data then support the development of a simple analytical model based on oxide rarefication and extrusion of bare metal through gaps in the oxide layer as driving the development of bonding. As a result, the model reproduces the experiments and provides guidance on optimization of bonding as a function of material and process parameters.

Analytical methods↗
DOE OSTI2024

SPH Modeling of Biomass Granular Flow: Engineering Application in Hoppers and Augers

Numerical modeling of granular biomass material flow in handling operations is indispensable to decipher flow upsets, commonly manifested as clogging and jamming in hoppers and augers. With a computational tool developed based on smoothed particle hydrodynamics (SPH), we simulated the hopper flow and auger feeding of six granular biomass materials. The good agreement between the experimental and numerical flow rates demonstrated the capability of the developed SPH solver in modeling the complicated flow of granular biomass materials. Further, the impact of physical and numerical parameters is investigated, and the major results show that the hopper flow pattern is controlled by shear band evolution; hopper clogging is collectively influenced by the opening size, wall friction, and material packing, and the impact of material compressibility on the auger feed rate is minimal. These parametric studies validate the solver’s robustness in simulating biomass flow in handling equipment and demonstrate that the SPH computational tool can provide insights about granular flow mechanics to facilitate handling equipment design and optimize handling operations.

09 BIOMASS FUELS↗
DOE OSTI2024

The impact of low-mode symmetry on inertial fusion energy output in the burning plasma state

Indirect Drive Inertial Confinement Fusion Experiments on the National Ignition Facility (NIF) have achieved a burning plasma state with neutron yields exceeding 170 kJ, roughly 3 times the prior record and a necessary stage for igniting plasmas. The results are achieved despite multiple sources of degradations that lead to high variability in performance. Results shown here, for the first time, include an empirical correction factor for mode-2 asymmetry in the burning plasma regime in addition to previously determined corrections for radiative mix and mode-1. Analysis shows that including these three corrections alone accounts for the measured fusion performance variability in the two highest performing experimental campaigns on the NIF to within error. Here we quantify the performance sensitivity to mode-2 symmetry in the burning plasma regime and apply the results, in the form of an empirical correction to a 1D performance model. Furthermore, we find the sensitivity to mode-2 determined through a series of integrated 2D radiation hydrodynamic simulations to be consistent with the experimentally determined sensitivity only when including alpha-heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Plasma decay of nanosecond pulsed laser-produced Ar and Ar–H 2 O sparks at atmospheric pressure

Time-resolved diagnostics were applied to investigate free-electron properties in nanosecond laser-produced discharges generated in atmospheric pressure Ar and in Ar–3%H 2 O. The discharges were generated using 23 ns, 1064 nm laser pulses. Broadband plasma imaging and laser Thomson scattering were combined with optical emission spectroscopy, with particular emphasis on the Stark broadening of the H α and H β lines. The plasma exhibited a bright emission that persists for up to 30–40 µs after breakdown. Plasma emission was then followed by a very weak glow emission that persisted for up to 19 ms after breakdown. Peak electron number density of ∼2 × 10 17 cm −3 and electron temperature of ∼7 eV were measured. An excellent agreement between both techniques was obtained regarding absolute electron number densities. The inferred free-electron temporal decay dynamics are consistent with processes dominated by hydrodynamic expansion and two- and three-body electron–ion recombination. These results provide benchmark data for modeling nanosecond laser discharges and demonstrate the reliability of combining Thomson scattering with Stark broadening in atmospheric laser sparks.

Thomson scattering↗
DOE OSTI2020

Joint galaxy–galaxy lensing and clustering constraints on galaxy formation

ABSTRACT We compare predictions for galaxy–galaxy lensing profiles and clustering from the Henriques et al. public version of the Munich semi-analytical model (SAM) of galaxy formation and the IllustrisTNG suite, primarily TNG300, with observations from KiDS + GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, and stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast \gt 10^{11}\, \mathrm{ M}_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode active galactic nucleus accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high-mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only (χ2 = 1.81 compared to χ2 = 7.79 for the SAM) and locally brightest galaxy samples (χ2 = 3.80 compared to χ2 = 5.01). With added dust corrections to the colours it matches the SDSS clustering signal well for red low-mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\, {{\ \rm per\ cent}}$ excessive lensing signals for intermediate-mass red galaxies with 10.2 < log10M*[M⊙] < 11.2 at $r \approx 0.6\, h^{-1}\, \text{Mpc}$, which require further theoretical development.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf galaxy evolution

ABSTRACT We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR), respectively, to study 78 radio active galactic nuclei (AGN) in nearby (z < 0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages, and accretion rates are ∼1035 W, ∼5 Myr, and ∼10−3.4 M⊙ yr−1, respectively. The median mechanical energy output of these sources is ∼100 times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability, and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf galaxy evolution, as is the case in massive galaxies.

Davis, F. (ORCID:0000000285224980)↗
DOE OSTI2024

Self-assembled reconfigurable pump architectures via magnetic colloidal swarms

Self-assembled swarms of interactive active units, which are adaptive and dynamically reconfigurable to accommodate different functionalities, represent a promising platform for the development of next-generation robotics. Here, we utilize the emergent collective behavior of active magnetic colloids confined in quasi-two-dimensional arrays of overlapping wells to demonstrate the self-organization of a colloidal swarm into a dynamic pump architecture capable of controlled transport of passive cargo particles. This dynamic architecture provides a global unidirectional looping flow pattern along the entire length of the system. We show that the flow direction of the dynamic swarm-based pump can be externally controlled by a phase shift of a driving magnetic field energizing the swarm. The experimental observations are supported by computational modeling based on phenomenological coarse-grained particle dynamics coupled to shallow-water Navier-Stokes hydrodynamics. In conclusion, our findings demonstrate how the emergent collective behavior of a swarm can be orchestrated into a desired functionality by exploiting the interplay between activity and confinement potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Experimental measurement of the effective contact angle for solvent/packing interactions in a structured packed column for CO 2 capture

The contact angle is a critical factor for determining the effective mass transfer area for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Many widely used commercial packings involve microscale features (perforation, corrugation, etc.) that may also affect the wetting behavior. This study proposes a systematic method of using a modified Wilhelmy plate to measure the effective contact angle to characterize the solvent and featured packing interaction. In lieu of computational efforts relying on assumptions guided by semi-informed correlation, the proposed method directly measures the effective contact angle as a function of the solvent and packing thermophysical and hydrodynamic properties. The characterization of the effective contact angle is then integrated in the computational fluid dynamics modeling to reduce uncertainty in the prediction of effective mas transfer area. Experiments were conducted for stainless steel coupons using water and aqueous sodium hydroxide (NaOH) solvent for verification of the proposed experimental protocol. The surface tension of the aqueous NaOH solvent was altered using surfactant and antifoam. The effective contact angle increases with the increased value of the surface tension for flat stainless steel sheets. On the other hand, effective contact angles do not vary in Mellapak coupons for aqueous monoethanolamine (MEA) and NaOH solvent. In this case, surface textures and sheet design play a dominant role in the surface tension of solvents. Furthermore, CO 2 loading has a significant effect on the contact angle for Mellapak coupons. As expected, the contact angle decreases with the increasing temperature of CO 2 capture solvents (MEA, EEMPA). The effective contact angle measurement using the Wilhelmy plate method can provide more accurate and efficient solvents for characterizing the solvent-packing surface interactions. Subsequently, the present method enhances the accuracy of the prediction of the effective mass transfer area in carbon capture by solvent absorption.

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

Numerical analysis of a two-dimensional nonsteady detonations

In the present work a system of two-dimensional nonsteady hydrodynamic and chemical kinetic equations was numerically integrated for an exothermic system. Assumed two-step reaction model simulates practically an oxyhydrogen mixture. The calculation starts from a plane Chapman-Jouguet detonation as an initial condition. Two-dimensional disturbances are generated by artificially placing nonuniformities ahead of the detonation front. Regardless of the difference of the given initial disturbances, a fixed number of triple shock waves were produced for a fixed combination of mixture model and geometry when the transition period was over. This shows that for a given detonation tube geometry any exothermic system has its own characteristic multidimensional structure. The obtained number of triple shock waves contained in the detonation front was in agreement with existing experimental observations under the same condition.

Taki, S.↗
NASA NTRS1997

Gamma-rays from supernova remnants and the signatures of diffusive shock acceleration

A nonlinear shock acceleration model which generates non-thermal proton distributions and includes a self-consistent determination of shock hydrodynamics, is considered. Gamma ray spectra are obtained for supernova remnants, allowing for the cessation of acceleration to high energies due to the finite ages and the sizes of the remnants. Gamma ray spectral cutoffs can be observed in the TeV range for reasonable remnant parameters and deviations from power law behavior are found at all energies from 1 MeV to cutoff. Correlated observations by the International Gamma Ray Astrophysics Laboratory, Whipple and other instruments may provide stringent constraints to understanding supernova remnants.

Baring, Matthew G.↗
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↗