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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Screening Diffusion of Small Molecules in Flexible Zeolitic Imidazolate Frameworks Using a DFT-Parameterized Force Field

Quantitatively modeling adsorbate diffusion through zeolitic imidazolate frameworks (ZIFs) must account for the inherent flexibility of these materials. Here, the lack of a transferable intramolecular ZIF force field (FF) for use in classical simulations has previously made an accurate simulation of adsorbate diffusion in many ZIFs impossible. We resolve this problem by introducing a density functional theory parameterized force field (FF) for ZIFs named the intraZIF-FF, which includes perturbations to the class I force fields previously used to model ZIFs. This FF outperforms ad hoc force fields at predicting ab initio relative energies and atomic forces taken from fully periodic ab initio molecular dynamics simulations of SALEM-2, ZIF-7, ZIF-8, and ZIF-90. We use the intraZIF-FF to predict the infinite dilution self-diffusion coefficients of 30 adsorbates with molecular diameters ranging from 2.66 to 7.0 Å in these 4 ZIFs. These results greatly increase the number of adsorbates for which accurate information about molecular diffusion in ZIFs is available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Dual-pace transient heat conduction in vertically aligned carbon nanotube arrays induced by structure separation

Vertically aligned arrays of carbon nanotubes (CNTs) are attractive for a wide range of macroscopic applications which can exploit the remarkable properties of individual nanotubes. In this work, an abnormal behavior of CNT bundles in heat conduction is discovered under the transient electro-thermal characterization. The measured voltage change over the sample shows a dual-pace thermal response (DTR), which could not be fitted using a single thermal diffusivity heat transfer model. Instead, two thermal diffusivities are determined from the DTR phenomenon. Three rounds of cryogenic experiments are conducted to investigate the physics behind DTR phenomenon. It starts to emerge when the temperature is reduced to a certain level. After two rounds of cryogenic experiments, the DTR phenomenon becomes permanent from 295 K to 10 K. The nano-scale structure separation induced by the increased thermal strain leads to macroscale structure separation, which results in two parallel heat conduction paths responsible for the DTR phenomenon. By building a new parallel heat transfer model taking both the transient and steady-state electrical and thermal response into consideration, the area ratio of separated CNTs is determined. Furthermore, this result uncovers the existence of coiled morphology and nano-structure evolution of CNTs under cryogenic state and its effect on thermal and electrical transport, especially on their transient behaviors. Knowledge of the thermal transport in these CNTs arrays is important considering the fact that transient thermal response strongly affects their mechanical, optical, and electrical behaviors.

42 ENGINEERING↗
DOE OSTI2021

Hygro-thermo-mechanical modeling of partially saturated air-entrained concrete containing dissolved salt and exposed to freeze-thaw cycles

Highlights: • A model for cyclic freezing of air-entrained concrete containing salt is proposed. • The freeze-thaw hysteresis, also including the effect of dissolved salt, is modeled. • The diffusive and convective flux of salt is coupled to the freeze-thaw processes. • The freeze-thaw response of air-entrained concrete containing salt is reproduced. • Liquid absorption from reservoirs increases with an increased concentration of salt. In cold regions, understanding the freeze-thaw behavior of air-entrained concrete is important for designing durable structures and assessing the remaining service life of existing structures. This study presents a hygro-thermo-mechanical multiphase model that describes the cyclic freeze-thaw behavior of partially saturated air-entrained concrete containing dissolved salt. An equilibrium and a non-equilibrium approach are adopted to model the ice formation, including the freeze-thaw hysteresis, inside the porous network. The model also considers the diffusive and convective transport of the dissolved salt coupled to the freeze-thaw processes. Two examples are presented to verify and highlight the capabilities of the model. The first example shows that the model is capable of reproducing the experimentally observed mechanical response of specimens containing NaC1-solutions of different concentrations. In the second example, a larger absorption of liquid from an external reservoir is obtained with an increasing salt concentration in the reservoir, which is consistent with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Theoretical limit of spatial resolution in diffuse optical tomography using a perturbation model

We have assessed the limit of spatial resolution of timedomain diffuse optical tomography (DOT) based on a perturbation reconstruction model. From the viewpoint of the structure reconstruction accuracy, three different approaches to solving the inverse DOT problem are compared. The first approach involves reconstruction of diffuse tomograms from straight lines, the second – from average curvilinear trajectories of photons and the third – from total banana-shaped distributions of photon trajectories. In order to obtain estimates of resolution, we have derived analytical expressions for the point spread function and modulation transfer function, as well as have performed a numerical experiment on reconstruction of rectangular scattering objects with circular absorbing inhomogeneities. It is shown that in passing from reconstruction from straight lines to reconstruction using distributions of photon trajectories we can improve resolution by almost an order of magnitude and exceed the accuracy of reconstruction of multi-step algorithms used in DOT. (optical tomography)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Diffuse-Interface Capturing Methods for Compressible Two-Phase Flows

Simulation of compressible flows became a routine activity with the appearance of shock-/contact-capturing methods. These methods can determine all waves, particularly discontinuous ones. However, additional difficulties may appear in two-phase and multimaterial flows due to the abrupt variation of thermodynamic properties across the interfacial region, with discontinuous thermodynamical representations at the interfaces. To overcome this difficulty, researchers have developed augmented systems of governing equations to extend the capturing strategy. These extended systems, reviewed here, are termed diffuse-interface models, because they are designed to compute flow variables correctly in numerically diffused zones surrounding interfaces. In particular, they facilitate coupling the dynamics on both sides of the (diffuse) interfaces and tend to the proper pure fluid–governing equations far from the interfaces. This strategy has become efficient for contact interfaces separating fluids that are governed by different equations of state, in the presence or absence of capillary effects, and with phase change. As a result, more sophisticated materials than fluids (e.g., elastic–plastic materials) have been considered as well.

42 ENGINEERING↗
DOE OSTI2015

Evaluation of Technetium Getters to Improve the Performance of Cast Stone

Cast Stone has been selected as the preferred waste form for solidification of aqueous secondary liquid effluents from the Hanford Tank Waste Treatment and Immobilization Plant (WTP) process condensates and low-activity waste (LAW) melter off-gas caustic scrubber effluents. Cast Stone is also being evaluated as a supplemental immobilization technology to provide the necessary LAW treatment capacity to complete the Hanford tank waste cleanup mission in a timely and cost effective manner. One of the major radionuclides that Cast Stone has the potential to immobilize is technetium (Tc). The mechanism for immobilization is through the reduction of the highly mobile Tc(VII) species to the less mobile Tc(IV) species by the blast furnace slag (BFS) used in the Cast Stone formulation. Technetium immobilization through this method would be beneficial because Tc is one of the most difficult contaminants to address at the U.S. Department of Energy (DOE) Hanford Site due to its complex chemical behavior in tank waste, limited incorporation in mid- to high-temperature immobilization processes (vitrification, steam reformation, etc.), and high mobility in subsurface environments. In fact, the Tank Closure and Waste Management Environmental Impact Statement for the Hanford Site, Richland, Washington (TC&WM EIS) identifies technetium-99 ( 99 Tc) as one of the radioactive tank waste components contributing the most to the environmental impact associated with the cleanup of the Hanford Site. The TC&WM EIS, along with an earlier supplemental waste-form risk assessment, used a diffusion-limited release model to estimate the release of different contaminants from the WTP process waste forms. In both of these predictive modeling exercises, where effective diffusivities based on grout performance data available at the time, groundwater at the 100-m down-gradient well exceeded the allowable maximum permissible concentrations for 99 Tc. (900 pCi/L). Recent relatively short-term (63 day) leach tests conducted on both LAW and secondary waste Cast Stone monoliths indicated that 99 Tc diffusivities were at or near diffusivities where the groundwater at the 100-m down-gradient well would exceed the allowable maximum permissible 99 Tc concentrations. There is, therefore, a need and an opportunity to improve the retention of Tc in the Cast Stone waste form. One method to improve the performance of the Cast Stone waste form is through the addition of “getters” that selectively sequester Tc inside Cast Stone.

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

Characterization and Modeling of Grain Boundary Chemistry Evolution in Ferritic Steels under Irradiation

Ferritic/martensitic (FM) steels such as HT-9, T-91 and NF12 with chromium concentrations in the range of 9-12 at.% Cr and high Cr ferritic steels (oxide dispersion strengthened steels with 12-18% Cr) are receiving increasing attention for advanced nuclear applications, e.g. cladding and duct materials for sodium fast reactors, pressure vessels in Generation IV reactors and first wall structures in fusion reactors, thanks to their advantages over austenitic alloys. Predicting the behavior of these alloys under radiation is an essential step towards the use of these alloys. Several radiation-induced phenomena need to be taken into account, including phase separation, solute clustering, and radiation-induced segregation or depletion (RIS) to point defect sinks. RIS at grain boundaries has raised significant interest because of its role in irradiation assisted stress corrosion cracking (IASCC) and corrosion of structural materials. Numerous observations of RIS have been reported on austenitic stainless steels where it is generally found that Cr depletes at grain boundaries, consistently with Cr atoms being oversized in the fcc Fe matrix. While FM and ferritic steels are also subject to RIS at grain boundaries, unlike austenitic steels, the behavior of Cr is less clear with significant scatter and no clear dependency on irradiation condition or alloy type. In addition to the lack of conclusive experimental evidence regarding RIS in F-M alloys, there have been relatively few efforts at modeling RIS behavior in these alloys. The need for predictability of materials behavior and mitigation routes for IASCC requires elucidating the origin of the variable Cr behavior. A systematic detailed high-resolution structural and chemical characterization approach was applied to ion-implanted and neutron-irradiated model Fe-Cr alloys containing from 3 to 18 at.% Cr. Atom probe tomography analyses of the microstructures revealed slight Cr clustering and segregation to dislocations and grain boundaries in the ion-irradiated alloys. More significant segregation was observed in the neutron irradiated alloys. For the more concentrated alloys, irradiation did not affect existing Cr segregation to grain boundaries and segregation to dislocation loops was not observed perhaps due to a change in the dislocation loop structure with increasing Cr concentration. Precipitation of α’ was observed in the neutron irradiated alloys containing over 9 at.% Cr. However ion irradiation appears to suppress the precipitation process. Initial low dose ion irradiation experiments strongly suggest a cascade recoil effect. The systematic analysis of grain boundary orientation on Cr segregation was significantly challenged by carbon contamination during ion irradiation or by existing carbon and therefore carbide formation at grain boundaries (sensitization). The combination of the proposed systematic experimental approach with atomistic modeling of diffusion processes directly addresses the programmatic need for developing and benchmarking predictive models for material degradation taking into account atomistic kinetics parameters

36 MATERIALS SCIENCE↗
NASA NTRS1981

Transport, charge exchange and loss of energetic heavy ions in the earth's radiation belts - Applicability and limitations of theory

Computer simulations of processes which control the relative abundances of ions in the trapping regions of geospace are compared with observations from discriminating ion detectors. Energy losses due to Coulomb collisions between ions and exospheric neutrals are considered, along with charge exchange losses and internal charge exchanges. The time evolution of energetic ion fluxes of equatorially mirroring ions under radial diffusion is modelled to include geomagnetic and geoelectric fluctutations. Limits to the validity of diffusion transport theory are discussed, and the simulation is noted to contain provisions for six ionic charge states and the source effect on the radiation belt oxygen ion distributions. Comparisons are made with ion flux data gathered on Explorer 45 and ISEE-1 spacecraft and results indicate that internal charge exchanges cause the radiation belt ion charge state to be independent of source charge rate characteristics, and relative charge state distribution is independent of the radially diffusive transport rate below the charge state redistribution zone.

Spjeldvik, W. N.↗
NASA NTRS1987

Effect of double layers on magnetosphere-ionosphere coupling

The dynamic aspects of auroral current structures are reviewed with emphasis on consequences for models of microscopic turbulence (MT). A number of models of MT are introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. The effect of a double layer (DL) electric field which scales with the plasma temperature and the Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is shown that the DL model is less diffusive than the resistive model, indicating the possibility of narrow intense current structures.

Lysak, Robert L.↗
NASA NTRS1999

Ocean Turbulence: Two-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 3

In papers 1 and 2 we have presented the results of the most updated 1-point closure model for the turbulent vertical diffusivities of momentum, heat and salt, K(sub m,h,s). In this paper, we derive the analytic expressions for K(sub m,h,s) using a new 2-point closure model that has recently been developed and successfully tested against some approx. 80 turbulence statistics for different flows. The new model has no free parameters. The expressions for K(sub m, h. s) are analytical functions of two stability parameters: the Turner number R(sub rho) (salinity gradient/temperature gradient) and the Richardson number R(sub i) (temperature gradient/shear). The turbulent kinetic energy K and its rate of dissipation may be taken local or non-local (K-epsilon model). Contrary to all previous models that to describe turbulent mixing below the mixed layer (ML) have adopted three adjustable "background diffusivities" for momentum. heat and salt, we propose a model that avoids such adjustable diffusivities. We assume that below the ML, K(sub m,h,s) have the same functional dependence on R(sub i) and R(sub rho) derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking measured by Gargett et al. (1981). The procedure frees the model from adjustable background diffusivities and indeed we use the same model throughout the entire vertical extent of the ocean. Using the new K(sub m,h, s), we run an O-GCM and present a variety of results that we compare with Levitus and the KPP model. Since the traditional 1-point (used in papers 1 and 2) and the new 2-point closure models used here represent different modeling philosophies and procedures, testing them in an O-GCM is indispensable. The basic motivation is to show that the new 2-point closure model gives results that are overall superior to the 1-point closure in spite of the fact that the latter rely on several adjustable parameters while the new 2-point closure has none. After the extensive comparisons presented in papers 1 and 2, we conclude that the new model presented here is overall superior for it not only is parameter free but also 2 because is part of a more general turbulence model that has been previously successfully tested on a wide variety of other types of turbulent flows.

Canuto, V. M.↗
DOE OSTI2023

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

The modeling of gas phase permeation through iron and nickel membranes

The gas phase permeation of hydrogen through metal membranes encompasses many kinetic processes. This paper reviews a permeation model which incorporates second order gas-surface reaction kinetics with simple bulk diffusion. The model is used to investigate the effect of this particular surface reaction of steady-state permeation. The dependence of the steady-state permeation flux on temperature, pressure, and thickness of the membrane has been calculated. The model predicts that the bulk controlled steady-state flux will change to a surface limited steady-state flux as either the temperature or thickness of the membrane is reduced. Finally, using independently derived parameters, the model is compared with permeation measurements on iron and nickel membranes.

Kuhn, David K.↗