Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion models”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Perpendicular Diffusion of Solar Energetic Particles: Model Results and Implications for Electrons

The processes responsible for the effective longitudinal transport of solar energetic particles (SEPs) are still not completely understood. We address this issue by simulating SEP electron propagation using a spatially 2D transport model that includes perpendicular diffusion. By implementing, as far as possible, the most reasonable estimates of the transport (diffusion) coefficients, we compare our results, in a qualitative manner, to recent observations at energies of 55–105 keV, focusing on the longitudinal distribution of the peak intensity, the maximum anisotropy, and the onset time. By using transport coefficients that are derived from first principles, we limit the number of free parameters in the model to (i) the probability of SEPs following diffusing magnetic field lines, quantified by a∈[0,1], and (ii) the broadness of the Gaussian injection function. It is found that the model solutions are extremely sensitive to the magnitude of the perpendicular diffusion coefficient and relatively insensitive to the form of the injection function as long as a reasonable value of a = 0.2 is used. We illustrate the effects of perpendicular diffusion on the model solutions and discuss the viability of this process as a dominant mechanism by which SEPs are transported in longitude. Lastly, we try to quantity the effectiveness of perpendicular diffusion as an interplay between the magnitude of the relevant diffusion coefficient and the SEP intensity gradient driving the diffusion process. It follows that perpendicular diffusion is extremely effective early in an SEP event when large intensity gradients are present, while the effectiveness quickly decreases with time thereafter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A semi-analytical method for simulating matrix diffusion in numerical transport models

We report a semi-analytical approximation for transient matrix diffusion is developed for use in numerical contaminant transport simulators. This method is an adaptation and extension of the heat conduction method of Vinsome and Westerveld (1980) used to simulate heat losses during thermally enhanced oil recovery. The semi-analytical method is used in place of discretization of the low permeability materials, and it represents the concentration profile in the low permeability materials with a fitting function that is adjusted in each element at each time-step. The resulting matrix diffusion fluxes are added to the numerical model as linear concentration-dependent source/sink terms. Since only the high permeability zones need to be discretized, the numerical formulation is extremely efficient compared to traditional approaches that require discretization of both the high and low permeability zones. The semi-analytical method compares favorably with the analytical solution for transient one-dimensional diffusion with first order decay, with a two-layer aquifer/aquitard solution, with the solution for transport in a fracture with matrix diffusion and decay, and with a fully numerical solution for transport in a thin sand zone bounded by clay with variable decay rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Diffuse interface approach to modeling crystal plasticity with accommodation of grain boundary sliding

In diffuse-interface (phase field) models of polycrystalline solids, the grain boundaries naturally possess a narrow region of finite thickness. A Ginzburg-Landau type kinetic equation for plastic deformation of a solid is derived that combines crystal plasticity and J 2 plasticity. Within such a framework for a polycrystal the grain interior can be defined by crystal plasticity, and the grain boundary region can be assigned J 2 type plasticity. The two are connected smoothly. This enables the modeling framework to accommodate grain boundary sliding (GBS), an important deformation mechanism for creep. The relaxation of elastic modulus of a polycrystalline solid and the power-law creep (both accommodated by GBS) are reported on by 2D and 3D simulations, respectively, and the results are compared to the 2D finite-element simulations of Ghahremani and Crossman-Ashby. A strong grain shape dependence and orientation dependence (for non-equiaxed grains) of the GBS effect are predicted by this model.

36 MATERIALS SCIENCE↗
DOE OSTI2015

GRIZZLY Model of Multi-Reactive Species Diffusion, Moisture/Heat Transfer and Alkali-Silica Reaction for Simulating Concrete Aging and Degradation

Concrete is widely used in the construction of nuclear facilities because of its structural strength and its ability to shield radiation. The use of concrete in nuclear power plants for containment and shielding of radiation and radioactive materials has made its performance crucial for the safe operation of the facility. As such, when life extension is considered for nuclear power plants, it is critical to have accurate and reliable predictive tools to address concerns related to various aging processes of concrete structures and the capacity of structures subjected to age-related degradation. The goal of this report is to document the progress of the development and implementation of a fully coupled thermo-hydro-mechanical-chemical model in GRIZZLY code with the ultimate goal to reliably simulate and predict long-term performance and response of aged NPP concrete structures subjected to a number of aging mechanisms including external chemical attacks and volume-changing chemical reactions within concrete structures induced by alkali-silica reactions and long-term exposure to irradiation. Based on a number of survey reports of concrete aging mechanisms relevant to nuclear power plants and recommendations from researchers in concrete community, we’ve implemented three modules during FY15 in GRIZZLY code, (1) multi-species reactive diffusion model within cement materials; (2) coupled moisture and heat transfer model in concrete; and (3) anisotropic, stress-dependent, alkali-silica reaction induced swelling model. The multi-species reactive diffusion model was implemented with the objective to model aging of concrete structures subjected to aggressive external chemical attacks (e.g., chloride attack, sulfate attack, etc.). It considers multiple processes relevant to external chemical attacks such as diffusion of ions in aqueous phase within pore spaces, equilibrium chemical speciation reactions and kinetic mineral dissolution/precipitation. The moisture/heat transfer module was implemented to simulate long-term spatial and temporal evolutions of the moisture and temperature fields within concrete structures at both room and elevated temperatures. The ASR swelling model implemented in GRIZZLY code can simulate anisotropic expansions of ASR gel under either uniaxial, biaxial and triaxial stress states, and can be run simultaneously with the moisture/heat transfer model and coupled with various elastic/inelastic solid mechanics models that were implemented in GRIZZLY code previously. This report provides detailed descriptions of the governing equations, constitutive equations and numerical algorithms of the three modules implemented in GRIZZLY during FY15, simulation results of example problems and model validation results by comparing simulations with available experimental data reported in the literature. The close match between the experiments and simulations clearly demonstrate the potential of GRIZZLY code for reliable evaluation and prediction of long-term performance and response of aged concrete structures in nuclear power plants.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1980

The longitudinal galactic cosmic ray intensity modulation in a diffusive and a scatter-free model of the inner heliosphere

Measurements of energetic solar flare electrons, Jovian electrons, and low-energy solar protons have led to the suggestion that the inner heliosphere is essentially scatter free. Consideration is given to a model in which scattering is only important within the corotating interaction regions (CIR's), while the cavity formed by the CIR's in the inner heliosphere is scatter free, as compared to an interplanetary diffusion model wherein scattering is determined by the solar wind and magnetic field microstructure everywhere in the heliosphere. The two models are compared in terms of the solar wind stream associated Forbush decreases at 1 AU and the longitudinal cosmic ray nucleonic intensity variations at larger heliocentric distances (Pioneer 11 measurements at about 3 AU). It is shown that the diffusion model is able to explain consistently all longitudinal galactic cosmic ray intensity modulation measurements.

Morfill, G.↗
DOE OSTI2018

Thermal resistance of pressed contacts of aluminum and niobium at liquid helium temperatures

Here, we examine the resistance to heat flow across contacts of mechanically pressed aluminum and niobium near liquid helium temperatures for designing a thermally conducting joint of aluminum and superconducting niobium. Measurements in the temperature range of 3.5 K to 5.5 K show the thermal contact resistance to grow as a near-cubic function of decreasing temperature, indicating phonons to be the primary heat carriers across the interface. In the 4 kN to 14 kN range of pressing force the contact resistance shows linear drop with the increasing force, in agreement with the model of micro-asperity plastic deformation at pressed contacts. Several thermal contact resistance models as well as the phonon diffuse mismatch model of interface thermal resistance are compared with the experimental data. The diffuse mismatch model shows closest agreement. The joints are further augmented with thin foil of indium, which lowers the joint resistance by an order of magnitude. The developed joint has nearly 1 K*cm2/W of thermal resistance at 4.2 K, is demountable, and free of the thermally resistive interfacial alloy layer that typically exists at welded, casted, or soldered joints of dissimilar metals.

43 PARTICLE ACCELERATORS↗
NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

This paper will include a detailed comparison of heat transfer models that rely upon the thermal diffusivity. The goals are to inform users of the development history of the various models and the resulting differences in model formulations, as well as to evaluate the models on a variety of validation cases so that users might better understand which models are more broadly applicable.

turbulence models↗
NASA NTRS1994

Nonadiabatic nonradial p-mode frequencies of the standard solar model, with and without helium diffusion

The nonadiabatic frequencies of a standard solar model and a solar model that includes helium diffusion are discussed. The nonadiabatic pulsation calculation includes physics that describes the losses and gains due to radiation. Radiative gains and losses are modeled in both the diffusion approximation, which is only valid in optically thick regions, and the Eddington approximation, which is valid in both optically thin and thick regions. The calculated pulsation frequencies for modes with l less than or equal to 1320 are compared to the observed spectrum of the Sun. Compared to a strictly adiabatic calculation, the nonadiabatic calculation of p-mode frequencies improves the agreement between model and observation. When helium diffusion is included in the model the frequencies of the modes that are sensitive to regions near the base of the convection zone are improved (i.e., brought into closer agreement with observation), but the agreement is made worse for other modes. Cyclic variations in the frequency spacings of the Sun as a function of frequency of n are presented as evidence for a discontinuity in the structure of the Sun, possibly located near the base of the convection zone.

Guenther, D. B.↗
NASA NTRS1973

Solute diffusion in liquid metals

A gas model of diffusion in liquid metals is presented. In this model, ions of liquid metals are assumed to behave like the molecules in a dense gas. Diffusion coefficient of solute is discussed with reference to its mass, ionic size, and pair potential. The model is applied to the case of solute diffusion in liquid silver. An attempt was made to predict diffusion coefficients of solutes with reasonable accuracy.

Bhat, B. N.↗
NASA NTRS1944

An Experimental Investigation of the Flow of Air in a Flat Broadening Channel

The wide use of diffusers, in various fields of technology, has resulted in several experimental projects to study the action and design of diffusers. Most of the projects dealt with steam (steam turbine nozzles). But diffusers have other applications - that is, ventilators, smoke ducts, air coolers, refrigeration, drying, and so forth. At present there is another application for diffusers in wind-tunnel design. Because of higher requirements and increased power of such installations more attention must be paid to the correctness of work and the decrease in losses due to every section of the tunnel. A diffuser, being one of the component parts of a tunnel , can in the event of faulty construction introduce considerable losses. Therefore, in the design of the new CAHI wind tunnel, it was suggested that an experimental study of diffusers be made, with a view to applying the results to wind tunnels. The experiments conducted by K. K. Baulin in the laboratories of CAHI upon models of diffusers of different cross sections, lengths, and angles of divergence, were a valuable source of experimental data. They were of no help, however, in reaching any conclusion regarding the optimum shape because of the complexity and diversity of the factors which all appeared simultaneously, thereby precluding the.study of the effects of any one factor separately. On the suggestion of the director of the CAHI,Prof. B. N. Ureff, it was decided to experiment on a two-dimensional diffuser model and determine the effect, of the angle of divergence. The author is acquainted with two experimental projects of like nature: the first was conducted with water, the other with air. The first of these works, although containing a wealth of experimental data, does not indicate the nature of flow or its relation to the angle of divergence. The second work is limited to four angles - that is, 12 deg, 24 deg, 45 deg, 90 deg. The study of this diffuser did not supply any information about the effect of smaller angles which, because of their advantages, are more commonly used, The author was able to acquaint himself with the second work only after the experiments were started. For these reasons, as well as because on the basis of those works no conclusion can be reached regarding the nature of flow distribution, of eddies, and so forth, experimental work was continued. The need for determining flow patterns follows from the fact that from them are determined methods of measurement - that is, the determination of velocities by means of the pitot tube, which, as is well known, gives correct indications only when placed with its axis parallel to the axis of flow. The data contained in this report were obtained from experiments conducted by the Aerodynamical Laboratories of the CAHI. The solutions to some. of the mathematical problems connected with the experiments are due to Prof. S. A, Chapligin.

Vedernikoff, A. N.↗
DOE OSTI2016

Investigating the Effects of Anisotropic Mass Transport on Dendrite Growth in High Energy Density Lithium Batteries

Dendrite formation on the electrode surface of high energy density lithium (Li) batteries causes safety problems and limits their applications. Suppressing dendrite growth could significantly improve Li battery performance. Dendrite growth and morphology is a function of the mixing in the electrolyte near the anode interface. Most research into dendrites in batteries focuses on dendrite formation in isotropic electrolytes (i.e., electrolytes with isotropic diffusion coefficient). In this work, an anisotropic diffusion reaction model is developed to study the anisotropic mixing effect on dendrite growth in Li batteries. The model uses a Lagrangian particle-based method to model dendrite growth in an anisotropic electrolyte solution. The model is verified by comparing the numerical simulation results with analytical solutions, and its accuracy is shown to be better than previous particle-based anisotropic diffusion models. Several parametric studies of dendrite growth in an anisotropic electrolyte are performed and the results demonstrate the effects of anisotropic transport on dendrite growth and morphology, and show the possible advantages of anisotropic electrolytes for dendrite suppression.

SPH, high energy density lithium batteries, precip↗
DOE OSTI2015

Albedo and Diffuse POA Measurements to Evaluate Transposition Model Uncertainty

Albedo and diffuse plane of array (DPOA) measurements are used in addition to standard global horizontal irradiance (GHI), direct normal irradiance (DNI), diffuse horizontal irradiance (DNI), and plane of array irradiance (POA) measurements to determine the impact of albedo on transposition model performance. Albedo measurements averaged 0.214. Daily albedo values ranged from 0.148 to 0.236 and were found to be correlated to daily relative humidity. DPOA measurements were compared to calculated DPOA values (from POA and DNI), and helped identify a suspected deviation from due south in the azimuth of the POA measurement. Since the measured albedo average was close to the typical fixed $^{albedo = 0.2}$ assumption, little difference was seen between using measured and fixed albedo (~0.15% differences in mean bias difference (MBD) and root mean squared difference (RMSD)). However, evaluation of transposition models at other fixed albedos showed an albedo change of 0.1 is found to increase total modeled insolation by approximately 1%. Thus, for locations with different ground surfaces (e.g., persistent snow cover of black surfaces), the impact of using measured albedo instead of the fixed $^{albedo = 0.2}$ assumption may be greater. Measurement deviations resulted in up to 2% changes in MBD and RIVISD when switching between interrelated measurements (e.g., GHI and DHI as inputs to transposition models versus DNI and DHI as inputs). Variation among transposition models was also up to 2% MBD and RMSD. Thus, for this data set, measurement deviation and transposition model selection are found to have more impact than using measured albedo instead fixed albedo.

42 ENGINEERING↗
DOE OSTI2019

Reshaping and sintering of 3D fcc metal nanoclusters: Stochastic atomistic modeling with realistic surface diffusion kinetics

Far-from-equilibrium evolution of metallic nanocluster shapes is highly sensitive to the atomistic-level details of surface diffusion for diverse local surface configurations. A stochastic model was developed incorporating realistic values for the multiple diffusion barriers (contrasting previous unrealistic generic prescriptions) based upon insights from homoepitaxial film growth. Kinetic Monte Carlo simulation then elucidates the conversion of Ag nanocubes to Wulff polyhedra mediated by nucleation of new {100} facets, the pinch-off of sufficiently elongated Ag nanorods, and key aspects of sintering for orientationally aligned Ag and Au nanoclusters. In conclusion, the time scale for sintering of Au nanoclusters observed in high-resolution transmission electron microscopy studies was also recovered.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction

Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.

Wang, Ziyan [Energy Geosciences Division Lawrence ↗
DOE OSTI2025

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗
DOE OSTI2022

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗