Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modelling Capabilities”

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 901 records · Page 50

A coupled hydrodynamic (HEC-RAS 2D) and water quality model (WASP) for simulating flood-induced soil, sediment, and contaminant transport

Increased intensity and frequency of floods raise concerns about the release and transport of contaminated soil and sediment to and from rivers and streams. To model these processes during flooding events, we developed an External Coupler in Python to link the Hydrologic Engineering Center-River Analysis System (HEC-RAS) 2D hydrodynamic model to the Water Quality Analysis Simulation Program (WASP). Accurate data transfer from a hydrodynamic model to a water quality model is critical. Our test results showed the External Coupler successfully linked HEC-RAS and WASP and addressed technical challenges in aggregating flow data and conserving mass during the flood event. We ran the coupled models for a 100-year flood event to calculate flood-induced transport of sediment-associated arsenic in Woodbridge Creek, NJ. Change in surface sediment and arsenic at the end of 48-h flood simulation ranged from a net loss of 13.5 cm to a net gain of 11.6 cm, and 16.2 to 2.9 mg/kg, respectively, per model segment, which demonstrates the capability of the coupled model for simulating sediment and contaminant transport in flood.

54 ENVIRONMENTAL SCIENCES↗

Leveraging Artificial Intelligence to Predict Novel Eutectic Alloys

The goal of this project was to train an artificial neural network (ANN) to predict the fractional composition and melting point of eutectic alloys using fundamental atomic properties as inputs. The fundamental properties considered include atomic number, atomic weight, atomic radius, valence electron concentration, electronegativity, and electron affinity. The project involved several phases, starting with data preparation, where phase diagram data was harvested from the ASM International database. Approximately 1300 binary eutectics were collected and cleaned to ensure relevance and accuracy. A regression model was selected for training, utilizing a rectified linear unit as the activation function. Various model configurations were evaluated for predictive accuracy, with validation techniques employed to ensure robustness. The model demonstrated predictive capabilities above random guessing and was able to achieve up to 11% accuracy under certain conditions. An ablative test identified atomic radius and valence electron concentration as critical inputs for model performance. Incorporating the melting point of atomic constituents improved accuracy significantly, although ultimately the model’s predictive capability still fell short of the 80% target. This report details the methodology, results, and implications of the research, contributing to the understanding of employing artificial intelligence to predict the phase transition behavior of eutectic alloys.

36 MATERIALS SCIENCE↗

REopt Lite User Manual

REopt Lite evaluates the economic viability of grid-connected solar photovoltaics, wind, combined heat and power, and electric and thermal storage at commercial and small industrial sites. It allows building owners to identify the system sizes and dispatch strategies that minimize the site’s life cycle cost of energy. REopt Lite also estimates the amount of time on-site generation and storage can sustain the site's critical load during a grid outage and allows the user the choice of optimizing for energy resilience. It is primarily used to inform project development decisions and to support research on the factors that drive project feasibility for market development and policy analysis. It is available through a web interface, application programing interface, and open-source code. This user manual provides an overview of the model, including its capabilities and typical applications; inputs and outputs; economic calculations; technology descriptions; and model parameters, variables, and equations. The model is highly flexible and is continually evolving to meet the needs of each analysis. Therefore, this report is not an exhaustive description of all capabilities, but rather a summary of the core components of the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗

First-Principles Evaluation of Proton Hopping in Tetrahedral Oxide Motifs

Proton-conducting oxides (PCOs) are important materials used as ionic conductors for energy conversion technologies. Existing research efforts on PCO optimization and discovery generally focus on complex perovskite-based oxides that require doping and alloying to engineer oxygen deficiency and high proton conductivity. However, the variety of chemical compositions and coordination environments in oxides poses challenges for efficient materials design. In this computational study, we construct a database of simplified motifs to elucidate the relationship between fundamental materials chemistry and proton kinetics. Specifically, we focus on the zincblende crystal structure as a proxy for tetrahedral metal–oxide (M–O) coordination environments. We systematically quantified the effects of cation type, oxidation states, and M–O bond lengths on the proton hopping barrier, and found that strong M–O bonds and metal cations with large and variable oxidation states (e.g., Mo 6+ , V 5+ ) lead to smaller proton hopping barriers. By mapping the candidate cations and their preferred bond geometries onto materials databases such as the Inorganic Crystal Structure Database (ICSD) and Materials Project, we identified real materials containing the corresponding metal–oxide units. In general, we observed good agreement between the calculated proton hopping barriers obtained in real crystal structures and those predicted by our motif database. We also discuss the limitations of our model and possible future extensions to improve its predictive capabilities. Overall, our model provides a first step for the rational design and quick screening of energy-efficient PCOs.

organic↗

Static and dynamic magnetic properties of honeycomb lattice antiferromagnets Na 2 M 2 TeO 6 , M = Co and Ni

The magnetic structures and spin dynamics of Na 2 CO 2 TeO 6 and Na 2 Ni 2 TeO 6 are investigated by means of elastic and inelastic neutron scattering measurements and the results are discussed in the context of a generalized Kitaev-Heisenberg model on a honeycomb lattice with strong spin-orbit coupling. The large number of parameters involved in the Hamiltonian model is evaluated by using an iterative optimization algorithm capable of extracting model solutions and simultaneously estimating their uncertainty. The analyses indicate that both Co 2+ (d 7 ) and Ni 2+ (d 8 ) antiferromagnets realize bond-dependent anisotropic nearest-neighbor interactions and support the theoretical predictions for the realization of Kitaev physics in 3 d electron systems with effective spins S = 1/2 and 1. By studying the Na-doped system Na 2.4 Ni 2 TeO 6 , we show that the control of Na content can provide an effective route for fine-tuning the magnetic lattice dimensionality, as well as to controlling the bond-dependent anisotropic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of Rapid Permafrost Thaw on CO2 and CH4 Fluxes in a Warmer and Wetter Future (Final Technical Report)

When ice-rich permafrost thaws, the ground subsides, creating thermokarst landscapes with dramatically different soil conditions and carbon fluxes than the original ecosystem. Roughly 20% of the northern permafrost region is susceptible to thermokarst formation (Olefeldt et al., 2016). Thermokarst formation is often rapid; tens of meters of permafrost can thaw within a few years (Schuur et al., 2015). In topographically low areas, thermokarst thaw converts boreal forest or tundra dry shrub ecosystems into sedge or Sphagnum moss wetlands (Olefeldt et al., 2016). While this type of landscape transformation releases carbon stored in permafrost into the atmosphere, on longer time scales, it facilitates sequestration of atmospheric carbon in plant biomass because permafrost thaw releases plant-available nutrients and wetlands are highly productive (M. C. Jones et al., 2017). However, wetlands also generate methane, which is a potent greenhouse gas. Methane emissions from thermokarst wetlands can cause these carbon-sequestering systems to have a positive global warming potential (Johansson et al., 2006; Turetsky et al., 2007). Our project objective was to improve Earth System and environmental predictability by advancing understanding of how CO 2 and CH 4 flux in permafrost thaw-induced wetlands (thermokarst) will change in the future as temperatures and climate conditions shift. Northern latitudes are expected to get warmer and wetter (IPCC 2013), and initiation and expansion of thermokarst thaw is expected to increase (Jorgenson et al. 2006; Zhang et al. 2017). Given these expected changes, our work sought to address three broad questions: Q1) How will northern latitude CO2 and CH4 emissions respond to warming temperatures? Q2) What is the impact of precipitation on permafrost thaw and carbon emissions? Q3) How do CO2 and CH4 emissions change as wetlands age after permafrost thaw? To answer these questions, we took both a modeling and measurement approach. Modeling work was conducted with DOE’s Energy Exascale Earth System Model (E3SM) land model (ELM). Empirical work took place in two primary locations. The first was a thermokarst site near Fairbanks, AK. The site is part of the Bonanza Creek Long Term Ecological Research program and is well instrumented (Neumann et al., 2019). The second was an isolated thawing permafrost wetland located on Kenai Peninsula — Brown’s Lake bog (B. M. Jones et al., 2016) — where the current climate is representative of what is expected at higher latitudes in the future. This site provides an ideal opportunity to test our hypotheses about how thermokarst wetland dynamics will respond to future environmental conditions. In addition, the project collaborated with researchers asking similar questions who were collecting measurements in high latitude post-glacial lakes located in collaborators tackling similar questions in Sweden (Stordalen Mire) (Emerson et al., 2021). Project efforts directly aligned with the stated goal of the funding opportunity announcement (FOA), which was “to improve the understanding and representation of terrestrial ecosystems in ways that advance Earth system model parameterizations and capabilities... thereby improving the quality of Earth and environmental model projections and providing the scientific foundation needed to support DOE’s science and energy missions.” Specifically, the project improved sophistication and accuracy of the Energy Exascale Earth System Model (E3SM), which is being developed primarily at DOE National Laboratories to support scientific research and decision-making.

54 ENVIRONMENTAL SCIENCES↗

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗

Non-Nuclear Test Facility Development to Support Verification and Validation Needs

This study focuses on developing a non-nuclear test facility to help meet the verification and validation (V&V) requirements of a water-cooled small modular reactor (SMR) system. The scaling analysis necessary for reactor system experimentation and assessment model V&V considers various reactor system accident scenarios such as main steam line break (MSLB), steam generator tube rupture (SGTR), and containment condensation and heat transfer. The goal is to identify key physics based phenomena of interest (POI), specific figures of merit (FOMs), and the state of knowledge (SOK) by employing phenomena identification and ranking table (PIRT) studies. Important physics phenomena such as heat transfer, fluid dynamics, structural integrity, and material selection must be addressed when scaling a thermal-hydraulics system. Alongside these considerations, modifications to the facility may be needed to align the scaled system with the available footprint (height and space) of the facility. These modifications may include infrastructure upgrades, changes to the piping layout, implementation of safety measures, and installation of new instrumentation and control systems. The test facility aims to provide a scaled-down version of prototypic facilities for comprehensive V&V activities by overcoming these challenges and making the necessary modifications. This approach minimizes scaling distortion, simulates prototypic plant conditions, and accurately models accident event progressions. Thus, the facility will be capable of adequately modeling a scaled SMR system to support the design, development and licensing needs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

COUPLED MULTIPHYSICS PRIMARY LOOP SIMULATIONS OF THE Mk1-FHR IN THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. A coupled multiphysics model of the Mk1 pebble-bed fluoride-salt-cooled high-temperature reactor (PB-FHR) is presented here. The analysis leverages NEAMS tools (Griffin), SAM, Pronghorn, and the MOOSE heat conduction module) for core neutronics, thermal hydraulics of the core and the primary loop, and multiscale fuel performance simulations. The analysis was entirely created by coupling standalone simulations of the reactor previously available on the VTB. All input files and documentation developed for this example are available on the VTB website: \url{mooseframework.inl.gov/virtual_test_bed/}}. This model was featured in the NRIC Tech Talk presented in December 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-Temperature Latent Heat Storage System using Transportable Heat Pipes for Versatile Integration with Emerging Microreactors

With the growing interest in the development of micro nuclear reactors and their application to microgrid, high-temperature thermal energy storage (TES) has also been receiving great attention as a technology to reinforce the energy utilization and profitability of this class of reactors. Given the importance and research needs of the TES technology, a research project is underway at Idaho National Laboratory (INL) to develop high-temperature TES that can be easily interfaced with emerging micro nuclear reactors. Currently under development is the Heat pipe-Integrated Thermal Battery (HITB), a novel latent heat storage design that enables flexible and robust integration with any design of currently considered micro nuclear reactors. The development of detailed design based on the proposed TES concept, HITB, requires not only development and validation of modeling and simulation capabilities, but also research that spans thermal performance analysis, component- and system-scale design optimization, and experimental demonstrations including high-temperature feasibility testing. The current research efforts (and plans) are divided into the following four areas: (i) Thermal and structural modeling and analysis at both system- and component-scales for evaluating performance and conducting design optimization, (ii) Small-scale experimental demonstrations for the key components of the proposed HITB design, (iii) Modelica-based dynamic model development for developing a control strategy under dynamic operation within integrated energy systems configuration, (iv) Integrated experimental testing and demonstration. This paper introduces the overall research efforts that have been and are being conducted to develop and demonstrate the novel design of the high-temperature latent heat storage device, HITB. The key design features, ongoing research efforts including the development of modeling and simulation capabilities, small-scale demonstrations of the HITB’s key functionalities prior to scaled-up experimental demonstration, are discussed. Then, this paper concludes with a remark on the future research plans for this project.

25 ENERGY STORAGE↗

Coupled Decay Heat and Thermal Hydraulic Capability for Loss-of-Coolant Accident Simulations

As the nuclear energy industry considers ways to achieve improved economics in the current fleet of light-water reactors (LWRs), one possible approach is to operate each cycle for longer durations. This causes a greater portion of the fuel to be burned and reduces the frequency of outages, which ultimately reduces the cost to operate the reactor. However, this also leads to higher burnup fuels than has traditionally been allowed in these reactors. Thus, there are concerns about integrity of high-burnup (HBu) fuel, especially during accident conditions such as loss-of-coolant accidents (LOCAs), as shown by Capps et al.. To investigate these concerns, advanced modeling and simulation capabilities are being leveraged to determine the susceptibility of HBu fuel to fuel fragmentation, relocation, and dispersion (FFRD). Improvements have previously been made to fuel performance capabilities to more accurately model these phenomena; multiphysics simulations have also been conducted to determine the power and burnup histories of the HBu fuel, which are needed as inputs for the fuel performance calculations. Most recently, new statistical approaches have been developed to identify a subset of fuel rods that have greater FFRD susceptibility, reducing the total number of fuel performance simulations required. Prior LOCA simulations have relied on the TRACE systems code, which can model the core and primary loop during accident conditions. TRACE includes many models for various aspects of the primary loop, but two sets of models are important for this report. First, TRACE uses a lumped-fuel approach for modeling the core. This approximates the ~50,000 fuel rods in the core with a much smaller number of rods. The rods can be lumped in various ways as determined by the user. For example, one lumped rod may be used to represent all rods in an assembly, sometimes with an additional rod representing the hottest fuel rod. However, due to runtime constraints and complexity of modeling, a more common approach is to group several assemblies or larger regions of the core into single lumped rods. These lumping schemes apply not only to fuel rods but to flow channels as well. Second, TRACE has several different models for treating decay heat, ranging from pregenerated decay heat curves based on an ANSI/ANS-5.1 standard (hereinafter abbreviated simply as ANSI) to explicit time-dependent heat inputs from the user. None of these models account for differences in isotopics between different rods, which is an approximation the work in this report seeks to eliminate. This report focuses on the implementation of coupled decay heat capabilities in the Virtual Environment for Reactor Applications (VERA) code suite to address a gap identified in previous LOCA simulations. This constitutes an improvement for both the lumped-fuel and decay heat models in TRACE. VERA has been developed to perform high-fidelity, whole-core multiphysics simulations for LWRs. Previously, during the Consortium for Advanced Simulation of LWRs (CASL) program, the emphasis was on providing accurate steady-state analysis—with a secondary focus on reactivity insertion accident (RIA) analysis—to address operational challenges in the nuclear energy industry. Under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, these capabilities are being extended to a broader range of transient analyses with the goal of quantifying the risk of fuel failures such as FFRD. To properly model such conditions with VERA, decay heat calculations have been integrated with the multiphysics to enable rod-by-rod thermal hydraulic (TH) conditions to be driven by the decay heat in long-running accidents such as LOCAs.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗

High Fidelity Interface Modeling to Enable Enhanced Electromagnetics ModSim Capabilities

Years of work by 1350 and others has shown that the phenomenology behind EM penetration of joints and seams is a major driver in the shielding effectiveness of ND systems. Via analysis of a canonical cylindrical geometry and comparison against experimental data, we will provide evidence supporting the theory that proper treatment of contact phenomenology including joint deformation, asperity-induced contact impedance, and appropriate treatment of machining tolerance values is required to match electromagnetics modeling and simulation results to experimental data.

36 MATERIALS SCIENCE↗

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗

Robust Machine Learning Inference from X-ray Absorption Near Edge Spectra through Featurization

X-ray absorption spectroscopy (XAS) is a commonly employed technique for characterizing functional materials. In particular, X-ray absorption near edge spectra (XANES) encode local coordination and electronic information, and machine learning approaches to extract this information are of significant interest. To date, most ML approaches for XANES have primarily focused on using the raw spectral intensities as input, overlooking the potential benefits of incorporating spectral transformations and dimensionality reduction techniques into ML predictions. Here, in this work, we focused on systematically comparing the impact of different featurization methods on the performance of ML models for XAS analysis. We evaluated the classification and regression capabilities of these models on computed data sets and validated their performance on previously unseen experimental data sets. Our analysis revealed an intriguing discovery: the cumulative distribution function feature achieves both high prediction accuracy and exceptional transferability. This remarkably robust performance can be attributed to its tolerance to horizontal shifts in the spectra, which is crucial when validating models using experimental data. While this work exclusively focuses on XANES analysis, we anticipate that the methodology presented here will hold promise as a versatile asset to the broader spectroscopy community.

36 MATERIALS SCIENCE↗

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗