Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary quality”

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 OSTI2026

Development of a Test-Bed for Testing and Refining EarthEn’s Supercritical CO 2 Based Energy Storage System

EarthEn’s energy storage concept leverages supercritical carbon dioxide (sCO 2 ) as a working fluid and relies on compact, high-performance components operating at elevated pressures and temperatures. To accelerate component development and reduce technical risk prior to larger-scale demonstrations, Oak Ridge National Laboratory (ORNL) developed a 100 kW-scale sCO 2 test-bed under a Cooperative Research and Development Agreement with EarthEn (CRADA NO. NFE-24-10050). The objective of the work was to design and construct a flexible experimental facility capable of reproducing key thermodynamic state points and heat-transfer conditions relevant to EarthEn’s thermal energy storage (TES) cycle, with particular emphasis on enabling development and evaluation of next-generation heat exchangers and TES concepts. The test-bed consists of a closed-loop sCO 2 circulation system housed within an open-topped enclosure. In its as-installed configuration, dense-phase sCO 2 is recirculated through a printed circuit recuperator, an electrically heated section, a throttling device used to simulate turbine expansion, and a water-cooled printed circuit heat exchanger that rejects heat to the building chilled-water system before returning to the pump. The pump is driven by a variable frequency drive, enabling controlled adjustment of flow and operating point. A comprehensive instrumentation suite was integrated to support both safe operation and high-quality data collection. Installed sensors include Coriolis flow meters for sCO 2 flow rate and density, resistance temperature detectors and thermocouples distributed throughout the loop (including the heated section and key heat exchanger ports), and pressure transducers for absolute and differential pressure measurements. The facility was designed to support high-pressure (19 MPa nominal) and high-temperature (575°C nominal) operation with credited overpressure protection provided by a rupture disk. Nominal operating conditions were selected to support 100 kW-class testing while maintaining flexibility for non-heated and heated shakedown, control development, and future integration of advanced TES test sections. In parallel with facility development, a system-level thermal-hydraulic model was created using Modelica-based tools to support component sizing, anticipate performance over targeted test conditions, and establish a framework for future model calibration against experimental data. At the conclusion of the project performance period, the facility was in final assembly, and the pressure boundary was nearly completed. However, several practical challenges associated with high-pressure/high-temperature systems and specialized component procurement impacted schedule and prevented initial pump-driven operation and full commissioning within the available resources. This report documents the as-built design, operating capabilities, and instrumentation, and it summarizes key lessons learned related to heater fabrication and testing, first-of-a-kind assembly factors, specialty flange supply constraints, and fill pump corrective actions. Finally, it outlines a phased plan for future commissioning and experimental campaigns, including control and instrumentation shakedown, heater characterization, model calibration, and testing at state points representative of EarthEn’s TES cycle.

25 ENERGY STORAGE↗
DOE OSTI2022

Dry Deposition Methods Based on Turbulence Kinetic Energy: 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single‐Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thallium Bromide Semiconductor Radiation Materials and Detectors Characterization Studies

Two overarching requirements of crucial importance for the commercial and scientific establishment of semiconductor radiation detectors are: (1) exceptionally high-purity crystals with impurity concentrations less than 1 part-per-billion and (2) single crystals that are relatively free from subgrain boundaries, secondary phases, dislocations, and other electrically-active defects. We investigated the properties of thallium bromide (TlBr) material that affect its performance with the goal of increasing the material's commercial viability for radiation detection applications. For this purpose, we used beamlines at the National Synchrotron Light Source (NSLS) at Brookhaven National Laboratory (BNL), performing Micron-scale X-ray Detector Mapping, White Beam X-ray Diffraction Topography, and Micron-scale X-ray Fluorescence. These characterization methods improve the industry's understanding of TlBr and lead to the production of improved instrumentation. In this particular study we report electro-migration measurements of positive-ion Cu, Ag, and Au impurities in TlBr detectors under electric field strengths typically used for device operation. BNL improved TlBr detectors with an electrode design that corrects the response non-uniformities caused by crystal defects. This design can achieve improved energy resolution while using typical-grade, commercial crystals with relaxed quality requirements, thus reducing the overall cost of detectors. Additional characterization will be necessary to fully understand the structure and performance of TlBr with the ultimate goal to achieve the highest energy resolution. Our findings from past work will be presented along with our recommendations for additional investigation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-supervised physics-informed generative networks for phase retrieval from a single X-ray hologram

X-ray phase contrast imaging significantly improves the visualization of structures with weak or uniform absorption, broadening its applications across a wide range of scientific disciplines. Propagation-based phase contrast is particularly suitable for time- or dose-critical in vivo/in situ/operando (tomography) experiments because it requires only a single intensity measurement. However, the phase information of the wave field is lost during the measurement and must be recovered. Conventional algebraic and iterative methods often rely on specific approximations or boundary conditions that may not be met by many samples or experimental setups. In addition, they require manual tuning of reconstruction parameters by experts, making them less adaptable for complex or variable conditions. Here we present a self-learning approach for solving the inverse problem of phase retrieval in the near-field regime of Fresnel theory using a single intensity measurement (hologram). A physics-informed generative adversarial network is employed to reconstruct both the phase and absorbance of the unpropagated wave field in the sample plane from a single hologram. Unlike most state-of-the-art deep learning approaches for phase retrieval, our approach does not require paired, unpaired, or simulated training data. This significantly broadens the applicability of our approach, as acquiring or generating suitable training data remains a major challenge due to the wide variability in sample types and experimental configurations. The algorithm demonstrates robust and consistent performance across diverse imaging conditions and sample types, delivering quantitative, high-quality reconstructions for both simulated data and experimental datasets acquired at beamline P05 at PETRA III (DESY, Hamburg), operated by Helmholtz-Zentrum Hereon. Furthermore, it enables the simultaneous retrieval of both phase and absorption information.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deep-Learning-Based Segmentation of Keyhole in In-Situ X-ray Imaging of Laser Powder Bed Fusion

In laser powder bed fusion processes, keyholes are the gaseous cavities formed where laser interacts with metal, and their morphologies play an important role in defect formation and the final product quality. The in-situ X-ray imaging technique can monitor the keyhole dynamics from the side and capture keyhole shapes in the X-ray image stream. Keyhole shapes in X-ray images are then often labeled by humans for analysis, which increasingly involves attempting to correlate keyhole shapes with defects using machine learning. However, such labeling is tedious, time-consuming, error-prone, and cannot be scaled to large data sets. To use keyhole shapes more readily as the input to machine learning methods, an automatic tool to identify keyhole regions is desirable. In this paper, a deep-learning-based computer vision tool that can automatically segment keyhole shapes out of X-ray images is presented. The pipeline contains a filtering method and an implementation of the BASNet deep learning model to semantically segment the keyhole morphologies out of X-ray images. The presented tool shows promising average accuracy of 91.24% for keyhole area, and 92.81% for boundary shape, for a range of test dataset conditions in Al6061 (and one AliSi10Mg) alloys, with 300 training images/labels and 100 testing images for each trial. Prospective users may apply the presently trained tool or a retrained version following the approach used here to automatically label keyhole shapes in large image sets.

36 MATERIALS SCIENCE↗
DOE OSTI2026

In-situ sensor monitoring of multi-class gas porosity formation in laser powder bed fusion using convolutional neural network

In-situ monitoring of defect formation remains a significant challenge in the laser powder bed fusion (LPBF) process. Recent advances have enabled real-time defect detection with machine learning and in-situ sensing technologies; however, most studies focus on binary classification of keyhole pores, limiting nuanced multi-class pore differentiation and formation mechanisms. This work introduces a multi-class pore detection framework (no pore, small pores < 15 µm, and large pores > 15 µm) by leveraging photodiode sensor data alongside high-fidelity synchrotron X-ray imaging. The 15 µm threshold is selected to distinguish between two fundamentally different defect mechanisms, following the physical size-mechanism boundary established by prior high-resolution synchrotron X-ray characterization of Al6061 LPBF. Distinguishing these classes is critical because large keyhole pores are structurally detrimental, whereas small gas pores are often benign, requiring different process control strategies. Thermal emission monitoring data collected simultaneously with high-speed X-ray imaging at the Stanford Synchrotron Radiation Lightsource (SSRL), are correlated with subsurface melt pool dynamics to establish ground truth. Continuous Wavelet Transform (CWT) with optimized parameters converts the photodiode time-series signals into time–frequency images, facilitating feature extraction. Convolutional Neural Networks (CNN) are then applied for real-time multi-class pore classification in an average inference time of 1 ms per signal window. It achieves 79% accuracy and an Area Under the Receiver Operating Characteristic curve (AUC ROC) score of 0.89 with five-fold cross-validation. The results demonstrate that coupling CWT-based feature engineering with CNN architecture enables reliable multi-class pore detection in Al6061 builds using affordable in-situ sensors. This approach advances scalable and affordable quality assurance in additive manufacturing by moving beyond binary defect detection toward more nuanced classification of porosity mechanisms with in-situ sensors and machine learning.

Laser powder bed fusion, Multi-class pores, In-sit↗
DOE OSTI2025

Replace Human Intelligence with Fast and Smart Geometric Reasoning and Graph Neural Network to Accelerate Next Gen ModSim Workflows

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Closing the Loop between In Situ Stress Complexity and EGS Fracture Complexity

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

42 ENGINEERING↗
DOE OSTI2020

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2022

Techno-economic assessment for the production of algal fuels and value-added products: opportunities for high-protein microalgae conversion

Abstract Background Microalgae possess numerous advantages for use as a feedstock in producing renewable fuels and products, with techno-economic analysis (TEA) frequently used to highlight the economic potential and technical challenges of utilizing this biomass in a biorefinery context. However, many historical TEA studies have focused on the conversion of biomass with elevated levels of carbohydrates and lipids and lower levels of protein, incurring substantial burdens on the ability to achieve high cultivation productivity rates relative to nutrient-replete, high-protein biomass. Given a strong dependence of algal biomass production costs on cultivation productivity, further TEA assessment is needed to understand the economic potential for utilizing potentially lower-cost but lower-quality, high-protein microalgae for biorefinery conversion. Results In this work, we conduct rigorous TEA modeling to assess the economic viability of two conceptual technology pathways for processing proteinaceous algae into a suite of fuels and products. One approach, termed mild oxidative treatment and upgrading (MOTU), makes use of a series of thermo-catalytic operations to upgrade solubilized proteins and carbohydrates to hydrocarbon fuels, while another alternative focuses on the biological conversion of those substrates to oxygenated fuels in the form of mixed alcohols (MA). Both pathways rely on the production of polyurethanes from unsaturated fatty acids and valorization of unconverted solids for use as a material for synthesizing bioplastics. The assessment found similar, albeit slightly higher fuel yields and lower costs for the MA pathway, translating to a residual solids selling price of $899/ton for MA versus $1033/ton for MOTU as would be required to support a $2.50/gallon gasoline equivalent (GGE) fuel selling price. A variation of the MA pathway including subsequent upgrading of the mixed alcohols to hydrocarbon fuels (MAU) reflected a required solids selling price of $975/ton. Conclusion The slight advantages observed for the MA pathway are partially attributed to a boundary that stops at oxygenated fuels versus fungible drop-in hydrocarbon fuels through a more complex MOTU configuration, with more comparable results obtained for the MAU scenario. In either case, it was shown that an integrated algal biorefinery can be economical through optimal strategies to utilize and valorize all fractions of the biomass.

09 BIOMASS FUELS↗
DOE OSTI2023

Simulating organic aerosol in Delhi with WRF-Chem using the volatility-basis-set approach: exploring model uncertainty with a Gaussian process emulator

The nature and origin of organic aerosol in the atmosphere remain unclear. The gas–particle partitioning of semi-volatile organic compounds (SVOCs) that constitute primary organic aerosols (POAs) and the multigenerational chemical aging of SVOCs are particularly poorly understood. The volatility basis set (VBS) approach, implemented in air quality models such as WRF-Chem (Weather Research and Forecasting model with Chemistry), can be a useful tool to describe emissions of POA and its chemical evolution. However, the evaluation of model uncertainty and the optimal model parameterization may be expensive to probe using only WRF-Chem simulations. Gaussian process emulators, trained on simulations from relatively few WRF-Chem simulations, are capable of reproducing model results and estimating the sources of model uncertainty within a defined range of model parameters. In this study, a WRF-Chem VBS parameterization is proposed; we then generate a perturbed parameter ensemble of 111 model runs, perturbing 10 parameters of the WRF-Chem model relating to organic aerosol emissions and the VBS oxidation reactions. This allowed us to cover the model's uncertainty space and to compare outputs from each run to aerosol mass spectrometer observations of organic aerosol concentrations and O:C ratios measured in New Delhi, India. The simulations spanned the organic aerosol concentrations measured with the aerosol mass spectrometer (AMS). However, they also highlighted potential structural errors in the model that may be related to unsuitable diurnal cycles in the emissions and/or failure to adequately represent the dynamics of the planetary boundary layer. While the structural errors prevented us from clearly identifying an optimized VBS approach in WRF-Chem, we were able to apply the emulator in the following two periods: the full period (1–29 May) and a subperiod period of 14:00–16:00 h LT (local time) on 1–29 May. The combination of emulator analysis and model evaluation metrics allowed us to identify plausible parameter combinations for the analyzed periods. We demonstrate that the methodology presented in this study can be used to determine the model uncertainty and to identify the appropriate parameter combination for the VBS approach and hence to provide valuable information to improve our understanding of OA production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2023

COMBLE Radar b1 Processing: Corrections, Calibrations, and Processing Report

The U.S. Department of Energy’s (DOE) Atmospheric Radiation Measurement (ARM) user facility recently concluded its Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE), with its campaign emphasis on marine boundary-layer clouds and mixed-phase clouds during cold-air outbreaks. The COMBLE campaign featured the deployment of the first ARM Mobile Facility (AMF1) to northern Scandinavia (Andenes, Norway), including its standard complement of ARM cloud radars. In keeping with user demands stemming from the previous AMF Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign in Argentina, a post-campaign radar mentor effort was initiated for COMBLE. This activity was intended to improve the usability of the ARM cloud radar data sets in response to overall demands for calibrated, corrected data sets for downstream studies and retrieval applications. In addition to the terrain complexities previously important to CACTI data sets (e.g., clutter designation and/or removal), COMBLE presented new challenges to existing ARM radar mentor capabilities. These included the extension of existing methodologies (e.g., relative calibration adjustment [RCA] target techniques) to frozen environments and the potential issues with their applicability when considering mixed-phase precipitation conditions. As in the previous CACTI documentation (Hardin et al. 2020), the overall calibration and conditioning process in ARM nomenclature is referred to as generating a “b1” datastream. For the radars, these “b1” standards refer to a datastream that has been calibrated (and cross-calibrated), with effort to deliver the highest-quality (well-characterized) data possible. The “b1” radar mentor reporting (this current document) is intended to detail (i) the status/quality of the original “a1” (raw) data sets during the COMBLE AMF campaign, (ii) the corrections and calibrations that are applied to generate the b1 datastreams available on ARM’s Data Discovery, and (iii) the details of the applied methods, e.g., how radar offset/calibration numbers were determined.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2026

Data Center Cybersecurity, Supply Chain Risk Management, and Emerging Regulation Cohort Summary: Takeaways and Action Plans

This report summarizes the outcomes of the Data Center Cohort under the Department of Energy’s Technical Assistance for Digital Assurance (TADA) initiative, aimed at enhancing grid resilience through cybersecurity, supply chain risk management (SCRM), and Cyber-Informed Engineering (CIE). The cohort engaged 17 organizations across utilities, data center operators, vendors, and technology providers in three sessions combining presentations, discussions, and exercises. Key topics included AI-driven load behavior, cybersecurity vulnerabilities in UPS/BESS and cooling systems, governance gaps at utility–data center boundaries, and supply chain integrity. Five cross-cutting themes emerged: interconnection architecture vulnerabilities, fragmented governance, AI-driven stability risks, lack of regulatory frameworks, and long-term supply chain concerns. Actionable recommendations were developed, including implementing DMZ segmentation, formalizing vendor access agreements, designing AI workload limits, and advancing standards through NERC and state-level programs. These strategies aim to strengthen resilience, clarify responsibilities, and ensure secure integration of data centers into the grid.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Colossal Grains for a New Structural Paradigm in Thin-Film Photovoltaics

II-VI based thin-film photovoltaics (PV) have experienced substantial commercial success on a foundation of fast large-scale deposition and processing. This type of processing, regardless of material type, nearly always results in small grain (nm to several um diameters) polycrystalline material that, if not for effective grain boundary (GB) and defect passivation, would severely limit PV performance. In the case of CdTe based PV technology, the nearly ubiquitous CdCl2 treatment is used, in part, to passivate the polycrystalline p-type absorber even though Cl is known to (and used to) compensate p-type II-VI materials. Given that VOC is still the most limited performance metric and it can only be improved through higher p-type absorber concentrations, CdCl2 treatments may be holding the technology back. In this talk, we will discuss results related to a process by which small grain CdSeTe material is rapidly transformed into very large grain material (100-1000um diameters, Fig. 1a). This explosive recrystallization, which we have termed colossal grain growth (CGG), occurs by a Cl free process with thin films (<10um) deposited using the same fast deposition techniques the industry has relied on for more than a decade. Our modeling suggests that GB recombination becomes significantly less important at these grain sizes and thus may eliminate the need for a CdCl2 treatment. This result is complimented by previous results of epitaxial growth at high rates on single crystal seed layers. Thin CdSeTe CGG seed layers can be used for rapid epitaxy of CdTe to complete a standard device stack (Fig.1b). This is an optimistic view and the current limitations and challenges related to this process and these materials will also be discussed, including alloy and substrate roughness sensitivity, intragrain material quality, and complex defect chemistry. These challenges are not insignificant but overcoming them may enable a new structural paradigm for thin-film PV.

CdTe↗