Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Isopycnal mixing by mesoscale eddies significantly impacts oceanic anthropogenic carbon uptake

Abstract Anthropogenic carbon dioxide uptake varies across Earth System Models for reasons that have remained obscure. When varied within a single model, the lateral eddy mixing coefficient A Redi produces a range of uptake similar to the modeled range. The highest uptake, resulting from a simulation with a constant A Redi of 2400 m 2 /s, simulates 15% more historical carbon uptake than a model with A Redi = 400 m 2 /s. A sudden doubling in carbon dioxide produces a 21% range in carbon uptake across the models. Two spatially dependent representations of A Redi produce uptake that lies in the middle of the range of constant values despite predicting very large values in the subtropical gyres. One‐dimensional diffusive models of the type used for integrated assessments can be fit to the simulations, with A Redi accounting for a substantial fraction of the effective vertical diffusion. Such models, however, mask significant regional changes in stratification and biological carbon storage.

Gnanadesikan, Anand↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2021

Modeling fission product diffusion in TRISO fuel particles with BISON

Diffusion of fission products in intact TRISO particles depends on particle geometry, fission product source rates, time, temperature, and temperature-dependent diffusion coefficients. Simulating this diffusion process requires models for source rates and diffusion coefficients, plus computation of the temperature field if not prescribed. In addition, simulation quality depends on discretization of the geometry, appropriate time stepping, and the accuracy of the solution method. In this paper, we explore the simulation of fission product diffusion in TRISO fuel particles using the finite element method via the fuel performance code Bison. Recent material model development has occurred in Bison for each material present in tri-structural isotropic (TRISO) fuel particles: the buffer, inner pyrolytic carbon, silicon carbide, and outer pyrolytic carbon layers, as well as the fuel kernel. Also, new mesh generation and fission product release fraction capabilities have been added. Diffusion capabilities are shown to converge to the correct solution via formal verification tests. A large number of code benchmarking problems are also given, with good results, showing that Bison’s computed release fractions closely match those of other software tools. Finally, a significant validation effort is detailed in which fission product release, measured as part of the AGR-1 capsule experiments, is compared to Bison outputs. Bison outputs compare very well to the experimental data and to PARFUME results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the Flow of a Cement Suspension: The Effects of Nano-Silica and Fly Ash Particles

Additives such as nano-silica and fly ash are widely used in cement and concrete materials to improve the rheology of fresh cement and concrete and the performance of hardened materials and increase the sustainability of the cement and concrete industry by reducing the usage of Portland cement. Therefore, it is important to study the effect of these additives on the rheological behavior of fresh cement. In this paper, we study the pulsating Poiseuille flow of fresh cement in a horizontal pipe by considering two different additives and when they are combined (nano-silica, fly ash, combined nano-silica, and fly ash). To model the fresh cement suspension, we used a modified form of the power-law model to demonstrate the dependency of the cement viscosity on the shear rate and volume fraction of cement and the additive particles. The convection–diffusion equation was used to solve for the volume fraction. After solving the equations in the dimensionless forms, we conducted a parametric study to analyze the effects of nano-silica, fly ash, and combined nano-silica and fly ash additives on the velocity and volume fraction profiles of the cement suspension. According to the parametric study presented here, larger nano-silica content results in lower centerline velocity of the cement suspension and larger non-uniformity of the volume fraction. Compared to nano-silica, fly ash exhibits an opposite effect on the velocity. Larger fly ash content results in higher centerline velocity, while the effect of the fly ash on the volume fraction is not obvious. For cement suspension containing combined nano-silica and fly ash additives, nano-silica plays a dominant role in the flow behavior of the suspension. The findings of the study can help the design and operation of the pulsating flow of fresh cement mortars and concrete in the 3D printing industry.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Structure of Laminar Permanently Blue, Opposed-Jet Ethylene-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently blue) ethylene-fueled diffusion flames at various strain rates were studied both experimentally and computationally using an opposed-jet configuration. Measurements of gas velocities, temperatures, and compositions were carried out along the stagnation stream line. Corresponding predictions of flame structure were obtained, based on numerical simulations using several contemporary reaction mechanisms for methane oxidation. Flame conditions studied included ethylene-fueled opposed-jet diffusion flames having stoichiometric mixture fractions of 0.7 with measurements involving strain rates of 60-240/s and predictions involving strain rates of 0-1140/s at normal temperature and pressure. It was found that measured major gas species concentrations and temperature distributions were in reasonably good agreement with predictions using mechanisms due to GRI-Mech and Peters and that effects of preferential diffusion significantly influence flame structure even when reactant mass diffusivities are similar. Oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet, decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures were found to exist over a broad range of strain rates, providing potential for significant computational simplifications for modeling purposes in some instances.

Lin, K.-C.↗
NASA NTRS2000

Structure of Laminar Permanently Blue, Opposed-Jet Ethylene-Fueled Diffusion Flames

The structure and state relationships of laminar soot-free (permanently blue) ethylene-fueled diffusion flames at various strain rates were studied both experimentally and computationally using an opposed-jet configuration. Measurements of gas velocities, temperatures, and compositions were carried out along the stagnation stream line. Corresponding predictions of flame structure were obtained, based on numerical simulations using several contemporary reaction mechanisms for methane oxidation. Flame conditions studied included ethylene-fueled opposed-jet diffusion flames having stoichiometric mixture fractions of 0.7 with measurements involving strain rates of 60-240/s and predictions involving strain rates of 0-1140/s at normal temperature and pressure. It was found that measured major gas species concentrations and temperature distributions were in reasonably good agreement with predictions using mechanisms due to GRI-Mech and Peters and that effects of preferential diffusion significantly influence flame structure even when reactant mass diffusivities are similar. Oxygen leakage to fuel-rich conditions and carbon monoxide leakage to fuel-lean conditions both increased as strain rates increased. Furthermore, increased strain rates caused increased fuel concentrations near the flame sheet, decreased peak gas temperatures, and decreased concentrations of carbon dioxide and water vapor throughout the flames. State relationships for major gas species and gas temperatures were found to exist over a broad range of strain rates, providing potential for significant computational simplifications for modeling purposes in some instances.

Lin, K.-C.↗
DOE OSTI2022

A window into atmospheric escape on early Mars provided by argon isotopes in a Martian meteorite

Novel analytical approaches to determining the trapped 38 Ar/ 36 Ar ratio of gases contained within Martian meteorites are presented and applied to the Martian regolith breccia Northwest Africa (NWA) 7034 and paired stone NWA 11220. The resulting data indicate that extensive mass-dependent fractionation of atmospheric Ar may have occurred within 150 million years of planetary formation, ostensibly as a result of diffusion-limited hydrodynamic escape. The inferred fractional loss of Ar and lighter atmospheric constituents exceeds 50%. In conclusion, these data suggest that volatiles derived from planetary outgassing and/or impactors may dominate the present abundance of atmospheric Ar.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

A Herschel/HIFI legacy survey of HF and H{sub 2}O IN the galaxy: Probing diffuse molecular cloud chemistry

We combine Herschel observations for a total of 12 sources to construct the most uniform survey of HF and H{sub 2}O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-infrared (HIFI) allows us to compare the HF and H{sub 2}O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H{sub 2}O velocity distributions follow each other almost perfectly and establish that HF and H{sub 2}O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H{sub 2} in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H{sub 2} in our sample, we find that the N(H{sub 2}O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H{sub 2}O might be used as a tracer of H{sub 2}—within a factor of 2.5—in the diffuse interstellar medium (ISM). We show that the measured factor of ∼2.5 variation around the median is driven by true local variations in the H{sub 2}O abundance relative to H{sub 2} throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H{sub 2}O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse ISM (n{sub H} ⩽ 1000 cm{sup −3}).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Assessing Jet-Induced Spatial Mixing in a Rich, Reacting Crossflow

In many advanced low NOx gas turbine combustion techniques, such as rich-burn/quick-mix/lean-burn (RQL), jet mixing in a reacting, hot, fuel-rich crossflow plays an important role in minimizing all pollutant emissions and maximizing combustion efficiency. Assessing the degree of mixing and predicting jet penetration is critical to the optimization of the jet injection design strategy. Different passive scalar quantities, including carbon, oxygen, and helium are compared to quantify mixing in an atmospheric RQL combustion rig under reacting conditions. The results show that the O2-based jet mixture fraction underpredicts the C-based mixture fraction due to jet dilution and combustion, whereas the He tracer overpredicts it possibly due to differences in density and diffusivity. The He-method also exhibits significant scatter in the mixture fraction data that can most likely be attributed to differences in gas density and turbulent diffusivity. The jet mixture fraction data were used to evaluate planar spatial unmixedness, which showed good agreement for all three scalars. This investigation suggests that, with further technique refinement, O2 or a He tracer could be used instead of C to determine the extent of reaction and mixing in an RQL combustor.

Demayo, T. N.↗
DOE OSTI2015

Theory of activated penetrant diffusion in viscous fluids and colloidal suspensions

We heuristically formulate a microscopic, force level, self-consistent nonlinear Langevin equation theory for activated barrier hopping and non-hydrodynamic diffusion of a hard sphere penetrant in very dense hard sphere fluid matrices. Penetrant dynamics is controlled by a rich competition between force relaxation due to penetrant self-motion and collective matrix structural (alpha) relaxation. In the absence of penetrant-matrix attraction, three activated dynamical regimes are predicted as a function of penetrant-matrix size ratio which are physically distinguished by penetrant jump distance and the nature of matrix motion required to facilitate its hopping. The penetrant diffusion constant decreases the fastest with size ratio for relatively small penetrants where the matrix effectively acts as a vibrating amorphous solid. Increasing penetrant-matrix attraction strength reduces penetrant diffusivity due to physical bonding. For size ratios approaching unity, a distinct dynamical regime emerges associated with strong slaving of penetrant hopping to matrix structural relaxation. A crossover regime at intermediate penetrant-matrix size ratio connects the two limiting behaviors for hard penetrants, but essentially disappears if there are strong attractions with the matrix. Activated penetrant diffusivity decreases strongly with matrix volume fraction in a manner that intensifies as the size ratio increases. We propose and implement a quasi-universal approach for activated diffusion of a rigid atomic/molecular penetrant in a supercooled liquid based on a mapping between the hard sphere system and thermal liquids. Calculations for specific systems agree reasonably well with experiments over a wide range of temperature, covering more than 10 orders of magnitude of variation of the penetrant diffusion constant.

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

The effects of complex chemistry on triple flames

The structure, ignition, and stabilization mechanisms for a methanol (CH3OH)-air triple flame are studied using Direct Numerical Simulations (DNS). The methanol (CH3OH)-air triple flame is found to burn with an asymmetric shape due to the different chemical and transport processes characterizing the mixture. The excess fuel, methanol (CH3OH), on the rich premixed flame branch is replaced by more stable fuels CO and H2, which burn at the diffusion flame. On the lean premixed flame side, a higher concentration of O2 leaks through to the diffusion flame. The general structure of the triple point features the contribution of both differential diffusion of radicals and heat. A mixture fraction-temperature phase plane description of the triple flame structure is proposed to highlight some interesting features in partially premixed combustion. The effects of differential diffusion at the triple point add to the contribution of hydrodynamic effects in the stabilization of the triple flame. Differential diffusion effects are measured using two methods: a direct computation using diffusion velocities and an indirect computation based on the difference between the normalized mixture fractions of C and H. The mixture fraction approach does not clearly identify the effects of differential diffusion, in particular at the curved triple point, because of ambiguities in the contribution of carbon and hydrogen atoms' carrying species.

Echekki, T.↗
DOE OSTI2023

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas – subsurface transport, turbulence, and anomalous materials – in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Examination of Dendritic Growth During Solidification of Ternary Alloys via a Novel Quantitative 3D Cellular Automaton Model

A three-dimensional (3D) quantitative cellular automaton (CA) model was developed to simulate dendritic growth during solidification processing of ternary alloys. A detailed method was proposed to solve solute diffusion and calculate the solid fraction for ternary alloys during solidification. The present model has been shown to accurately predict dendrite morphologies and solute distributions of both single equiaxed dendrite and columnar dendrites during solidification. The model was also used to study the influence of the concentration of a third component (Mg in the Al-Si-Mg system) on dendritic growth. With increasing Mg concentration, the steady-state dendrite tip growth velocity was shown to decrease, resulting in a shorter primary dendrite length. The 3D CA simulation results agree well with the prediction of the LGK theoretical model. Multi-columnar dendrite growth was simulated with different cooling rates. Primary and secondary dendrite morphology was shown to have little variation at low cooling rates (1.0 to 2.0 K/s) and a constant undercooling. The average secondary dendrite arm spacing was shown to decrease with the increase of cooling rate within this range. The 3D CA simulation results are in good agreement with the experimental directional solidification experiments of a commercial A356 (Al-7 wt pct Si-0.3 wt pct Mg).

36 MATERIALS SCIENCE↗
DOE OSTI2026

Refined Sensitivity Estimates for Single-Molecule Magnet Dark Matter Detectors

We revisit the sensitivity of Single Molecule Magnet (SMM) crystals as detectors for low-mass dark matter. In previous work, we established the concept of the ``magnetic bubble chamber'', where energy deposited by dark matter triggers a magnetic avalanche in a metastable crystal. The original sensitivity estimates relied on a conservative criterion requiring the spin relaxation time to be strictly shorter than the thermal diffusion time. Here, we demonstrate that this criterion effectively ignores the stochastic nature of spin relaxation. We derive a refined analytic estimate which accounts for the fraction of spins that relax even when diffusion is fast. We show that the Zeeman energy released by this fraction contributes to local heating, significantly lowering the energy threshold for avalanche formation. We present simulation results confirming this effect and report on experimental verification of the assumed low-temperature thermal properties of two representative SMM crystals, Mn$_{12}$-acetate and Mn$_{32}$. Together, these efforts extend this pathfinder program toward the realization of SMM-based detectors with controlled material properties and enhanced dark matter sensitivity.

Eberhardt, Andrew [Tokyo U., IPMU]↗
NASA NTRS1987

Oxygen-18 tracer study of the passive thermal oxidation of silicon

This work focuses on the thermal oxidation of silicon near 1273 K using the double-tracer oxidation method. The results confirm that oxidation occurs by the transport of electrically neutral nonnetwork oxygen through the interstitial space of the vitreous silica (v-SiO2) scale. Simultaneously, self- (or isotopic-) diffusion occurs in the network, resulting in characteristic isotopic fraction distributions near the gas-scale interface. The self-diffusion coefficients calculated from these profiles agree with those reported for tracer diffusion in v-SiO2, and the diffusion coefficient calculated from the scale growth is consistent with reported O2 permeation data. An important parameter that describes the double-oxidation behavior is a ratio relating the scale thickness grown during the second oxidation, the network self-diffusion coefficient for oxygen, and the time of the second oxidation.

Cawley, J. D.↗