Search NASA⌕ Search

SEARCH · Search NASA

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,549 records · Page 86

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

An iterative radiative transfer code for ocean-atmosphere systems

The details of an iterative radiative transfer code for computing the intensity and degree of polarization of diffuse radiation in models of the ocean-atmosphere system are described. The present code neglects the upwelling radiation from below the ocean surface and as such can be applied to the part of the spectrum where the absorption by water is strong. To establish the reliability of the numerical scheme and the computer code, the results are compared with those of Fraser and Walker (1968), Dave (1972), and Mullamaa (1964); they are found to be in excellent agreement. The computations also show that both the intensity and the degree of polarization of the upwelling diffuse radiation at the top of the atmosphere vary significantly when the rough ocean at the base of the atmosphere is replaced by a Lambertian surface that reflects the same energy as the rough ocean.

Ahmad, Z.↗

Gravitational effects on chemically reacting laminar boundary layer flows over a horizontal flat plate

A theoretical study of the chemically reacting laminar boundary layer flow over a horizontal flat plate with gravitationally induced buoyant force is presented. A diffusion flame sheet model was used to describe the combustion process. The effects of gravity on the purely force convection flow can be characterized by a dimensionless coordinate quantity, which is involved in the generation of the governing equations. A numerical solution of the zero and first order governing equations subject to the appropriate physical boundary conditions was obtained. It is shown that the cross stream buoyancy induced body force acts effectively to produce a streamwise pressure gradient in the fluid adjacent to the plate surface. It is concluded that buoyancy plays an important role in boundary layer diffusion flames.

Lavid, M.↗

Bioethanol enrichment using zeolite membranes: Molecular modeling, conceptual process design and techno-economic analysis

Lignocellulosic biomass-derived ethanol offers a potential alternative to fossil-derived fuels. However, the energy-intensive nature of its recovery process, i.e., a sequence of two distillation columns – beer and rectification columns - enriching from 5 to 37, and 37–93 wt%, respectively, limits its environmental and economic benefits. Here, in this paper, we assess the potential of using zeolite membranes for bioethanol recovery. Atomistic-level and molecular dynamic simulations are performed to determine adsorption and diffusion properties of the ethanol-water mixture in zeolite membranes. These properties are modeled using the real adsorption solution theory and the Maxwell-Stefan equations to describe permeation through the zeolite membranes. A comparison of steady state permeance and selectivity predicted by the model with the values from experiments suggests that the simulated membranes are more permeable and selective than the real membranes. This is attributed to the presence of structural non-idealities and hydrophilic defects in real membranes, while the adsorption and diffusion properties obtained using molecular simulations reflect the behavior in ideal crystals. Thus, a reduced diffusivity model with an empirical relation for enhanced water adsorption is used to capture similar performance as that obtained by real membranes. This model is further used in developing conceptual process designs to assess the viability of zeolite membranes for bioethanol enrichment in industry. Both hydrophobic and hydrophilic zeolite membranes are considered. Hydrophobic zeolite membranes show potential for energy savings but lack in separation performance. On the other hand, hydrophilic zeolite membranes can achieve the separation target but result in no energy savings. Thus, a configuration that uses a combination of hydrophobic and hydrophilic membranes is proposed. It can achieve the separation target and results in 15% energy savings over distillation. Techno-economic analysis suggests that ~ 10-fold improvements in permeation or equivalent cost reductions are required for economic viability of this scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reaction rates for reaction-diffusion kinetics on unstructured meshes

The reaction-diffusion master equation is a stochastic model often utilized in the study of biochemical reaction networks in living cells. It is applied when the spatial distribution of molecules is important to the dynamics of the system. A viable approach to resolve the complex geometry of cells accurately is to discretize space with an unstructured mesh. Diffusion is modeled as discrete jumps between nodes on the mesh, and the diffusion jump rates can be obtained through a discretization of the diffusion equation on the mesh. Reactions can occur when molecules occupy the same voxel. Here in this paper, we develop a method for computing accurate reaction rates between molecules occupying the same voxel in an unstructured mesh. For large voxels, these rates are known to be well approximated by the reaction rates derived by Collins and Kimball, but as the mesh is refined, no analytical expression for the rates exists. We reduce the problem of computing accurate reaction rates to a pure preprocessing step, depending only on the mesh and not on the model parameters, and we devise an efficient numerical scheme to estimate them to high accuracy. We show in several numerical examples that as we refine the mesh, the results obtained with the reaction-diffusion master equation approach those of a more fine-grained Smoluchowski particle-tracking model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulated interannual variability of the Greenland Sea deep water formation and its connection to surface forcing

A fully prognostic Arctic ice-ocean model is used to study the interannual variability of deepwater formation in the Greenland Sea Gyre based on the simulations for the Arctic ice-ocean system for the period 1955 and 1960 - 1985. The model uses monthly climatology for thermodynamic forcing components (such as air temperature and cloudiness), together with constant annual net precipitation and river runoff. The daily wind forcing is derived from analyzed sea level air pressures from the National Center for Atmospheric Research (NCAR). In summary, the model shows that the occurence of deep convection in the Greenland Sea Gyre is controlled by the extensive Fram Strait ice export and/or local wind conditions in the Greenland Sea. In the latter case the weakening of the local wind curl allows the Polar Front to move eastward. The movement of the Polar Front causes adverse ice conditions, often together with much larger than normal ice export from the Arctic, such as in 1968, which can block convection in the gyre. The density difference between upper and lower layers is investigated as an indication of water mass formation through convection, occurring as strong diffusion in the model. The model-simulated density difference between the average top 100 m and deep levels reveals that the period 1960 - 1985 had only a few distinct years with weak stratification, and, especially, the model predicts no deep convection since the nid-1970s. The common factor for the years of the weakest decrease of the model-predicted heat content of the upper 2000 m which can, to a high degree, be explained by local heat loss.

Haekkinen, Sirpa↗

Quantitative characterization of soil micro-aggregates: New opportunities from sub-micron resolution synchrotron X-ray microtomography

Soil microaggregates are the fundamental building block, at the micron scale, of the highly hierarchical structure of soils, and can exert a significant control on the local biological metabolism and microbial community partitioning. In this study we propose an analysis protocol for the morphometric characterization of complete soil microaggregates based on sub-micron resolution synchrotron X-ray microtomography. A comprehensive characterization of the aggregate morphology is the first step towards a complete characterization of the soil microaggregates, when trying to correlate morphometric parameters with physical and/or biological properties, or when building models (e.g., effective diffusivity, microbial distribution, etc.). In this study, we demonstrate our characterization approach on two single microaggregate samples from dramatically different soil environments: one from Kansas, primarily composed by inorganic particles, and one from Barrow (Alaska) dominated by plant fragments. A series of state-of-the-art morphometric analysis techniques have been employed providing quantitative results highlighting specific differences of the two samples. The role of the microstructure in a scenario microbial population development has been discussed and it has been found that the Barrow microaggregate seems to be more favorable, from a purely geometrical point of view, as also confirmed by a simple model presented in this work. Finally, the potential of this approach, when coupled with chemical and biological analysis for a fully comprehensive characterization of soil aggregates in the larger picture of enhanced biological activity, is evident.

54 ENVIRONMENTAL SCIENCES↗

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗

The Importance of a Moving Boundary Approach for Modeling the SEI Layer Growth to Predict Capacity Fade

One of the contributing factors to the aging of lithium-ion batteries is the growth of the solid-electrolyte interphase (SEI) layer. The growth of the SEI layer leads to the irreversible loss of lithium available for cycling and increases the resistance of the battery. Physics-based models in literature model the kinetically limited or solvent diffusion-limited growth. In such models, the interface resistance is a constant, and the contribution to the overpotential of the intercalation reaction from the SEI layer is considered to be ohmic. In this study, we propose a model that describes the growth of the SEI layer on the electrode surface as a moving interface. The transport of lithium ions and the solvent in the electrolyte are affected by this moving interface. The equations that govern the species transport and the potential drop across the SEI layer are derived from dilute solution theory and solved by transforming the coordinates of the moving boundary. The ion transport induces changes in the conductivity across the SEI layer, which affects the potential drop that arises due to its growth. The effects of this potential on capacity fade are studied over cycling the battery.

25 ENERGY STORAGE↗

Monte Carlo simulation of reentry flows with ionization

The Direct Simulation Monte Carlo method is applied to a rarefied, weakly ionized, hypersonic flow over a blunt axisymmetric body. An ionization model based on the concept of ambipolar diffusion is used and a model for the sheath is presented. The effects of the new modeling techniques are investigated for flow over the Project Fire II configuration at 11.37 km/s at an altitude of 84.6 km. The calculated results are presented and compared with both experimental data and solutions where ionization effects were not included. In general, the calculated results overpredict the experimental values by about 15-20 percent.

Taylor, Jeff C.↗

Geographic variations in the thermal and diffusive stability of ground ice on Mars

Time-dependent models of the thermal and diffusive behavior of the Martian regolith are developed to investigate the stability of ground ice within the top several meters of the regolith. The geographic distribution of thermal inertia and albedo as well as the latitudinal variation in insolation are included in calculations of surface and subsurface temperatures between +/- 60 deg latitude. Ground ice is found to be stable where the annual mean surface and subsurface temperatures were below the atmospheric frost point. Calculations of temperatures at high and low obliquity suggest that ground ice would be stable globally at high obliquity and would not be stable between +/- 60 deg latitude at low obliquity. The time scales for condensation of ground ice are found to be comparable to that of orbital oscillations, suggesting that the present geographic distribution of ground ice may depend on the orbital history of Mars.

Mellon, Michael T.↗

Comparison of PDF and Moment Closure Methods in the Modeling of Turbulent Reacting Flows

In modeling turbulent reactive flows, Probability Density Function (PDF) methods have an advantage over the more traditional moment closure schemes in that the PDF formulation treats the chemical reaction source terms exactly, while moment closure methods are required to model the mean reaction rate. The common model used is the laminar chemistry approximation, where the effects of turbulence on the reaction are assumed negligible. For flows with low turbulence levels and fast chemistry, the difference between the two methods can be expected to be small. However for flows with finite rate chemistry and high turbulence levels, significant errors can be expected in the moment closure method. In this paper, the ability of the PDF method and the moment closure scheme to accurately model a turbulent reacting flow is tested. To accomplish this, both schemes were used to model a CO/H2/N2- air piloted diffusion flame near extinction. Identical thermochemistry, turbulence models, initial conditions and boundary conditions are employed to ensure a consistent comparison can be made. The results of the two methods are compared to experimental data as well as to each other. The comparison reveals that the PDF method provides good agreement with the experimental data, while the moment closure scheme incorrectly shows a broad, laminar-like flame structure.

Norris, Andrew T.↗

The large scale gas and dust distribution in the galaxy: Implications for star formation

Infrared Astronomy Observations are presented for the diffuse infrared (IR) emissions from the galactic plane at wavelengths of 60 and 100 microns and the total far infrared intensity and its longitudinal variations in the disk were derived. Using available CO, 5 GHz radio-continuum, and HI data, the IR luminosity per hydrogen mass and the ingrared excess (IRE) ratio in the Galaxy were derived. The longitudinal profiles of the 60 and 100 micron emission were linearly decomposed into three components that are associated with molecular (H2), neutral (HI), and ionized (HII) phases in the interstellar medium (ISM), and the relevant dust properties were derived in each phase. Implications of the findings for various models of the diffuse IR emisison and for star formation in the galactic disk are discussed.

Sodroski, T. J.↗

A study of plasmaspheric density distributions for diffusive equilibrium conditions

The plasmaspheric density distribution has been modeled for a range of solar cycle, seasonal and diurnal conditions with a magnetic flux tube dependent diffusive equilibrium model by using experimentally determined values of ionospheric parameters at 675 km as boundary conditions. Data is presented in terms of plasmaspheric H(+) and He(+) density contours, total flux tube content, and equatorial plasma density for a range of L-values from 1.15 to 3.0. The variation of equatorial density with L-value shows good agreement with the 1/L exp 4 dependence observed experimentally. The results show that the model predicts larger solar cycle and diurnal variation in equatorial plasma density than observed using whistler techniques. However, the whistler method requires a model to deduce the equatorial density and is therefore open to interpretation. Seasonal variations are rather artificial since in this general model no attempt has been made to match equatorial densities for flux tubes emanating from the winter and summer hemispheres.

Li, W.↗

Analytic Couple Modeling Introducing Device Design Factor, Fin Factor, Thermal Diffusivity Factor, and Inductance Factor

A set of convenient thermoelectric device solutions have been derived in order to capture a number of factors which are previously only resolved with numerical techniques. The concise conversion efficiency equations derived from governing equations provide intuitive and straight-forward design guidelines. These guidelines allow for better device design without requiring detailed numerical modeling. The analytical modeling accounts for factors such as i) variable temperature boundary conditions, ii) lateral heat transfer, iii) temperature variable material properties, and iv) transient operation. New dimensionless parameters, similar to the figure of merit, are introduced including the device design factor, fin factor, thermal diffusivity factor, and inductance factor. These new device factors allow for the straight-forward description of phenomenon generally only captured with numerical work otherwise. As an example a device design factor of 0.38, which accounts for thermal resistance of the hot and cold shoes, can be used to calculate a conversion efficiency of 2.28 while the ideal conversion efficiency based on figure of merit alone would be 6.15. Likewise an ideal couple with efficiency of 6.15 will be reduced to 5.33 when lateral heat is accounted for with a fin factor of 1.0.

Seebeck effect↗

Modeling the stellar contribution to the Galactic component of the diffuse soft X-ray background. I - Background fluxes and number counts

A combination of a stellar Galaxy model based on optical data, stellar X-ray luminosity functions derived from the full Einstein base, and a model for X-ray absorption derived from hydrogen column densities are used to estimate the contribution to the diffuse soft X-ray background flux from the main-sequence A, F, G, K, and M stars, and RS CVn stars, at various energies ranging from 0.1 to about 5 keV. Previous discrepancies between earlier estimates of the stellar contribution to the diffuse soft X-ray background are resolved; this stellar contribution to the diffuse soft X-ray background is found to be under 3 percent for photon energies less than energy I and J bands (about 0.8-2 keV), at a threshold sensitivity for point source detection about 10 exp -10 ergs/s sq cm. At low latitudes, stellar contribution estimates are less than 3 percent below 0.3 keV, 7-40 percent in the medium-energy bands, and 27-70 percent in the I and J bands. It is shown that while dM stars are the major contributors to the diffuse stellar flux, other stellar types contribute as much as 40 percent of this flux at the higher energies in the passband studied.

Kashyap, V.↗

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗