SEARCH · Engineering Papers

Results for “Diffusion models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗

DOE OSTI2016

Collisional transport across the magnetic field in drift-fluid models

Drift ordered fluid models are widely applied in studies of low-frequency turbulence in the edge and scrape-off layer regions of magnetically confined plasmas. Here, we show how collisional transport across the magnetic field is self-consistently incorporated into drift-fluid models without altering the drift-fluid energy integral. We demonstrate that the inclusion of collisional transport in drift-fluid models gives rise to diffusion of particle density, momentum, and pressures in drift-fluid turbulence models and, thereby, obviates the customary use of artificial diffusion in turbulence simulations. We further derive a computationally efficient, two-dimensional model, which can be time integrated for several turbulence de-correlation times using only limited computational resources. The model describes interchange turbulence in a two-dimensional plane perpendicular to the magnetic field located at the outboard midplane of a tokamak. The model domain has two regions modeling open and closed field lines. The model employs a computational expedient model for collisional transport. Numerical simulations show good agreement between the full and the simplified model for collisional transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data

Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.

MPPD↗

DOE OSTI2021

Efficient Space–Time Reduced Order Model for Linear Dynamical Systems in Python Using Less than 120 Lines of Code

A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space–time ROM for linear dynamical problems has been developed [Choi et al., Space–tume reduced order model for large-scale linear dynamical systems with application to Boltzmann transport problems, Journal of Computational Physics, 2020], which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 10−5 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space–time least-squares Petrov–Galerkin (LSPG) projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space–time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space–time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our LSPG method with the traditional Galerkin method and show that the space–time ROMs can achieve O(10−3) to O(10−4) relative errors for these problems. Depending on parameter–separability, online speed-ups may or may not be achieved. For the FOMs with parameter–separability, the space–time ROMs can achieve O(10) online speed-ups. Finally, we present an error analysis for the space–time LSPG projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

A Lagrangian Simulation of Subsonic Aircraft Exhaust Emissions

To estimate the effect of subsonic and supersonic aircraft exhaust on the stratospheric concentration of NO(y), we employ a trajectory model initialized with air parcels based on the standard release scenarios. The supersonic exhaust simulations are in good agreement with 2D and 3D model results and show a perturbation of about 1-2 ppbv of NO(y) in the stratosphere. The subsonic simulations show that subsonic emissions are almost entirely trapped below the 380 K potential temperature surface. Our subsonic results contradict results from most other models, which show exhaust products penetrating above 380 K, as summarized. The disagreement can likely be attributed to an excessive vertical diffusion in most models of the strong vertical gradient in NO(y) that forms at the boundary between the emission zone and the stratosphere above 380 K. Our results suggest that previous assessments of the impact of subsonic exhaust emission on the stratospheric region above 380 K should be considered to be an upper bound.

Schoeberl, M. R.↗

NASA NTRS1990

Modeling of the SiC chemical vapor deposition process and comparison with experimental results

A model of the NASA Lewis SiC CVD process is described, and the results presented. A key feature of the model is the direct coupling of gas-phase chemical kinetics with diffusion. A deposition chemistry model using reactive sticking coefficients for each gas-phase species provides the diffusion boundary conditions. SiC deposition rates predicted by the model agree reasonably well with experimentally observed rates.

Annen, K. D.↗

DOE OSTI2018

ISD3: a particokinetic model for predicting the combined effects of particle sedimentation, diffusion and dissolution on cellular dosimetry for in vitro systems

The development of particokinetic models describing the delivery of insoluble or poorly soluble nanoparticles to cells in liquid cell culture systems has improved the basis for dose-response analysis, hazard ranking from high-throughput systems, and now allows for translation of exposures across in vitro and in vivo test systems. Complimentary particokinetic models that address processes controlling delivery of both particles and released ions to cells, and the influence of particle size changes from dissolution on particle delivery for cell-culture systems would help advance our understanding of the role of particles ion dosimetry on cellular toxicology. We developed ISD 3 , an extension of our previously published model for insoluble particles, by deriving a specific formulation of the Population Balance Equation for soluble particles. ISD3 describes the time, concentration and particle size dependent dissolution of particles, their delivery to cells, and the delivery and uptake of ions to cells in in vitro liquid test systems. The model is modular, and can be adapted by application of any empirical model of dissolution, alternative approaches to calculating sedimentation rates, and cellular uptake or treatment of boundary conditions. We apply the model to calculate the particle and ion dosimetry of nanosilver and silver ions in vitro after calibration of two empirical models, one for particle dissolution and one for ion uptake. The results demonstrate utility and accuracy of the ISD 3 framework for dosimetry in these systems. Total media ion concentration, particle concentration and total cell-associated silver time-courses were well described by the model, across 2 concentrations of 20 and 110 nm particles. ISD 3 was calibrated to dissolution data for 20 nm particles as a function of serum protein concentration, but successfully described the media and cell dosimetry time-course for both particles at all concentrations and time points. We also report the finding that protein content in media has effects both on the initial rate of dissolution and the resulting near-steady state ion concentration in solution.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1989

Observational constraints on an embedded cloud model for the soft X-ray diffuse background

The reason for the observed anticorrelation between soft X-ray diffuse background (SXRB) and neutral hydrogen column densities are considered. Using the B and C band data from an all-sky survey of the SXRB, as well as data in the Be band, it is shown that the Jakobsen and Kahn (1986) embedded cloud model, in which several statistical properties of the SXRB in the B and C bands were used to constrain their model parameters, cannot fit the observed band ratios, except in the limit in which virtually all of the emission originates in the near side of the absorbing gas. It was also shown that model parameters corresponding to appreciably mixed emission and absorption, as proposed by Bunner et al. (1969) and Davidsen et al. (1972), do not fit the data. Model parameters corresponding to the absorption models of Bowyer et al. (1968), Bunner et al. (1969), Davidsen et al. (1972), and Yentis et al. (1972), cannot fit the data unless the absorbing material is more strongly clumped than is permitted by 21-cm observations.

Burrows, D. N.↗

DOE OSTI2015

Study on hydrogen isotopes permeation in fluidized state of liquid lithium-lead

Lithium-lead (Li-Pb) is one of the most promising candidate materials for the liquid blanket of fusion reactors. Hydrogen transfer under a fluidized condition of Li-Pb is investigated experimentally to design a Li-Pb blanket system. Li-Pb eutectic alloy flows inside a Ni tube in the experimental system, where H{sub 2} permeates into and out of the forced Li-Pb flow. The overall H{sub 2} permeation rate is analyzed using a mass balance model. Hydrogen atoms diffuse in Ni and Li-Pb. The steady-state H{sub 2} permeation rate obtained by this experiment is smaller than the result of the calculation model. A resistance factor is introduced to the present analysis in order to evaluate the influence of other H{sub 2} transfer mechanisms, such as diffusion in Li-Pb and dissolution reaction between Ni and Li-Pb. The contribution of the resistance to the overall H{sub 2} permeation rate becomes large when the flow rate of Li-Pb is low. This is because the boundary layer thickness between Ni and Li-Pb affects the overall H{sub 2} permeation rate. When the flow velocity of Li-Pb increases, the thickness of the boundary layer becomes thin, and the driving force of H{sub 2} permeation through the Ni wall becomes bigger. (authors)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Cation–polymer interactions and local heterogeneity determine the relative order of alkali cation diffusion coefficients in PEGDA hydrogels

Current research efforts are focused on endowing polymer membranes with ion–ion selectivity by incorporating ion–polymer interactions into materials to bias the selective partitioning and or diffusivity of one species over another. However, little is known about the impact of such interactions on the mechanisms of ion transport. In this study, we probe the influence of cation–polymer interactions on cation, anion, and salt diffusivity in a model membrane material, poly(ethylene glycol) diacrylate (PEGDA) by modeling concentrated polyethylene oxide solutions via molecular dynamics simulations. These results are compared to published experimental data for LiCl, NaCl, and KCl diffusion in PEGDA. Experimentally, the order of salt and cation diffusion coefficients for LiCl, NaCl, and KCl deviate from the order in aqueous solutions. Here, simulations identify these deviations to arise from cation–polymer coordination in the membrane. Both the fraction of bound cations and the average binding lifetime increases with decreasing cation hydration free energy (moving down the alkali series), leading to different diffusivity trends in the membrane compared to solution. However, to recover the experimentally observed order of diffusivities cations and salt in our simulations, we needed to incorporate membrane heterogeneity explicitly via a polymer charge scaling procedure. Together, our results indicate that cation–polymer interactions, as well as spatial heterogeneity within the membrane, play a critical role in dictating the observed order of alkali cation and salt diffusion coefficients in membranes.

36 MATERIALS SCIENCE↗

NASA NTRS1979

Characteristics of acoustic-gravity waves in a diffusively separated atmosphere

The effects of wave-induced diffusion on the characteristics of plane acoustic gravity waves in the earth's thermosphere are investigated by means of a two fluid model atmosphere in diffusive equilibrium. At wave periods greater than the shortest characteristic diffusion time for the minor gas, diffusion induced by the wave is considered to act to eliminate amplitude and phase differences between the wave-induced density fluctuations of individual gases and restore the perturbed fluid to diffusive equilibrium. Vertical diffusion is most important at large scale sizes, but horizontal diffusion dominates for horizontal wavelengths of less than several hundred kilometers. As a result of wave-induced diffusion, AE-C satellite measurements of neutral density fluctuations of thermospheric constituents are only compatible with relatively small scale internal gravity waves. Velocity and temperature differences for acoustic-gravity waves are greatest at periods near the mean collision time and diffusion time and lead to the dissipation of wave energy. The predictions of the theory are consistent with a high latitude, low altitude source for the medium scale waves and a more localized random source for the waves observed by AE-C.

Del Genio, A. D.↗

NASA NTRS1985

The 1973-1984 Solar Modulation of Cosmic Ray Nuclei

As a continuation of the program of solar modulation studies, new measurements were carried out with the cosmic ray telescope on the Earth satellite IMP-8, of the intensity time variations and the energy spectra of galactic cosmic ray protons, helium, carbon and oxygen from 1980 through 1984 including the recent solar maximum. In order to test the applicability of a steady state model of solar modulation during a period which includes times of rapidly changing modulation, these fluxes were equated with the predictions of a conventional model of solar modulation which assumes equilibrium between modulation mechanisms. It is found that for a reasonable range of variations of the diffusion coefficient the model predictions can be made to agree with the measurements at essentially all times during the studied period. The model can account also for the observed hysteresis effects between cosmic rays of different rigidities.

Garcia-Munoz, M.↗

DOE OSTI2015

Background model systematics for the Fermi GeV excess

The possible gamma-ray excess in the inner Galaxy and the Galactic center (GC) suggested by Fermi-LAT observations has triggered a large number of studies. It has been interpreted as a variety of different phenomena such as a signal from WIMP dark matter annihilation, gamma-ray emission from a population of millisecond pulsars, or emission from cosmic rays injected in a sequence of burst-like events or continuously at the GC. We present the first comprehensive study of model systematics coming from the Galactic diffuse emission in the inner part of our Galaxy and their impact on the inferred properties of the excess emission at Galactic latitudes 2° < |b| < 20° and 300 MeV to 500 GeV. We study both theoretical and empirical model systematics, which we deduce from a large range of Galactic diffuse emission models and a principal component analysis of residuals in numerous test regions along the Galactic plane. We show that the hypothesis of an extended spherical excess emission with a uniform energy spectrum is compatible with the Fermi-LAT data in our region of interest at 95% CL. Assuming that this excess is the extended counterpart of the one seen in the inner few degrees of the Galaxy, we derive a lower limit of 10.0° (95% CL) on its extension away from the GC. We show that, in light of the large correlated uncertainties that affect the subtraction of the Galactic diffuse emission in the relevant regions, the energy spectrum of the excess is equally compatible with both a simple broken power-law of break energy E{sub break} = 2.1 ± 0.2 GeV, and with spectra predicted by the self-annihilation of dark matter, implying in the case of b-bar b final states a dark matter mass of m{sub χ}=49{sup +6.4}{sub −5.4} GeV.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

36 MATERIALS SCIENCE↗

DOE OSTI2015

HD/H{sub 2} AS A PROBE OF THE ROLES OF GAS, DUST, LIGHT, METALLICITY, AND COSMIC RAYS IN PROMOTING THE GROWTH OF MOLECULAR HYDROGEN IN THE DIFFUSE INTERSTELLAR MEDIUM

We modeled recent observations of UV absorption of HD and H{sub 2} in the Milky Way and toward damped/subdamped Lyα systems at z = 0.18 and z >1.7. N(HD)/N(H{sub 2}) ratios reflect the separate self-shieldings of HD and H{sub 2} and the coupling introduced by deuteration chemistry. Locally, observations are explained by diffuse molecular gas with 16 cm{sup –3} ≲ n(H) ≲ 128 cm{sup –3} if the cosmic-ray ionization rate per H nucleus ζ {sub H} =2 × 10{sup –16} s{sup –1}, as inferred from H{sub 3} {sup +} and OH{sup +}. The dominant influence on N(HD)/N(H{sub 2}) is the cosmic-ray ionization rate with a much weaker downward dependence on n(H) at solar metallicity, but dust extinction can drive N(HD) higher as with N(H{sub 2}). At z > 1.7, N(HD) is comparable to the Galaxy but with 10 times smaller N(H{sub 2}) and somewhat smaller N(H{sub 2})/N(H I). Comparison of our Galaxy with the Magellanic Clouds shows that smaller H{sub 2}/H is expected at subsolar metallicity, and we show by modeling that HD/H{sub 2} increases with density at low metallicity, opposite to the Milky Way. Observations of HD would be explained with higher n(H) at low metallicity, but high-z systems have high HD/H{sub 2} at metallicity 0.04 ≲ Z ≲ 2 solar. In parallel, we trace dust extinction and self-shielding effects. The abrupt H{sub 2} transition to H{sub 2}/H ≈ 1%-10% occurs mostly from self-shielding, although it is assisted by extinction for n(H) ≲ 16 cm{sup –3}. Interior H{sub 2} fractions are substantially increased by dust extinction below ≲ 32 cm{sup –3}. At smaller n(H), ζ {sub H}, small increases in H{sub 2} triggered by dust extinction can trigger abrupt increases in N(HD)

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Influence of Atomic Diffusion on Stellar Ages and Chemical Tagging

In the era of large stellar spectroscopic surveys, there is an emphasis on deriving not only stellar abundances but also the ages for millions of stars. In the context of Galactic archeology, stellar ages provide a direct probe of the formation history of the Galaxy. We use the stellar evolution code MESA to compute models with atomic diffusion—with and without radiative acceleration—and extra mixing in the surface layers. The extra mixing consists of both density-dependent turbulent mixing and envelope overshoot mixing. Based on these models we argue that it is important to distinguish between initial, bulk abundances (parameters) and current, surface abundances (variables) in the analysis of individual stellar ages. In stars that maintain radiative regions on evolutionary timescales, atomic diffusion modifies the surface abundances. We show that when initial, bulk metallicity is equated with current, surface metallicity in isochrone age analysis, the resulting stellar ages can be systematically overestimated by up to 20%. The change of surface abundances with evolutionary phase also complicates chemical tagging, which is the concept that dispersed star clusters can be identified through unique, high-dimensional chemical signatures. Stars from the same cluster, but in different evolutionary phases, will show different surface abundances. We speculate that calibration of stellar models may allow us to estimate not only stellar ages but also initial abundances for individual stars. In the meantime, analyzing the chemical properties of stars in similar evolutionary phases is essential to minimize the effects of atomic diffusion in the context of chemical tagging.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1986

Radial diffusion and ion partitioning in the Io torus

A model is presented for radial diffusion and charge state partitioning of sulfur and oxygen ions in the Io torus, including effects of electron impact and charge exchange. When applied to Voyager 1 radial profiles of total ion flux tube content, the model shows that the ion residence time in the torus, tau(D), as defined in spectroscopic studies of ion partitioning, is related to the radial diffusion coefficient, D(LL), at L = 7 by tau(D) approximately 8/D(LL)(7). This result appears to bring spectroscopic estimates of the ion residence time (tau/D/ greater than about 60 to 100 days) into reasonable agreement with estimates of D(LL) from magnetospheric diffusion studies, D(LL) equals approximately 10 to the -6th/s.

Cheng, A. F.↗

DOE OSTI2016

UO 2 Grain Growth: Developing Phase Field Models for Pore Dragging, Solute Dragging and Anisotropic Grain Boundary Energies

A detailed phase field model for the effect of pore drag on grain growth kinetics was implemented in MARMOT. The model takes into consideration both the curvature-driven grain boundary motion and pore migration by surface diffusion. As such, the model accounts for the interaction between pore and grain boundary kinetics, which tends to retard the grain growth process. Our 2D and 3D simulations demonstrate that the model capture all possible pore-grain boundary interactions proposed in theoretical models. For high enough surface mobility, the pores move along with the migrating boundary as a quasi-rigid-body, albeit hindering its migration rate compared to the pore-free case. For less mobile pores, the migrating boundary can separate from the pores. For the pore-controlled grain growth kinetics, the model predicts a strong dependence of the growth rate on the number of pores, pore size, and surface diffusivity in agreement with theroretical models. An evolution equation for the grain size that includes these parameters was derived and showed to agree well with numerical solution. It shows a smooth transition from boundary-controlled kinetics to pore-controlled kinetics as the surface diffusivity decreases or the number of pores or their size increases. This equation can be utilized in BISON to give accurate estimate for the grain size evolution. This will be accomplished in the near future. The effect of solute drag and anisotropy of grain boundary on grain growth will be investigated in future studies.

36 MATERIALS SCIENCE↗