SEARCH · Engineering Papers

Results for “Monte Carlo integration”

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 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A Laboratory-driven Multiscale Investigation of X-Ray Induced Mass Loss and Photochemical Evolution in Cosmic Carbon and Silicate Dust

We present the results of an integrated laboratory and modeling investigation into the impact of stellar X-rays on cosmic dust. Carbonaceous grains were prepared in a cooled (< 200 K) supersonic expansion from aromatic molecular precursors, and were later irradiated with 970 eV X-rays. Silicate (enstatite) grains were prepared via laser ablation, thermally annealed, and later irradiated with 500 eV X-rays. Infrared spectra of the 3.4 μm band of the carbon sample prepared with benzene revealed 84% ± 5% band area loss for an X-ray dose of 5.2 x 10 23 eV.cm -2 . Infrared spectra of the 8 -12 μm Si-O band of the silicate sample revealed band area loss up to 63% ± 5% for doses of 2.3 x 10 23 eV.cm -2 . A hybrid Monte Carlo particle trajectory approach was used to model the impact of X-rays and ensuing photoelectrons, Auger and collisionally ionized electrons through the bulk. As a result of x-ray ionization and ensuing Coulomb explosions on surface molecules, the calculated mass loss is 60% for the carbonaceous sample and 46% for the silicate sample, within a factor of 2 of the IR band loss, supporting an X-ray induced mass-loss mechanism. We apply the laboratory X-ray destruction rates to estimate the lifetimes of dust grains in protoplanetary disks surrounding 1 M ⊙ and 0.1 M ⊙ G and M stars. In both cases, X-ray destruction timescales are short (a few million years) at the disk surface, but are found to be much longer than typical disk lifetimes (≳10 Myr) over the disk bulk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗

DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

ATF-2 Ramp Conceptual Design Report

The following report documents the conceptual design for the ATF-2 Ramp experiment. The experiment aims to fill the in-pile irradiation testing gap to conduct integral ramp testing, which was created by the closure of R2, Osiris, and Halden test reactors. The concept involves the simultaneous ramping of three fuel pins using a power axial locator mechanism in the Loop-2A testing facility in the center flux trap of the advanced test reactor (ATR). The three fuel pins are ramped in individual coolant channels containing a prototypic PWR environment. Rod failure will be detected using a fuel rod elongation sensor attached to the fuel pin upper end cap. An axial stack of concentric hafnium and zirconium shrouds will be used to shape the flux around the test pins to create different power levels in each pin and to ensure the peak power location of each pin remains in the center of the test rod. Monte Carlo simulations are used to demonstrate the viability of this design concept.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Study of quark and gluon jet substructure in Z+jet and dijet events from pp collisions

Measurements of jet substructure describing the composition of quark- and gluon-initiated jets are presented. Proton-proton (pp) collision data at √s = 13 TeV collected with the CMS detector are used, corresponding to an integrated luminosity of 35.9 fb -1 . Generalized angularities are measured that characterize the jet substructure and distinguish quark- and gluon-initiated jets. These observables are sensitive to the distributions of transverse momenta and angular distances within a jet. The analysis is performed using a data sample of dijet events enriched in gluon-initiated jets, and, for the first time, a Z+jet event sample enriched in quark-initiated jets. The observables are measured in bins of jet transverse momentum, and as a function of the jet radius parameter. Additionally, each measurement is repeated applying a “soft drop” grooming procedure that removes soft and large angle radiation from the jet. Using these measurements, the ability of various models to describe jet substructure is assessed, showing a clear need for improvements in Monte Carlo generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Combining Simulation and Optimization to Derive Operating Policies for a Concentrating Solar Power Plant

Optimizing short-term decisions over a rolling horizon and/or using deterministic penalties to capture system stochasticity can lead to myopic policies that fail to consider unplanned events and their long-term adverse effects. We present a methodology that integrates an off-line optimization model with a simulation procedure to determine the profitability of different operating strategies; specifically, the latter is used to generate additional constraints for the former when failures occur according to system component operating lifetimes that (i) are subject to exogenous uncertainty, and (ii) may degrade more quickly under specific operating conditions. We use the feedback provided by the simulation model in a parametric analysis to obtain penalties that can be used in short-term operations scheduling to maximize the long-term revenues obtained by the optimization model. We apply this research to a concentrating solar power plant; our results show that the methodology can be used to choose an operating policy that balances maximizing profit while accounting for maintenance costs. Integrating the optimization model with a simulation procedure reveals that aggressive prices for cycling yield about 55% fewer startups and 30% fewer failures compared to using a more typical start-up operating strategy, and can save hundreds of thousands to millions of dollars in repair costs over the lifetime of the plant.

concentrating solar power↗

DOE OSTI2024

The ICRP, MELODI, and ALLIANCE workshop on effects of ionizing radiation exposure in offspring and next generations: a summary of discussions

Task Group 121 – Effects of ionizing radiation exposure in offspring and next generations – is a task group under the Committee 1 of the International Commission on Radiological Protection (ICRP), approved by the Main Commission on 18th November 2021. The main goals of Task Group 121 are to (1) review and update the scientific literature of relevance to radiation-related effects in the offspring of parent(s) exposed to ionizing radiation in both human and non-human biota; (2) to assess preconceptional and intrauterine effects of radiation exposure and related morbidity and mortality; and, (3) to provide advice about the level of evidence and how to consider these preconceptional and postconceptional effects in the system of radiological protection for humans and non-human biota. The Task Group is reviewing relevant literature since Publication 90 ‘Biological effects after prenatal irradiation (embryo and fetus)’ (2003) and will include radiation-related effects on future generations in humans, animals, and plants. This review will be conducted to account for the health effects on offspring and subsequent generations in the current system of radiological protection. Radiation detriment calculation will also be reviewed. Finally, preliminary recommendations will be made to update the integration of health effects in offspring and next generations in the system of radiological protection. A Workshop, jointly organized by ICRP Task Group 121 and European Radiation Protection Research Platforms MELODI and ALLIANCE was held in Budapest, Hungary, from 31st May to 2nd June 2022. Participants discussed four important topics: (1) hereditary and epigenetic effects due to exposure of the germ cell line (preconceptional exposure), (2) effects arising from exposure of the embryo and fetus (intrauterine exposure), (3) transgenerational effects on biota, and (4) its potential impact on the system of radiological protection. Based on the discussions and presentations during the breakout sessions, newer publications, and gaps on the current scientific literature were identified. For instance, there are some ongoing systematic reviews and radiation epidemiology reviews of intrauterine effects. There are newer methods of Monte Carlo simulation for fetal dosimetry, and advances in radiation genetics, epigenetics, and radiobiology studies. While the current impact of hereditary effects on the global detriment was reported as small, the questions surrounding the effects of radiation exposure on offspring and the next generation are crucial, recurring, and with a major focus on exposed populations. Here, this article summarizes the workshop discussions, presentations, and conclusions of each topic and introduces the special issue of the International Journal of Radiation Biology resulting from the discussions of the meeting.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2020

Origin of mean-field behavior in an elastic Ising model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Precise measurement of the tt¯ production cross-section and lepton differential distributions in eμ dilepton events from s=13TeVpp collisions with the ATLAS detector

The inclusive top quark pair (tt¯$$t\bar{t}$$) cross-section σtt¯$$\sigma _{t\bar{t}}$$ has been measured in proton–proton collisions at s=13TeV$$\sqrt{s}=13\,\text {TeV}$$, using 140fb-1$$140\,{\text {fb}^{-1}} $$ of data collected by the ATLAS experiment at the Large Hadron Collider. Using events with an opposite-charge eμ$$e\mu $$ pair and b-tagged jets, the cross-section is measured to be: σtt¯=829.3±1.3(stat)±8.0(syst)±7.3(lumi)±1.9(beam)pb,$$\begin{aligned} \sigma _{t\bar{t}} & = 829.3\pm 1.3\,\mathrm {(stat)}\ \pm 8.0\,\mathrm {(syst)}\ \pm 7.3\,\mathrm {(lumi)}\ \\ & \quad \pm 1.9\,\mathrm {(beam)}\,\textrm{pb}, \end{aligned}$$where the uncertainties reflect the limited size of the data sample, experimental and theoretical systematic effects, the integrated luminosity, and the proton beam energy, giving a total uncertainty of 1.3%. The result is used to determine the top quark pole mass via the dependence of the predicted cross-section on mtpole$${m_{t}^\textrm{pole}}$$, giving mtpole=172.8-1.7+1.5$${m_{t}^\textrm{pole}}=172.8^{+1.5}_{-1.7}$$ GeV$$\text {GeV}$$. The same event sample is used to measure absolute and normalised differential cross-sections for the tt¯→eμνν¯bb¯$$t\bar{t} \rightarrow e\mu u \bar{ u }b\bar{b} $$ process as a function of single-lepton and dilepton kinematic variables. Complementary measurements of eμbb¯$$e\mu b\bar{b} $$ production, treating both tt¯$$t\bar{t}$$ and Wt events as signal, are also provided. Both sets of differential cross-sections are compared to the predictions of various Monte Carlo event generators, demonstrating that the state-of-the-art generators Powheg MiNNLO and Powhegbb4l describe the data better than Powheghvq.

Aad, G↗

DOE OSTI2024

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

DOE OSTI2022

Application of Monte Carlo code SHIFT for MSBR dose rate calculations

This paper presents a study that applied the Monte Carlo code SHIFT to calculate the radiological environment within a detailed Molten Salt Breeder Reactor (MSBR) model. It represents one of the first applications of the SHIFT code in stand-alone simulations for non-light water reactors to help demonstrate the code's potential uses in the design, licensing, and operation of advanced reactors. The radiological conditions of the MSBR were modeled when the reactor is at two different operation modes: normal full power and a drained state. The Forward Weighted-Consistent Adjoint Driven Importance Sampling (FW-CADIS) hybrid method in SHIFT was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the SHIFT simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

Physics-Informed Gaussian Process Regression for States Estimation and Forecasting in Power Grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method. We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power. For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗

DOE OSTI2020

CALPHAD Uncertainty Quantification and TDBX

CALPHAD uncertainty quantification (UQ) is the foundation of materials design with quantified confidence. We report a framework and software packages to enable CALPHAD UQ assessment and calculation using commercial CALPHAD software (Thermo-Calc). This Bayesian inference framework is coupled with a Markov chain Monte Carlo algorithm to establish uncertainty traces with a given thermodynamic database file (TDB) and corresponding experimental data points. This general framework is demonstrated with the Ni–Cr binary system. The algorithm is firstly validated on synthetic data with known ground truth. Then it is applied to real experimental data to generate posterior traces. We develop a file format named TDBX, which provides a single source of truth by combining the original TDB content and the traces for each assessed Gibbs energy parameter. CALPHAD UQ calculations are performed based on the TDBX file, from which uncertainties for phase boundaries, enthalpy curves, and solidification range are collected as examples of basic design parameters. This TDBX file with corresponding scripts are made open-source. Finally, the combination of CALPHAD UQ assessments and calculations connected by TDBX supports uncertainty-assisted modeling, enabling the integrated application of modern design with uncertainty methodologies to computational materials design.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗