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,621 records · Page 90

A model for the compositions of non-stoichiometric intermediate phases formed by diffusion reactions, and its application to Nb 3 Sn superconductors

In this work we explore the compositions of non-stoichiometric intermediate phases formed by diffusion reactions: a mathematical framework is developed and tested against the specific case of Nb 3 Sn superconductors. In the first part, the governing equations for the bulk diffusion and interphase interface reactions during the growth of a compound are derived, numerical solutions to which give both the composition profile and growth rate of the compound layer. The analytic solutions are obtained with certain approximations made. In the second part, we explain an effect that the composition characteristics of compounds can be quite different depending on whether it is the bulk diffusion or grain boundary diffusion that dominates in the compounds, and that “frozen” bulk diffusion leads to unique composition characteristics that the bulk composition of a compound layer remains unchanged after its initial formation instead of varying with the diffusion reaction system; here the model is modified for the case of grain boundary diffusion. Lastly, we apply this model to the Nb 3 Sn superconductors and propose approaches to control their compositions.

36 MATERIALS SCIENCE↗

The Hot Horizontal-Branch Stars in NGC288 - Effects of Diffusion and Stratification on Their Atmospheric Parameters*

Context. NGC288 is a globular cluster with a well developed blue horizontal branch covering the so-called u-jump which indicates the onset of diffusion. It is therefore well suited to study the effects of diffusion in blue horizontal branch (HB) stars. Aims. We compare observed abundances to predictions from stellar evolution models calculated with diffusion and from stratified atmospheric models. We verify the effect of using stratified model spectra to derive atmospheric parameters. In addition we investigate the nature of the overluminous blue HB stars around the u-jump. Methods. We define a new photometric index sz from uvby measurements that is gravity sensitive between 8 000K and 12 000 K. Using medium-resolution spectra and Stroemgren photometry we determine atmospheric parameters (Teff, logg) and abundances for the blue HB stars. We use both homogeneous and stratified model spectra for our spectroscopic analyses. Results. The atmospheric parameters and masses of the hot HB stars in NGC288 show a behaviour seen also in other clusters for temperatures between 9 000K and 14 000 K. Outside this temperature range, however, they follow rather the results found for such stars in (omega)Cen. The abundances derived from our observations are for most elements (except He and P) within the abundance range expected from evolutionary models that include the effects of atomic diffusion and assume a surface mixed mass of 10(exp −7) M. The abundances predicted by stratified model atmospheres are generally significantly more extreme than observed, except for Mg. The use of stratified model spectra to determine effective temperatures, surface gravities and masses moves the hotter stars to a closer agreement with canonical evolutionary predictions. Conclusions. Our results show definite promise towards solving the long-standing issue of surface gravity and mass discrepancies for hot HB stars, but there is still much work needed to arrive at a self-consistent solution.

horizontal branch (HB)↗

On the distribution of He+ in the plasmasphere from observations of resonantly scattered He II 304-A radiation

The intensity of the He II 304-A radiation in the night sky was measured from a sounding rocket launched from White Sands Missile Range on June 9, 1972. The data obtained on this flight are compared with predictions of theoretical models of the He(+) altitude distribution within the plasmasphere. The observations are not consistent with an isothermal equilibrium model of the helium ion in the plasmasphere, confirming earlier results. The data are best fit by a diffussive equilibrium model of exospheric He(+) having a plasma temperature gradient of 0.5 K per km and a number density at 1600 km of 200,000 He(+) ions per cu cm. In this model the ion densities are calculated without allowance for thermal diffusion. A similar model that allows for thermal diffusion and that has a He(+) number density at 1600 km of 1500 per cu cm and a temperature gradient of 1.2 K per km yields a less satisfactory fit to the observations but is more consistent with in situ satellite observations.

Paresce, F.↗

Numerical modeling of HgCdTe solidification: Effects of phase diagram, double-diffusion convection and microgravity level

A numerical model of HgCdTe solidification was implemented using finite the element code FIDAP. Model verification was done using both experimental data and numerical test problems. The model was used to evaluate possible effects of double-diffusion convection in molten material, and microgravity level on concentration distribution in the solidified HgCdTe. Particular attention was paid to incorporation of HgCdTe phase diagram. It was found, that below a critical microgravity amplitude, the maximum convective velocity in the melt appears virtually independent on the microgravity vector orientation. Good agreement between predicted interface shape and an interface obtained experimentally by quenching was achieved. The results of numerical modeling are presented in the form of video film.

Bune, Andris V.↗

Exploring Li-Air Batteries for High Specific Energy and High Power Applications: A Simulation Study

Commercialization of lithium-air batteries face many challenges, such as electrolyte decomposition, short cycle life, low energy efficiency, low power density, etc. However, commercialization of Li-air batteries for mass sensitive applications such as electric vehicles, portable power source, and drones is more challenging due to additional constraints of safety, electrolyte evaporation, high specific energy requirements, and reliable discharge times. In this presentation, we will present our finite element simulation results comparing Li-O2 and Li-air batteries using power density, energy density, specific power, and discharge times as metrics to evaluate different electrolytes and electrode geometry to reduce total mass and maximize discharge current. We use a finite element model and a discharge product model developed in which is based on porous electrode and concentrated electrolyte theories and the discharge product is modeled using quantum tunneling; for reaction kinetics and oxygen diffusion, an improved model was used. The electrolyte properties such as ion conductivity and ion diffusion were obtained from Molecular Dynamics (MD) simulations while the other parameters for the finite element model were calibrated to match experiments at high discharge current densities (>1.5 mA/cm2). The mass densities of different electrolytes were computed using MD simulations as well. For this presentation, we examine the practical electrochemical mass of a system at different current ratings, the sensitivity of mass to the use of ambient air as compared to pure oxygen as well as the electrolyte, which affects maximum current density and the total mass associated with the electrolyte (which includes the mass of additional components), and, the optimization of battery geometry for total discharge time, average discharge voltage, maximum discharge current, and minimum electrochemical mass.

Mehta, M.↗

The diffusion approximation versus the telegraph equation for modeling solar energetic particle transport with adiabatic focusing. I. Isotropic pitch-angle scattering

The diffusion approximation to the Fokker-Planck equation is commonly used to model the transport of solar energetic particles in interplanetary space. In this study, we present exact analytical predictions of a higher order telegraph approximation for particle transport and compare them with the corresponding predictions of the diffusion approximation and numerical solutions of the full Fokker-Planck equation. We specifically investigate the role of the adiabatic focusing effect of a spatially varying magnetic field on an evolving particle distribution. Comparison of the analytical and numerical results shows that the telegraph approximation reproduces the particle intensity profiles much more accurately than does the diffusion approximation, especially when the focusing is strong. However, the telegraph approximation appears to offer no significant advantage over the diffusion approximation for calculating the particle anisotropy. The telegraph approximation can be a useful tool for describing both diffusive and wave-like aspects of the cosmic-ray transport.

79 ASTRONOMY AND ASTROPHYSICS↗

Shell Model for Atomistic Simulation of Lithium Diffusion in Mixed Mn/Ti Oxides

Mixed Mn/Ti oxides present attractive physicochemical properties such as their ability to accommodate Li for application in Li-ion batteries. In this work, atomic parameters for Mn were developed to extend an existing shell model of the Li-Ti-O system and allow simulations of pure and lithiated Mn and mixed Mn/Ti oxide polymorphs. The shell model yielded good agreement with experimentally-derived structures (i.e. lattice parameters and inter-atomic distances) and represented an improvement over existing potential models. The shell model was employed in molecular dynamics (MD) simulations of Li diffusion in the 1×1 c direction channels of LixMn1 yTiyO2 with the rutile structure, where 0 ≤ x ≤ 0.25 and 0 ≤ y ≤ 1. In the infinite dilution limit, the arrangement of Mn and Ti ions in the lattice was found to have a significant effect on the activation energy for Li diffusion in the c channels due to the destabilization of half of the interstitial octahedral sites. Anomalous diffusion was demonstrated for Li concentrations as low as x = 0.125, with a single Li ion positioned in every other c channel. Further increase in Li concentration showed not only the substantial effect of Li-Li repulsive interactions on Li mobility but also their influence on the time dependence of Li diffusion. The results of the MD simulations can inform intrinsic structure-property relationships for the rational design of improved electrode materials for Li-ion batteries.

Environmental Molecular Sciences Laboratory↗

Coupled social and infrastructure approaches for enhancing solar energy adoption. Final Report

The goal of this project was to work with rural electric cooperatives to facilitate the diffusion of solar energy adoption in households located in the rural and semi-urban areas of Virginia by identifying social and behavioral factors that might be unique to rural regions; and develop a model to calculate the solar adoption propensity score for household based on their demographics, social and behavioral characteristics which would provide an objective metric to cooperatives that can be further used to do targeted marketing of rooftop solar panels. This was achieved through the following tasks: (1) Conducted a survey of the members of Virginia electric cooperatives to identify demographic, social, financial and behavioral attributes of individuals who are likely to adopt rooftop solar panels. (2) Developed a highly detailed, data-driven, agent-based model of the population of Virginia, focusing on the rural regions. (3) Developed diffusion models that use social, behavioral, and demographic factors, and peer effects to study their impact on solar adoption in rural areas. (4) Built a prototype tool based on the diffusion model to help study market segmentation in rural areas and made it available to National Rural Electric Cooperative Association (NRECA). (5) Results and recommendations derived from the model were provided to NRECA to be shared with participating cooperatives. (6) Results were published in peer reviewed journals, conference proceedings and book chapters, and ideas disseminated through presentations and newsletters. There were several important methodological contributions made under this project which are detailed in the published papers, including: (1) Built a decision-adjusted model for predicting adoptors with imbalanced training data; (2) designed seeding strategies to maximize adoption given a fixed budget; (3) built a methodology to compare different agent based models; (4) created models to identify important factors that influence decision to adopt solar panels; and (5) built a methodology for building household profiles of solar generation to study the duck curve phenomenon. The team included members from the University of Virginia (lead), National Rural Electric Cooperative Association (NRECA), Arizona State University, Virginia Tech and Sandia National Laboratory. Note that no individual entity or stakeholder has incentive to promote solar in rural regions. Most of the research and work focuses around urban regions where the potential for growth in solar adoption is higher due to higher population density. This puts rural areas at a disadvantage. By improving the diffusion of solar adoption in rural parts of the country, we can not only provide clean energy to rural areas but also promote job growth and improves energy independence.

14 SOLAR ENERGY↗

Effects of Hydrogen Bonding on Nuclear Data Development of Liquid Anhydrous HF

Anhydrous Hydrogen Fluoride (HF) at high temperatures and pressures is used to process and manufacture nuclear fuel. As HF is often used directly with uranium, correct neutron thermal scattering cross sections are crucial to criticality safety applications. Classical molecular dynamics (CMD) simulation of the flexible HF system was used to create the thermal scattering law (TSL) and cross sections. The initial 2-site model is used in LAMMPS, and it can not capture the H-bond. To correctly represent the H-bond effects, a second, 3-site model was constructed in GROMACS. The 3-site model handled H-bonds by connecting a massless charge to the molecule. Key model parameters were compared to experimental data to verify the approach and models. To get the normalized VACF, the model was compared using hydrogen and fluorine bond length, density, potential energy, and diffusion coefficient. The phonon DOSs for both models were derived from the normalized VACF. DOSs were used to estimate the TSL ( S ( α, β )) and neutron thermal scattering cross sections for hydrogen in HF. The TSLs were evaluated using the FLASSH code with the Schofield diffusion model. It was observed that the representation of the hydrogen bonding changes the TSL's diffusional contributions. This is represented in the low energy scattering cross section, where intermolecular binding effects shift the cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

Kinetics of pack aluminization of nickel

The kinetics of pack aluminization of unalloyed nickel in packs of varying aluminum activity with various halide activators were studied. Surface compositions of the coatings as functions of time, temperature, and pack composition were obtained in order to establish the boundary conditions for diffusion in the system. The structure of the packs was also examined in order to clarify the mechanism of aluminum transport. The results indicate that the kinetics of pack aluminization are controlled jointly by gas diffusion in the pack and solid diffusion in the coating. Levine and Caves' model for gas diffusion was combined with calculations of rates of diffusion in the solid to formulate a more complete theory for the kinetics of pack aluminization.

Seigle, L. L.↗

Diffusion in intact secondary cell wall models of plants at different equilibrium moisture content

Secondary plant cell walls are composed of carbohydrate and lignin polymers, and collectively represent a significant renewable resource. Leveraging these resources depends in part on a mechanistic understanding for diffusive processes within plant cell walls. Common wood protection treatments and biomass conversion processes to create biorefinery feedstocks feature ion or solvent diffusion within the cell wall. X-ray fluorescence microscopy experiments have determined that ionic diffusion rates are dependent on cell wall hydration as well as the ionic species through non-linear relationships. In this work, we use classical molecular dynamics simulations to map the diffusion behavior of different plant cell wall components (cellulose, hemicellulose, lignin), ions (Na + , K + , Cu2 + , Cl - ) and water within a model for an intact plant cell wall at various hydration states (3-30 wt% water). From these simulations, we analyze the contacts between different plant cell wall components with each other and their interaction with the ions. Generally, diffusion increases with increasing hydration, with lignin and hemicellulose components increasing diffusion by an order of magnitude over the tested hydration range. Ion diffusion depends on charge. Positively charged cations preferentially interact with hemicellulose components, which include negatively charged carboxylates. As a result, positive ions diffuse more slowly than negatively charged ions. Measured diffusion coefficients are largely observed to best fit piecewise linear trends, with an inflection point between 10 and 15% hydration. These observations shed light onto the molecular mechanisms for diffusive processes within secondary plant cell walls at atomic resolution.

09 BIOMASS FUELS↗

Modeling of column apparatuses: A review

This paper presents a review of the modeling method on the base of the physical approximations of the mechanics of continua, which have been developed for processes in column apparatuses. This method includes diffusion type of model for modeling of mass transfer with chemical reaction in column apparatuses with and without circulation zones. The diffusion type of model is used for modeling of scale effect in column apparatuses too. The study concluded that the proposal method is possibility for investigation the influence of radial non uniformity of the velocity distribution on the process efficiency, influence of zones breadths on the mass transfer efficiency in the column. The method of the column apparatuses modeling can be used for modeling of physical and chemical absorption, chemical adsorption, homogeneous and heterogeneous (catalytic) chemical reactions, airlift reactors for chemical and photochemical reactions.

42 ENGINEERING↗

Modeling Radionuclide Releases from TRISO Particles by Simultaneous Diffusion Through and Corrosion of the Silicon Carbide Barrier Layer

TRISO nuclear fuel particles that are less than 1 mm in diameter are designed with multiple barrier layers to retain fission products both during reactor operations and for long-term geological disposal. The primary barrier is a 35 μm thick silicon carbide (SiC-a highly impermeable semi-metal) layer for which data are available on the diffusion of short-lived fission products at high temperatures (> 1000 °C). However, for a geological repository, this layer may contact brine and hence corrode even at ambient temperatures. As an initial approach to assess the effectiveness of the SiC barrier for geological repositories, ranges of fission product diffusivities and corrosion rates for SiC are modeled concurrently with the simultaneous effect of radioactive decay. Using measured corrosion rates of SiC, if the diffusivity is more than about 10 -20 m 2 /s, fission product releases may occur before the SiC barrier has corroded to the point of breach. For diffusivities less than about 10 -21 m 2 /s there may not be significant diffusional releases prior to SiC barrier removal/breach by corrosion. This work shows the importance of estimating diffusivities in SiC at geological repository temperatures, and highlights the relevance of evaluating the porosity/permeability evolution of the SiC layer in a geologic environment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Leveraging inter-firm influence in the diffusion of energy efficiency technologies: An agent-based model

Energy efficiency technologies (EETs) are crucial for saving energy and reducing carbon dioxide emissions. However, the diffusion of EETs in small and medium-sized enterprises is rather slow. Literature shows the interactions between innovation adopters and potential adopters have significant impacts on innovation diffusion. Enterprises lack the motivation to share information, and EETs usually lack observability, which suppress the inter-firm influence. Thus, an information platform, together with proper policies encouraging or forcing enterprises to disclose EET-related information, should help harness inter-firm influence to accelerate EETs' diffusion. To explore whether and how such an information platform affects EETs' diffusion in small and medium-sized enterprises, this report builds an agent-based model to mimic EET diffusion processes. Based on a series of controlled numerical experiments, some counter-intuitive phenomena are discovered and explained. The results show that the information platform is a double-edged sword that notably accelerates EETs' diffusion by approximately 47% but may also boost negative information to diffuse even faster and delay massive adoption of EETs. Increasing network density and the intensity of inter-firm influence are effective to speed EET diffusion, but their impacts diminish drastically after reaching some critical values (0.05 and 0.15 respectively) and eventually harm the stability of the system. Ultimately, the findings implicate that EET suppliers should carefully launch their promising but immature products; policies that can reduce the perceived risk by enterprises and the effort to maintain an informative rather than judgmental information platform can prominently mitigate the negative side effects brought by high fluidity of information.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A laserball calibration device for the SNO+ scintillator phase

Located 2 km underground in SNOLAB, Sudbury, Canada, SNO+ is a large scale liquid scintillator experiment that primarily aims to search for neutrinoless double beta decay. Whilst SNO+ has light and radioactive calibration sources external to the inner volume, an internally deployed optical source is necessary for the full characterization of the detector model. A laser diffuser ball developed for SNO has previously demonstrated to be an effective optical calibration device for both SNO and SNO+ water phase. Since the introduction of liquid scintillator for SNO+, the material compatibility, cleanliness, and radiopurity requirements of any materials in contact with the internal medium have increased. Improving on the original SNO laserball design, a new laserball calibration device has been developed for the SNO+ scintillator phase with the goal of measuring the optical properties of the detector and performing routine PMT gain and timing calibrations. Simulations have been written to model the diffusion properties to optimise optical and temporal performance for calibration. Prototype laserballs have been built and characterised, demonstrating sub-ns timing resolution and a quasi-isotropic light distribution.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effects of detailed droplet heating models on turbulent sprays vaporization behavior

The effects of three different droplet heating models on the vaporization history and internal structure of turbulent liquid fuel sprays were investigated. The models considered were the infinite-diffusion, diffusion-limit, and effective-conductivity models. A numerical solution for the models was developed and implemented in the KIVA-II computer code. Low temperature and relatively high temperature numerical studies were conducted. The low temperature computations were compared with existing experimental data. The comparisons showed that while the infinite-diffusion and diffusion-limit models respectively overpredicted and underpredicted the fuel vapor peak concentration and distribution in the combustor, the effective-conductivity model gave results that were in better agreement with measurements. A limited study for the high temperature case was performed due to lack of experimental data and predictions using the three models were compared with each other.

Mawid, M. A.↗

Hierarchical model reduction driven by a proper orthogonal decomposition for parametrized advection-diffusion-reaction problems

This work combines the Hierarchical Model (HiMod) reduction technique with a standard Proper Orthogonal Decomposition (POD) to solve parametrized partial differential equations for the modeling of advection-diffusion-reaction phenomena in elongated domains (e.g., pipes). This combination leads to what we define as HiPOD model reduction, which merges the reliability of HiMod reduction with the computational efficiency of POD. Two HiPOD techniques are presented and assessed by an extensive numerical verification.

97 MATHEMATICS AND COMPUTING↗