Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Boundary layer blockage in expansion tube nozzles

The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Sudnitsin, Olga↗
NASA NTRS1990

COS-B gamma ray sources beyond the predicted diffuse emission

COS-B data were reanalyzed using for background subtraction the modeled galactic diffuse gamma-ray emission based on HI- and CO-line surveys and the gamma-ray data itself. A methodology was developed for this purpose with the following three features: automatic generation of source catalogs using correlation analysis, simulation of trials to derive significance thresholds for source detection, and bootstrap sampling to drive error boxes and confidence intervals for source parameters. The analysis shows that about half of the 2CG sources are explained by concentrations in the distribution of molecular hydrogen. Indication for a few weak new sources is also obtained.

Mayer-Hasselwander, H. A.↗
NASA NTRS1992

Cosmic-ray-modified supernova remnant shocks

The evolution of cosmic-ray-modified SNRs during the early Sedov phase is investigated using the time-dependent, two-fluid model for diffusive shock acceleration. Consideration is given to the sensitivity of net acceleration efficiency to model assumptions regarding density structure of the external medium as well as time dependence in the diffusion coefficient and the cosmic-ray adiabatic index. Of the model assumptions explored, the greatest sensitivity was found to time variations in the specific heat ratio to cosmic rays. It was found that if tau is greater than or approximately equal to few times 100, dynamically significant cosmic-ray pressures are produced early in the Sedov inefficient cosmic-ray accelerators. For the SNR models under consideration, the total energy channeled into the cosmic rays can be of order 10 percent of the initial blast energy for both the uniform density and inverse-square density interstellar medium models.

Jones, T. W.↗
DOE OSTI2013

A simplified model for thermal-wave cavity self-consistent measurement of thermal diffusivity

A simplified theoretical model was developed for the thermal-wave cavity (TWC) technique in this study. This model takes thermal radiation into account and can be employed for absolute measurements of the thermal diffusivity of gas and liquid samples without any knowledge of geometrical and thermal parameters of the components of the TWC. Using this model and cavity-length scans, thermal diffusivities of air and distilled water were accurately and precisely measured as (2.191 ± 0.004) × 10{sup −5} and (1.427 ± 0.009) × 10{sup −7} m{sup 2} s{sup −1}, respectively, in very good agreement with accepted literature values.

42 ENGINEERING↗
DOE OSTI2021

Improved Prediction of Cold-Air Pools in the Weather Research and Forecasting Model Using a Truly Horizontal Diffusion Scheme for Potential Temperature

The terrain-following vertical coordinate system used by many atmospheric models, including the Weather Research and Forecasting (WRF) Model, is prone to errors in regions of complex terrain. These errors stem, in part, from the calculation of horizontal gradients within the diffusion term of the momentum or scalar evolution equations. In WRF, such gradients can be calculated along coordinate surfaces, or using metric terms that help account for grid skewness. However, neither of these options ensures a truly horizontal gradient calculation, especially if a grid cell is skewed enough that the heights of the neighboring grid points used in the calculation fall outside the vertical range of the cell. In this work, an improved scheme that uses Taylor series approximations to vertically interpolate variables to the level necessary for a truly horizontal gradient calculation is implemented in WRF for the diffusion of potential temperature. The scheme is validated using an atmosphere-at-rest configuration, in which spurious flows develop only as a result of numerical errors and can thus be used as a proxy for model performance. Following validation, the method is applied to the simulation of cold-air pools (CAPs), which occur in regions of complex terrain and are characterized by strong near-surface temperature gradients. Using the truly horizontal scheme, idealized simulations demonstrate reduced numerical mixing in a quiescent CAP, and a realistic case study in the Columbia River basin shows a reduction in positive wind speed bias by up to roughly 20% compared to observations from the Second Wind Forecast Improvement Project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Mass dependence of overshooting beneath convective envelopes

State of the art stellar evolution codes [Paxton et al., 2010, Demarque et al., 2004, Weiss and Schlattl, 2008, Siess et al., 2013, Christensen-Dalsgaard, 2008] evoke a diffusive process to model convective overshooting. The diffusion coefficient [Freytag et al., 1996, Pratt et al., 2017] used in the model can be set to change with the classification of the convective zone as non-burning, H-burning, He-burning, or metal-burning. It can also be set to change at defined evolutionary points, such as the bottom of the asymptotic giant branch, or during the third dredge up [Herwig, 2000, Lugaro et al., 2003]. Aside from these abrupt changes, the diffusion coefficient is typically locked to a percentage of the pressure scale height measured at the convective boundary. However, there is no theoretical reason for convective overshooting, or indeed other convective properties, to change in the same way that the pressure scale height changes as a star evolves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The GeV Excess Shining Through: Background Systematics for the Inner Galaxy Analysis

Recently, a spatially extended excess of gamma rays collected by the Fermi-LAT from the inner region of the Milky Way has been detected by different groups and with increasingly sophisticated techniques. Yet, any final conclusion about the morphology and spectral properties of such an extended diffuse emission are subject to a number of potentially critical uncertainties, related to the high density of cosmic rays, gas, magnetic fields and abundance of point sources. We will present a thorough study of the systematic uncertainties related to the modelling of diffuse background and to the propagation of cosmic rays in the inner part of our Galaxy. We will test a large set of models for the Galactic diffuse emission, generated by varying the propagation parameters within extreme conditions. By using those models in the fit of Fermi-LAT data as Galactic foreground, we will show that the gamma-ray excess survives and we will quantify the uncertainties on the excess emission morphology and energy spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

ICAT: A numerical scheme to minimize numerical diffusion in advection-dispersion modeling and its application in identifying flow channeling

The advection-dispersion equation for scalar transport is essential for the numerical modeling of many fluid dynamics problems. However, solutions from numerical schemes always suffer from numerical diffusion and/or oscillation. In this study, we develop an Intra-Cell Advection Tracking (ICAT) scheme to minimize numerical diffusion and preserve monotonicity for advection-dispersion modeling. The key idea is to introduce “queues” in each discretized cell, and using a sequential transport rule and a flow distribution mechanism to track the scalar transport in these queues temporally and spatially. The capability and limitations of ICAT are first investigated through three test cases. Compared with the results obtained from other numerical schemes, the results from ICAT show substantially reduced numerical diffusion and agree better with analytical solutions. We also employ ICAT to simulate the transport process of a conservative tracer in a fracture with a highly heterogeneous aperture distribution. We report discrete flow channels in the fracture are better discerned by ICAT than by other numerical schemes, indicating the suitability of ICAT for modeling tracer transport in channelized flow fields.

15 GEOTHERMAL ENERGY↗
NASA NTRS1994

Radial diffusion of relativistic electrons in Neptune's magnetosphere

Measurements of greater than or = 1 MeV electrons in Neptune's magnetosphere from L approximately = 4 to 14 are well represented by solutions of a model radial diffusion equation that includes losses due to absorption by the satellite Proteus. The model provides estimates of the radial diffusion coefficient, which is 5 x 10(exp -8) (L/5)(exp 8)/s, and the outer boundary energy spectrum, which is an exponential with an e-folding energy of 0.1 MeV. The diffusion coefficient is consistent with theoretical estimates based on the assumption that diffusion is driven by atmospheric neutral winds.

Selesnick, R. S.↗
DOE OSTI2026

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗
DOE OSTI2013

On the asymptotic preserving property of the unified gas kinetic scheme for the diffusion limit of linear kinetic models

The unified gas kinetic scheme (UGKS) of K. Xu et al. (2010) [37], originally developed for multiscale gas dynamics problems, is applied in this paper to a linear kinetic model of radiative transfer theory. While such problems exhibit purely diffusive behavior in the optically thick (or small Knudsen) regime, we prove that UGKS is still asymptotic preserving (AP) in this regime, but for the free transport regime as well. Moreover, this scheme is modified to include a time implicit discretization of the limit diffusion equation, and to correctly capture the solution in case of boundary layers. Contrary to many AP schemes, this method is based on a standard finite volume approach, it does neither use any decomposition of the solution, nor staggered grids. Several numerical tests demonstrate the properties of the scheme.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2019

MODIS Aqua Reflective Solar Bands Calibration for NASA's R2018 Ocean Color Products

Remote sensing ocean color products have stringent requirements on radiometric calibration stability. To address a calibration deficiency in MODIS Aqua in recent years, the NASA Ocean Biological Processing Group (OBPG) developed a new calibration for the reflective solar bands. Prior to the reprocessing of NASA's ocean color products for 2018 (R2018), the OBPG MODIS products had been based on calibration provided by the MODIS Calibration Support Team (MCST). Several modifications were made to the MCST calibration approach to improve the calibration accuracy for ocean color products. These include: 1) apply 936nm detector normalization to solar diffuser stability monitor (SDSM) data to reduce coherent noise; 2) modeling solar diffuser (SD) degradation wavelength dependency to determine SD degradation in near-infrared and shortwave infrared wavelength; 3) detector gains are computed using SD screen-closed data to better match ocean radiance levels in all bands; 4) a simple atmospheric correction is performed to reduce bidirectional reflectance distribution function (BRDF) effects in desert trends; 5) modulated relative spectral response (RSR) impact on ocean data is estimated and used to adjust the calibration coefficients. 6) smoothing is used to characterize the temporal change in calibration, and response versus scan angle (RVS) changes are characterized using 2nd order polynomials to improve spatial/temporal calibration stability. Relative to the previous R2014 ocean color products, the R2018 calibration removed the suspect late-mission global trends in blue-band water-leaving reflectance, and some anomalously large short-term variability (spikes) in the temporal trend of chlorophyll concentration. This paper will describe the OBPG calibration with a focus on the differences between the MCST and OBPG approaches.

calibration↗
NASA NTRS1993

Alveolar ventilation to perfusion heterogeneity and diffusion impairment in a mathematical model of gas exchange

This study describes a two-compartment model of pulmonary gas exchange in which alveolar ventilation to perfusion (VA/Q) heterogeneity and impairment of pulmonary diffusing capacity (D) are simultaneously taken into account. The mathematical model uses as input data measurements usually obtained in the lung function laboratory. It consists of two compartments and an anatomical shunt. Each compartment receives fractions of alveolar ventilation and blood flow. Mass balance equations and integration of Fick's law of diffusion are used to compute alveolar and blood O2 and CO2 values compatible with input O2 uptake and CO2 elimination. Two applications are presented. The first is a method to partition O2 and CO2 alveolar-arterial gradients into VA/Q and D components. The technique is evaluated in data of patients with chronic obstructive pulmonary disease (COPD). The second is a theoretical analysis of the effects of blood flow variation in alveolar and blood O2 partial pressures. The results show the importance of simultaneous consideration of D to estimate VA/Q heterogeneity in patients with diffusion impairment. This factor plays an increasing role in gas alveolar-arterial gradients as severity of COPD increases. Association of VA/Q heterogeneity and D may produce an increase of O2 arterial pressure with decreasing QT which would not be observed if only D were considered. We conclude that the presented computer model is a useful tool for description and interpretation of data from COPD patients and for performing theoretical analysis of variables involved in the gas exchange process.

NASA Discipline Environmental Health↗
NASA NTRS1983

Diffusion processes in Al2O3 scales - Void growth, grain growth, and scale growth

The internal microstructure and growth kinetics of Al2O3 scales on Ni-15Cr-13Al (wt percent) are investigated by TEM and analyzed in relation to models of diffusivity. Polished arc-melted specimens were oxidized in 1-atm air at 1100 C for 0.1, 1.0, and 20 hours and ion-thinned for TEM at 100 kV. The frequency distribution of void size and grain size is determined for different oxidation times and scale depths. The kinetics of microvoid growth and of grain and scale growth are plotted and related via simplified models to lattice and grain-boundary oxygen diffusivity, respectively. Good agreement is found between model predictions and data obtained by Oishi and Kingery (1960) on oxygen diffusion in bulk Al2O3. The further implications and limitations of these findings are discssed.

Smialek, J. L.↗
DOE OSTI2018

CFD Modeling of Boron Removal from Liquid Silicon with Cold Gases and Plasma

The present study focuses on a specific step of the metallurgical path of purification to provide solar-grade silicon: the removal of boron through the injection of H{sub 2}O(g)-H{sub 2}(g)-Ar(g) (cold gas process) or of Ar-H{sub 2}-O{sub 2} plasma (plasma process) on stirred liquid silicon. We propose a way to predict silicon and boron flows from the liquid silicon surface by using a CFD model (©Ansys Fluent) combined with some results on one-dimensional diffusive-reactive models to consider the formation of silica aerosols in a layer above the liquid silicon. The comparison of the model with experimental results on cold gas processes provided satisfying results for cases with low and high concentrations of oxidants. This confirms that the choices of thermodynamic data of HBO(g) and the activity coefficient of boron in liquid silicon are suitable and that the hypotheses regarding similar diffusion mechanisms at the surface for HBO(g) and SiO(g) are appropriate. The reasons for similar diffusion mechanisms need further enquiry. We also studied the effect of pressure and geometric variations in the cold gas process. For some cases with high injection flows, the model slightly overestimates the boron extraction rate, and the overestimation increases with increasing injection flow. A single plasma experiment from SIMaP (France) was modeled, and the model results fit the experimental data on purification if we suppose that aerosols form, but it is not enough to draw conclusions about the formation of aerosols for plasma experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2015

ANALYSIS OF MAGNETOROTATIONAL INSTABILITY WITH THE EFFECT OF COSMIC-RAY DIFFUSION

We present the results obtained from the linear stability analysis and 2.5 dimensional magnetohydrodynamic (MHD) simulations of magnetorotational instability (MRI), including the effects of cosmic rays (CRs). We took into account the CR diffusion along the magnetic field but neglected the cross-field-line diffusion. Two models are considered in this paper: the shearing box model and differentially rotating cylinder model. We studied how MRI is affected by the initial CR pressure (i.e., energy) distribution. In the shearing box model, the initial state is uniform distribution. Linear analysis shows that the growth rate of MRI does not depend on the value of the CR diffusion coefficient. In the differentially rotating cylinder model, the initial state is a constant angular momentum polytropic disk threaded by a weak uniform vertical magnetic field. Linear analysis shows that the growth rate of MRI becomes larger if the CR diffusion coefficient is larger. Both results are confirmed by MHD simulations. The MHD simulation results show that the outward movement of matter by the growth of MRI is not impeded by the CR pressure gradient, and the centrifugal force that acts on the concentrated matter becomes larger. Consequently, the growth rate of MRI is increased. On the other hand, if the initial CR pressure is uniform, then the growth rate of the MRI barely depends on the value of the CR diffusion coefficient.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Atoms in carbon cages as a source of interstellar diffuse lines

A model to describe the resonance absorption lines of various atoms trapped in closed carbon cages is presented. These systems may be responsible for some of the as yet unexplained diffuse interstellar bands. Model potentials for possible atom-C60 systems are obtained and used to calculate the resonance lines. The trapped atoms considered are O, N, Si, Mg, Al, Na, and S, and in all cases the resonance lines are shifted toward the red as compared to the isolated atoms. The calculated wavelengths are compared to the range of wavelengths observed for the diffuse interstellar bands, and good agreement is found for Mg and Si resonance lines. Other lines may be caused by other than resonance transitions or by trapped molecules. The oscillator strengths and the abundances are evaluated and compared with observation. Mechanisms to explain the observed band width of the lines and the existence of certain correlated pairs of lines are discussed.

Ballester, J. L.↗