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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2014

PREPARATION AND TESTING OF CORROSIONAND SPALLATION-RESISTANT COATINGS

This Energy & Environmental Research Center (EERC) project is designed to determine if plating APMT®, a specific highly oxidation-resistant oxide dispersion-strengthened FeCrAl alloy made by Kanthal, onto nickel-based superalloy turbine parts is a viable method for substantially improving the lifetimes and maximum use temperatures of the parts. The method for joining the APMT plate to the superalloys is called evaporative metal bonding and involves placing a thin foil of zinc (Zn) between the plate and the superalloy, clamping them together, and heating in an atmosphere-controlled furnace. Upon heating, the Zn melts and dissolves the oxide skins of the alloys at the bond line, allowing the two alloys to diffuse into each other. The Zn then diffuses through the alloys and evaporates from their surfaces. Laboratory testing to determine the diffusion rate of Zn through the alloys has been completed. However, an analytical solution does not exist to model the diffusion of zinc through the alloys. For this reason, a finite difference algorithm using MATLAB was developed. It makes use of the hopscotch algorithm. The model allows the user to specify the dimensions of the metal parts, the Zn concentration at the bondline, the mesh size, time step, and Zn diffusivity. The experimentally measured values of diffusivity for Zn in APMT and Rene 80/CM 247LC are approximately 2.7 × 10-12 and 4 × 10-14 m2/s, respectively. While the qualitative behavior of the model appears correct, a comparison of the diffusion predictions with the experimental results from earlier in the project indicates that the expected Zn concentration is significantly higher than that measured experimentally. The difference depends on the assumed initial concentration, which is difficult to quantify exactly under experimental conditions for t = 0. In addition to the diffusion work, the coefficients of thermal expansions were determined for each of the alloys as a function of temperature. This information has been entered into a finite element model using ANSYS so that appropriate force-applying structures can be designed for use in joining structures composed of APMT and the nickel alloys. Finite element modeling has been performed to finalize the fabrication geometry for the corrosion-testing phase. The addition of another bolt increases stress uniformity away from the region where the clamping is applied. It appears that a bolt spacing of approximately 25 mm in each jig is appropriate. This will allow the fabrication of 50-mm-wide sections of joints for the corrosion-testing task. Gasifier sampling activities continue to determine what types of trace contaminants may occur in cleaned syngas that could lead to corrosion or deposition in turbines firing coal syngas. The EERC has several pilot-scale gasifiers that are continually used in a variety of test configurations as determined by the needs of the projects that are funding the tests. We are sampling both noncombusted and combusted syngas produced during some of the pilot-scale gasifier tests. This year sampling was performed of both syngas and combusted syngas while the entrained-flow gasifier (EFG) was firing subbituminous coal from the Antelope Mine in Wyoming. Results of scanning electron microscope analyses of the syngas before combustion showed no submicron particles, only flakes of iron oxide that had likely formed on steel surfaces inside the combustor. As shown in the 2013 annual report, soot was also collected from the syngas when the much-lower-temperature fluid-bed gasifier (FBG) was fired, indicating that the much higher temperature of the EFG prevented soot formation. However, particles collected from the combusted syngas consist almost entirely of submicron soot, and little to no vaporized metals made it past the warm-gas filters and scrubbers in the high-temperature EFG system which could then deposit in a turbine system burning a higher hydrogen syngas. These results are consistent with the analyses of the particulates collected from combusted syngas when the lowertemperature FBG system is used.

Electric↗