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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Apparent permeability prediction of organic shale with generalized lattice Boltzmann model considering surface diffusion effect

Gas flow in shale is associated with both organic matter (OM) and inorganic matter (IOM) which contain nano-pores ranging in size from a few to hundreds of nano-meters. In addition to the non-continuum effect which leads to an apparent permeability of gas higher than the intrinsic permeability, the surface diffusion of adsorbed gas in organic pores also can influence the apparent permeability through its own transport mechanism. In this study, a generalized lattice Boltzmann model (GLBM) is employed for gas flow through the reconstructed shale matrix consisting of OM and IOM. The Expectation–Maximization (EM) algorithm is used to assign the pore size distribution to each component, and the dusty gas model (DGM) and generalized Maxwell–Stefan model (GMS) are adopted to calculate the apparent permeability accounting for multiple transport mechanisms including viscous flow, Knudsen diffusion and surface diffusion. Effects of pore radius and pressure on permeability of both IOM and OM as well as effects of Langmuir parameters on OM are investigated. The effect of total organic content and distribution on the apparent permeability of the reconstructed shale matrix at different surface diffusivity is also studied. It is found that the influence of pore size and pressure on the apparent permeability of organic matter is affected by the surface diffusion of adsorbed gas. Furthermore, surface diffusion plays a significant role in determining apparent permeability and the velocity distribution of shale matrix.

02 PETROLEUM↗
DOE OSTI2012

Diffusion Dominant Solute Transport Modelling in Fractured Media Under Deep Geological Environment - 12211

Deep geologic disposal of high activity and long-lived radioactive waste is gaining increasing support in many countries, where suitable low permeability geological formation in combination with engineered barriers are used to provide long term waste contaminant and minimise the impacts to the environment and risk to the biosphere. This modelling study examines the solute transport in fractured media under low flow velocities that are relevant to a deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes. The effects of water velocity in the fracture, matrix porosity and diffusion on solute transport are investigated and discussed. Some illustrative modelled results are presented to demonstrate the use of the model to examine the effects of media degradation on solute transport, under the influences of hydrogeological (diffusion dominant) and microbially mediated chemical processes. The challenges facing the prediction of long term degradation such as cracks evolution, interaction and coalescence are highlighted. The potential of a novel microstructure informed modelling approach to account for these effects is discussed, particularly with respect to investigating multiple phenomena impact on material performance. The GRM code is used to examine the effects of media degradation for a geological waste disposal package, under the combined hydrogeological (diffusion dominant) and chemical effects in low groundwater flow conditions that are typical of deep geological disposal systems. An illustrative reactive transport modelling application demonstrates the use of the code to examine the interplay of kinetic controlled biogeochemical reactive processes with advective and diffusive transport, under the influence of media degradation. The initial model results are encouraging which show the disposal system to evolve in a physically realistic manner. In the example presented the reactive-transport coupling develops chemically reducing zones, which limit the transport of uranium. This illustrates the potential significance of media degradation and chemical effect on the transport of radionuclides which would need to be taken into account when examining the long-term behaviour and containment properties of the geological disposal system. Microstructure-informed modelling and its potential linkage with continuum flow modelling is a subject of ongoing studies. The approach of microstructure-informed modelling is discussed to provide insight and a mechanistic understanding of macroscopic parameters and their evolution. The proposed theoretical and methodological basis for microstructure-informed modelling of porous quasi-brittle media has the potential to develop into an explanatory and predictive tool for deriving mechanism-based, as opposed to phenomenological, evolution laws for macroscopic properties. These concepts in micro-scale modelling are likely to be applicable to the diffusion process, in addition to advective transport illustrated here for porous media. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Modelling of monovacancy diffusion in W over wide temperature range

The diffusion of monovacancies in tungsten is studied computationally over a wide temperature range from 1300 K until the melting point of the material. Our modelling is based on Molecular Dynamics technique and Density Functional Theory. The monovacancy migration barriers are calculated using nudged elastic band method for nearest and next-nearest neighbour monovacancy jumps. The diffusion pre-exponential factor for monovacancy diffusion is found to be two to three orders of magnitude higher than commonly used in computational studies, resulting in attempt frequency of the order 10{sup 15} Hz. Multiple nearest neighbour jumps of monovacancy are found to play an important role in the contribution to the total diffusion coefficient, especially at temperatures above 2/3 of T{sub m}, resulting in an upward curvature of the Arrhenius diagram. The probabilities for different nearest neighbour jumps for monovacancy in W are calculated at different temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Modeling and Simulating Diffused Aeration for Hydrogen Removal from Expansion Tanks of Parabolic Trough Solar Thermal Power Plants

High-temperature decomposition of organic heat transfer fluids (HTFs) can cause hydrogen build-up in the receiver annulus of parabolic trough power plants. This build-up increases the receiver thermal losses and results in a decline in power output. Prior work has shown that removal of hydrogen from the expansion tank can be an effective mitigation strategy. This paper presents a modeling methodology and simulation results for passive and diffused aeration removal of hydrogen from the expansion tank. Focus is on estimating the mass transfer of hydrogen from the HTF into the headspace gas across the surface and bubble liquid/gas interfaces. Simulations for the operation conditions and expansion tank geometry of the Nevada Solar One power plant (located near Las Vegas, Nevada, USA) show that passive removal (surface mass transfer only) may be sufficient to obtain hydrogen removal rates of 1.7 x 10-4 mol/s (minimal removal rate necessary). A diffused aeration system (surface plus bubble mass transfer) with gas injection volume rates up to 100 L/s may be necessary to obtain a hydrogen removal rate of 2.0 x 10-4 mol/s (optimal removal rate).

14 SOLAR ENERGY↗
DOE OSTI2014

Interface scattering in polycrystalline thermoelectrics

We study the effect of electron and phonon interface scattering on the thermoelectric properties of disordered, polycrystalline materials (with grain sizes larger than electron and phonons' mean free path). Interface scattering of electrons is treated with a Landauer approach, while that of phonons is treated with the diffuse mismatch model. The interface scattering is embedded within a diffusive model of bulk transport, and we show that, for randomly arranged interfaces, the overall system is well described by effective medium theory. Using bulk parameters similar to those of PbTe and a square barrier potential for the interface electron scattering, we identify the interface scattering parameters for which the figure of merit ZT is increased. We find the electronic scattering is generally detrimental due to a reduction in electrical conductivity; however, for sufficiently weak electronic interface scattering, ZT is enhanced due to phonon interface scattering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum

Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2013

Investigation of Fission Product Transport into Zeolite-A for Pyroprocessing Waste Minimization

Methods to improve fission product salt sorption into zeolite-A have been investigated in an effort to reduce waste associated with the electrochemical treatment of spent nuclear fuel. It was demonstrated that individual fission product chloride salts were absorbed by zeolite-A in a solid-state process. As a result, recycling of LiCl-KCl appears feasible via adding a zone-freezing technique to the current treatment process. Ternary salt molten-state experiments showed the limiting kinetics of CsCl and SrCl2 sorption into the zeolite. CsCl sorption occurred rapidly relative to SrCl2 with no observed dependence on zeolite particle size, while SrCl2 sorption was highly dependent on particle size. The application of experimental data to a developed reaction-diffusion-based sorption model yielded diffusivities of 8.04 × 10-6 and 4.04 × 10-7 cm2 /s for CsCl and SrCl2, respectively. Additionally, the chemical reaction term in the developed model was found to be insignificant compared to the diffusion term.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Effect of Diffusion Length in Modeling of Equiaxed Dendritic Solidification under Buoyancy Flow in a Configuration of Hebditch–Hunt Experiment

Modeling of equiaxed solidification is vital for understanding the solidification process of metallic alloys. In this work, an extended literature review is given for the models currently used for equiaxed solidification simulations. Based on this analysis, we present a three-phase multiscale equiaxed solidification model in which some approximations regarding solute transport at microscopic scale are put together in a new way and incorporated into macroscopic transport equations. For the latter, a term relating to the momentum exchange between the two phases is revised, and a modification for the grain packing algorithm is proposed. A modernized model for equiaxed dendrite growth is tested using a case of solidification of Sn-5 wt pct Pb alloy in a parallelepiped cavity that mimics the Hebditch–Hunt experiment. The results obtained using two approaches to calculate diffusion length are presented and compared both with each other and with numerical results from elsewhere. It is demonstrated that diffusion length has a crucial effect on the final segregation pattern.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Quasilinear model for energetic particle diffusion in radial and velocity space

A quasilinear model for passive energetic particle (EP) turbulent diffusion in radial and velocity space is fitted and tested against nonlinear gyrokinetic tokamak simulations with the GYRO code [J. Candy and R. E. Waltz, Phys. Rev. Lett. 91, 045001 (2003)]. Off diagonal elements of a symmetric positive definite 2 Multiplication-Sign 2 EP diffusion matrix account for fluxes up radial (energy) gradients driven by energy (radial) gradients of the EP velocity space distribution function. The quasilinear ratio kernel of the model is provided by a simple analytic formula for the EP radial and velocity space EP diffusivity relative to radial thermal ion energy diffusivity at each linear mode of the turbulence driven by the thermal plasma. The TGLF [G. M. Staebler, J. E. Kinsey, and R. E. Waltz, Phys. Plasmas 14, 0055909 (2007); ibid. 15, 0055908 (2008)] tokamak transport model provides the linear mode frequency and growth rates to the kernel as well as the nonlinear spectral weight for each mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Compact Models for Defect Diffusivity in Semiconductor Alloys.

Predicting transient effects caused by short - pulse neutron irradiation of electronic devices is an important part of Sandia's mission. For example , predicting the diffusion of radiation - induced point defects is needed with in Sandia's Qualification Alternative to the Sandia Pulsed Reactor (QASPR) pro gram since defect diffusion mediates transient gain recovery in QASPR electronic devices. Recently, the semiconductors used to fabricate radiation - hard electronic devices have begun to shift from silicon to III - V compounds such as GaAs, InAs , GaP and InP . An advantage of this shift is that it allows engineers to optimize the radiation hardness of electronic devices by using alloy s such as InGaAs and InGaP . However, the computer codes currently being used to simulate transient radiation effects in QASP R devices will need to be modified since they presume that defect properties (charge states, energy levels, and diffusivities) in these alloys do not change with time. This is not realistic since the energy and properties of a defect depend on the types of atoms near it and , therefore, on its location in the alloy. In particular, radiation - induced defects are created at nearly random locations in an alloy and the distribution of their local environments - and thus their energies and properties - evolves with time as the defects diffuse through the alloy . To incorporate these consequential effects into computer codes used to simulate transient radiation effects, we have developed procedures to accurately compute the time dependence of defect energies and properties and then formulate them within compact models that can be employed in these computer codes. In this document, we demonstrate these procedures for the case of the highly mobile P interstitial (I P ) in an InGaP alloy. Further dissemination only as authorized to U.S. Government agencies and their contractors; other requests shall be approved by the originating facility or higher DOE programmatic authority.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Density-Functional-Theory Modeling of Cation Diffusion in Bulk La 1 - x Sr x MnO 3 ± δ ( x = 0.0 – 0.25 ) for Solid-Oxide Fuel-Cell Cathodes

In this work, the A - and B -site cation migration pathways involving defect complexes in bulk La 1-x Sr x MnO 3±δ (LSM) at x = 0.0-0.25 are investigated based on density-functional-theory modeling for solid-oxide fuel-cell (SOFC) cathode applications. We propose a dominant A -site cation migration mechanism which involves an A -site cation (e.g., La$x\atop{A}$) V A "' of a V A "' -V B "' cluster, where La$x\atop{A}$, V A "' and V B "' are La 3+ , A-site vacancy, and B-site vacancy in bulk LSM, respectively, and V A "' -V B "' is the first nearest-neighbor V A "' and V B "' pair. This hop exhibits an approximately 1.6-eV migration barrier as compared to approximately 2.9 eV of the La$x\atop{A}$ hop into a V A "'. This decrease in the cation migration barrier is attributed to the presence of the V B "' relieving the electrostatic repulsion and steric constraints to the migrating A-site cations in the transition-state image configurations.

03 NATURAL GAS↗
DOE OSTI2012

THE FERMI BUBBLES. II. THE POTENTIAL ROLES OF VISCOSITY AND COSMIC-RAY DIFFUSION IN JET MODELS

The origin of the Fermi bubbles recently detected by the Fermi Gamma-ray Space Telescope in the inner Galaxy is mysterious. In the companion paper Guo and Mathews (Paper I), we use hydrodynamic simulations to show that they could be produced by a recent powerful active galactic nucleus (AGN) jet event. Here, we further explore this scenario to study the potential roles of shear viscosity and cosmic-ray (CR) diffusion on the morphology and CR distribution of the bubbles. We show that even a relatively low level of viscosity ({mu}{sub visc} {approx}> 3 g cm{sup -1} s{sup -1}, or {approx}0.1%-1% of Braginskii viscosity in this context) could effectively suppress the development of Kelvin-Helmholtz instabilities at the bubble surface, resulting in smooth bubble edges as observed. Furthermore, viscosity reduces circulating motions within the bubbles, which would otherwise mix the CR-carrying jet backflow near bubble edges with the bubble interior. Thus viscosity naturally produces an edge-favored CR distribution, an important ingredient to produce the observed flat gamma-ray surface brightness distribution. Generically, such a CR distribution often produces a limb-brightened gamma-ray intensity distribution. However, we show that by incorporating CR diffusion that is strongly suppressed across the bubble surface (as inferred from sharp bubble edges) but is close to canonical values in the bubble interior, we obtain a reasonably flat gamma-ray intensity profile. The similarity of the resulting CR bubble with the observed Fermi bubbles strengthens our previous result in Paper I that the Fermi bubbles were produced by a recent AGN jet event. Studies of the nearby Fermi bubbles may provide a unique opportunity to study the potential roles of plasma viscosity and CR diffusion on the evolution of AGN jets and bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Developing a laser shockwave model for characterizing diffusion bonded interfaces

The US National Nuclear Security Agency has a Global Threat Reduction Initiative (GTRI) with the goal of reducing the worldwide use of high-enriched uranium (HEU). A salient component of that initiative is the conversion of research reactors from HEU to low enriched uranium (LEU) fuels. An innovative fuel is being developed to replace HEU in high-power research reactors. The new LEU fuel is a monolithic fuel made from a U-Mo alloy foil encapsulated in Al-6061 cladding. In order to support the fuel qualification process, the Laser Shockwave Technique (LST) is being developed to characterize the clad-clad and fuel-clad interface strengths in fresh and irradiated fuel plates. LST is a non-contact method that uses lasers for the generation and detection of large amplitude acoustic waves to characterize interfaces in nuclear fuel plates. However, because the deposition of laser energy into the containment layer on a specimen's surface is intractably complex, the shock wave energy is inferred from the surface velocity measured on the backside of the fuel plate and the depth of the impression left on the surface by the high pressure plasma pulse created by the shock laser. To help quantify the stresses generated at the interfaces, a finite element method (FEM) model is being utilized. This paper will report on initial efforts to develop and validate the model by comparing numerical and experimental results for back surface velocities and front surface depressions in a single aluminum plate representative of the fuel cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Atomistic modeling of intrinsic and radiation-enhanced fission gas (Xe) diffusion in UO2 +/- x: Implications for nuclear fuel performance modeling

Based on density functional theory (DFT) and empirical potential calculations, the diffusivity of fission gas atoms (Xe) in UO2 nuclear fuel has been calculated for a range of non-stoichiometry (i.e. UO2x), under both out-of-pile (no irradiation) and in-pile (irradiation) conditions. This was achieved by first deriving expressions for the activation energy that account for the type of trap site that the fission gas atoms occupy, which includes the corresponding type of mobile cluster, the charge state of these defects and the chemistry acting as boundary condition. In the next step DFT calculations were used to estimate migration barriers and internal energy contributions to the thermodynamic properties and calculations based on empirical potentials were used to estimate defect formation and migration entropies (i.e. pre-exponentials). The diffusivities calculated for out-of-pile conditions as function of the UO2x nonstoichiometrywere used to validate the accuracy of the diffusion models and the DFT calculations against available experimental data. The Xe diffusivity is predicted to depend strongly on the UO2x non-stoichiometry due to a combination of changes in the preferred Xe trap site and in the concentration of uranium vacancies enabling Xe diffusion, which is consistent with experiments. After establishing the validity of the modeling approach, it was used for studying Xe diffusion under in-pile conditions, for which experimental data is very scarce. The radiation-enhanced Xe diffusivity is compared to existing empirical models. Finally, the predicted fission gas diffusion rates were implemented in the BISON fuel performance code and fission gas release from a Risø fuel rod irradiation experiment was simulated. 2014 Elsevier B.V. All rights

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2013

Diffusion Dominant Solute Transport Modelling In Deep Repository Under The Effect of Emplacement Media Degradation - 13285

Deep geologic disposal of high activity and long-lived radioactive waste is being actively considered and pursued in many countries, where low permeability geological formations are used to provide long term waste contaminant with minimum impact to the environment and risk to the biosphere. A multi-barrier approach that makes use of both engineered and natural barriers (i.e. geological formations) is often used to further enhance the containment performance of the repository. As the deep repository system subjects to a variety of thermo-hydro-chemo-mechanical (THCM) effects over its long 'operational' lifespan (e.g. 0.1 to 1.0 million years, the integrity of the barrier system will decrease over time (e.g. fracturing in rock or clay)). This is broadly referred as media degradation in the present study. This modelling study examines the effects of media degradation on diffusion dominant solute transport in fractured media that are typical of deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes, while the effects of degradation is studied using a pore network model that considers the media diffusivity and network changes. Model results are presented to demonstrate the use of a 3D pore-network model, using a novel architecture, to calculate macroscopic properties of the medium such as diffusivity, subject to pore space changes as the media degrade. Results from a reactive transport model of a representative geological waste disposal package are also presented to demonstrate the effect of media property change on the solute migration behaviour, illustrating the complex interplay between kinetic biogeochemical processes and diffusion dominant transport. The initial modelling results demonstrate the feasibility of a coupled modelling approach (using pore-network model and reactive transport model) to examine the long term behaviour of deep geological repositories with media property change under complex geochemical conditions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Grassmannian diffusion maps based surrogate modeling via geometric harmonics

Abstract A novel surrogate model based on the Grassmannian diffusion maps (GDMaps) and utilizing geometric harmonics (GH) is developed for predicting the response of complex physical phenomena. The method utilizes GDMaps to obtain a low‐dimensional representation of the underlying behavior of physical/mathematical systems with respect to uncertain input parameters. Using this representation, GH, an out‐of‐sample extension technique, is employed to create a global map from the input parameter space to a Grassmannian diffusion manifold. GH is further employed to locally map points on the diffusion manifold onto the tangent space of a Grassmann manifold. The exponential map is then used to project the points in the tangent space onto the Grassmann manifold, where reconstruction of the full solution is performed. The performance of the proposed surrogate model is verified with three examples. The first problem is a toy example used to illustrate the technique. In the second example, errors associated with the various mappings are assessed by studying response predictions of the electric potential of a dielectric cylinder in a homogeneous electric field. The last example applies the method for uncertainty prediction in the strain field evolution in a model amorphous material using the shear transformation zone theory of plasticity.

42 ENGINEERING↗
DOE OSTI2023

Galactic Bar Resonances with Diffusion: An Analytic Model with Implications for Bar–Dark Matter Halo Dynamical Friction

Abstract The secular evolution of disk galaxies is largely driven by resonances between the orbits of “particles” (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari, and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies—effects that have, with few exceptions, been previously explored only with complex N -body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles’ slow actions. We solve this equation for various values of the dimensionless diffusion strength Δ, and then apply our theory to the calculation of bar–halo dynamical friction. For Δ = 0, we recover the classic result of Tremaine and Weinberg that friction ultimately vanishes, owing to the phase mixing of resonant orbits. However, for Δ > 0, we find that diffusion suppresses phase mixing, leading to a finite torque. Our results suggest that stochasticity—be it physical or numerical—tends to increase bar–halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time.

79 ASTRONOMY AND ASTROPHYSICS↗