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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Long-time integration methods for mesoscopic models of pattern-forming systems

Spectral methods for simulation of a mesoscopic diffusion model of surface pattern formation are evaluated for long simulation times. Backwards-differencing time-integration, coupled with an underlying Newton-Krylov nonlinear solver (SUNDIALS-CVODE), is found to substantially accelerate simulations, without the typical requirement of preconditioning. Quasi-equilibrium simulations of patterned phases predicted by the model are shown to agree well with linear stability analysis. Simulation results of the effect of repulsive particle-particle interactions on pattern relaxation time and short/long-range order are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Report on simulation of fission gas and fission product diffusion in UO 2

In UO 2 nuclear fuel, the retention and release of fission gas atoms such as xenon (Xe) are important for nuclear fuel performance by, for example, reducing the fuel thermal conductivity, causing fuel swelling that leads to mechanical interaction with the clad, increasing the plenum pressure and reducing the fuel–clad gap thermal conductivity. We use multi-­scale simulations to determine fission gas diffusion mechanisms as well as the corresponding rates in UO 2 under both intrinsic and irradiation conditions. In addition to Xe and Kr, the fission products Zr, Ru, Ce, Y, La, Sr and Ba have been investigated. Density functional theory (DFT) calculations are used to study formation, binding and migration energies of small clusters of Xe atoms and vacancies. Empirical potential calculations enable us to determine the corresponding entropies and attempt frequencies for migration as well as investigate the properties of large clusters or small fission gas bubbles. A continuum reaction-­diffusion model is developed for Xe and point defects based on the mechanisms and rates obtained from atomistic simulations. Effective fission gas diffusivities are then obtained by solving this set of equations for different chemical and irradiation conditions using the MARMOT phase field code. The predictions are compared to available experimental data. The importance of the large Xe U3O cluster (a Xe atom in a uranium + oxygen vacancy trap site with two bound uranium vacancies) is emphasized, which is a consequence of its high mobility and high binding energy. We find that the Xe U3O cluster gives Xe diffusion coefficients that are higher for intrinsic conditions than under irradiation over a wide range of temperatures. Under irradiation the fast-­moving Xe U3O cluster recombines quickly with irradiation-induced interstitial U ions, while this mechanism is less important for intrinsic conditions. The net result is higher concentration of the Xe U3O cluster for intrinsic conditions than under irradiation. We speculate that differences in the irradiation conditions and their impact on the Xe U3O cluster can explain the wide range of diffusivities reported in experimental studies. However, all vacancy-­mediated mechanisms underestimate the Xe diffusivity compared to the empirical radiation-­enhanced rate used in most fission gas release models. We investigate the possibility that diffusion of small fission gas bubbles or extended Xe-­vacancy clusters may give rise to the observed radiation-­enhanced diffusion coefficient. These studies highlight the importance of U divacancies and an octahedron coordination of uranium vacancies encompassing a Xe fission gas atom. The latter cluster can migrate via a multistep mechanism with a rather low effective barrier, which together with irradiation-induced clusters of uranium vacancies, gives rise to the irradiation-enhanced diffusion coefficient observed in experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Synthesis and Modeling of Hollow Intermetallic Ni-Zn Nanoparticles Formed by the Kirkendall Effect

Intermetallic Ni–Zn nanoparticles (NPs) were synthesized via the chemical conversion of nickel NPs using a zerovalent organometallic zinc precursor. After the injection of a diethylzinc solution, Ni NPs progressively transformed from a solid to a hollow Ni–Zn intermetallic structure with time. During the transformation of Ni NPs to intermetallic structures, they retained their overall spherical morphology. The growth mechanism for the solid-to-hollow nanoparticle transformation is ascribed to the nanoscale Kirkendall effect due to unequal diffusion rates of Ni and Zn. We develop a diffusion model for nonreactive, homogeneous, diffusion-controlled intermetallic hollow NP formation including moving boundaries at the interfaces of void–solid and solid–bulk solutions. Apparent diffusion coefficients for both metals and vacancy were evaluated from modeling the time-dependent growth of the void. The apparent diffusion coefficients obtained in this system compared favorably with results from measurement at grain boundaries in bulk Ni–Zn. This study represents the first combined experimental modeling of the formation of hollow nanostructures by the nanoscale Kirkendall effect.

Jana, Subhra↗
DOE OSTI2014

Determination of charge-carrier diffusion length in the photosensing layer of HgCdTe n-on-p photovoltaic infrared focal plane array detectors

In the present paper, we propose a method for evaluating the bulk diffusion length of minority charge carriers in the photosensing layer of photovoltaic focal plane array (FPA) photodetectors. The method is based on scanning a strip-shaped illumination spot with one of the detector diodes at a low level of photocurrents j{sub ph} being registered; such scanning provides data for subsequent analysis of measured spot-scan profiles within a simple diffusion model. The asymptotic behavior of the effective (at j{sub ph} ≠ 0) charge-carrier diffusion length l{sub d} {sub eff} as a function of j{sub ph} for j{sub ph} → 0 inferred from our experimental data proved to be consistent with the behavior of l{sub d} {sub eff} vs j{sub ph} as predicted by the model, while the obtained values of the bulk diffusion length of minority carriers (electrons) in the p-HgCdTe film of investigated HgCdTe n-on-p FPA photodetectors were found to be in a good agreement with the previously reported carrier diffusion-length values for HgCdTe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Controlled plasma–droplet interactions: a quantitative study of OH transfer in plasma–liquid interaction

Plasmas in contact with liquids are a rich source of OH radicals and have been extensively studied in the last decade to leverage the ability to generate chemically reactive species in gas phase plasmas to decompose organics. Multiphase transfer of OH radicals is highly transport limited and to overcome transport limits, the plasma activation of aerosols, small liquid droplets, interspersed in the plasma has been proposed. In this work, we report a combined experimental and modeling study of a controlled plasma–droplet interaction experiment using a diffuse RF glow discharge in He + 0.2% H 2 O with detailed plasma diagnostics, ex situ analysis of the plasma-induced chemistry in the droplet containing formate, droplet trajectory and size measurements. This enables a quantitative study of the reactivity transfer of OH from the gas phase plasma to the liquid phase and how its diffusion limitations impact formate decomposition in the water droplet. For a droplet with a diameter of 36 μm, we observed 50% reduction in formate concentration in the droplets after plasma treatment for droplet residence times in the plasma of ~10 ms. These short droplet residence times in the plasma allow in some cases for droplet size reductions of ~5% in spite gas temperatures of 360 K. A one-dimensional reaction–diffusion model was used to calculate the OH transport and formate oxidation inside the droplet and was able to predict the conversion of formate by plasma in a droplet without any fitting parameters. The model further shows that formate conversion is dominated by near-interfacial reactions with OH radicals and is limited by diffusion of formate in the droplet. Here, the results show that a controlled plasma–micro-droplet reactor as reported in this study might be an excellent tool for detailed quantitative plasma–liquid interaction studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Simulating the density of organic species in the atmosphere of Titan with a coupled ion-neutral photochemical model

We introduce a one-dimensional coupled ion-neutral photochemical kinetics and diffusion model to study the atmospheric composition of Titan in light of new theoretical kinetics calculations and scientific findings from the Cassini–Huygens mission. The model extends from the surface to the exobase. The atmospheric background, boundary conditions, vertical transport and aerosol opacity are all constrained by the Cassini–Huygens observations. The chemical network includes reactions between hydrocarbons, nitrogen and oxygen bearing species. It takes into account neutrals and both positive and negative ions with masses extending up to 116 and 74 u, respectively. We incorporate high-resolution isotopic photoabsorption and photodissociation cross sections for N 2 as well as new photodissociation branching ratios for CH 4 and C 2 H 2 . Ab initio transition state theory calculations are performed in order to estimate the rate coefficients and products for critical reactions. Critical reactions of production and loss for neutrals and ions are quantitatively assessed and thoroughly discussed. The vertical distributions of neutrals and ions predicted by the model generally reproduce observational data, suggesting that for the small species most chemical processes in Titan’s atmosphere and ionosphere are adequately described and understood; some differences are highlighted. Notable remaining issues include (i) the total positive ion density (essentially HCNH + ) in the upper ionosphere, (ii) the low mass negative ion densities (CN - , C 3 N - /C4H - ) in the upper atmosphere, and (iii) the minor oxygen-bearing species (CO 2 , H 2 O) density in the stratosphere. Pathways towards complex molecules and the impact of aerosols (UV shielding, atomic and molecular hydrogen budget, nitriles heterogeneous chemistry and condensation) are evaluated in the model, along with lifetimes and solar cycle variations.

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

Evaluation of light collection from highly scattering media using wavelength-shifting fibers

Opaque scintillators are designed to have a short scattering length such that scintillation photons are probabilistically confined to a small region of space about their origin. The benefit of this feature is that information on the interaction event topology can be recorded with greater fidelity than in traditional highly transparent media with sensors at large distances from the light production region. Opaque scintillator detectors rely on wavelength-shifting fibers to extract the scintillation light; however, the efficiency of light collection has not yet been directly measured in experiment. We measured the efficiency of light collection as a function of the optical parameters of an opaque liquid and the distance from the origin of the light to the fiber. We use the experimental data to validate a Monte Carlo model of light transport and collection and discuss a simple diffusion model that reproduces the results of Monte Carlo simulation with high fidelity. Furthermore, this combination of validated models has the potential for use in predictions of performance in various designs of future opaque scintillator detectors such as LiquidO.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Interpreting inorganic compositional depth profiles to understand the rate-limiting step in vapor phase infiltration processes

Vapor phase infiltration (VPI) is a post-polymerization modification technique that infuses inorganics into polymers to create organic–inorganic hybrid materials with new properties. Much is yet to be understood about the chemical kinetics underlying the VPI process. The aim of this study is to create a greater understanding of the process kinetics that govern the infiltration of trimethyl aluminum (TMA) and TiCl 4 into PMMA to form inorganic-PMMA hybrid materials. To gain insight, this paper initially examines the predicted results for the spatiotemporal concentrations of inorganics computed from a recently posited reaction–diffusion model for VPI. This model provides insight on how the Damköhler number (reaction versus diffusion rates) and non-Fickian diffusional processes (hindering) that result from the material transforming from a polymer to a hybrid can affect the evolution of inorganic concentration depth profiles with time. Subsequently, experimental XPS depth profiles are collected for TMA and TiCl 4 infiltrated PMMA films at 90 °C and 135 °C. The functional behavior of these depth profiles at varying infiltration times are qualitatively compared to various computed predictions and conclusions are drawn about the mechanisms of each of these processes. TMA infiltration into PMMA appears to transition from a diffusion-limited process at low temperatures (90 °C) to a reaction-limited process at high temperatures (135 °C) for the film thicknesses investigated here (200 nm). While TMA appears to fully infiltrate these 200 nm PMMA films within a few hours, TiCl 4 infiltration into PMMA is considerably slower, with full saturation not occurring even after 2 days of precursor exposure. Infiltration at 90 °C is so slow that no clear conclusions about mechanism can be drawn; however, at 135 °C, the TiCl 4 infiltration into PMMA is clearly a reaction-limited process, with TiCl 4 permeating the entire thickness (at low concentrations) within only a few minutes, but inorganic loading continuously increasing in a uniform manner over a course of 2 days. Near-surface deviations from the uniform-loading expected for a reaction-limited process also suggest that diffusional hindering is high for TiCl 4 infiltration into PMMA. In conclusion, these results demonstrate a new, ex situ analysis approach for investigating the rate-limiting process mechanisms for vapor phase infiltration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗
DOE OSTI2021

Neutron Tunneling: A New Mechanism to Power Explosive Phenomena in Neutron Stars, Magnetars, and Neutron Star Mergers

Neutron tunneling between neutron-rich nuclei in inhomogeneous dense matter encountered in neutron star crusts can release enormous energy on a short timescale to power explosive phenomena in neutron stars. In this work, we clarify aspects of this process that can occur in the outer regions of neutron stars when oscillations or cataclysmic events increase the ambient density. We use a time-dependent Hartree–Fock–Bogoliubov formalism to determine the rate of neutron diffusion and find that large amounts of energy can be released rapidly. The roles of nuclear binding, two-body interaction, and pairing in neutron diffusion times are investigated. We consider a one-dimensional quantum diffusion model and extend our analysis to study the impact of diffusion in three dimensions. Here, we find that these novel neutron transfer reactions can generate energy in the amount of ≃ 10 40 –10 44 erg under suitable conditions and assumptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Flow Rate Dependent 1D Model for Thermally Stratified Hot-Water Energy Storage

Stratified tank models are used to simulate thermal storage in applications such as residential or commercial hot-water storage tanks, chilled-water storage tanks, and solar thermal systems. The energy efficiency of these applications relates to the system components and the level of stratification maintained during various flow events in the tank. One-dimensional (1D) models are used in building energy simulations because of the short computation time but often do not include flow-rate dependent mixing. The accuracy of 1D models for plug flow, plug flow with axial conduction, and two convection eddy-diffusivity models were compared with experimental data sets for discharging a 50-gal residential tank and recharging the tank with hot water from an external hot-water source. A minimum and maximum relationship for the eddy diffusivity factor were found at Re <2100 and >10,000 for recirculation of hot water to the top of the tank and vertical tubes inletting cold water at the bottom. The root mean square error decreased from >4 °C to near 2 °C when considering flow-based mixing models during heating, while the exponential decay of the eddy diffusion results in a root mean square error reduction of 1 °C for cone-shaped diffusers that begin to relaminarize flow at the inlet.

1D flow model↗
DOE OSTI2016

THE ROLE OF COSMIC-RAY PRESSURE IN ACCELERATING GALACTIC OUTFLOWS

We study the formation of galactic outflows from supernova (SN) explosions with the moving-mesh code AREPO in a stratified column of gas with a surface density similar to the Milky Way disk at the solar circle. We compare different simulation models for SN placement and energy feedback, including cosmic rays (CRs), and find that models that place SNe in dense gas and account for CR diffusion are able to drive outflows with similar mass loading as obtained from a random placement of SNe with no CRs. Despite this similarity, CR-driven outflows differ in several other key properties including their overall clumpiness and velocity. Moreover, the forces driving these outflows originate in different sources of pressure, with the CR diffusion model relying on non-thermal pressure gradients to create an outflow driven by internal pressure and the random-placement model depending on kinetic pressure gradients to propel a ballistic outflow. CRs therefore appear to be non-negligible physics in the formation of outflows from the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A semi-automated algorithm for designing stellarator divertor and limiter plates and application to HSX

We present a semi-automated algorithm for designing three-dimensional divertor or limiter plates targeting low heat loads. The algorithm designs the plates in two stages: firstly, the parallel heat flux distribution is caught on vertically-inclined plates at one or several toroidal locations. Secondly, the power per unit area is reduced by stretching, tilting and bending the plates toroidally. Heat transport is modelled using the EMC3-Lite code, which uses an anisotropic diffusion model. We apply this scheme to HSX, a medium-sized stellarator located at the University of Wisconsin–Madison. Starting from the current machine with an extended vessel wall, we construct plates which are able to effectively catch and spread the heat for three different magnetic configurations. The scheme has a computational cost in the order of tens of CPU-minutes, making it a powerful tool for semi-automated plasma-facing component design in three-dimensional environments.

anisotropic diffusion↗
DOE OSTI2022

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗
DOE OSTI2017

Kinetic and Mechanistic Evaluation of Inorganic Arsenic Species Adsorption onto Humic Acid Grafted Magnetite Nanoparticles

In this study, humic acid coated magnetic iron oxide nanoparticles (HA-MNP) were synthesized, characterized, and studied for the removal of toxic inorganic arsenic species from aqueous media. The adsorption of As(III) and As(V) followed pseudo-second order kinetics and the observed data were accurately modeled employing Freundlich adsorption isotherm. Application of the Weber and Morris intraparticle diffusion model to the observed kinetic data suggest that the adsorption occurs in three distinct stages, fast, intermediate, and slow steps. We propose the initial stage is governed by surface association, followed by intraparticle diffusion of arsenic through the HA matrix and finally, chemical reaction or bonding between the arsenic species and HA functionality. The HAMNP nano-adsorbent absorbs > 95% of the inorganic arsenic species with an absorption capacity of 12.2-12.6 mg/g from aqueous media and is effective under a variety of conditions. Post arsenic adsorption characterization of the nanoparticles suggests that As(III) binds with the carboxylate group of HA through a proposed ester type linkage while electrophilic reactions can occur between the nucleophilic functional groups present in HA and the electrophilic arsenic atom in As(V). Finally, the results obtained demonstrated that HA-MNP are robust and have promise for effective As (III) and As(V) remediation.

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

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Measurements of Transverse Beam Diffusion Rates in the Fermilab Tevatron Collider

The transverse beam diffusion rate vs. particle oscillation amplitude was measured in the Tevatron using collimator scans. All collimator jaws except one were retracted. As the jaw of interest was moved in small steps, the local shower rates were recorded as a function of time. By using a diffusion model, the time evolution of losses could be related to the diffusion rate at the collimator position. Preliminary results of these measurements are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

SOLAR MODELS WITH NEW LOW METAL ABUNDANCES

In the past decade, the photospheric abundances of the Sun had been revised several times by many observers. The standard solar models constructed with the new low-metal abundances disagree with helioseismic results and detected neutrino fluxes. The solar model problem has puzzled some stellar physicists for more than 10 years. Rotation, enhanced diffusion, convection overshoot, and magnetic fields are used to reconcile the new abundances with helioseismology. The too low helium subsurface abundance in enhanced diffusion models can be improved by the mixing caused by rotation and magnetic fields. The problem of the depth of the convective zone in rotating models can be resolved by convection overshoot. Consequently, the Asplund–Grevesse–Sauval rotation model including overshooting (AGSR) reproduces the seismically inferred sound-speed and density profiles and the convection zone depth as well as the Grevesse and Sauval model computed before. But this model fails to reproduce the surface helium abundance, which is 0.2393 (2.6 σ away from the seismic value), and neutrino fluxes. The magnetic model called AGSM keeps the agreement of the AGSR and improves the prediction of the surface helium abundance. The observed separation ratios r {sub 02} and r {sub 13} are reasonably reproduced by AGSM. Moreover, neutrino fluxes calculated by this model are not far from the detected neutrino fluxes and the predictions of previous works.

79 ASTRONOMY AND ASTROPHYSICS↗