Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Biomass pyrolysis devolatilization kinetics of herbaceous and woody feedstocks

Devolatilization kinetics were determined using a modified micropyrolyzer reactor for several biomass feedstocks: switchgrass, corn stover, red oak, and pine. Here, the micropyrolyzer was directly coupled to a flame ionization detector (FID) to track the release of volatiles from the biomass. Time series data from these experiments was analyzed to determine apparent devolatilization rates. Care was taken to assure the experiments were isothermal and kinetically limited calculating a Biot number less than 0.1 and pyrolysis numbers greater than 10, which simplifies the derivation of devolatilization rates. A single, first order reaction was able to model devolatilization rates at temperatures up to 500 °C. No correlation was found between the inorganic content of the biomass and its rate of devolatilization. Apparent activation energies were in the range of 54.9–88.4 kJ mol –1 . The rate coefficient at 500 °C was calculated as 1.90–5.14 s –1 for the four feedstocks.

09 BIOMASS FUELS↗
DOE OSTI2024

Substantial oxygen loss and chemical expansion in lithium-rich layered oxides at moderate delithiation

Delithiation of layered oxide electrodes triggers irreversible oxygen loss, one of the primary degradation modes in lithium-ion batteries. However, the delithiation-dependent mechanisms of oxygen loss remain poorly understood. Here we investigate the oxygen non-stoichiometry in Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ electrodes as a function of Li content by using cycling protocols with long open-circuit voltage steps at varying states of charge. Surprisingly, we observe substantial oxygen loss even at moderate delithiation, corresponding to 2.5, 4.0 and 7.6 ml O 2 per gram of Li 1.18–x Ni 0.21 Mn 0.53 Co 0.08 O 2–δ after resting at upper capacity cut-offs of 135, 200 and 265 mAh g −1 for 100 h. Our observations suggest an intrinsic oxygen instability consistent with predictions of high oxygen activity at intermediate potentials versus Li/Li + . In addition, we observe a large chemical expansion coefficient with respect to oxygen non-stoichiometry, which is about three times greater than those of classical oxygen-deficient materials such as fluorite and perovskite oxides. Furthermore, our work challenges the conventional wisdom that deep delithiation is a necessary condition for oxygen loss in layered oxide electrodes and highlights the importance of calendar ageing for investigating oxygen stability.

Atomistic models↗
DOE OSTI2021

Application of Machine Learning Techniques to an Agent-Based Model of Pantoea

Agent-based modeling (ABM) is a powerful simulation technique which describes a complex dynamic system based on its interacting constituent entities. While the flexibility of ABM enables broad application, the complexity of real-world models demands intensive computing resources and computational time; however, a metamodel may be constructed to gain insight at less computational expense. Here, we developed a model in NetLogo to describe the growth of a microbial population consisting of Pantoea . We applied 13 parameters that defined the model and actively changed seven of the parameters to modulate the evolution of the population curve in response to these changes. We efficiently performed more than 3,000 simulations using a Python wrapper, NL4Py . Upon evaluation of the correlation between the active parameters and outputs by random forest regression, we found that the parameters which define the depth of medium and glucose concentration affect the population curves significantly. Subsequently, we constructed a metamodel, a dense neural network, to predict the simulation outputs from the active parameters and found that it achieves high prediction accuracy, reaching an R 2 coefficient of determination value up to 0.92. Our approach of using a combination of ABM with random forest regression and neural network reduces the number of required ABM simulations. The simplified and refined metamodels may provide insights into the complex dynamic system before their transition to more sophisticated models that run on high-performance computing systems. The ultimate goal is to build a bridge between simulation and experiment, allowing model validation by comparing the simulated data to experimental data in microbiology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Oxidation kinetics of biochar from woody and herbaceous biomass

The goal of this study was to determine oxidation kinetics for biochar produced from fast pyrolysis of various biomass feedstocks. In particular, the role of inherent ash content on the oxidation rate was evaluated. Thermogravimetric analysis (TGA) and fluidized bed combustion experiments were used to explore oxidation kinetics of six fast pyrolysis produced biochars with diverse ash content. Here, the reaction rates varied by a factor of three under chemical kinetic-limited conditions, demonstrating inorganic content impacts oxidation rate. Chemical kinetic rate coefficients were proposed as a function of compositional parameters to determine overall fit and impact. Potassium content was found to have a positive correlation, best describing the differences in the oxidation kinetic rate coefficients. Additionally, feedstocks were subjected to a 1 M HCl acid wash mitigating the catalytic activity of the metals. Acid washed biochars had lower oxidation kinetic rates compared to their unwashed counterparts, indicating the removal of catalytically active metals reduced oxidation rate. Gas composition (CO/CO 2 ) was measured during fluidized bed experiments for both acid-washed and unwashed biochars, which varied between the six biochars. Formation of CO 2 was greatly affected by catalytic metals, finding potassium content to correlate well with a higher percentage of CO 2 formation as compared to CO. Comparison of oxidation rates were made between the two experimental apparatuses to measure the effect of attrition on biochar oxidation.

09 BIOMASS FUELS↗
DOE OSTI2021

A biokinetic model for systemic sodium

This paper describes an updated biokinetic model for systemic sodium (Na), developed for use in a series of reports by the International Commission on Radiological Protection (ICRP) on occupational intake of radionuclides. In contrast to the ICRP's previous model for intake of radio-sodium by workers, the updated model depicts realistic directions of movement of Na in the body including recycling of activity between blood and tissues. The updated model structure facilitates extension of the baseline transfer coefficients for adults to different age groups and to special exposure scenarios such as transfer of radio-sodium from the mother to the foetus or the nursing infant. Dose coefficients for 22Na and 24Na based on the updated model generally do not differ greatly from those based on the ICRP's previous Na model when both models are connected to the ICRP's latest dosimetry system. The main exception is that the updated model yields roughly twofold higher dose coefficients for endosteal bone surface than does the previous model due to the dosimetrically cautious assumption in the updated model that exchangeable Na in bone resides on bone surface.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Reactions of NO 3 with aromatic aldehydes: gas-phase kinetics and insights into the mechanism of the reaction

Rate coefficients for the reaction of NO3 radicals with a series of aromatic aldehydes were measured in a 7300 L simulation chamber at ambient temperature and pressure by relative and absolute methods. The rate coefficients for benzaldehyde (BA), ortho-tolualdehyde (O-TA), meta-tolualdehyde (M-TA), para-tolualdehyde (P-TA), 2,4-dimethyl benzaldehyde (2,4-DMBA), 2,5-dimethyl benzaldehyde (2,5-DMBA) and 3,5-dimethyl benzaldehyde (3,5-DMBA) were k 1 = 2.6 ± 0.3, k 2 = 8.7 ± 0.8, k 3 = 4.9 ± 0.5, k 4 = 4.9 ± 0.4, k 5 = 15.1 ± 1.3, k 6 = 12.8 ± 1.2 and k 7 = 6.2 ± 0.6, respectively, in the units of 10 -15 cm 3 molec. -1 s -1 at 298 ± 2 K. The rate coefficient k 13 for the reaction of the NO 3 radical with deuterated benzaldehyde (benzaldehyde-d1) was found to be half that of k 1 . The end product of the reaction in an excess of NO 2 was measured to be C 6 H 5 C(O)O 2 NO 2 . Furthermore, theoretical calculations of aldehydic bond energies and reaction pathways indicate that the NO 3 radical reacts primarily with aromatic aldehydes through the abstraction of an aldehydic hydrogen atom. The atmospheric implications of the measured rate coefficients are briefly discussed.

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

Optimization of La 2 NiO 4+δ Electrolysis Cell Oxygen Electrode through Surfactant-Enabled LaCoO 3±δ Nanocatalyst Deposition

Lanthanum nickelate (LNO) has shown promise as a Cr-resistant air electrode material for SOECs but has suboptimal surface oxygen exchange properties. Nanocoating of the LNO surface with lanthanum cobaltite (LCO) was chosen to improve cell performance as a surface oxygen conductor. The work focused on the implementation of a two-step nano-LCO film deposition utilizing catechol molecules in a porous LNO electrode. The subgoals of the work were to maintain nanosized LCO particles/ grains to increase active surface area and to control the regularity/ homogeneity of the coating across the microstructure. To achieve these goals, a novel surfactant-enhanced liquid infiltration method was utilized, where nucleation sites were spread across the electrode structure to control the location and size of LCO particles. Various catechol surfactant compositions were evaluated for their ability to control the kinetics of nanoparticle deposition and the homogeneity of the coating. Chelated LCO was characterized by X-ray diffraction (XRD), which found a substantial improvement in LCO formation with surfactant addition and determined polymerized norepinephrine to be the best-performing surfactant, with 88.4% pure LCO formed at low temperature. X-ray photoelectron spectroscopy (XPS) confirmed LCO nanostructures formed by the two-step infiltration process, showing no impurities and a stable perovskite structure. Deposition kinetics were analyzed using atomic force microscopy (AFM), correlating infiltration times and solution molarity to nanoparticle size and distribution, the results of which were confirmed in symmetrical cell samples by scanning electron microscopy (SEM). Electrochemical impedance spectroscopy (EIS) testing demonstrated substantial improvements in polarization resistance, where the nanocoating reduced the resistance by ∼55% to 0.152 Ω·cm 2 at 700 °C and 0.039 Ω·cm 2 at 800 °C. Electrical conductivity relaxation (ECR) at this temperature confirmed an improved surface oxygen exchange coefficient of the LCO + LNO heterostructure predicted by the Bode data from EIS, alongside a reduction in activation energy by about 30%.

Deposition↗
DOE OSTI2025

Visualizing energy transfer between redox-active colloids

Redox-active colloids (RACs) represent a novel class of energy carriers that exchange electrical energy upon contact. Understanding contact-mediated electron transfer dynamics in RACs offers insights into physical contact events in colloidal suspensions and enables quantification of electrical energy transport in nonconjugated polymers. Redox-based electron transport was directly observed in monolayers of micron-sized RACs containing ethyl-viologen side groups via fluorescence microscopy through an unexpected nonlinear electrofluorochromism that is quantitatively coupled to the redox state of the colloid. Via imaging studies, using this electrofluorochromism, the apparent charge transfer diffusion coefficient D CT of the RAC was easily determined. The visualization of energy transport within suspensions of redox-active colloids was also demonstrated. Our work elucidates fundamental mechanisms of energy transport in colloidal systems, informs the development of next-generation redox flow batteries, and may inspire new designs of smart active soft matter including conductive polymers for applications ranging from electrochemical sensors and organic electronics to colloidal robotics.

Qu, Alan Subing [Univ. of Illinois at Urbana-Champ↗
DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗
DOE OSTI2020

Thermodynamically consistent forward and reverse degrees of rate control in reversible reactions

The net rate of a composite reaction is the difference between the forward and reverse reaction rates, which are kinetically distinct despite sharing elementary reaction steps and therefore have different rate-controlling transition states and species. Thus, degrees of rate control defined to identify rate-controlling transition states and species for the net rate confound contributions from the forward and reverse reactions. Herein, the forward and reverse degrees of rate control are defined to independently quantify the extent that species and transition states control the forward and reverse rates in reversible reactions. These degrees of rate control are defined as the relative change in the forward and reverse reaction rates per k B T decrease in the standard-state molecular free energies of transition states and species, and they are related to experimentally measurable quantities such as steady-state fractional coverages, reaction orders, and energies and entropies of activation of the forward and reverse reaction rates. The forward/reverse degrees of rate control represent stoichiometric coefficients for species and transition states in an equilibrium between the apparent initial states and the apparent transition states of the apparent rate-controlling steps of forward and reverse overall reactions. At equilibrium, the apparent transition states for the forward reaction and reverse reaction converge, and thus the forward and reverse rate-controlling steps combine to form a single apparent rate-controlling step at equilibrium. Here, this apparent rate-controlling step is comprised of an apparent initial state, transition state, and final state, where the apparent final state of the forward reaction is the apparent initial state of the reverse reaction. The apparent rate-controlling step behaves identically to an elementary step reaction at equilibrium with a pseudo-mass-action rate function given by the transition-state-theory (TST) form rate function (Foley and Bhan, 2020) with a stoichiometric number equal to the affinity-averaged stoichiometric number, $\bar{σ}$ The ratio of the forward and reverse TST-form rate functions is identical to the overall thermodynamic equilibrium relation, consistent with the principles of microscopic reversibility and detailed balance at equilibrium.

42 ENGINEERING↗
DOE OSTI2021

Predicting Thermal Excursions During in situ Oxidative Regeneration of Packed Bed Catalytic Fast Pyrolysis Catalyst

Ex situ catalytic fast pyrolysis (CFP) uses a secondary reactor to upgrade biomass pyrolysis vapors to stabilized CFP oils with reduced oxygen content. In one configuration, the secondary reactor is operated as a packed-bed swing reactor system which allows coke-deactivated beds to be decarbonized in situ while other beds remain online for vapor upgrading. In situ decarbonization must be done carefully to avoid irreversible deactivation and/or physical degradation of catalyst pellets. Given that packed bed reactors are well known to have poor heat transfer characteristics, this is a critical issue impacting scaleability and commercial viability of the technology. To predict thermal excursions during regeneration, finite element computational models have been built to assist in scaling up oxidative decarbonization of a Pt/TiO2 CFP catalyst (0.5 mm spheres) from a bench scale packed bed with 100 g of catalyst to a pilot scale packed bed with 2 kg of catalyst and internal cooling tubes. Based on transient measurements of outlet temperature and effluent CO2 concentration, and using an assumed coke profile and activation energy, this paper demonstrates that specific combinations of effective thermal conductivity and wall heat transfer coefficient can fit bench scale oxidative regeneration data equally well. For the upscaled 2 kg bed, four bench-scale “best fit” parameter pairs give different predictions for location and magnitude of thermal excursions, with the maximum computed bed temperature gradients ranging from 30 °C cm-1 to as high as 3000 °C cm-1. The larger the fraction of heat removal by conduction through the cooling tubes, the greater the differences between the parameter pairs. The modelling results presented in this paper cast doubt on the industrial viability of the proposed combination of catalyst, bed and regeneration process, and point to the need for alternate reactor designs. However, there is considerable uncertainly in some of the key model parameters. The reliability of model predictions can be increased by adding more temperature measurements at key bed locations, testing additional variations in process conditions, performing careful bed dissections to determine the true coke profile, and perhaps most importantly, directly measuring the effective thermal conductivity of the catalyst pellets.

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

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗
DOE OSTI2026

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Guideline to Mitigate Interfacial Degradation Processes in Solid‐State Batteries Caused by Cross Diffusion

Abstract Diffusion of transition metals across the cathode–electrolyte interface is identified as a key challenge for the practical realization of solid‐state batteries. This is related to the formation of highly resistive interphases impeding the charge transport across the materials. Herein, the hypothesis that formation of interphases is associated with the incorporation of Co into the Li 7 La 3 Zr 2 O 12 lattice representing the starting point of a cascade of degradation processes is investigated. It is shown that Co incorporates into the garnet structure preferably four‐fold coordinated as Co 2+ or Co 3+ depending on oxygen fugacity. The solubility limit of Co is determined to be around 0.16 per formula unit, whereby concentrations beyond this limit causes a cubic‐to‐tetragonal phase transition. Moreover, the temperature‐dependent Co diffusion coefficient is determined, for example, D 700 °C = 9.46 × 10 −14 cm 2 s −1 and an activation energy E a = 1.65 eV, suggesting that detrimental cross diffusion will take place at any relevant process condition. Additionally, the optimal protective Al 2 O 3 coating thickness for relevant temperatures is studied, which allows to create a process diagram to mitigate any degradation with a minimum compromise on electrochemical performance. This study provides a tool to optimize processing conditions toward developing high energy density solid‐state batteries.

Chemistry↗
DOE OSTI2021

Strongly temperature dependent ferroelectric switching in AlN, Al 1-x Sc x N, and Al 1-x B x N thin films

This manuscript reports the temperature dependence of ferroelectric switching in Al 0.84 Sc 0.16 N, Al 0.93 B 0.07 N, and AlN thin films. Polarization reversal is demonstrated in all compositions and is strongly temperature dependent. Between room temperature and 300 °C, the coercive field drops by almost 50% in all samples, while there was very small temperature dependence of the remanent polarization value. Furthermore, over this same temperature range, the relative permittivity increased between 5% and 10%. Polarization reversal was confirmed by piezoelectric coefficient analysis and chemical etching. Applying intrinsic/homogeneous switching models produces nonphysical fits, while models based on thermal activation suggest that switching is regulated by a distribution of pinning sites or nucleation barriers with an average activation energy near 28 meV.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Numerical analysis of gas exhaust in Wendelstein 7-X using the direct simulation Monte Carlo method

The present work is focused on a 3D numerical assessment of the Wendelstein 7-X (W7-X) particle exhaust. For all the numerical simulations the direct simulation Monte Carlo solver of the DIVGAS workflow, has been employed. The complex 3D geometry of the sub-divertor region includes the pumping gap panel, supporting structures, cooling pipes as well as the cryo-vacuum pump. All the considered flow simulations correspond to the Standard magnetic configuration of W7-X. The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; The coupling between EMC3-EIRENE and DIVGAS, which considers the fact that the incoming neutral particle flux at the sub-divertor is based on realistic plasma background, has been demonstrated. Three plasma scenarios have been considered, for which is clearly seen that by increasing the heating power, the neutral pressure as well as the resulting pumping efficiency is increased. The obtained numerical results of the neutral pressure in the sub-divertor lie within a more general scan matrix, which assumes a wider range of incoming particle flux, namely 10 19 –10 24 (s −1 ). It has been observed that, the sub-divertor neutral pressure is proportional to the incoming neutral particle flux, with the effective pumping speed to be a constant of proportionality. The influence of switching off the cryo-vacuum pump on the sub-divertor pressure is rather modest and a weak increase of the neutral pressure in the sub-divertor is expected. Correlations of the sub-divertor pressure with the total incoming particle flux as well as the individual pumped flux at each of the AEH and AEP sections have been deduced. Moreover, it has been demonstrated that the influence of the incoming neutral particle flux on the albedo coefficient at the AEH and AEP pumping gaps is rather weak. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of future experimental campaigns.

DSMC method↗
DOE OSTI2021

Active Tuning of Optical Constants in the Visible–UV: Praseodymium-Doped Ceria—a Model Mixed Ionic–Electronic Conductor

Mixed ionic-electronic conductors offer chemical and electrical means for active tuning of their optical constants, e.g. with variations in oxygen non-stoichiometry in Pr 0.1 Ce 0.9 O 2-δ , enabling implementation of adaptive thin film optical devices. Here, we demonstrate in-situ chemo-tuning of the extinction coefficient in Pr 0.1 Ce 0.9 O 2-δ at elevated temperatures and provide a tuning model that treats the interdependence of mobile oxygen vacancies and small polarons coupled to variations in optically active praseodymium ions. Further, a new means for electro-tuning of the optical constants of mixed ionic-electronic conductors was demonstrated experimentally and modeled for Pr 0.1 Ce 0.9 O 2-δ thin films deposited on grid-like electrode structures. Modeling of non-steady state optical transmittance modulations in the latter allows for estimation of oxygen vacancy mobility that determines the switching speed of the device. Quenched-in values of nr and k to room temperature become nonvolatile, providing a modulation range in the extinction coefficient of Δk ~0.1 (change of ~800%) and in the refractive index of Δn r ~0.1 (relative to initial n r of ~2.35). Key figures of merit including, transmission optical modulation of ~0.04 per 1 mV nm -1 , switching energy per area of 670 pJ μm -2 and switching times of seconds, were demonstrated, with further improvements possible.

36 MATERIALS SCIENCE↗