SEARCH · Search NASA

Results for “parametric optimization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2024

LENS: Learning Enabled Network Synthesis

RTRC and UMD have developed novel machine learning based methods under the ARPA-E DIFFERENTIATE program for rapid acceleration of hypothesis generation in complex architecture design spaces involving both discrete choices of component inclusion and interconnection and continuous parametric decisions. The project named Learning Enabled Network Synthesis (LENS) further demonstrated the developed methods on challenging electrical power converter design problems by identifying the most suitable circuit topologies and simultaneously selecting the most appropriate components to achieve optimized design of power converter with improved performances. We demonstrated that LENS could enable exploration of very large design space of circuit topologies and components by addressing the limitations of conventional design process in non-linear, high switching speed, multi-dimensional power converter design and optimization. The key innovation developed in LENS is the seamless integration of statistical learning and logical reasoning techniques and building on the individual strengths of these techniques for rapid hypothesis discovery. The main component of LENS comprises of: 1) Graph Reasoning Engine (GRE) to enforce composition rules that rapidly reject all discrete architectures that are composed incorrectly and generates an adaptive database of feasible designs which can be used by ML modules, 2) Graph Generative Learning module which is a deep neural network based generative model for graph architectures which can enable design space exploration beyond the dataset generated by the GRE, 3) Graph Reduced Order Model (ROM) for graph domains for accelerating computation of output metrics, and 4) Active learning and Rule Discovery module for sample efficient learning and extracting logical rules from the learned ML models which will be integrated in the GRE to enhance the filtering effectiveness. LENS approach can be applied to any design domains where designs can be represented as multi-attribute graphs. The LENS team integrated the various technical innovations listed above into an optimization pipeline and exercised the optimization pipeline on the converter design problem. The LENS project demonstrated that the developed AI/ML technologies can be used to generate novel converter circuits >45x faster than experts on chosen use-cases. This can enable faster design space exploration and identification of new designs which are not considered by experts due to the increasing design space complexity. This has significant potential impact on the public and energy needs of the country. It is currently estimated that 30% of all electrical powers generated passes through power converters. The future estimate is that 80% of all power generated would be passing through converters. LENS fills a critical gap in this space since by accelerating the design process the designers would be able to generate more efficient converters which can lead to significant energy savings for the country.

42 ENGINEERING↗

DOE OSTI2021

HAIMOS Ensemble Forecasts for Intra-day and Day- Ahead GHI, DNI and Ramps

The objective of this research is to develop a hybrid physics-based/data-driven forecast model to improve direct normal and global horizontal irradiance (DNI and GHI) prediction for horizons ranging from 1 to 72 hours. Project objectives also address key gaps in state-of-the-art solar forecasting: accurate probabilistic solar forecasts and the forecasting of large irradiance ramps (ramp onset and magnitude). The proposed model ensembles Numerical Weather Prediction (NWP) forecasts, determinist physics-based algorithms, and new-generation cloud cover products (high-resolution rapid refresh satellite images and Large Eddy Simulations). The result is the Hybrid Adaptive Input Model Objective Selection (HAIMOS) ensemble model. HAIMOS blends state of the art machine learning methodologies with physics-based models for cloud cover and cloud optical depth forecasts. The technical activities followed a two-pronged strategy. First, the preprocessing of data, the selection of inputs to the nonlinear approximators, the type of approximator and objective functions, and post-processing ensembling techniques included in HAIMOS were all optimized adaptively to find the best model for a specific goal (reduce DNI/GHI forecast error, improve the prediction of ramp onset, etc.). Second, a large effort was put in improving cloud identification and the forecast of cloud cover and cloud optical depth. To this end, new-generation cloud parametrization products were developed in this work. These include improved algorithms to assist in cloud identification, cloud classification and cloud parametrization from satellite images – three key factors in the accuracy of 1 to 6-hours irradiance forecasts and prediction of ramp onset. Furthermore, we also included cloud information extracted from high resolution rapid refresh satellite images (GOES-16) and Large Eddy Simulations (LES). LES was used to model the atmosphere in detail over locations of interest and produce cloud optical depth forecasts. Once these data streams were validated, they were used as input data to the HAIMOS forecast. The model was developed using data from several climatologically distinct locations with potential for high solar penetration. In the last year of the project, we conducted a validation campaign according to the guidelines stipulated by the Topic Area 1 project as described in the FOA. This effort brings, for the first time, proven machine-learning methodologies for generating state-of-the-art solar forecasts interweaved with detailed physics-based models for cloud detection, and cloud optical depth forecasts. HAIMOS will generate accurate irradiance probabilistic forecast to assist in reducing solar generation prediction error. Globally optimized solar forecast models are more likely to impact solar energy stakeholders. The goal of this project was to increase the state-of-the-art forecast skill from their present values of 10 to 35%. At the end of the project, we achieved between 30% and 50% forecast skill across a wide range of horizons for both GHI and DNI.

14 SOLAR ENERGY↗

DOE OSTI2024

The hyperplane of early-type galaxies: using stellar population properties to increase the precision and accuracy of the fundamental plane as a distance indicator

ABSTRACT We use deep spectroscopy from the SAMI (Sydney-AAO Multi-object Integral) Galaxy Survey to explore the precision of the fundamental plane (FP) of early-type galaxies as a distance indicator for future single-fibre spectroscopy surveys. We study the optimal trade-off between sample size and signal-to-noise ratio (SNR), and investigate which additional observables can be used to construct hyperplanes with smaller intrinsic scatter than the FP. We add increasing levels of random noise (parametrized as effective exposure time) to the SAMI spectra to study the effect of increasing measurement uncertainties on the FP- and hyperplane-inferred distances. We find that, using direct-fit methods, the values of the FP and hyperplane best-fitting coefficients depend on the spectral SNR, and reach asymptotic values for a mean $\langle \mathrm{ SNR} \rangle =40\, \mathrm{\mathring{\rm A}}^{-1}$. As additional variables for the FP we consider three stellar-population observables: light-weighted age, stellar mass-to-light ratio, and a novel combination of Lick indices ($I_\mathrm{age}$). For an $\langle \mathrm{ SNR} \rangle =45~\mathrm{\mathring{\rm A}}^{-1}$ (equivalent to 1-h exposure on a 4-m telescope), all three hyperplanes outperform the FP as distance indicators. Being an empirical spectral index, $I_\mathrm{age}$ avoids the model-dependent uncertainties and bias underlying age and mass-to-light ratio measurements, yet yields a 10 per cent reduction of the median distance uncertainty compared to the FP. We also find that, as a by-product, the $I_\mathrm{age}$ hyperplane removes most of the reported environment bias of the FP. After accounting for the different SNR, these conclusions also apply to a 50 times larger sample from SDSS-III (Sloan Digital Sky Survey). However, in this case, only $\mathrm{ age}$ removes the environment bias.

D’Eugenio, Francesco (ORCID:0000000323888172)↗

DOE OSTI2024

CMB-S4: Iterative Internal Delensing and r Constraints

Abstract The tightest constraints on the tensor-to-scalar ratio r can only be obtained after removing a substantial fraction of the lensing B -mode sample variance. The planned Cosmic Microwave Background (CMB)-S4 experiment ( cmb-s4.org ) will remove the lensing B -mode signal internally by reconstructing the gravitational lenses from high-resolution observations. We document here a first lensing reconstruction pipeline able to achieve this optimally for arbitrary sky coverage. We make it part of a map-based framework to test CMB-S4 delensing performance and its constraining power on r , including inhomogeneous noise and two non-Gaussian Galactic polarized foreground models. The framework performs component separation of the high-resolution maps, followed by the construction of lensing B -mode templates, which are then included in a parametric small-aperture map cross-spectra-based likelihood for r . We find that the lensing reconstruction and framework achieve the expected performance, compatible with the target σ ( r ) ≃ 5 · 10 −4 in the absence of a tensor signal, after an effective removal of 92%–93% of the lensing B -mode variance, depending on the simulation set. The code for the lensing reconstruction can also be used for cross-correlation studies with large-scale structures, lensing spectrum reconstruction, cluster lensing, or other CMB lensing-related purposes. As part of our tests, we also demonstrate the joint optimal reconstruction of the lensing potential with the lensing curl potential mode at second order in the density fluctuations.

Astronomy & Astrophysics↗

DOE OSTI2025

Helical Photonic Metamaterials for Encrypted Chiral Holograms

Helical structures are among the most quintessential three-dimensional (3D) forms that exhibit mirror asymmetry, a hallmark of chirality. Various structural parameters of helices directly linked to chiroptical properties highlight their importance as essential optical metamaterials for polarization-resolved sensors, imaging, and spectroscopies. However, such function-defining properties remain incompletely understood due to fabrication challenges and the lack of a relationship between structure and optical properties. Here, helical structures are analyzed parametrically, and correlations are established that are applicable to the design of chiral helical optical metamaterials. By systematically varying independent parameters—such as from single-turn to five-turn helices and from small major radii to larger ones optimized to fit the unit cell—the underlying relationships with ellipticty are revealed. In addition to theoretical modeling, the findings are experimentally validated using 3D printing and terahertz spectroscopy. The results demonstrate that optimized helical structures are mechanically tunable and exhibit unprecedented optical properties, including broadband and high-magnitude ellipticity spectra. Being embedded in soft elastomers, helical arrays can serve as soft, stretchable optical-mechanical sensors and holograms containing encoded information, such as barcodes and quick response (QR) codes. Chiral QR codes are realized using pixelated single helices with different handedness, demonstrating their potential as advanced encryption/decryption systems for security applications and chiral metaholograms.

Encrypted QR codes↗

DOE OSTI2026

Poisson Log-Normal Process for Count Data Prediction

Modeling count data is important in physics and other scientific disciplines, where measurements often involve discrete, non-negative quantities such as photon or neutrino detection events. Traditional parametric approaches can be trained to generate integer-count predictions but may struggle with capturing complex, non-linear dependencies often observed in the data. Gaussian process (GP) regression provides a robust non-parametric alternative to modeling continuous data; however, it cannot generate integer outputs. We propose the Poisson Log-Normal (PoLoN) process, a framework that employs GP to model Poisson log-rates. As in GP regression, our approach relies on the correlations between data points captured via GP kernel structure rather than explicit functional parameterizations. We demonstrate that the PoLoN predictive distribution is Poisson-LogNormal and provide an algorithm for optimizing kernel hyperparameters. Furthermore, we adapt the PoLoN approach to the problem of detecting weak localized signals superimposed on a smoothly varying background - a task of considerable interest in many areas of science and engineering. Our framework allows us to predict the strength, location and width of the detected signals. We evaluate PoLoN's performance using both synthetic and real-world datasets, including the open dataset from CERN which was used to detect the Higgs boson at the Large Hadron Collider. Our results indicate that the PoLoN process can be used as a non-parametric alternative for analyzing, predicting, and extracting signals from integer-valued data.

Saha, Anushka [Rutgers U., Piscataway]↗

DOE OSTI2020

GMRES with embedded ensemble propagation for the efficient solution of parametric linear systems in uncertainty quantification of computational models

In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.

BLAS↗

DOE OSTI2022

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗

DOE OSTI2024

Photovoltaic and Cost Analysis for Winston-Salem, North Carolina

This study assesses the feasibility of solar installations at various sites in Winston-Salem, focusing on factors such as solar resource availability, electricity costs, and rooftop area for photovoltaic systems. The analysis begins by estimating daily energy requirements based on annual electricity use, followed by adjusting for seasonal and operational fluctuations. Using regional solar data from the National Solar Radiation Database (NSRDB), we calculate average peak sun hours to determine effective system sizing. A parametric approach using the System Advisor Model (SAM) refines this sizing process, incorporating a safety margin of 1.2 to address demand peaks. Each PV system is designed to fit available rooftop space, and a coverage threshold of 70% is identified as optimal for maximizing cost savings and sustainability. This threshold allows installations to meet substantial energy demands, supporting energy resilience and enhancing economic returns.

14 SOLAR ENERGY↗

DOE OSTI2020

EVALUATION OF CONTROLLED END GAS AUTO IGNITION WITH EXHAUST GAS RECIRCULATION IN A STOICHIOMETRIC, SPARK IGNITED, NATURAL GAS ENGINE

Many stationary and heavy-duty on-road natural gas fueled engines today operate under stoichiometric conditions with a three-way catalyst. The disadvantage of stoichiometric natural gas engines compared to lean-burn natural gas and diesel engines is lower efficiency, resulting primarily from lower power density and compression ratio. Exhaust gas recirculation (EGR) coupled with advanced combustion controls can enable operation with higher compression ratio and power density, which yields higher efficiency. This also results in engine operation between the limits of knock and misfire. Operation between these limits has been named controlled end gas auto-ignition (C-EGAI) and can be used to improve the brake efficiency of the engine. Various methods of cylinder pressure-based knock quantification were explored to implement C-EGAI. The indicated quantification methods are used for the implementation of a control scheme for C-EGAI with a relation to the fractional heat release due to auto-ignition. A custom EGR system was built and the effect of EGR on the performance of a stochiometric, spark ignited, natural gas engine is evaluated. C-EGAI is implemented and the optimal parameters are determined for peak performance under EGR and C-EGAI conditions. In this study, knock detection is used for the recognition, magnitude, and location of the auto-ignition events. Cylinder pressure-based knock detection was the primary method for determining the occurrence and location of knock but was also used for implementing the ignition control scheme for controlled end gas auto-ignition. The combustion intensity metric (CIM) enabled parametric ignition timing control which allowed for the creation of a relationship between fractional heat release due to auto-ignition and CIM. Many stationary and heavy-duty on-road natural gas fueled engines today operate under stoichiometric conditions with a three-way catalyst. The disadvantage of stoichiometric natural gas engines compared to lean-burn natural gas and diesel engines is lower efficiency, resulting primarily from lower power density and compression ratio. Exhaust gas recirculation (EGR) coupled with advanced combustion controls can enable operation with higher compression ratio and power density, which yields higher efficiency. This also results in engine operation between the limits of knock and misfire. Operation between these limits has been named controlled end gas auto-ignition (C-EGAI) and can be used to improve the brake efficiency of the engine. Various methods of cylinder pressure-based knock quantification were explored to implement C-EGAI. The indicated quantification methods are used for the implementation of a control scheme for C-EGAI with a relation to the fractional heat release due to auto-ignition. A custom EGR system was built and the effect of EGR on the performance of a stochiometric, spark ignited, natural gas engine is evaluated. C-EGAI is implemented and the optimal parameters are determined for peak performance under EGR and C-EGAI conditions. In this study, knock detection is used for the recognition, magnitude, and location of the auto-ignition events. Cylinder pressure-based knock detection was the primary method for determining the occurrence and location of knock but was also used for implementing the ignition control scheme for controlled end gas auto-ignition. The combustion intensity metric (CIM) enabled parametric ignition timing control which allowed for the creation of a relationship between fractional heat release due to auto-ignition and CIM. Both exhaust gas recirculation and controlled end gas auto-ignition were analyzed with a cooperative fuel research (CFR) engine modified for boosted fuel/air intake. The data was interpreted to provide a proper evaluation of unique analytical methods to quantify the results of C_EGAI and characterize the live auto-ignition events. The control variables for this method of C-EGAI were optimized with EGR conditions to generate the point of peak performance on the CFR engine under stoichiometric, spark ignited, natural gas conditions.

Bayliff, Scott Michael↗

DOE OSTI2022

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

A Design Method for Low-Pressure Venturi Nozzles

The purpose of this work is to provide empirical design models for low-pressure, subsonic Venturi nozzles. Experimentally validated simulations were used to determine the effect of nozzle geometry and operating conditions on the suction ratio (ratio of suction mass flow rate to motive mass flow rate) of low-pressure, subsonic Venturi nozzles, over a wide range of geometries and operating conditions, through a parametric study. The results of the parametric study were used to develop seven empirical models, each with a different range of applicability or calculating a different indicator of nozzle performance (i.e., suction ratio, momentum ratio, or dynamic pressure ratio), of the Venturi nozzles using a constrained multi-variable global optimization method. Of the seven empirical models, the best models were found to be those for low- (less than one) and high-suction ratios (greater than one), with mean absolute percentage errors of 5% and 18%, respectively. These empirical models provide a design tool for subsonic, low-pressure Venturi nozzles that is more than an order of magnitude more accurate than a governing equation approach or conventional flow head calculations. These newly-developed empirical models can be applied for initial nozzle design when precise suction ratios are required.

O’Hern, Hannah↗

DOE OSTI2021

Learning from many collider events at once

There have been a number of recent proposals to enhance the performance of machine learning strategies for collider physics by combining many distinct events into a single ensemble feature. To evaluate the efficacy of these proposals, we study the connection between single-event classifiers and multievent classifiers under the assumption that collider events are independent and identically distributed. We show how one can build optimal multievent classifiers from single-event classifiers, and we also show how to construct multievent classifiers such that they produce optimal single-event classifiers. This is illustrated for a Gaussian example as well as for classification tasks relevant for searches and measurements at the Large Hadron Collider. We extend our discussion to regression tasks by showing how they can be phrased in terms of parametrized classifiers. Empirically, we find that training a single-event (per-instance) classifier is more effective than training a multievent (per-ensemble) classifier, as least for the cases we studied, and we relate this fact to properties of the loss function gradient in the two cases. While we did not identify a clear benefit from using multievent classifiers in the collider context, we speculate on the potential value of these methods in cases involving only approximate independence, as relevant for jet substructure studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

System Analysis of an Internal Combustion Engine (ICE) - Solid Oxide Fuel Cell (SOFC) Hybrid Cycle

The variability of renewable energy sources poses challenges for reliable grid operation. Conventional thermal power sources, though reliable, often lack operational flexibility. Hybrid energy systems that integrate Solid Oxide Fuel cells (SOFC) with Internal Combustion Engines (ICE) offer a promising solution by achieving relatively higher efficiency and grid-following capability. This study investigates performance of a 100-kW pressurized SOFC-ICE hybrid cycle. In this configuration, unutilized SOFC fuel is used to drive the engine, with a turbocharger providing air supply and an external reformer generating syngas. System components were numerically modeled using MATLAB/SIMULINK for the SOFC and reformer, and EBSILON® for the ICE and balance of plant. Parametric studies varied fuel utilization (70-90%), reformer temperature (600 – 1000K), anode off-gas recirculation (0 – 70%), and current density (0.2 – 0.55 A/cm2). Results show that the SOFC and ICE operate as thermally independent topping and bottoming cycles, achieving peak efficiency of 62% under optimized conditions.

hybrid↗

DOE OSTI2021

MTW-OPAL: a technology development platform for ultra-intense optical parametric chirped-pulse amplification systems

Optical parametric chirped-pulse amplification implemented using multikilojoule Nd:glass pump lasers is a promising approach for producing ultra-intense pulses (>10 23 W/cm 2 ). We report on the MTW-OPAL Laser System, an optical parametric amplifier line (OPAL) pumped by the Nd:doped portion of the multi-terawatt (MTW) laser. This midscale prototype was designed to produce 0.5-PW pulses with technologies scalable to tens of petawatts. Technology choices made for MTW-OPAL were guided by the longer-term goal of two full-scale OPALs pumped by the OMEGA EP to produce 2 × 25-PW beams that would be co-located with kilojoule–nanosecond ultraviolet beams. Several MTW-OPAL campaigns that have been completed since “first light” in March 2020 show that the laser design is fundamentally sound, and optimization continues as we prepare for “first-focus” campaigns later this year.

42 ENGINEERING↗

DOE OSTI2026

Reparameterizing a Lipid Force Field Using Small-Angle X-ray Scattering to Improve Predictions of Multicomponent Membranes under Organic Solvent Stress

Understanding and predicting properties of lipid-bilayer membranes are essential to elucidating organismal physiology and pathophysiology. Therefore, substantial efforts have been undertaken to derive accurate molecular-mechanics force fields (FF) to allow simulation of their properties; however, much of these past efforts focused on tuning force fields to accurately reproduce the properties of model single-component membranes and not multicomponent or stressed membranes. Here, we tuned the CHARMM36 FF for a simple 2-component model of a Gram-positive bacterium. This updated force field is parametrized against a single condensed-phase property, the small-angle X-ray scattering (SAXS) intensities of the 2-component lipid system, using a version of ForceBalance implemented previously (named ForceBalance-SAS) with SAXS and small-angle neutron scattering (SANS) intensities as optimization targets. After tuning, we observe improved agreement with the experimental SAXS of the membrane under 1-butanol and tetrahydrofuran solvent stresses, with a factor of 10.0 and 7.5 reduction, respectively, in χ 2 , the measure of the discrepancy between experimental and computed SAXS intensities. Furthermore, this reparametrized force field yielded improved agreement between experimental and simulated membrane thicknesses for the pure lipid systems. However, we note limitations in transferability and diagnose the source of such limitations, particularly the need for SANS in addition to SAXS.

Demerdash, Omar N. A. [Oak Ridge National Laborato↗

DOE OSTI2024

Sequential Design of Experiments for Pilot Testing of Novel Solvent System

The CCSI2 program is supporting a six-month test campaign at the National Carbon Capture Center (NCCC) for evaluation of a novel water-lean solvent. This presentation describes CCSI2’s efforts in process modeling of the solvent system for both coal and natural gas-based flue gas sources and initial uncertainty quantification (UQ) work to estimate parametric uncertainty in key sub-models of interest (e.g., thermodynamics, mass transfer, reaction kinetics). Moreover, perspective is provided on how UQ and sequential design of experiments (SDoE) tools are used to assess the impact of model uncertainty on projected process performance, use this information to optimize data collection during the campaign, and refine process models through data collection. This framework is expected to reduce the overall model uncertainty, and thus risk associated with scale-up as the process moves towards commercialization.

Morgan, Joshua↗

DOE OSTI2025

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗