Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis of Knowledge”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Science and Solutions for a Planet Under Pressure: Synthetic Research, Analytics, and Community Engagement on Five Priority Theme Areas

Support from the Department of Energy proved instrumental in supporting the work of Global Council for Science and the Environment (GCSE) on synthetic research, analysis, and community engagement, oriented around the themes and discussions from the GCSE Drawdown 2021 Conference — Research to Action: Science and Solutions for a Planet Under Pressure. These efforts helped inform and enhance the societal relevance of the knowledge that has been generated through the conference and, at the same time, advance the scientific insights acquired through synthesis and scientific translation efforts. The work centered on three main objectives: 1) documentation of findings from the completed GCSE 2021, culminating in the publication of a synthesis report; 2) development and vetting of integrated use cases that illustrate synthesis methods and enable scientific outreach, communications, and translational discussions among universities and the broader research community; and 3) strengthening the scientific basis for energy education research by sharing findings in university-educator fora.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

First-Principles Simulation of Carrier Recombination Mechanisms in Halide Perovskites

In recent years, there have been remarkable developments in halide perovskites, which are used in highly efficient optoelectronic devices and exhibit intriguing materials physics. Detailed knowledge of carrier recombination mechanisms is essential for understanding their excellent performance and to further increase their efficiency. Obtaining such knowledge is challenging however, and different studies have reached divergent conclusions in some cases. This progress report outlines the critical developments in understanding the carrier recombination mechanisms in halide perovskites from a computational perspective. The primary focus is radiative and Auger recombination, since they have not been systematically assessed and discussed before, and a number of important issues have been actively debated. This comprehensive discussion of the carrier recombination mechanisms is aimed at establishing physically justified insights that can form the basis for better materials and devices design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

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

Codiscovering graphical structure and functional relationships within data: A Gaussian Process framework for connecting the dots

Most problems within and beyond the scientific domain can be framed into one of the following three levels of complexity of function approximation. Type 1: Approximate an unknown function given input/output data. Type 2: Consider a collection of variables and functions, some of which are unknown, indexed by the nodes and hyperedges of a hypergraph (a generalized graph where edges can connect more than two vertices). Given partial observations of the variables of the hypergraph (satisfying the functional dependencies imposed by its structure), approximate all the unobserved variables and unknown functions. Type 3: Expanding on Type 2, if the hypergraph structure itself is unknown, use partial observations of the variables of the hypergraph to discover its structure and approximate its unknown functions. These hypergraphs offer a natural platform for organizing, communicating, and processing computational knowledge. While most scientific problems can be framed as the data-driven discovery of unknown functions in a computational hypergraph whose structure is known (Type 2), many require the data-driven discovery of the structure (connectivity) of the hypergraph itself (Type 3). We introduce an interpretable Gaussian Process (GP) framework for such (Type 3) problems that does not require randomization of the data, access to or control over its sampling, or sparsity of the unknown functions in a known or learned basis. Its polynomial complexity, which contrasts sharply with the super-exponential complexity of causal inference methods, is enabled by the nonlinear ANOVA capabilities of GPs used as a sensing mechanism.

Science & Technology - Other Topics↗
DOE OSTI2019

A data-driven framework for sparsity-enhanced surrogates with arbitrary mutually dependent randomness

The challenge of quantifying uncertainty propagation in real-world systems is rooted in the high-dimensionality of the stochastic input and the frequent lack of explicit knowledge of its probability distribution. Traditional approaches show limitations for such problems, especially when the size of the training data is limited. To address these diculties, we have developed a general framework of constructing surrogate models on spaces of stochastic input with arbitrary probability measure irrespective of the mutual dependencies between individual components of the random inputs and the analytical form. The present Data-driven Sparsity-enhancing Rotation for Arbitrary Randomness (DSRAR) framework includes a data-driven construction of multivariate polynomial basis for arbitrary mutually dependent probability measure and a sparsity enhancement rotation procedure. This sparsity-enhancing rotation method was initially proposed in our previous work [1] for Gaussian density distributions, which may not be feasible for non-Gaussian distributions due to the loss of orthogonality after the rotation. To remedy such diculties, we developed a new data-driven approach to construct orthonormal polynomials for polynomials for arbitrary mutually dependent (amdP) randomness, ensuring the constructed basis maintains the orthogonality/near-orthogonality with respect to the density of the rotated random vector, where directly applying the regular polynomial chaos including arbitrary polynomial chaos (aPC) [2] shows limitations due to the assumption of the mutual independence between the components of the random inputs. The developed DSRAR framework leads to accurate recovery, with only limited training data, of a sparse representation of the target functions. The eectiveness of our method is demonstrated in challenging problems such as PDEs and realistic molecular systems within high-dimensional conformational space (O(10)) where the underlying density is implicitly represented by a large collection of sample data, as well as systems with explicitly given non-Gaussian probabilistic measures.

Lei, Huan↗
DOE OSTI2021

Sensitivity Analysis of Particle-In-Cell Modeling Parameters in Settling Bed, Bubbling Fluidized Bed and Circulating Fluidized Bed

The objective of the work presented is to perform a preliminary sensitivity analysis of particle-in-cell (PIC) model parameters when applied to settling bed, bubbling fluidized bed, and circulating fluidized bed simulations. These examples correspond to widely different flow conditions commonly seen in chemical engineering applications. Simulations were performed using the PIC method in the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). As part of the non-intrusive uncertainty quantification (UQ) analysis, simulation campaigns were generated using Nodeworks. Sampling locations or settings for PIC model parameters were determined using the Latin Hypercube method. Response surfaces were created using radial basis functions (RBF), and Sobol’ indices were estimated to quantify the influence of model parameters on the quantities of interest (QoI). This study marks a first step towards systematically determining optimal ranges for model parameters used in MFiX-PIC. Based on limited experience, it is expected that these values would depend strongly on flow conditions. Given the complexity of the multiphase flow systems under analysis, a non-intrusive UQ based approach is used to identify the most influential parameters in each case. This prior knowledge will help in proposing an effective design of experiments (DoE) and determine optimal parameters through techniques such as deterministic or Bayesian calibration, which will be pursued in the future.

42 ENGINEERING↗
DOE OSTI2019

Designing a Site Closure Strategy for Uranium in Groundwater and Surface Water - 19534

A site closure strategy was designed and mapped for residual uranium present in groundwater and surface water at the US DOE Office of Legacy Management (LM) Monticello Mill Tailings Site (MMTS) Operable Unit (OU) III, which is regulated under CERCLA. The MMTS is a former mill site that was used to process uranium and vanadium ore in southeastern Utah until 1960. OU III includes surface water in Montezuma Creek and groundwater present in a shallow alluvial aquifer along the creek valley for several kilometers downgradient of the former mill site. The current remedy consists of an operating pump-and-treat system, monitored natural attenuation (MNA), and institutional controls (ICs). Performance evaluations indicate that a prolonged restoration period is likely, as well as an asymptotic trend in uranium concentrations that still exceed the remediation goal. Multiple strategies and scenarios were evaluated before selecting a preferred site closure strategy. Recommendations were developed to support advancing the site towards closure over the next few years and guide the process of regulatory interactions and decision making at the MMTS OU III. Site closure is an important milestone that reflects successful remedial progress at LM sites. Post-closure site management would be maintained to ensure conditions are protective of human health and the environment. The site closure strategy was considered in the context of site remedial progress. The primary sources of groundwater contamination have been removed, including mill tailings and contaminated soils, sediment, and debris. Five-year reviews have regularly concluded that the remedy is protective of human health and environment and that ICs are effective in preventing exposure to contaminated groundwater. The affected aquifer is low-yielding and has no history of domestic use. In 2004, the final remedy selected in the Record of Decision (ROD) was MNA with ICs. However, a pump-and-treat system was installed and continues to operate as a contingency remedy in the area with the highest residual uranium concentrations. The need for a contingency pump-and-treat system was triggered by a departure of the long-term monitoring data from model projections that were used as the basis for remedy performance metrics specified in the ROD. In 2015, the pump-and-treat system was expanded. Average uranium concentrations have declined locally over time. However, monitoring data and column studies indicate that the system effectiveness will eventually become ineffective in reducing uranium concentrations, likely due to the reequilibration of uranium present in the solid matrix in various states (i.e., sorbed, precipitated). Currently, there are no agreed-upon metrics for discontinuing the contingency pump-and-treat system. Three closure strategies effectively applied at other CERCLA sites were considered at MMTS OU III: (1) restoring natural attenuation conditions and conducting MNA; (2) applying alternate concentration limits (ACLs); and (3) waiving a remedial goal based on technical impracticability (TI). Each closure strategy was evaluated in the following manner: (1) identifying the criteria and requirements associated with current regulations, policy, and guidance, (2) assessing the ability to meet those criteria based on the current site knowledge and the conceptual site model, and (3) identifying data gaps that need to be filled to further assess the closure strategy. Following this evaluation, a preferred closure strategy was selected: to transition to a remedy based on MNA and maintenance of ICs, without active remediation. Several scenarios for transitioning to the preferred closure strategy were identified. Each scenario consisted of 'if/then logic' associated with a larger decision diagram. The decision diagram considered the current conceptual site model (CSM) and potential outcomes of studies that will be conducted to fill data gaps, as well as the performance of the existing treatment system and long-term monitoring data. Two scenarios for transitioning to the preferred closure strategy were highlighted for future planning and were described in detail. For both scenarios, the preferred closure strategy includes demonstrating criteria for MNA acceptance, pump-and-treat termination, long-term monitoring, and eventual site delisting. Recommended actions over the next few years were also developed, including updating the CSM, conducting numerical modeling to predict remedial performance and time frames, conducting geochemical analyses to assess uranium fate and transport, and summarizing the lines of evidence and basis for MNA. This detailed example of the development of a site-specific closure strategy provides the thought process and decision logic used to assess uncertainty using an if/then framework; identify common elements and recommended actions to advance the site towards closure; and navigate the inherent uncertainty in the timing of future events, outcome of regulatory discussions, and remedial progress over time that is common to many complex US DOE facilities. In addition, the national regulations, policy, guidance, and state of practice regarding commonly used site closure strategies was summarized for consideration, including MNA, ACLs, and TI. (authors)

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

Experimental investigation of abnormal transverse flow enhancement of α particles in heavy-ion collisions

The mass dependence of the transverse flow for Z = 1–5 fragments from the collisions of 40 Ar + 27 Al, 40 Ar + 48 Ti, and 40 Ar + 58 Ni at 47 MeV/nucleon is investigated experimentally in this article. The transverse flow values are determined using the in-plane components of the fragment transverse momenta, where three conventional methods, i.e., the kinetic flow tensor method, the transverse momentum analysis method, and the azimuthal correlation method, are applied to reconstruct the reaction plane in an event-by-event basis. It is demonstrated from the comparison of the present experimental mass dependent flow measurements and the model simulations using an improved antisymmetrized molecular dynamics model that the experimentally observed abnormal α transverse flow enhancement is closely related to the reaction plane reconstruction procedure in the flow extraction. We further investigate the physical existence of the abnormal α flow behavior using a two-particle azimuthal correlation method, which allows us to provide the relative flow magnitude information with an identification of fragment charge number without the knowledge of the reaction plane differing from the three conventional methods. It is found that the relative flow magnitudes deduced from the two-particle azimuthal correlation functions with an identification of Z, with the correction for the recoil effect imposed by the momentum conservation, show a monotonically increasing trend as a function of fragment charge number, with no exception of the α flow enhancement. These results, in addition to those from the improved antisymmetrized molecular dynamics model simulations, definitely provide experimental evidences for the inexistence of the abnormal α flow behavior in the heavy-ion collisions at the present incident energy region in nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Harmonizing model organism data in the Alliance of Genome Resources

The Alliance of Genome Resources (the Alliance) is a combined effort of 7 knowledgebase projects: Saccharomyces Genome Database, WormBase, FlyBase, Mouse Genome Database, the Zebrafish Information Network, Rat Genome Database, and the Gene Ontology Resource. The Alliance seeks to provide several benefits: better service to the various communities served by these projects; a harmonized view of data for all biomedical researchers, bioinformaticians, clinicians, and students; and a more sustainable infrastructure. The Alliance has harmonized cross-organism data to provide useful comparative views of gene function, gene expression, and human disease relevance. The basis of the comparative views is shared calls of orthology relationships and the use of common ontologies. The key types of data are alleles and variants, gene function based on gene ontology annotations, phenotypes, association to human disease, gene expression, protein–protein and genetic interactions, and participation in pathways. The information is presented on uniform gene pages that allow facile summarization of information about each gene in each of the 7 organisms covered (budding yeast, roundworm Caenorhabditis elegans, fruit fly, house mouse, zebrafish, brown rat, and human). The harmonized knowledge is freely available on the alliancegenome.org portal, as downloadable files, and by APIs. We expect other existing and emerging knowledge bases to join in the effort to provide the union of useful data and features that each knowledge base currently provides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

NGEE Arctic Integrated Modeling (IM2): Improved subgrid hillslope hydrologic connectivity

This data product represents the integration of new code capability for arctic tundra hillslope hydrologic processes into the Energy Exascale Earth System Model (E3SM), through the E3SM Land Model (ELM) component. This code integration is the result of collaborative effort between the NGEE Arctic project and the E3SM project. The current ELM represents water movement primarily through vertical processes, such as precipitation, canopy interception, evaporation, infiltration, and soil water movement. Lateral water movement—such as surface runoff, subsurface flow, and river transport—plays a significant role in the hydrological cycle, especially in regions with varied topography. While E3SM includes a runoff routing component representing water transport in the river network, the lateral transport of water at the subgrid scale within the land model has previously not been taken into account. With the recent development of topographic units within the ELM subgrid data structure, there is an opportunity to simulate hillslope hydrologic connectivity by introducing water transport along topographic gradients. We expect that more realistic representation of hillslope hydrologic processes will lead to improved predictions of both soil water content and river network flows. Lateral transport of water at and near the surface is represented as a sub-grid process in this new code development. Water is tracked as it moves from higher to lower elevations within a gridcell. This capability uses the nested hierarchical sub-grid scheme within ELM to connect water fluxes from sub-grid elements with higher elevation to those with lower elevation. This data record consists of a single document (pdf format) that describes the theoretical basis for the hillslope hydrology processes added to ELM, and describes the modifications made to the ELM code. The Methods section of this metadata record includes a link to the public E3SM code repository where the exact code modifications as integrated in E3SM can be accessed. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

Thornton, Peter E [ORNL] (ORCID:0000000247595158)↗
DOE OSTI2020

Experimental Investigation of Low Cost, Low Thermal Conductivity Thermal Barrier Coating on HCCI Combustion, Efficiency, and Emissions

In-cylinder surface temperature is of heightened importance for Homogeneous Charge Compression Ignition (HCCI) combustion since the combustion mechanism is thermo-kinetically driven. Thermal Barrier Coatings (TBCs) selectively manipulate the in-cylinder surface temperature, providing an avenue for improving thermal and combustion efficiency. A surface temperature swing during combustion/expansion reduces heat transfer losses, leading to more complete combustion and reduced emissions. At the same time, achieving a highly dynamic response sidesteps preheating of charge during intake and eliminates the volumetric efficiency penalty. The magnitude and temporal profile of the dynamic surface temperature swing is affected by the TBC material properties, thickness, morphology, engine speed, and heat flux from the combustion process. This study follows prior work of authors with Yttria Stabilized Zirconia, which systematically engineered coatings for HCCI combustion. Herein, a modeling study was used to assess the impacts of various TBC material properties, e.g. density, thickness, and thermal conductivity on the temperature swing effect. Reducing conductivity emerged as a most promising avenue, rather than reducing both the density and effective conductivity by increasing porosity, the current work emphasizes a material with natively low conductivity. A novel TBC formulation was developed, leverages a class of materials that, to the author’s best knowledge, have not been used as a thermal barrier coating previously. Here, a systematic experimental investigation was carried out using single-cylinder research engine. Experimental engine studies utilizing the novel ‘glassy’ low-K coating exhibit advanced ignition phasing and reduced combustion duration relative to the baseline engine. Heat transfer measurements indicate a net reduction in heat flux over the entire cycle, although main effect is felt during expansion, and the reduced heat transfer losses manifest in higher gross indicated thermal efficiency by approximately 5-6% on a relative basis.

42 ENGINEERING↗
DOE OSTI2020

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗
DOE OSTI2021

A multi-dimensional parametric study of variability in multi-phase flow dynamics during geologic CO 2 sequestration accelerated with machine learning

Successful geologic CO 2 storage projects depend on numerical simulations to predict reservoir performance during site selection, injection verification, and post-injection monitoring phases of the project. These numerical simulations solve non-linear sets of coupled partial differential equations, while accounting for multi-phase fluid dynamics on the basis of constitutive equations that are embedded into the solution scheme. As a consequence, individual simulations often require tens to hundreds of hours to complete on high-performance computing clusters. Moreover, laboratory experiments reveal that parametric functions for capillary pressure and relative permeability exhibit substantial variability, even within the same rock type. This combination of computational expense and wide-ranging parametric variability means that there remains substantial uncertainty in the behavior of multi-phase CO 2 -water systems, particularly in the context of feedbacks between relative permeability and capillary pressure. To bridge this knowledge gap, here we develop a novel workflow that utilizes physics-based numerical simulation to train an artificial neural network (ANN) emulator for interrogating the multivariate parameter space that governs both capillary pressure and relative permeability. With this approach, the ANN is trained to emulate both fluid pressure distribution and CO 2 saturation, which are then interrogated quantitatively to generate parametric response surface mappings with high-fidelity resolution. Results from this study initially show that capillary entry pressure is the dominant control on both CO 2 plume geometry and fluid pressure propagation when considering the combined effects of capillary pressure and relative permeability, particularly when phase interference is low and residual CO 2 saturation is high. Moreover, the ANN emulator provides tremendous computational speed-up by computing 2691 individual simulations in several minutes; whereas, the same simulation ensemble would have required ~3 years of simulation time using only physics-based simulation methods (25,000 times speed up).

58 GEOSCIENCES↗
DOE OSTI2020

Modernization of Technical Requirements for Licensing of Advanced Non-Light Water Reactors: Selection and Evaluation of Licensing Basis Events

This report, “Modernization of Technical Requirements for Licensing of Advanced Non-Light Water Reactors: Selection and Evaluation of Licensing Basis Events,” represents a key element in the development of a methodology for the efficient licensing of advanced non-light water reactors (non-LWRs). It is the result of a Licensing Modernization Project (LMP) led by Southern Company and cost-shared by the U.S. Department of Energy (DOE). The LMP will result in detailed proposals for establishing licensing technical requirements to facilitate risk-informed and performance-based design and licensing of advanced non-LWRs. Such a methodology acknowledges enhancements in safety achievable with advanced designs and reflects more recent states of knowledge regarding safety and design innovation, creating an opportunity for reduced regulatory complexity with increased levels of safety. The project builds on best practices as well as previous activities through DOE and industry-sponsored advanced reactor licensing initiatives.

42 ENGINEERING↗
DOE OSTI2022

IRSN Work in Support of NCSP: FY2021 and Plan for Future [Slides]

IRSN and NCSP collaboration has been very fruitful. Despite challenges presented by COVID-19, work has been performed regularly and successfully by video conference. This includes monthly basis meeting with LLNL and LANL and specifics meetings with ORNL and SNL when needed. After two years without any IRSN participation to experiments, ISRN is looking forward to fully participating in the three proposed integral experiments scheduled in 2022. IRSN is working on several new projects including TEX-MOX and TEX-low temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scaling High-Resolution Soil Organic Matter Composition to Improve Predictions of Potential Soil Respiration Across the Continental United States

Despite the importance of microbial soil organic matter (SOM) respiration in regulating the flux of carbon between soils and the atmosphere, soil carbon cycling models remain primarily based on climate and soil properties, leading to large uncertainty in predictions. To address this knowledge gap, we analyzed high-resolution water-extractable SOM profiles from soil cores collected across the United States by the 1,000 Soils Pilot of the Molecular Observation Network. Our innovation lies in using machine learning to distill thousands of SOM formula into tractable units; and it enables integrating data from molecular measurements into soil respiration models. In surface soils, SOM chemistry provided better estimates of potential soil respiration than soil physicochemistry, and using them combined yielded the best prediction. Overall, we identify specific subsets of organic molecules that may improve predictions of global soil respiration and create a strong basis for developing new representations in process-based models.

54 ENVIRONMENTAL SCIENCES↗