Search NASA⌕ Search

SEARCH · Search NASA

Results for “activity coefficient model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Volumetrically Absorbing Thermal Insulator (VATI) for High-temperature Receivers

This project seeks to exploit volumetric absorption of concentrated solar irradiation in a thin, high-temperature, thermally conductive medium. High efficiency thermal conversion is achieved by volumetric absorption of both incoming solar irradiation (short wavelength) and emitted irradiation (long wavelength), demonstrating a volumetrically absorbing thermally insulating (VATI) effect. The thermalized energy is conducted to the back wall, enclosing the working fluid (molten salt or supercritical CO 2 ). Reliance on conductive thermal transport requires a thermally conductive medium, in the absence of which large temperature gradients drive losses due to emission. The project’s outcomes are important to realize a cost-effective approach to reduce optical and thermal losses from CSP receivers at high temperatures (720°C). High-temperature stable and commercially available porous SiC structures (open-cell foams) were explored as a volumetrically absorbing and thermally insulating layer. To the best of our knowledge, the current state-of-the-art for industrially deployed coatings is Pyromark. However, Pyromark suffers degradation at 700+°C and diurnal temperature changes. This necessitates periodic recoating, and the downtime increases the overall levelized cost of energy (LCOE) production. On the other hand, contemporary research activities have generated significant advances in the development of selective emitter coatings, which present challenges with costs, scalability, and stability. The pursued approach alleviates these concerns by developing a receiver that utilizes the inherent structure of high-temperature stable porous materials to enable robust and cost-effective receivers which require no periodic maintenance downtimes. The overall goal of the project is to experimentally demonstrate a figure of merit (FOM) of 0.92 at a temperature of 720°C and solar irradiation of 1000x concentration with porous receivers. The relevant crystallographic (phase) optical and thermal properties of porous SiC were first characterized. Second, the 3-D geometry of the porous SiC was analyzed using micro-X-ray tomography and converted to CAD data using image processing analysis. Utilizing this 3-D geometry and relevant optical/thermal properties, Monte Carlo- Ray Tracing (MCRT) analysis was performed to extract important parameters governing solar-thermal energy conversion such as extinction coefficient (β, 1/m), scattering albedo (ω) and the scattering phase function (Φ). These properties were then integrated into an in-house radiative and conduction transport model to solve for temperature and transport fluxes characterizing the solar-thermal energy conversion. This model was utilized to predict the FOM for various SiC porous geometries and to identify the highest possible FOM. Testing of the optimized porous structures will be accomplished with a custom-built high-accuracy (< ±4%) FOM measurement test-stand and a 1000x solar concentrator. When neglecting convective losses and resistance at the open boundary and the back wall, the optimized SiC foam leads to a FOM of 0.84. This FOM does not surpass the performance of Pyromark 2500. Yet, conversely to Pyromark 2500, SiC is stable at high temperatures and does not degrade over time. Also, it is possible to boost the FOM of SiC by engineering its effective thermal conductivity and its scattering albedo. A FOM of 0.92 is predicted for an optimized foam by accounting for convective losses and resistance at the open boundary and the back wall. This largely exceeds the FOM of Pyromark 2500 (~0.86) predicted by neglecting convective losses and resistance. Therefore, further engineering of the foam could lead to unprecedented FOMs.

14 SOLAR ENERGY↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2023

Cooling performance of an active-passive hybrid composite phase change material (HcPCM) finned heat sink: Constant operating mode

Here, the present study explores a hybrid thermal management technology based on air cooling and hybrid composite phase change material (HcPCM) filled finned heat sink for cooling performance of lower to medium heat flux dissipating electronic devices. Two-dimensional numerical simulations are conducted to study the conjugate heat transfer effects of three types of finned heat sink: air-cooled finned heat sink, HcPCM-cooled finned heat sink, and hybrid (air-HcPCM) cooled finned heat sink. A heat sink with a constant volume faction of plate-fins is designed in all cases and simultaneous effects of hybrid nanoparticles and air are investigated to keep the heat sink temperature at safe operating conditions between 40–60°C. The effect of air is incorporated into the heat sink by applying the convective heat transfer coefficient of h c = 10–100W /m 2 .K which tends to create the natural convection and forced convection heat transfer characteristics. The heat flux is varied from 25–40kW /m 2 in the current study. The hybrid nanoparticles of carbon additives (GO and MWCNTs) are dispersed into the RT-35HC, used as a PCM, with a volume fraction of 0% to 6%. Transient simulations are carried out using COMSOL Multiphysics to solve the governing equations for PCM based conjugate heat transfer model. The results showed that forced convection heat transfer improved the cooling performance of the hybrid heat sink compared to natural convection heat transfer. The addition of nanoparticles further enhanced thermal enhancement and uniform melting distribution of PCM inside the finned heat sink. The h c between 30 to 50W/m 2 .K shows optimized values for forced convection heat transfer operating conditions. The volume fraction of 2% of GO+MWCNTs nanoparticles in recommended or optimum concentration for uniform melting of PCM inside the finned heat sink.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Model for the Softening Factor within Stages of Work Hardening

A formulation for a dimensionless coefficient c b is derived that represents a scale of microstructural softening for alloys which follow Kocks-Mecking (K-M) work hardening behavior. The variation of the true plastic strain ε p between the proportional limit σ y and the strength σ u at the instability is determined using the Considère criterion. Parameterization of the model is limited to variables expressly measured within tensile experiments. Further development is now made for the softening factor c bi through the individual and sequential stages 3 and 4 of plastic deformation during K-M work hardening Q behavior. Application is shown for tensile test results of Ti-6Al-4V made by different additively manufactured (AM) processes. It is found that the variation in plastic strain as a function of c bi produces a continuous curve representative of the alloy system. The results of data analysis indicate that Θ o3 and c b3 increase, while Θ o4 and c b4 decrease, as the total plastic strain ε p increases. Furthermore, formulations derived for K-M stages 3 and 4 enable the evaluation of other material parameters such as the activation volume ν* for the onset of plastic deformation. This activation volume is found to be near constant, at a ν*-value of 0.353±0.036 nm 3 as computed using a strain-rate sensitivity of strength exponent m of 0.014, irrespective of the AM method used to produce the Ti-6Al-4V alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗
DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-Temperature 2D Optical Relaxation Visualizes Enhanced Oxygen Exchange Kinetics at Metal-Mixed Conducting Oxide Interfaces

Solid-state heterointerfaces are of interest for emergent local behavior that is distinct from either bulk parent compound. One technologically relevant example is the case of mixed ionic/electronic conductor (MIEC)–metal interfaces, which play an important role in electrochemistry. Metal–MIEC composite electrodes can demonstrate improved catalytic activity vs single-phase MIECs, improving fuel cell efficiency. Similarly, MIEC surface reaction kinetics are often evaluated using techniques that place metal current collectors in contact with the surface under evaluation, potentially altering the response vs the native surface. Here, techniques enabling direct and local in situ observation of the behavior at and around such heterointerfaces are needed. Here, we develop a spatially resolved optical transmission relaxation (2D-OTR) method providing continuous evaluation of local, high-temperature, controlled atmosphere defect kinetics across a ~1 cm2 sample area simultaneously in a contact-free manner. We apply it to observe the spatial variance of oxygen incorporation and evolution rates at ~525–620 °C, in response to step changes in oxygen partial pressure, on MIEC SrTi 0.65 Fe 0.35 O 3–x films as a function of distance from porous Pt and Au layers. Using this model geometry, we find significant enhancements in kinetics adjacent to the metals that decay over a few millimeter distance. To extract kinetic parameters, we fit the short-term optical data (initial portion of relaxations) with an exponential decay function appropriate for surface-exchange-limited kinetics, yielding apparent surface exchange coefficients (k chem ) with spatial resolution, decreasing with distance from the metal. To understand the kinetic processes governing the complete (long-term) optical relaxations, we performed COMSOL simulations, which demonstrated that a combination of laterally varying k chem and in-plane diffusion controls the observed kinetics over the full time range. Further support for spatially varying kchem comes from demonstrations of changing surface and bulk chemistry vs distance from the metal–MIEC interface, by X-ray photoelectron and optical absorption spectroscopies, respectively. Although microporous Pt and Au are not excellent electrodes in isolation, both metals exert a synergistic effect on the oxygen surface exchange rate in the presence of the mixed conducting film.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of a Full-Scale Connected U-Net for Reflectivity Inpainting in Spaceborne Radar Blind Zones

CloudSat’s Cloud Profiling Radar is a valuable tool for remotely monitoring high-latitude snowfall, but its ability to observe hydrometeor activity near the Earth’s surface is limited by a radar blind zone caused by ground clutter contamination. This study presents the development of a deeply supervised U-Net-style convolutional neural network to predict cold season reflectivity profiles within the blind zone at two Arctic locations. The network learns to predict the presence and intensity of near-surface hydrometeors by coupling latent features encoded in blind zone-aloft clouds with additional context from collocated atmospheric state variables (i.e., temperature, specific humidity, and wind speed). Results show that the U-Net predictions outperform traditional linear extrapolation methods, with low mean absolute error, a 38% higher Sørensen–Dice coefficient, and vertical reflectivity distributions 60% closer to observed values. The U-Net is also able to detect the presence of near-surface cloud with a critical success index (CSI) of 72% and cases of shallow cumuliform snowfall and virga with 18% higher CSI values compared to linear methods. An explainability analysis shows that reflectivity information throughout the scene, especially at cloud edges and at the 1.2-km blind zone threshold, along with atmospheric state variables near the tropopause, are the most significant contributors to model skill. This surface-trained generative inpainting technique has the potential to enhance current and future remote sensing precipitation missions by providing a better understanding of the nonlinear relationship between blind zone reflectivity values and the surrounding atmospheric state.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Considerations for AMI-Based Operations for Distribution Feeders

More than $5 billion in investments in advanced metering infrastructure (AMI) technologies, AMI deployments, as pervasive secondary network voltage monitoring systems, provide opportunities for utility operations and controls. This paper focuses on the considerations for AMI-based tools and techniques as the industry moves toward operationalizing such large data sets. Phase identification is a first such tool. Numerous distribution network analysis, monitoring, and control applications - including volt/volt-ampere reactive control, state estimation, and distribution automation - require accurate phase connectivity information in the system models. The phase connectivity database maintained by utilities is inaccurate because of a significant amount of missing data, restoration activities, and network reconfiguration. Existing phase identification techniques that estimate phase connectivity work well in distribution feeders that have low or no photovoltaic (PV) generation; however, they fail to identify the phases accurately when considerable PV generation is present. This work addresses the phase identification problem in the presence of high PV generation using statistical analysis methods. Further, insights into the AMI data requirements for this application in terms of data window length and resolution are provided using sensitivity analysis performed on an actual distribution feeder model of San Diego Gas & Electric Company. The results of this study show that the phase connectivity, even in the presence of high PV generation, can be accurately identified using statistical analysis of AMI data of 1 day.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Experimental and analytical study of the hydrodynamic and single and two-phase convective heat transfer performance of flexible PDMS microchannels with micropillar arrays

Various copper and silicon based thermal management systems are used in the cooling of electronics. However, the rigid nature of these materials along with their high thermal and electrical conductivity pose a difficulty in developing direct contact embedded flexible cooling systems that can offer robust cooling performance. The low density, thermal stability, chemical inertness, and electrical insulation of Polydimethylsiloxane (PDMS) make it an ideal material to develop lightweight direct contact thermal management systems for electronics. Its ease of fabrication with tunable flexibility provides the opportunity to go beyond traditional electronics and develop advanced active and passive thermal management systems for a wide–range of applications in foldable and wearable electronics, liquid cooling garments, microgravity, and electric motors. In this study, a flexible PDMS based microchannel with micropillar arrays, which enhance the thermal performance of the device through capillary-assisted flow, has been developed. The hydrodynamic and convection heat transfer performance of three PDMS wick pillar geometries, ranging from a porosity of 0.8–0.91, are investigated and compared under single-phase and two-phase conditions. Dielectric coolant FC-3283 is employed and permeability measurements are made for mass fluxes ranging from 53 kg/m 2 s to 369 kg/m 2 s. Given its conformability, the device demonstrates a deviation from Darcy’s Law, within the laminar regime, with an increasing permeability with mass flux at the rate of ~0.5–0.8 Darcy/(kg/m 2 s). A semi-analytical model has been developed and reported to quantify the conformability of the device. The heat transfer performance is experimentally evaluated using the same dielectric fluid for mass flux ranging from 105 kg/m 2 s to 420 kg/m 2 s with heat fluxes ranging from 1.5 W/cm 2 to 16 W/cm 2 . Heat transfer coefficients of up to 7000 W/m 2 K are observed, which are comparable to copper and silicon microchannels. The effect of porosity on the single phase thermal performance has been evaluated against the pumping power to provide a basis for thermal management system design. Finally, high-speed imaging is performed to study the two-phase flow characteristics to provide insight into the vapor formation and removal.

42 ENGINEERING↗
DOE OSTI2021

CFD modeling of pre-spark heat release in a boosted direct-injection spark-ignition engine

Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.

42 ENGINEERING↗
DOE OSTI2024

Revealing the role of redox reaction selectivity and mass transfer in current–voltage predictions for ensembles of photocatalysts

Photocatalysts are conceptually simple reaction units where nanoscale semiconductors integrated with catalysts drive a pair of redox reactions on illumination. However, the proximity of reaction sites performing cathodic and anodic reactions poses dire challenges to realize large light-to-fuel conversion efficiencies. In this study, a powerful, yet straightforward, equivalent-circuit detail-balance modeling framework is developed and applied to evaluate the performance of photocatalytic systems featuring multiple light absorbers. Specifically, low bandgap iridium-doped strontium titanate is modeled as a Z-scheme photocatalyst to achieve desirable hydrogen evolution and iron-based redox shuttle oxidation reactions. Our model has unique capabilities to simulate competing redox reactions and address mass-transfer limitations. In a significant departure from state-of-the-art circuit models, our study develops tools to perform load-line analyses by incorporating a net electrochemical load curve that includes both desired and competing redox reactions. Consequently, reaction selectivity is predicted from equivalent circuit models for photocatalytic and photoelectrochemical systems. Our investigation into ensembles comprised of multiple, semi-transparent light absorbers reveals their potential to outperform a single, optically thick light absorber, particularly when operated under mass-transfer-limited conditions. However, this outcome hinges on minimizing mass-transfer rates of select redox species to prevent undesired reactions of hydrogen oxidation and/or redox shuttle reduction. Our findings demonstrate that reaction selectivity can be achieved by tuning asymmetry in redox species mass-transfer even with perfectly symmetric electrocatalytic charge-transfer coefficients. The influences of various kinetic, mass-transfer, and thermodynamic parameters are explored to offer crucial insights for synthesis of the next-generation of photocatalysts and selective coatings, and reactor designs.

25 ENERGY STORAGE↗
DOE OSTI2023

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Verification of RESRAD-OFFSITE Code (V.4)

This report documents the verification of RESRAD-OFFSITE Version 4.0 and describes, where necessary, the verification of the following: • The data comprising the standard dose and risk coefficient libraries in the RESRAD database files Master_dcf_ICRP07.mdb and Master_dcf_2k.mdb. • The extraction and transfer of the data from the selected database file to the computational code by the RESRAD-OFFSITE 4.0 interface, ResOWin.exe. • The different processes that are modeled by the main computational code in RESRAD OFFSITE 4.0, ResOMain.exe. • The data displayed in the graphical and text reports. Many verifications were performed as part of the quality assurance quality control program associated with the development and release of RESRAD-OFFSITE 4.0, namely: • developer testing, • internal independent testing, and • release testing. Some were also performed in response to questions from users regarding the performance of the code. The main text of the report focuses on summarizing a subset of those tests, both independent and developer tests that verified the computations performed by the code. The verifications included in this report served as the basis for the development of the release tests of the computational executables and provided the quantitative results to be compared with the code output. The input and output interfaces and the data transfers between the various executables of the code were tested while performing the verification testing. They were tested intentionally during release testing. This report also provides some basic information to help in understanding the activities that were verified. The report: • outlines the components of RESRAD-OFFSITE 4.0 and the interconnections between these components, • outlines the processes modeled by the computational code, • provides summary figures and tables to offer confirmation of the verification of the computational components of the code, • reproduces the verifiers’ reports, if available, in individual appendices, • refers to the previous verification report (Yu et al. 2011) for more details about some of the verifications, and • reproduces the test cases and the testers’ reports from the release testing in individual appendices, when possible.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multi-site assessment of the potential of fine resolution red-edge vegetation indices for estimating gross primary production

Gross primary production (GPP) models driven by fine resolution remote sensing data characterize the spatial and temporal heterogeneities in plant photosynthesis, which is largely dependent on biome-specific maximum photosynthetic capacity. The red-edge reflectance, sensitive to leaf chlorophyll content, is a good proxy of maximum photosynthetic capacity. More importantly, studies show that the red-edge reflectance-related chlorophyll content index (CIr) multiplied by the incident photosynthetic active radiation (PAR in ) strongly correlates to GPP estimated at carbon flux towers (GPP flux ). Yet, to the best of our knowledge, there is no systematic study investigating the general relationship between fine spatial resolution CIr and GPP among biomes and the relationship between CIr and maximum photosynthetic capacity in GPP models. To provide an overview on incorporating space-borne CIr into a GPP model, we applied fine resolution Sentinel-2-derived CIr and GPP flux over 57 flux sites representative of 10 biomes. We investigated the relationship between CIr and GPP flux , and the spatio-temporal relationship between CIr and ecosystem maximum photosynthetic capacity indicated by the potential ecosystem light use efficiency (LUE pot ). We also evaluated the relationship between other five vegetation indices (VIs) and GPP flux . Results showed that the CIr multiplied by PARin has a higher agreement (R 2 > 0.5) with GPP than other VIs. A universal relationship exists between the CIr multiplied by PARin and GPP, except for forest biomes. The CIr also strongly (R 2 > 0.5) relates to the LUE pot during the peak of the growing season. The CIr has a low spatial variance (CV = 0.25) among biomes, highlighting that CIr can be a proxy of maximum photosynthetic capacity in GPP models that do not require biome-dependent coefficients. This study provides insight for incorporating CIr into GPP models and better quantifying global terrestrial photosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Development of an Integrated Multidiagnostic to Assess the High-Z Impurity Fluxes in the Metallic Environment of WEST Using IMAS

WEST is an actively cooled, long-pulse tokamak with nearly all plasma-facing components (PFC) made of tungsten. One of the aims of WEST is to study plasma operations with tungsten PFCs in preparation for long-pulse operations on high-Z divertor devices, such as ITER. For long-pulse operation, the high-Z impurity content and transport to the core plasma are critical concerns that require further measurement and interpretation to improve plasma performance and PFC durability. This work details the impurity influxes in WEST during a series of discharges in which the lower hybrid (LH) injected power was incrementally increased. An analysis has been performed of measurements collected from an array of edge diagnostics. Visible spectroscopy was utilized to measure the spectral radiances generated by fuel particles (D) and impurities (W, O, and C) at the divertor and at the antennas with a newly developed spectral peak-fitting tool used to analyze the data in WEST. The scrape-off layer (SOL) plasma conditions (n e and T e ) measured at the divertor target with flush Langmuir probes and near the outer mid-plane (OMP) using reciprocating Langmuir probes (RCPs) are used to evaluate the number of ionizations per photons (S/XB) coefficients required to estimate the impurity fluxes obtained with the collisional-radiative model ColRadPy. Here, the array of edge diagnostics discussed in this work, coupled with SOL plasma modeling tools, represents a multidiagnostic interpretative modeling workflow that will continue to be applied to upcoming experimental campaigns on the WEST experiment to assess high-Z impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Advancing the quantification of aerosol-cloud interactions with the CALIPSO-CloudSat-Aqua/MODIS record

Aerosol-cloud-precipitation interactions are assessed over the non-polar ocean using more than 11 years of combined Aqua-MODIS, CALIPSO-CALIOP, and CloudSat products. The analysis first shows the benefit of incorporating vertically resolved aerosol extinction coefficient (σext) in aerosol-cloud interactions (ACI) assessments, demonstrating that: σext vertically collocated with the cloud layer () correlates best with cloud droplet number concentration (Nd), column-integrated aerosol optical depth (AOD) cannot explain the Nd variability in the extratropics, and the S-shape of the AOD-Nd relationship reported in previous studies is not replicated when using instead of AOD, with a Nd- linearity more consistent with in-situ studies over the ocean. ACI metric, estimated as the log-scale regression between CALIOP and MODIS Nd reveals that the eastern Pacific is the region with the strongest ACI, followed by the Southern Ocean. The susceptibility of clouds to changes in their liquid water path (LWP) and frequency of precipitation followed a 2-step calculation by combining the Nd- regression (ACI) with the regression between these macrophysical variables and Nd. LWP susceptibility is negative (LWP decreases with aerosol loading) and statistically significant over the eastern Pacific, eastern Atlantic, and extratropics. In contrast, vast areas of the tropical and subtropical ocean feature negligible changes in LWP with aerosol. Precipitation frequency susceptibility is negative, but the values are only significant over the coastal eastern Pacific and Atlantic. The findings suggest that previous modeling assessments relying on AOD may need to be revisited by taking advantage of the synergy between passive and active sensors.

Li, Zhujun↗