Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional 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 901 records · Page 50

Developing the Science Basis for Understanding Polymer Encapsulant Degradation Mechanisms: DuraMAT 2.0 Final Project Report

Polymeric encapsulants are essential materials in photovoltaic modules, protecting sensitive electronics from the environment while providing mechanical integrity to the multilayered assembly. However, these polymeric materials are susceptible to degradation processes driven by the ingress of environmental species, ultraviolet radiation, thermal stresses, and mechanical loading. In this study, we employ a combined atomistic simulation and accelerated aging experimental approach to study the molecular-scale mechanisms of encapsulant degradation. Classical molecular dynamics simulations quantify the diffusion of environmental and degradation species through the polymer matrix, producing composition-specific diffusion coefficients. Reactive simulations characterize activation energy barriers and reaction rate constants for key chemical pathways. In parallel, thermal-desorption analyses coupled with mass spectrometry monitor the emergence and concentration profiles of degradation products under controlled stressor conditions. By integrating simulation and experiment, we establish quantitative correlations between polymer composition, species diffusivity, and chemical reactivity. We anticipate that these relations and quantitative values could serve as high-fidelity inputs to reaction-diffusion models, enabling physics-informed lifetime predictions and guiding the design of more durable encapsulant materials for solar energy applications.

36 MATERIALS SCIENCE↗

An efficient approach for treating composition-dependent diffusion within organic particles

Mounting evidence demonstrates that under certain conditions the rate of component partitioning between the gas and particle phase in atmospheric organic aerosol is limited by particle-phase diffusion. To date, however, particle-phase diffusion has not been incorporated into regional atmospheric models. An analytical rather than numerical solution to diffusion through organic particulate matter is desirable because of its comparatively small computational expense in regional models. Current analytical models assume diffusion to be independent of composition and therefore use a constant diffusion coefficient. To realistically model diffusion, however, it should be composition-dependent (e.g. due to the partitioning of components that plasticise, vitrify or solidify). This study assesses the modelling capability of an analytical solution to diffusion corrected to account for composition dependence against a numerical solution. Results show reasonable agreement when the gas-phase saturation ratio of a partitioning component is constant and particle-phase diffusion limits partitioning rate (<10% discrepancy in estimated radius change). However, when the saturation ratio of the partitioning component varies, a generally applicable correction cannot be found, indicating that existing methodologies are incapable of deriving a general solution. Until such time as a general solution is found, caution should be given to sensitivity studies that assume constant diffusivity. Furthermore, the correction was implemented in the polydisperse, multi-process Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) and is used to illustrate how the evolution of number size distribution may be accelerated by condensation of a plasticising component onto viscous organic particles.

54 ENVIRONMENTAL SCIENCES↗

Diffusion approximation for large absorption in radiative transfer

A new diffusion model is developed for radiative transfer in particulate media. It includes the effects of higher Legendre moments while avoiding the mathematical complexities of solving multiple coupled differential moment equations and satisfying higher-order boundary conditions. The method accurately extrapolates the conventional Eddington approximation to problems involving large absorption. Although the simplifying assumptions limit the model to nonbeam incidence and spatial homogeneity of the scattering material, they are nonrestrictive with regard to the size and shape of the medium, the character of background reflections, and the type of phase function.

Meador, W. E.↗

Coupled low-energy - ring current plasma diffusion in the Jovian magnetosphere

The outwardly diffusing Iogenic plasma and the simultaneously inwardly diffusing ring current plasma in the Jovian magnetosphere are described using a coupled diffusion model which incorporates the effects of the pressure gradient of the ring current into the cross-L diffusion coefficient. The coupled diffusion coefficient is derived by calculating the total energy available to drive the diffusion process. The condition is imposed that the diffusion coefficient takes on a local minimum value at some point in the region L = 7-8, at which point the gradient of the Io plasma density is specified as ramp value given by Siscoe et al. (1981). The hypothesis that the pressure gradient of the ring current causes the diminution of radial plasma transport is tested, and solution profiles for the Iogenic and ring current plasma densities are obtained which imply that the Io plasma ramp is caused by a high-density, low-energy component of the ring current hitherto unobserved directly.

Summers, D.↗

DSMC calculations for the double ellipse

The direct simulation Monte Carlo (DSMC) method involves the simultaneous computation of the trajectories of thousands of simulated molecules in simulated physical space. Rarefied flow about the double ellipse for test case 6.4.1 has been calculated with the DSMC method of Bird. The gas is assumed to be nonreacting nitrogen flowing at a 30 degree incidence with respect to the body axis, and for the surface boundary conditions, the wall is assumed to be diffuse with full thermal accommodation and at a constant wall temperature of 620 K. A parametric study is presented that considers the effect of variations of computational domain, gas model, cell size, and freestream density on surface quantities.

Moss, James N.↗

Three-dimensional modelling of horizontal chemical vapor deposition. I - MOCVD at atmospheric pressure

A systematic numerical study of the MOCVD of GaAs from trimethylgallium and arsine in hydrogen or nitrogen carrier gas at atmospheric pressure is reported. Three-dimensional effects are explored for CVD reactors with large and small cross-sectional aspect ratios, and the effects on growth rate uniformity of tilting the susceptor are investigated for various input flow rates. It is found that, for light carrier gases, thermal diffusion must be included in the model. Buoyancy-driven three-dimensional flow effects can greatly influence the growth rate distribution through the reactor. The importance of the proper design of the lateral thermal boundary conditions for obtaining layers of uniform thickness is emphasized.

Ouazzani, Jalil↗

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗

Characteristics of Non-Premixed Turbulent Flames in Microgravity

This project is concerned with the characteristics of turbulent hydrocarbon (primarily propane) gas-jet diffusion flames in microgravity. A microgravity environment provides the opportunity to study the structure of turbulent diffusion flames under momentum-dominated conditions (large Froude number) at moderate Reynolds number which is a combination not achievable in normal gravity. This paper summarizes progress made since the last workshop. Primarily, the features of flame radiation from microgravity turbulent jet diffusion flames in a reduced gravity environment are described. Tests were conducted for non-premixed, nitrogen diluted propane flames burning in quiescent air in the NASA Glenn 5.18 Second Zero Gravity Facility. Measured flame radiation from wedge-shaped, axial slices of the flame are compared for microgravity and normal gravity flames. Results from numerical computations of the flame using a k-e model for the turbulence are also presented to show the effects of flame radiation on the thermal field. Flame radiation is an important quantity that is impacted by buoyancy as has been shown in previous studies by the authors and also by Urban et al. It was found that jet diffusion flames burning under microgravity conditions have significantly higher radiative loss (about five to seven times higher) compared to their normal gravity counterparts because of larger flame size in microgravity and larger convective heat loss fraction from the flame in normal gravity. These studies, however, were confined to laminar flames. For the case of turbulent flames, the flame radiation is a function of time and both the time-averaged and time-dependent components are of interest. In this paper, attention is focused primarily on the time-averaged level of the radiation but the turbulent structure of the flame is also assessed from considerations of the radiation power spectra.

Hegde, U.↗

A quantitative model of charge injection by ruthenium chromophores connecting femtosecond to continuous irradiance conditions

A kinetic framework for the ultrafast photophysics of tris(2,2-bipyridine)ruthenium(II) phosphonated and methyl-phosphonated derivatives is used as a basis for modeling charge injection by ruthenium dyes into a semiconductor substrate. By including the effects of light scattering, dye diffusion, and adsorption kinetics during sample preparation and the optical response of oxidized dyes, quantitative agreement with multiple transient absorption datasets is achieved on timescales spanning femtoseconds to nanoseconds. In particular, quantitative agreement with important spectroscopic handles—the decay of an excited state absorption signal component associated with charge injection in the UV region of the spectrum and the dynamical redshift of a ~500 nm isosbestic point—validates our kinetic model. Pseudo-first-order rate coefficients for charge injection are estimated in this work, with an order of magnitude ranging from 10 11 to 10 12 s –1 . The model makes the minimalist assumption that all excited states of a particular dye have the same charge injection coefficient, an assumption that would benefit from additional theoretical and experimental exploration. Here, we have adapted this kinetic model to predict charge injection under continuous solar irradiation and find that as many as 68 electron transfer events per dye per second take place, significantly more than prior estimates in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport dynamics of self-consistent, near-marginal drift-wave turbulence. I. Investigation of the ability of external flows to tune the non-diffusive dynamics

The reduction of turbulent transport across sheared flow regions has been known for a long time in magnetically confined toroidal plasmas. However, details of the dynamics are still unclear, in particular, in what refers to the changes caused by the flow on the nature of radial transport itself. In Paper II, we have shown in a simplified model of drift wave turbulence that, when the background profile is allowed to evolve self-consistently with fluctuations, a variety of transport regimes ranging from superdiffusive to subdiffusive open up depending on the properties of the underlying turbulence. In this study, we show that externally applied sheared flows can, under the proper conditions, cause the transport dynamics to be diffusive or subdiffusive

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale model of metal alloy oxidation at grain boundaries

High temperature intergranular oxidation and corrosion of metal alloys is one of the primary causes of materials degradation in nuclear systems. In order to gain insights into grain boundary oxidation processes, a mesoscale metal alloy oxidation model is established by combining quantum Density Functional Theory (DFT) and mesoscopic Poisson-Nernst-Planck/classical DFT with predictions focused on Ni alloyed with either Cr or Al. Analysis of species and fluxes at steady-state conditions indicates that the oxidation process involves vacancy-mediated transport of Ni and the minor alloying element to the oxidation front and the formation of stable metal oxides. The simulations further demonstrate that the mechanism of oxidation for Ni-5Cr and Ni-4Al is qualitatively different. Intergranular oxidation of Ni-5Cr involves the selective oxidation of the minor element and not matrix Ni, due to slower diffusion of Ni relative to Cr in the alloy and due to the significantly smaller energy gain upon the formation of nickel oxide compared to that of Cr{sub 2}O{sub 3}. This essentially one-component oxidation process results in continuous oxide formation and a monotonic Cr vacancy distribution ahead of the oxidation front, peaking at alloy/oxide interface. In contrast, Ni and Al are both oxidized in Ni-4Al forming a mixed spinel NiAl{sub 2}O{sub 4}. Different diffusivities of Ni and Al give rise to a complex elemental distribution in the vicinity of the oxidation front. Slower diffusing Ni accumulates in the oxide and metal within 3 nm of the interface, while Al penetrates deeper into the oxide phase. Ni and Al are both depleted from the region 3–10 nm ahead of the oxidation front creating voids. The oxide microstructure is also different. Cr{sub 2}O{sub 3} has a plate-like structure with 1.2–1.7 nm wide pores running along the grain boundary, while NiAl{sub 2}O{sub 4} has 1.5 nm wide pores in the direction parallel to the grain boundary and 0.6 nm pores in the perpendicular direction providing an additional pathway for oxygen diffusion through the oxide. The proposed theoretical methodology provides a framework for modeling metal alloy oxidation processes from first principles and on the experimentally relevant length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations for solvent-based carbon capture. Part 2: Chemical absorption across a wetted wall column: Original Research Article: Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations

Part 1 of this paper presents a numerical model for non-reactive physical mass transfer across a wetted wall column (WWC). In Part 2, we improved the existing computational fluid dynamics (CFD) model to simulate chemical absorption occurring in a WWC as a bench-scale study of solvent-based carbon dioxide (CO 2 ) capture. In this study, to generate data for WWC model validation, CO 2 mass transfer across a monoethanolamine (MEA) solvent was first measured on a WWC experimental apparatus. The numerical model developed in this work can account for both chemical absorption and desorption of CO 2 in MEA. In addition, the overall mass transfer coefficient predicted using traditional/empirical correlations is conducted and compared with CFD prediction results for both steady and wavy falling films. A Bayesian statistical calibration algorithm is adopted to calibrate the reaction rate constants in chemical absorption/desorption of CO 2 across a falling film of MEA. The posterior distributions of the two transport properties, i.e., Henry's constant and gas diffusivity in the non-reacting nitrous oxide (N 2 O)/MEA system obtained from Part 1 of this study, serves as priors for the calibration of CO 2 reaction rate constants after using the N 2 O/CO 2 analogy method. Finally, the calibrated model can be used to predict the CO 2 mass transfer in a WWC for a wider range of operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the Role of Interfaces for Water Management in Platinum Group Metal-Free Electrodes in Polymer Electrolyte Fuel Cells

For this study, a systematic analysis, both experimental and model-assisted, has been performed over three main configurations of platinum group metal-free (PGM-free) electrodes in polymer electrolyte fuel cells (PEFCs): Catalyst-coated membrane CCM technology is being compared to the gas-diffusion electrode (GDE) method of electrode fabrication and juxtaposed to a hybrid/combined GDE-CCM method of membrane-electrode assembly (MEA) fabrication. The corresponding electrodes were evaluated for their electrochemical performance, modeled, and studied with in situ and operando X-ray computed tomography (X-ray CT). The study establishes that through-thickness inhomogeneities play the most important role in water withdrawal/water management and affect most significantly PGM-free PEFC performance. The catalyst integration technique results in formation of interfacial regions with increased porosity and surface roughness. These regions form critical interfaces de facto responsible for flooding type behavior of the PEFC as shown for a first time by operando X-ray CT. The computational model shows that the PEFC performance critically depends on liquid water formation and transport at cold and wet conditions.

25 ENERGY STORAGE↗

Physical and mechanical metallurgy of zirconium alloys for nuclear applications: a multi-scale computational study

In the post-Fukushima world, the stability of materials under extreme conditions is an important issue for the safety of nuclear reactors. Because the nuclear industry is going to continue using advanced zirconium cladding materials in the foreseeable future, it become critical to gain fundamental understanding of the several interconnected problems. First, what are the thermodynamic and kinetic factors affecting the oxidation and hydrogen pick-up by these materials at normal, off-normal conditions, and in long-term storage? Secondly, what protective coatings (if any) could be used in order to gain extremely valuable time at off-normal conditions, e.g., when temperature exceeds the critical value of 2200°F? Thirdly, the kinetics of oxidation of such protective coating or braiding needs to be quantified. Lastly, even if some degree of success is achieved along this path, it is absolutely critical to have automated inspection algorithms allowing identifying defects of cladding as soon as possible. This work strives to explore these interconnected factors from the most advanced computational perspective, utilizing such modern techniques as first-principles atomistic simulations, computational thermodynamics of materials, diffusion modeling, and the morphological algorithms of image processing for defect identification. Consequently, it consists of the four parts dealing with these four problem areas preceded by the introduction and formulation of the studied problems. In the 1st part an effort was made to employ computational thermodynamics and ab initio calculations to shed light upon the different stages of oxidation of ziraloys (2 and 4), the role of microstructure optimization in increasing their thermal stability, and the process of hydrogen pick-up, both in normal working conditions and in long-term storage. The 2nd part deals with the need to understand the influence and respective roles of the two different plasticity mechanisms in Zr nuclear alloys: twinning (at low T) and crystallographic slip (higher T’s). For that goal, a description of the advanced plasticity model is outlined featuring the non-associated flow rule in hcp materials including Zr. The 3rd part describes the kinetic theory of oxidation of the several materials considered to be perspective coating materials for Zr alloys: SiC and ZrSiO 4 . In the 4th part novel and advanced projectional algorithms for defect identification in zircaloy coatings are described. In so doing, the author capitalized on some 12 years of his applied industrial research in this area. Our conclusions and recommendations are presented in the 5th part of this work, along with the list of used literature and the scripts for atomistic, thermodynamic, kinetic, and morphological computations.

36 MATERIALS SCIENCE↗

Understanding and controlling water-organic co-transport in amorphous microporous materials

The transport of molecules in microporous material is a significant and active area of research in separation applications. The movement of vapor/liquid molecules in a high-loading condition in microporous spaces is especially challenging to interrogate. This project aims to understand and control the transport of complex water and organic solvent mixtures in varying structures of microporous carbon molecular sieve (CMS) and activated carbon membranes. Modeling of transport behavior in such carbonaceous samples using computational modeling is difficult due to their amorphous structure. Therefore, this research intends to study the transport mechanism of water-organic mixtures by experimentally revealing fundamental transport properties, such as guest sorption amounts as well as diffusion and permeation rates. Considering the different size and guest-host affinity of water and organic solvent molecules, the objective of this research is to understand the structural conditions within the ultramicropores and micropores that generate different types of molecular transport and selection mechanisms within such complex systems. Carbonaceous materials are developed using tailored pyrolysis techniques to pyrolyze polymeric precursors, including polyvinylidene chloride (PVDC),polyvinylidene fluoride (PVDF), polymer of intrinsic microporosity (PIM) -1, and fully aromatic polyamide. The initial stage of the research will focus on the microscopic diffusion and sorption studies of pure component water, xylene isomers, and n,n- dimethylformamide within various carbons. Additional structural investigation on CMS, such as gas physisorption and neutron scattering studies, will be executed to gain deeper insight into these structure-transport relationships. Scanning electron microscopy and X-ray photoelectron microscopy are also used to further characterize these materials. The research will progress to probe mutual diffusion of complex water-organic mixtures within various CMS microstructures. Competitive sorption and permeation studies of the binary mixtures in different microstructural CMS will be followed. The obtained binary mixture transport parameters were used to predict binary water-organic mixtures transport.

36 MATERIALS SCIENCE↗

Ultrafast Metal Electrodeposition Revealed by In Situ Optical Imaging and Theoretical Modeling towards Fast–Charging Zn Battery Chemistry

Metallic Zn is a preferred anode material for rechargeable aqueous batteries towards a smart grid and renewable energy storage. Importantly, understanding how the metal nucleates and grows at the aqueous Zn anode is a critical and challenging step to achieve full reversibility of Zn battery chemistry, especially under fast-charging conditions. Here, by combining in situ optical imaging and theoretical modeling, we uncover the critical parameters governing the electrodeposition stability of the metallic Zn electrode, that is, the competition among crystallographic thermodynamics, kinetics, and Zn 2+ -ion diffusion. Moreover, steady-state Zn metal plating/ stripping with Coulombic efficiency above 99 % is achieved at 10-100 mA cm -2 in a reasonably high concentration (3 M) ZnSO 4 electrolyte. Significantly, a long-term cycling-stable Zn metal electrode is realized with a depth of discharge of 66.7% under 50 mA cm -2 in both Zn || Zn symmetrical cells and MnO 2 || Zn full cells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of magnetic field line random walk in noisy reduced magnetohydrodynamic turbulence

When a magnetic field consists of a mean part and fluctuations, the stochastic wandering of its field lines is often treated as a diffusive process. Under suitable conditions, a stable value is found for the mean square transverse displacement per unit parallel displacement relative to the mean field. Here, we compute the associated field line diffusion coefficient for a highly anisotropic 'noisy' reduced magnetohydrodynamic model of the magnetic field, which is useful in describing low frequency turbulence in the presence of a strong applied DC mean magnetic field, as may be found, for example, in the solar corona, or in certain laboratory devices. Our approach is nonperturbative, based on Corrsin's independence hypothesis, and makes use of recent advances in understanding factors that control decorrelation over a range of parameters described by the Kubo number. Both Bohm and quasilinear regimes are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Flexible machine-learning interatomic potential for simulating structural disordering behavior of Li 7 La 3 Zr 2 O 12 solid electrolytes

Batteries based on solid-state electrolytes, including Li 7 La 3 Zr 2 O 12 (LLZO), promise improved safety and increased energy density; however, atomic disorder at grain boundaries and phase boundaries can severely deteriorate their performance. Machine-learning (ML) interatomic potentials offer a uniquely compelling solution for simulating chemical processes, rare events, and phase transitions associated with these complex interfaces by mixing high scalability with quantum-level accuracy, provided that they can be trained to properly address atomic disorder. To this end, we report the construction and validation of an ML potential that is specifically designed to simulate crystalline, disordered, and amorphous LLZO systems across a wide range of conditions. The ML model is based on a neural network algorithm and is trained using ab initio data. Performance tests prove that the developed ML potential can predict accurate structural and vibrational characteristics, elastic properties, and Li diffusivity of LLZO comparable to ab initio simulations. As a demonstration of its applicability to larger systems, we show that the potential can correctly capture grain boundary effects on diffusivity, as well as the thermal transition behavior of LLZO. Here these examples show that the ML potential enables simulations of transitions between well-defined and disordered structures with quantum-level accuracy at speeds thousands of times faster than ab initio methods.

25 ENERGY STORAGE↗