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 811 records · Page 45

Laser-induced generation of surface periodic structures in media with nonlinear diffusion

A model of fast formation of high-contrast periodic structure appearing on a semiconductor surface under action of laser radiation is proposed. The process of growing a surface structure due to the interaction surface plasmon- polaritons excited on nonequilibrium electrons with incident laser radiation are considered in the framework of a medium with nonlinear diffusion of nonequilibrium carriers (defects). A resonance effect of superfast pico- and subpicosecond amplification of the plasmon-polariton structure generated on the surface, the realization of which can result in a high-contrast defect lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Model-based analysis of the impact of diffuse radiation on CO 2 exchange in a temperate deciduous forest

Clouds and aerosols increase the fraction of global solar irradiance that is diffuse light. This phenomenon is known to increase the photosynthetic light use efficiency (LUE) of closed-canopy vegetation by redistributing photosynthetic photon flux density (400–700 nm) from saturated, sunlit leaves at the top of the canopy, to shaded leaves deeper in the canopy. We combined a process-based carbon cycle model with 10 years of eddy covariance carbon flux measurements and other ancillary data sets to assess 1) how this LUE enhancement influences interannual variation in carbon uptake, and 2) how errors in modeling diffuse fraction affect predictions of carbon uptake. Modeled annual gross primary productivity (GPP) increased by ≈0.94% when observed levels of diffuse fraction were increased by 0.01 (holding total irradiance constant). The sensitivity of GPP to increases in diffuse fraction was highest when the diffuse fraction was low to begin with, and lowest when the diffuse fraction was already high. Diffuse fraction also explained significantly more of the interannual variability of modeled net ecosystem exchange (NEE), than did total irradiance. Two tested radiation partitioning models yielded over- and underestimates of diffuse fraction at our site, which propagated to over- and underestimates of annual NEE, respectively. Our findings highlight the importance of incorporating LUE enhancement under diffuse light into models of global primary production, and improving models of diffuse fraction.

54 ENVIRONMENTAL SCIENCES↗

Modeling experimental results of diffusion of alkaline solutions through a compacted bentonite barrier

The interaction between concrete/cement and swelling clay (bentonite) has been modeled in the context of engineered barrier systems for deep geological disposal of high-level radioactive waste. The geochemical transformations observed in laboratory diffusion experiments at 60 and 90 {sup o}C between bentonite and different high-pH solutions (K-Na-OH and Ca(OH){sub 2}-saturated) were reconciled with the reactive transport code CrunchFlow. For K-Na-OH solutions (pH = 13.5 at 25 {sup o}C) partial dissolution of montmorillonite and precipitation of Mg-silicates (talc-like), hydrotalcite and brucite at the interface are predicted at 60 {sup o}C, while at 90 {sup o}C the alteration is wider. Alkaline cations diffused beyond the mineralogical alteration zone by means of exchange with Mg{sup 2+} in the interlayer region of montmorillonite. Very slow reactivity and minor alteration of the clay are predicted in the Ca(OH){sub 2}-bentonite system. The model is a reasonable description of the experiments but also demonstrates the difficulties in modeling processes operating at a small scale under a diffusive regime.

36 MATERIALS SCIENCE↗

Exits in order: How crowding affects particle lifetimes

Diffusive processes are often represented using stochastic random walk frameworks. The amount of time taken for an individual in a random walk to intersect with an absorbing boundary is a fundamental property that is often referred to as the particle lifetime, or the first passage time. The mean lifetime of particles in a random walk model of diffusion is related to the amount of time required for the diffusive process to reach a steady state. Mathematical analysis describing the mean lifetime of particles in a standard model of diffusion without crowding is well known. However, the lifetime of agents in a random walk with crowding has received much less attention. Since many applications of diffusion in biology and biophysics include crowding effects, here we study a discrete model of diffusion that incorporates crowding. Using simulations, we show that crowding has a dramatic effect on agent lifetimes, and we derive an approximate expression for the mean agent lifetime that includes crowding effects. Our expression matches simulation results very well, and highlights the importance of crowding effects that are sometimes overlooked.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗

Fuel Thermo-physical Characterization Project: Evaluation of Models to Calculate Thermal Diffusivity of Layered Composites

The Office of Material Management and Minimization Fuel Thermo-physical Characterization Project at Pacific Northwest National Laboratory (PNNL) is tasked with using PNNL facilities and processes to receive irradiated low enriched uranium-molybdenum fuel plate samples and perform analyses in support of the Office of Material Management and Minimization Reactor Conversion Program. This work is in support of the Fuel Development Pillar that is managed by Idaho National Laboratory. A key portion of the scope associated with this project was to measure the thermal properties of fuel segments harvested from plates that were irradiated in the Advanced Test Reactor. Thermal diffusivity of samples prepared from the fuel segments was measured using laser flash analysis. Two models, one developed by PNNL and the other developed by the Technische Universität München (TUM), were evaluated to extract the thermal diffusivity of the uranium-molybdenum alloy from measurements made on the irradiated, layered composites. The experimental data of the “TC” irradiated fuel segment was evaluated using both models considering a three-layer and five-layer system. Both models are in acceptable agreement with one another and indicate that the zirconium diffusion barrier has a minimal impact on the overall thermal diffusivity of the monolithic U-Mo fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

In Situ Observation of Calcium Aluminate Inclusions Dissolution into Steelmaking Slag

The dissolution rate of calcium aluminate inclusions in CaO-SiO{sub 2}-Al{sub 2}O{sub 3} slags has been studied using confocal scanning laser microscopy (CSLM) at elevated temperatures: 1773 K, 1823 K, and 1873 K (1500 °C, 1550 °C, and 1600 °C). The inclusion particles used in this experimental work were produced in our laboratory and their production technique is explained in detail. Even though the particles had irregular shapes, there was no rotation observed. Further, the total dissolution time decreased with increasing temperature and decreasing SiO{sub 2} content in the slag. The rate limiting steps are discussed in terms of shrinking core models and diffusion into a stagnant fluid model. It is shown that the rate limiting step for dissolution is mass transfer in the slag at 1823 K and 1873 K (1550 °C and 1600 °C). Further investigations are required to determine the dissolution mechanism at 1773 K (1500 °C). The calculated diffusion coefficients were inversely proportional to the slag viscosity and the obtained values for the systems studied ranged between 5.64 × 10{sup −12} and 5.8 × 10{sup −10} m{sup 2}/s.

36 MATERIALS SCIENCE↗

Transport of hydrogen in metals with occupancy dependent trap energies

Common diffusion trapping models for modeling hydrogen transport in metals are limited to traps with single de-trapping energies and a saturation occupancy of one. While they are successful in predicting typical mono isotopic ion implantation and thermal degassing experiments, they fail at describing recent experiments on isotope exchange at low temperatures. This paper presents a new modified diffusion trapping model with fill level dependent de-trapping energies that can also explain these new isotope exchange experiments. Density function theory (DFT) calculations predict that even mono vacancies can store between 6 and 12 H atoms with de-trapping energies that depend on the fill level of the mono vacancy. The new fill level dependent diffusion trapping model allows to test these DFT results by bridging the gap in length and time scale between DFT calculations and experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Realistic galaxy image simulation via score-based generative models

ABSTRACT We show that a denoising diffusion probabilistic model (DDPM), a class of score-based generative model, can be used to produce realistic mock images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument (DESI) grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real data set. We quantify the similarity by borrowing from the deep generative learning literature, using the ‘Fréchet inception distance’ to test for subjective and morphological similarity. We also introduce the ‘synthetic galaxy distance’ metric to compare the emergent physical properties (such as total magnitude, colour, and half-light radius) of a ground truth parent and synthesized child data set. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as adversarial networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate inpainting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we ‘DESI-fy’ cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-Linear Langevin and Fractional Fokker–Planck Equations for Anomalous Diffusion by Lévy Stable Processes

The numerical solutions to a non-linear Fractional Fokker–Planck (FFP) equation are studied estimating the generalized diffusion coefficients. The aim is to model anomalous diffusion using an FFP description with fractional velocity derivatives and Langevin dynamics where Lévy fluctuations are introduced to model the effect of non-local transport due to fractional diffusion in velocity space. Distribution functions are found using numerical means for varying degrees of fractionality of the stable Lévy distribution as solutions to the FFP equation. The statistical properties of the distribution functions are assessed by a generalized normalized expectation measure and entropy and modified transport coefficient. The transport coefficient significantly increases with decreasing fractality which is corroborated by analysis of experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling convective core overshoot and diffusion in Procyon constrained by asteroseismic data

We compare evolved stellar models, which match Procyon's mass and position in the HR diagram, to current ground-based asteroseismic observations. Diffusion of helium and metals along with two conventional core overshoot descriptions and the Kuhfuss nonlocal theory of convection are considered. We establish that one of the two published asteroseismic data reductions for Procyon, which mainly differ in their identification of even versus odd l values, is a significantly more probable and self-consistent match to our models than the other. The most probable models according to our Bayesian analysis have evolved to just short of turnoff, still retaining a hydrogen convective core. Our most probable models include Y and Z diffusion and have conventional core overshoot between 0.9 and 1.5 pressure scale heights, which increases the outer radius of the convective core by between 36% and 43%, respectively. We discuss the significance of this comparatively higher than expected core overshoot amount in terms of internal mixing during evolution. The parameters of our most probable models are similar regardless of whether adiabatic or nonadiabatic model p-mode frequencies are compared to the observations, although, the Bayesian probabilities are greater when the nonadiabatic model frequencies are used. All the most probable models (with or without core overshoot, adiabatic or nonadiabatic model frequencies, diffusion or no diffusion, including priors for the observed HRD location and mass or not) have masses that are within 1σ of the observed mass 1.497 ± 0.037 M {sub ☉}.

79 ASTRONOMY AND ASTROPHYSICS↗

Multilayer interface tracking model of pure tungsten oxidation

Here, we present a numerical model to predict oxide scale growth on tungsten surfaces under exposure to oxygen at high temperatures. The model captures the formation of four thermodynamically-compatible oxide sublayers, WO 2 , WO 2.72 , WO 2.9 , and WO 3 , on top of the metal substrate. Oxide layer growth is simulated by tracking the oxide/oxide and oxide/metal interfaces using a sharp-interface Stefan model coupled to diffusion kinetics. The model is parameterized using selected experimental measurements and electronic structure calculations of the diffusivities of all the oxide subphases involved. We simulate oxide growth at temperatures of 600°C and above, extracting the power law growth exponents in each case, which we find to deviate from classical parabolic growth in several cases. We conduct a comparison of the model predictions with an extensive experimental data set, with reasonable agreement at most temperatures. While many gaps in our understanding still exist, this work is a first attempt at embedding the thermodynamic and kinetic complexity of tungsten oxide growth into a comprehensive mesoscale kinetic model that attempts to capture the essential features of tungsten oxidation to fill existing knowledge gaps and guide and enhance future tungsten oxidation models.

36 MATERIALS SCIENCE↗

Numerical simulation of a thermodynamically consistent four‐species tumor growth model

We develop a thermodynamically consistent four‐species model of tumor growth on the basis of the continuum theory of mixtures. Unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction‐diffusion equation. The formulation involves systems of highly nonlinear partial differential equations of surface effects through diffuse‐interface models. A mixed finite element spatial discretization is developed and implemented to provide numerical results demonstrating the range of solutions this model can produce. A time‐stepping algorithm is then presented for this system, which is shown to be first order accurate and energy gradient stable. The results of an array of numerical experiments are presented, which demonstrate a wide range of solutions produced by various choices of model parameters.

97 MATHEMATICS AND COMPUTING↗

An adaptive tau-leaping method for stochastic simulations of reaction-diffusion systems

Stochastic modelling is critical for studying many biochemical processes in a cell, in particular when some reacting species have low population numbers. For many such cellular processes the spatial distribution of the molecular species plays a key role. The evolution of spatially heterogeneous biochemical systems with some species in low amounts is accurately described by the mesoscopic model of the Reaction-Diffusion Master Equation. The Inhomogeneous Stochastic Simulation Algorithm provides an exact strategy to numerically solve this model, but it is computationally very expensive on realistic applications. We propose a novel adaptive time-stepping scheme for the tau-leaping method for approximating the solution of the Reaction-Diffusion Master Equation. This technique combines effective strategies for variable time-stepping with path preservation to reduce the computational cost, while maintaining the desired accuracy. The numerical tests on various examples arising in applications show the improved efficiency achieved by the new adaptive method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Application of PDF mixing models to premixed flames with differential diffusion

Differential diffusion alters the balance of reaction and diffusion in turbulent premixed combustion, affecting the performance and emissions of combustion devices. Modelling combustion devices with Probability or Filtered Density Function (PDF or FDF) methods provides an exact treatment for the change in composition due to chemical reaction, while molecular mixing has to be modelled. Previous PDF molecular mixing models do not account for differential diffusion in a manner which satisfies realizability requirements. A new approach for treating differential diffusion, which ensures realizability, is proposed for pairwise-exchange mixing models in general, and applied in the Interaction by Exchange with the Mean (IEM) model of Dopazo [26], and in the Euclidean Minimum Spanning Tree (EMST) model of Subramaniam and Pope [5]. The new differential diffusion models are referred to as IEM-DD and EMST-DD respectively. Results from two and three-dimensional DNS of turbulent premixed methane–air combustion show that mixing rates and conditional statistics of species mass fractions depend on species diffusivities and the combustion regime. Zero-dimensional PDF model results obtained for the two-dimensional DNS case show that the EMST-DD model best reproduces the features that characterize differential diffusion in the DNS. The essential feature of the EMST-DD model, which accounts for its success in turbulent premixed combustion, is that differential mixing rates are imposed within a model which mixes locally in composition space.

Energy & Fuels↗

Diffusive flux in a model of stochastically gated oxygen transport in insect respiration

Oxygen delivery to insect tissues is controlled by transport through a branched tubular network that is connected to the atmosphere by valve-like gates, known as spiracles. In certain physiological regimes, the spiracles appear to be randomly switching between open and closed states. Quantitative analysis of this regime leads a reaction-diffusion problem with stochastically switching boundary condition. We derive an expression for the diffusive flux at long times in this problem. Our approach starts with the derivation of the passage probability for a single particle that diffuses between a stochastically gated boundary, which models the opening and closing spiracle, and the perfectly absorbing boundary, which models oxygen absorption by the tissue. This passage probability is then used to derive an expression giving the diffusive flux as a function of the geometric parameters of the tube and characteristic time scales of diffusion and gate dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modelling oxygen self-diffusion in UO 2 under pressure

Access to values for oxygen self-diffusion over a range of temperatures and pressures in UO 2 is important to nuclear fuel applications. Here, elastic and expansivity data are used in the framework of a thermodynamic model, the cBΩ model, to derive the oxygen self-diffusion coefficient in UO 2 over a range of pressures (0–10 GPa) and temperatures (300–1900 K). Furthermore, the significant reduction in oxygen self-diffusion as a function of increasing hydrostatic pressure, and the associated increase in activation energy, is identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗