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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Mechanism of charge recombination in meso-structured organic-inorganic hybrid perovskite solar cells: A macroscopic perspective

In the currently popular organic-inorganic hybrid perovskite solar cells, the slowness of the charge recombination processes is found to be a key factor for contributing to their high efficiencies and high open circuit voltages, but the underlying recombination mechanism remains unclear. In this work, we investigate the bimolecular recombination (BR) and the trap-assisted monomolecular recombination (MR) in meso-structured perovskite solar cells under steady state working condition, and try to reveal their roles on determining the device performance. Some interfacial effects such as the injection barriers at the selective contacts are examined as well. Based on the macroscopic device modeling, the recombination resistance-voltage (R{sub rec}−V) and the current density-voltage (J–V) curves are calculated to characterize the recombination mechanism and describe the device performance, respectively. Through comparison with the impedance spectroscopy extracted R{sub rec} data, it is found that under the typical BR reduction factor and deep trap densities observed in experiments, the MR dominates the charge recombination in the low voltage regime, while the BR dominates in the high voltage regime. The short circuit current and the fill factor could be reduced by the significant MR but the open circuit voltage is generally determined by the BR. The different electron injection barriers at the contact can change the BR rate and induce different patterns for the R{sub rec}–V characteristics. For the perovskites of increased band gaps, the R{sub rec}'s are significantly enhanced, corresponding to the high open circuit voltages. Finally, it is revealed that the reduced effective charge mobility due to the transport in electron and hole transporting material makes the R{sub rec} decrease slowly with the increasing voltage, which leads to increased open circuit voltage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2017

Development of Computational Tools for Modeling Thermal and Radiation Effects on Grain Boundary Segregation and Precipitation in Fe-Cr-Ni-based Alloys

This work aims at developing computational tools for modeling thermal and radiation effects on solute segregation at grain boundaries (GBs) and precipitation. This report described two major efforts. One is the development of computational tools on integrated modeling of thermal equilibrium segregation (TES) and radiation-induced segregation (RIS), from which synergistic effects of thermal and radiation, pre-existing GB segregation have been taken into consideration. This integrated modeling was used in describing the Cr and Ni segregation in the Fe-Cr-Ni alloys. The other effort is thermodynamic modeling on the Fe-Cr-Ni-Mo system which includes the major alloying elements in the investigated alloys in the Advanced Radiation Resistant Materials (ARRM) program. Through thermodynamic calculation, we provide baseline thermodynamic stability of the hardening phase Ni2(Cr,Mo) in selected Ni-based super alloys, and contribute knowledge on mechanistic understanding on the formation of Ni2(Cr,Mo) in the irradiated materials. The major outcomes from this work are listed in the following: 1) Under the simultaneous thermal and irradiation conditions, radiation-induced segregation played a dominant role in the GB segregation. The pre-existing GB segregation only affects the subsequent radiation-induced segregation in the short time. For the same element, the segregation tendency of Cr and Ni due to TES is opposite to it from RIS. The opposite tendency can lead to the formation of W-shape profile. These findings are consistent with literature observation of the transitory W-shape profile. 2) While TES only affects the distance of one or two atomic layers from GBs, the RIS can affect a broader distance from GB. Therefore, the W-shape due to pre-existing GB segregation is much narrower than that due to composition gradient formed during the transient state. Considering the measurement resolution of Auger or STEM analysis, the segregation tendency due to RIS should play a dominant role in the measured values. However, The GB segregation due to pre-existing GB segregation may affect the chemical potential of element at GB, and subsequently the corrosion resistance. 3) Based on the newly developed thermodynamic database of Fe-Cr-Ni-Mo, we predicted the Ni2(Cr,Mo) as a thermodynamically stable phase in all investigated low Fe-content Ni-based alloys. The calculated phase amount decreases with the increasing Fe content, being consistent with that observed in the irradiated materials. 4) The formation of the Ni2(Cr,Mo) phase in irradiated materials is due to irradiation enhanced diffusion. The calculated equilibrium Ni2(Cr,Mo) amount is more than that observed in the irradiated materials, suggesting that the amount of Ni2(Cr,Mo) is likely to increase more with further irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SISGR - Hydrogen Caged in Carbon-Exploration of Novel Carbon-Hydrogen Interactions

Hydrogen trapped in a carbon cage, captured through repulsive interactions, is a novel concept in hydrogen storage. Trapping hydrogen via repulsive interactions borrows an idea from macroscale hydrogen storage (i.e. compressed gas storage tanks) and reapplies these concepts on the nanoscale in specially designed molecular containers. Under extreme conditions of pressure, hydrogen solubility in carbon materials is expected to increase and carbon is expected to restructure to minimize volume via a mixed sp2/sp3 hydrogenated state. Thermodynamics dictate that pre-formed C-H structures will rearrange with increased pressure, yet the final carbon-hydrogen interactions may be dependent upon the mechanism by which hydrogen is introduced. Gas “trapping” is meant to denote gas present in a solid in a high density, adsorbed-like state, when the external pressure is much less than that necessary to provide a comparable fluid density. Trapping thus denotes a kinetically metastable state rather than thermodynamic equilibrium. This project probed mechanochemical means to polymerize select hydrocarbons in the presence of gases, in an attempt to form localized carbon cages that trap gases via repulsive interactions. Aromatic, polyaromatic, and hydroaromatic molecules expected to undergo cyclo-addition reactions were polymerized at high (~GPa) pressures to form extended hydrogenated amorphous carbon networks. Notably, aromatics with a pre-existing internal free volume (such as Triptycene) appeared to retain an internal porosity upon application of pressure. However, a high photoluminescence background after polymerization precluded in situ identification of trapped gases. No spectroscopic evidence was found after depressurization that would be indicative of pockets of trapped gases in a localized high-pressure environment. Control studies suggested this measurement may be insensitive to gases at low pressure. Similarly, no spectral fingerprint was found for gas-imbued spherical carbon nanoshells, even after chemical “capping” of the gas-imbued nanoshells to limit gas diffusivity. Subsequently, spectral probes of gas vibrational modes adsorbed in various carbon nanostructures (including activated carbons, single-wall carbon nanotubes, polymers of intrinsic microporosity (PIMs), and UV-irradiated PIMs with decreased pore size) were found only at high pressure. The vibrational mode of the adsorbed film became perturbed in high density films, and the perturbation was sensitive to surface functional groups, pore size, and pore dimension. Experimental results were corroborated with first-principle modeling using density functional theory. Development of semi-empirical correlations that relate the spectral features to pore dimension, geometry, and chemical potential of the adsorbed film are on-going.

08 HYDROGEN↗
DOE OSTI2017

First-Principles-Based Computational Study on Nucleation and Growth Mechanisms of U on Mo(110) Surface Solvated in an Eutectic LiCl-KCl Molten Salt

In nuclear engineering and industry electrochemical methods have been intensively employed in pyro-processing technology to recycle used nuclear fuels for substantial alleviation of the concern with the safe treatment of high-level radioactive nuclear wastes. Pyro-processing is based on sequential electrochemical reactions at high temperature (typically above 700 K at ambient pressure) to enhance the efficiency for decoupling nuclear fuels from unfavorable fission fragments. Using a molten salt in electrorefining process, at which actinides can be selectively deposited in solid cathode from ionic states, gives several advantages over an aqueous solution: higher thermal conductivity (it is liquid phase at high temperature, molten salt plays a role as medium at high temperature over 773 K in electrochemical reaction, lower kinetic barrier in ionic diffusion mobility and better stability against undesired side reactions with radioactive materials. Therefore, the volume of high-level nuclear wastes to be, otherwise, disposed into secured sites is dramatically reduced. Most of all it conforms to nuclear non-proliferation treaty. Electrochemical reactions in a molten salt medium are, however, much more complicated than in aqueous solution rendering the accurate understanding of thermodynamic and kinetic properties formidably challenging. It is represented by extensive researches in experiments and computations to identify underlying mechanisms on how U and other fission fragments behave in the electrorefining system. Regardless of the researches there is a general agreement that the initial stage of the nucleation and growth process of the U deposit is critical for determining the properties of final product and the electrorefining efficiency. In spite of the focused studies it still remains unclear. Theoretical and computational researches have been often carried out with various empirical parameter for interatomic potential energy, which are unknown variables until numerical fits to experimental measurements or just assumption is made. Thus, empirical simulations sometimes predicted outcomes are inconsistent depending on materials of interest with experimental measurements. This may be largely because the accuracy of an empirical interatomic potential energy models critically depends on the input parameters and thus, can be limited in transferability of the results. Furthermore, it may be devoid of fundamental physics or chemistry. This situation requires special caution when an electrorefining system is studied to recycle spent nuclear fuels because of the extreme difficulty in reliable acquisition of thermodynamic and kinetic data from direct experimental measurements. In this presentation, we extensively utilized first principles density functional theory (DFT) calculations and ab-initio molecular dynamic (AIMD) simulations to figure out underlying mechanisms of the initial stage of nucleation and growth behaviors of U on Mo(110) surface in a eutectic LiCl-KCl molten salt at T = 773 K. To rigorously investigate the effect of the molten salt electrolyte we compared the results with those calculated under vacuum condition.

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

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2015

Energy dynamics and current sheet structure in fluid and kinetic simulations of decaying magnetohydrodynamic turbulence

Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfvén waves, which interact and rapidly generate strong MHD turbulence. The total energy is conserved and the rate of turbulent energy decay is very similar in both codes, although the fluid code has numerical dissipation, whereas the kinetic code has kinetic dissipation. The inertial range power spectrum index is similar in both the codes. The fluid code shows a perpendicular wavenumber spectral slope of k{sub ⊥}{sup −1.3}. The kinetic code shows a spectral slope of k{sub ⊥}{sup −1.5} for smaller simulation domain, and k{sub ⊥}{sup −1.3} for larger domain. We estimate that collisionless damping mechanisms in the kinetic code can account for the dissipation of the observed nonlinear energy cascade. Current sheets are geometrically characterized. Their lengths and widths are in good agreement between the two codes. The length scales linearly with the driving scale of the turbulence. In the fluid code, their thickness is determined by the grid resolution as there is no explicit diffusivity. In the kinetic code, their thickness is very close to the skin-depth, irrespective of the grid resolution. This work shows that kinetic codes can reproduce the MHD inertial range dynamics at large scales, while at the same time capturing important kinetic physics at small scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The high-β{sub N} hybrid scenario for ITER and FNSF steady-state missions

New experiments on DIII-D have demonstrated the steady-state potential of the hybrid scenario, with 1 MA of plasma current driven fully non-inductively and β{sub N} up to 3.7 sustained for ∼3 s (∼1.5 current diffusion time, τ{sub R}, in DIII-D), providing the basis for an attractive option for steady-state operation in ITER and FNSF. Excellent confinement is achieved (H{sub 98y2} ∼ 1.6) without performance limiting tearing modes. The hybrid regime overcomes the need for off-axis current drive efficiency, taking advantage of poloidal magnetic flux pumping that is believed to be the result of a saturated 3/2 tearing mode. This allows for efficient current drive close to the axis, without deleterious sawtooth instabilities. In these experiments, the edge surface loop voltage is driven down to zero for >1 τ{sub R} when the poloidal β is increased above 1.9 at a plasma current of 1.0 MA and the ECH power is increased to 3.2 MW. Stationary operation of hybrid plasmas with all on-axis current drive is sustained at pressures slightly above the ideal no-wall limit, while the calculated ideal with-wall MHD limit is β{sub N} ∼ 4–4.5. Off-axis Neutral Beam Injection (NBI) power has been used to broaden the pressure and current profiles in this scenario, seeking to take advantage of higher predicted kink stability limits and lower values of the tearing stability index Δ′, as calculated by the DCON and PEST3 codes. Results based on measured profiles predict ideal limits at β{sub N} > 4.5, 10% higher than the cases with on-axis NBI. A 0-D model, based on the present confinement, β{sub N} and shape values of the DIII-D hybrid scenario, shows that these plasmas are consistent with the ITER 9 MA, Q = 5 mission and the FNSF 6.7 MA scenario with Q = 3.5. With collisionality and edge safety factor values comparable to those envisioned for ITER and FNSF, the high-β{sub N} hybrid represents an attractive high performance option for the steady-state missions of these devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗
DOE OSTI2016

A GLOBAL GALACTIC DYNAMO WITH A CORONA CONSTRAINED BY RELATIVE HELICITY

We present a model for a global axisymmetric turbulent dynamo operating in a galaxy with a corona that treats the parameters of turbulence driven by supernovae and by magneto-rotational instability under a common formalism. The nonlinear quenching of the dynamo is alleviated by the inclusion of small-scale advective and diffusive magnetic helicity fluxes, which allow the gauge-invariant magnetic helicity to be transferred outside the disk and consequently to build up a corona during the course of dynamo action. The time-dependent dynamo equations are expressed in a separable form and solved through an eigenvector expansion constructed using the steady-state solutions of the dynamo equation. The parametric evolution of the dynamo solution allows us to estimate the final structure of the global magnetic field and the saturated value of the turbulence parameter α{sub m}, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with α-quenching. We then solve these equations simultaneously to study the saturation of the large-scale magnetic field, its dependence on the small-scale magnetic helicity fluxes, and the corresponding evolution of the force-free field in the corona. The quadrupolar large-scale magnetic field in the disk is found to reach equipartition strength within a timescale of 1 Gyr. The large-scale magnetic field in the corona obtained is much weaker than the field inside the disk and has only a weak impact on the dynamo operation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Metastability and reliability of CdTe solar cells

Thin-film modules of all technologies often suffer from performance degradation over time. Some of the performance changes are reversible and some are not, which makes deployment, testing, and energy-yield prediction more challenging. Manufacturers devote significant empirical efforts to study these phenomena and to improve semiconductor device stability. Still, understanding the underlying reasons of these instabilities remains clouded due to the lack of ability to characterize materials at atomistic levels and the lack of interpretation from the most fundamental material science. The most commonly alleged causes of metastability in CdTe devices, such as ‘migration of Cu’, have been investigated rigorously over the past fifteen years. Still, the discussion often ended prematurely with stating observed correlations between stress conditions and changes in atomic profiles of impurities or CV doping concentration. Multiple hypotheses suggesting degradation of CdTe solar cell devices due to interaction and evolution of point defects and complexes were proposed, and none of them received strong theoretical or experimental confirmation. We believe it should be noted that atomic impurity profiles in CdTe provide very little intelligence on active doping concentrations. The same elements could form different energy states, which could be either donors or acceptors, depending on their position in crystalline lattice. Defects interact with other extrinsic and intrinsic defects; for example, changing the state of an impurity from an interstitial donor to a substitutional acceptor often is accompanied by generation of a compensating intrinsic interstitial donor defect. Moreover, all defects, intrinsic and extrinsic, interact with the electrical potential and free carriers so that charged defects may drift in the electric field and the local electrical potential affects the formation energy of the point defects. Such complexity of interactions in CdTe makes understanding of temporal changes in device performance even more challenging and a closed solution that can treat the entire system and its interactions is required.

14 SOLAR ENERGY↗
DOE OSTI2017

Final Report: In-Situ TEM Observations of Degradation Mechanisms in Next-Generation High-Energy Density Lithium-Ion Battery Systems

This project originally sought to characterize nanoscale processes associated with the degradation of next-generation high energy density lithium-ion battery electrodes via in-situ transmission electron microscopy (TEM). The research also focused on developing the experimental techniques necessary for investigating electrochemical systems by in-situ transmission electron microscopy and providing a framework for understanding and quantifying electron beam effects that influence experimental results. In line with these goals, this project demonstrated and published the first in-situ TEM cycling of a Li-ion battery electrode in commercially relevant non-aqueous electrolyte: a technique subsequently adopted by a number of researchers in the field. Experiments performed on Sn anodes, demonstrated porosity evolution during Li dealloying as a mechanism for strain accommodation during the first cycle, and observed partially reversible crystalline to amorphous phase transitions associated with Li insertion and extraction in subsequent cycles. Complementary study of in-situ dealloying in a simple model, Au-Cu, was also performed. We also published the earliest quantitative measurements of electron beam effects in aqueous and non-aqueous liquids. Specifically, we utilized a model system to quantify the fraction of incident electron beam that contributed to charge transfer mediated reactions in non-aqueous liquid, since such liquids are primarily dose sensitive. We also quantified the evolution of model aqueous chemistries as a function of dose rate and chemical potential to better understand the reaction mechanisms active in liquids exposed to electron beams. The original project goals included the use of controlled environments to simulate different battery operational conditions. To aid such experiments, we developed new experimental apparatuses and techniques to perform in-situ liquid cell TEM under controlled temperature, including heating and cooling. We demonstrated this capability in aqueous systems; specifically understanding phase selection and solid-liquid interface interactions with particles during water crystallization and hydrothermal precipitation kinetics of ZnO. We also published the first demonstration of in-situ TEM based high-field (plasma) experiments and in-situ TEM of aqueous photochemistry, which demonstrated an previously unrecognized early growth mode transition during Au sputtering and delocalization of gas bubble nucleation during water splitting, respectively. In addition to the scientific results, this portion of the project greatly expanded the range of systems that can be characterized in realistic environments at the nanoscale via in-situ TEM. Unfortunately, commercial battery electrolytes are quite sensitive to electron dose. This limits the time of acquisition and/or the spatial resolution of the experiments. The same problem has been encountered in all subsequent publications in the field. The major problem, relative to the project original goals, is that long term cycling is challenging to characterize and the desired high resolution or analytical capabilities are limited. In line with our original goals to characterize battery reaction and degradation mechanisms under controlled electrochemical cycling, we developed a completely new type of Li-ion electrochemical cell that has been useful for both structural and analytical characterization. To date, we have demonstrated this platform for in-situ cycling of Li-ion batteries during scanning electron microscopy (SEM), Auger electron spectroscopy (AES), and x-ray photoelectron spectroscopy (XPS). In-situ AES has been found to be a quite powerful tool to characterize the chemical evolution of Li-ion batteries cycled in-situ and much of our effort in the final project period has focused on charactering Li-ion cathodes using the technique. We utilized in-situ SEM to propose a new mechanism for Li whisker/dendrite formation and show the importance of diffusion mediated stress relaxation in fracture of Sn anode particles. A combination of in-situ XPS and AES was used to clarify the CuO anode reaction mechanism. In the final project period, we have gained new insights into surface reactions occurring on cathodes at high voltage using in-situ AES. Specifically, we have discovered that many of the SEI formation reactions typically attributed to electrolyte decomposition are actually more fundamental to the electrode material surface thermodynamics in the presence of an arbitrary carbon source. These insights greatly affect how we understand and control surface degradation induced capacity fade at cathodes. Finally, we have made some effort to apply this novel in-situ platform to TEM characterization and believe this is a fruitful avenue for future research.

25 ENERGY STORAGE↗
DOE OSTI2019

Revisiting the Analytical Solutions of Heat Transport in Fractured Reservoirs Using a Generalized Multirate Memory Function

Numerous analytical solutions have been developed for modeling thermal perturbations to underground formations caused by deep-well injection of fluids. Each solution has been derived for a specific boundary value problem and a simplified flow network with one set of parallel fractures. In this paper, new generalized solutions G*( x, s ) are developed using (existing) global transfer functions G$^*_0(x,s)$ and a new memory function g*( s ), where x and s are the space and Laplace variable. The memory function represents the solutions of conductive heat exchange between fractures and matrix blocks and between fractured aquifers and unfractured aquitards. The memory function is developed to account for multirate exchange induced by different shapes, sizes, properties, and volumetric fractions of matrix blocks bounded by multiple sets of orthogonal fractures with different spacing. The global transfer functions represent the fundamental solutions to convective, convective-conductive, and convective-dispersive heat transport in fractures (or aquifers) without exchange and are available for various (1-D linear, 1-D radial, 2-D dipole, and single-well injection-withdrawal) flow fields. The new solutions with exchange are developed using $G^*(x,s) = B^*(s)G^*_0{x, s[1 + Θg^*(s)]}$, thereby greatly simplifying solution development in a novel way, where Θ and B*( s ) are a fracture-matrix scaling factor and the boundary condition function. As a result, the new solutions are applied to several example problems, showing that heat transport in fractured aquifers is significantly impacted by (1) thermal dispersion in fractures that is rarely considered, (2) multirate heat exchange with a wide range of size and anisotropy of rectangular matrix blocks, and (3) heat exchange between aquifers and aquitards.

58 GEOSCIENCES↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2013

The Low-Recycling Lithium Boundary and Implications for Plasma Transport

Pumping of incident hydrogen and impurity ions by lithium enables control of the particle inventory and fueling profile in magnetic-confined plasmas, and may raise the plasma temperature near the wall. As a result, the density gradient is expected to contribute substantially to the free-energy, affecting particle and thermal transport from micro-turbulence which is typically the dominant transport mechanism in high-temperature fusion experiments. Transport in gyrokinetic simulations of density-gradient-dominated profiles is characterized by a small linear critical gradient, large particle flux, and preferential diffusion of cold particles. As a result, the heat flux is below 5/2 or even 3/2 times the particle flux, usually assumed to be the minimum for convection. While surprising, this result is consistent with increasing entropy. Coupled TEM-ITG (ion-temperature- gradient) simulations using ηe = ηi find η = ∇T /∇n∼0.8 maximizes the linear critical pressure gradient, which suggests that experiments operating near marginal ITG stability with larger η would increase the linear critical pressure gradient by transferring free-energy from the temperature gradient to the density gradient. Simulations were performed with profiles predicted for the Lithium Tokamak Experiment (LTX) if ion thermal transport was neoclassical, while electron thermal transport and particle transport were a fixed ratio above the neoclassical level. A robust TEM instability was found for the outer half radius, while the ITG was found to be driven unstable as well during gas puff fueling. This suggests that TEM transport will be an important transport mechanism in high-temperature low-recycling fusion experiments, and in the absence of stabilizing mechanisms, may dominate over neoclassical transport. A diagnostic suite has been developed to measure hydrogen and impurity emission in LTX in order to determine the lower bound on recycling that can be achieved in a small tokamak using solid lithium coatings, assess its dependence on the operating condition of the lithium surface, and evaluate its impact on the discharge. Coatings on the close-fitting stainless-steel substrate produce a significant reduction in recyling, so that the effective particle confinement times are as low as 1 ms. Measurements of particle inventory in the plasma and hydrogen Lyman-α emission indicate that hydrogen recycling at the surface increases as subsequent discharges are performed; nevertheless, strong pumping of hydrogen is observed even after almost double the cumulative fueling is applied that should saturate the lithium coating to the penetration depth of hydrogen ions. Probe measurements show that when external fueling is terminated, the scrape-off-layer of discharges with fresh coatings decays to lower density and rises to higher electron temperature than for discharges with a partially-passivated surface, consistent with reduced edge cooling from recycled particles. Near the end of the discharge, higher plasma current correlates with reduced τp * and hydrogen emission, suggesting that discharges with fresh coatings achieve higher electron temperature in the core. A novel approach using neutral modeling was developed for the inverse problem of determining the distribution of recycled particle flux from PFC surfaces given a large number of emission measurements, revealing that extremely low levels of recycling (Rcore∼0.6 and Rplate∼0.8) have been achieved with solid lithium coatings. Together with impurity emission measurements, modeling suggests that during periods of particularly low electron density, influx of impurities from the walls contributes substantially to the global particle balance.

magnetic confinement↗
DOE OSTI2025

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2013

PROTOSTELLAR DISK EVOLUTION OVER MILLION-YEAR TIMESCALES WITH A PRESCRIPTION FOR MAGNETIZED TURBULENCE

Magnetorotational instability (MRI) is the most promising mechanism behind accretion in low-mass protostellar disks. Here we present the first analysis of the global structure and evolution of non-ideal MRI-driven T-Tauri disks on million-year timescales. We accomplish this in a 1+1D simulation by calculating magnetic diffusivities and utilizing turbulence activity criteria to determine thermal structure and accretion rate without resorting to a three-dimensional magnetohydrodynamical (MHD) simulation. Our major findings are as follows. First, even for modest surface densities of just a few times the minimum-mass solar nebula, the dead zone encompasses the giant planet-forming region, preserving any compositional gradients. Second, the surface density of the active layer is nearly constant in time at roughly 10 g cm{sup -2}, which we use to derive a simple prescription for viscous heating in MRI-active disks for those who wish to avoid detailed MHD computations. Furthermore, unlike a standard disk with constant-{alpha} viscosity, the disk midplane does not cool off over time, though the surface cools as the star evolves along the Hayashi track. Instead, the MRI may pile material in the dead zone, causing it to heat up over time. The ice line is firmly in the terrestrial planet-forming region throughout disk evolution and can move either inward or outward with time, depending on whether pileups form near the star. Finally, steady-state mass transport is an extremely poor description of flow through an MRI-active disk, as we see both the turnaround in the accretion flow required by conservation of angular momentum and peaks in M-dot (R) bracketing each side of the dead zone. We caution that MRI activity is sensitive to many parameters, including stellar X-ray flux, grain size, gas/small grain mass ratio and magnetic field strength, and we have not performed an exhaustive parameter study here. Our 1+1D model also does not include azimuthal information, which prevents us from modeling the effects of Rossby waves.

79 ASTRONOMY AND ASTROPHYSICS↗