Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “energy use behavior”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Capturing plume behavior in complex terrain: an overview of the Nevada National Security Site Meteorological Experiment (METEX21)

METEX21 was an atmospheric tracer release experiment executed at the Department of Energy’s Nevada National Security Site (NNSS) in the southwestern U.S to study terrain-induced wind and thermodynamic conditions that influence local-scale (<5-km) plume transport under varying atmospheric forcing conditions. Meteorological observations were collected using 10-m tall meteorological towers, 2-m tall tripods with 3-d sonic anemometers, a 3-m tall eddy covariance flux tower, Doppler profiling lidars, Doppler scanning lidars, weather-balloon launched radiosondes, and a tethered balloon equipped with wind, temperature, and aerosol sensors at heights up to 800 m. A smoke tracer was released along three transects in the horizontal and vertical directions and observed with video cameras, aerosol sensors and lidars (via aerosol backscatter). The observations showed evidence of large-scale/synoptic transience as well as local-scale upslope and downslope flows, along-axis valley flows, recirculation eddies on leeward slopes, and periods of strong shear and veer aloft. The release days were classified as either synoptically-driven or locally-driven, and a single case day is presented in detail for each. Synoptically-forced days show relatively narrow smoke plumes traveling down the valley from north to south (with the predominant wind direction), with little deviation in transport direction regardless of the elevation or ground locations of the smoke releases, except near the presence of leeside recirculation eddies. Locally-forced days exhibit a wider range of plume behavior due to the combination of thermally-induced valley and slope flows, which are often flowing in different cardinal directions, and wind shear found aloft at higher altitudes and elevations. We saw evidence of smoke lofting on top of the mesas due to strong upslope flows on these days. A major finding of this experiment was the effectiveness of scanning lidars to measure 2-dimensional plume transport out to a 2–3 km distance; much farther than could be visibly observed. METEX21 was the first of three planned tracer experiments at NNSS, and future experiments will incorporate multiple tracers to improve individual plume identification so that finer resolution flow details can be attained from these measurements, as well as deploy a larger suite of meteorological instrumentation, including more temperature profiling data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Effects of Landau quantization on neutrino emission and absorption

Some neutron stars known as magnetars possess very strong magnetic fields, with surface fields as large as 10 15 G and internal fields that are possibly stronger. Recent observations of the radio pulsar GLEAM-X J1627 suggest it may have a surface field as strong as 10 16 G. In the presence of a strong magnetic field, Leinson and Perez showed and Baiko and Yakovlev confirmed that the energy levels of electrons and protons are quantized and the direct Urca process allows neutron stars to cool rapidly, even at low density. For the case of magnetic fields 𝐵 ≳ 10 16 G, we find features in the emissivity due to energy quantization that are not captured by the frequently employed quasiclassical approximation where energy levels are treated as nearly continuous. Resonances can result in amplification of the neutrino emissivity at specific densities compared to a calculation that neglects quantization, particularly at low temperature. These effects are not important for the thermal evolution of an entire neutron star, but may be relevant for phenomena that depend on behavior at specific densities. Here, we present a fully relativistic calculation of the direct Urca rate in a strong magnetic field using the standard V-A weak Lagrangian incorporating mean field nuclear effects and discuss approaches to the numerical challenge the modified wave functions present and a new semianalytic approximation. These tools are also applicable to calculating neutrino opacities in strong magnetic fields in the ejecta of binary neutron star mergers. We calculate the opacities for neutrinos capturing on free nucleons at subsaturation densities and temperatures exceeding an MeV. We find an enhancement to capture processes of the lowest energy neutrinos by an order of magnitude or more due to suppression of electron Pauli blocking in the case of capture on neutrons, and from the effect of the nucleon magnetic moments in the case of capture on protons.

astrophysical electromagnetic fields↗
DOE OSTI2021

Phase transition dynamics in one-dimensional halide perovskite crystals

Triiodide perovskites CsPbI 3 , CsSnI 3 , and FAPbI 3 (where FA is formamidinium) are highly promising materials for a range of optoelectronic applications in energy conversion. However, they are thermodynamically unstable at room temperature, preferring to form low-temperature (low-T) non-perovskite phases with one-dimensional anisotropic crystal structures. While such thermodynamic behavior represents a major obstacle toward realizing high-performance devices based on their high-temperature (high-T) perovskite phases, the underlying phase transition dynamics are still not well understood. In this work, we use in situ optical micro-spectroscopy to quantitatively study the transition from the low-T to high-T phases in individual CsSnI 3 and FAPbI 3 nanowires. We reveal a large blueshift in the photoluminescence (PL) peak (~38 meV) at the low-T/high-T two-phase interface of partially transitioned FAPbI 3 wire, which may result from the lattice distortion at the phase boundary. Compared to the experimentally derived activation energy of CsSnI 3 (~1.93 eV), the activation energy of FAPbI 3 is relatively small (~0.84 eV), indicating a lower kinetic energy barrier when transitioning from a face-sharing octahedral configuration to a corner-sharing one. Further, the phase propagation rate in CsSnI 3 is observed to be relatively high, which may be attributed to a high concentration of Sn vacancies. Our results could not only facilitate a deeper understanding of phase transition dynamics in halide perovskites with anisotropic crystal structures, but also enable controllable manipulation of optoelectronic properties via local phase engineering.

36 MATERIALS SCIENCE↗
NASA NTRS2020

An Improved Plastically Dilatant Unified Viscoplastic Constitutive Formulation for Multiscale Analysis of Polymer Matrix Composites Under High Strain Rate Loading

Polymer matrix composites are commonly used to fabricate energy-absorbing structures expected to experience impact loading. As such, a detailed understanding of the dynamic response of the constituent materials is necessary. Since the rate, temperature, and pressure dependence of carbon fiber reinforced polymer matrix composites are primarily manifestations of the rate, temperature, and pressure dependence of the polymer matrix, it is crucial that the constitutive behavior of the matrix be accurately characterized. In this work, an existing unified viscoplastic constitutive formulation is extended to ensure thermodynamic consistency and to more accurately account for the tension-compression asymmetry observed in the response of polymeric materials. A new plastic potential function is proposed, and elementary loading conditions are utilized to determine relations between model constants to ensure nonnegative plastic dissipation, a necessary thermodynamic requirement. Expressions for plastic Poisson’s ratios are derived and are bounded by enforcing nonnegative plastic dissipation. The model is calibrated against available experimental data from tests conducted over a range of strain rates, temperatures, and loading cases on a representative thermoset epoxy; good correlation between simulations and experimental data is obtained. Temperature rises due to the conversion of plastic work to heat are computed via the adiabatic heat energy equation. The viscoplastic polymer model is then used as a constitutive model in the generalized method of cells micromechanics theory to investigate the effects of matrix adiabatic heating on the high strain rate response of a unidirectional composite. The thermodynamic consistency of the model ensures plastic dissipation can only cause an increase in temperature. Simulation results indicate that significant thermal softening due to the conversion of plastic work to heat is observed in the composite for matrix dominated deformation modes.

Impact loading↗
DOE OSTI2021

Manipulating Pore Topology and Functionality to Promote Fluorocarbon-Based Adsorption Cooling

With the worldwide demand for refrigeration and cooling expected to triple, it is increasingly important to search for alternative energy resources to drive the refrigeration cycles with reduced electricity consumption. Recently, adsorption cooling has gained increased attention since energy reallocation in such systems is based on gas adsorption/desorption, which can be driven by waste/natural heat sources. Eco-friendly sorption-based cooling relies on the cyclic transfer of refrigerant gas from a high to low energy state by the pseudo-compression effect resulting from adsorption and desorption. The driving force for energy transfer relies on heat rather than electricity. The performance of a sorption chiller is primarily influenced by this cyclic sorption behavior, which is characterized as the working capacity of the porous sorbent. Thus increases in this working capacity directly translate to a more compact and efficient cooling system. How-ever, a lack of highly effective sorbent/refrigerant pairs lowers cooling performance, and therefore has limited applicability. To this end, synthetic metal-organic frameworks (MOFs) and covalent organic polymers (COPs) possess higher porosity and greater tunability leading to more substantial potential benefits for adsorption, compared to traditional sorbent materials. Similarly, hydrofluorocarbon refrigerants have more favorable applicability given ease of operation above atmospheric pressures due to suitable saturated vapor pressures and boiling points. For these reasons, our work focuses on an ongoing strategy to promote sorption cooling via improvements in the sorbent/refrigerant pair. Specifically, we target the interaction of hydrofluorocarbon refrigerants with MOF/COP materials at a molecular level by interpreting the host-guest chemistry and the role of framework pore topology. These molecular level differences translate to cooling performance, which is described herein. These strategies include engineering framework porosity (i.e., pore size, pore volume) by using elongated organic linkers and stereochemistry control during synthesis; manipulating the sorb-ate/sorbent interaction by introducing functional moieties or unsaturated metal centers to enhance working capacities in narrow pressure ranges; varying pore topology/morphology to impact adsorption isotherm behavior; and leveraging defective sites within the frameworks to further enhance adsorption capability. Here this atomic level understanding of sorb-ate-sorbent interactions is conducted using various in situ experimental techniques such as synchrotron-based X-ray diffraction, X-ray absorption spectroscopy, in situ Fourier transform infrared spectroscopy, and direct sorption energies determinization with calorimetry. Moreover, the experimentally studied interactions and the corresponding adsorption mechanism are corroborated by computational studies using density functional theory (DFT) and grand canonical Monte Carlo (GCMC) simulations. Using this approach, we envision the capability to engineer designed frameworks with precise molecular control to target refrigerant molecules and thereby enhance the performance of desired working pairs for sorption-based cooling.

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

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗
NASA NTRS1976

Work function determination of promising electrode materials for thermionic energy converters

The work function determinations of candidate materials for low temperature (1400 K) thermionics through vacuum emission tests are discussed. Two systems, a vacuum emission test vehicle and a thermionic emission microscope are used for emission measurements. Some nickel and cobalt based super alloys were preliminarily examined. High temperature physical properties and corrosion behavior of some super alloy candidates are presented. The corrosion behavior of sodium is of particular interest since topping cycles might use sodium heat transfer loops. A Marchuk tube was designed for plasma discharge studies with the carbide and possibly some super alloy samples. A series of metal carbides and other alloys were fabricated and tested in a special high temperature mass spectrometer. This information coupled with work function determinations was evaluated in an attempt to learn how electron bonding occurs in transition alloys.

Jacobson, D.↗
NASA NTRS2016

Li-Doped Ionic Liquid Electrolytes: From Bulk Phase to Interfacial Behavior

Ionic liquids have been proposed as candidate electrolytes for high-energy density, rechargeable batteries. We present an extensive computational analysis supported by experimental comparisons of the bulk and interfacial properties of a representative set of these electrolytes as a function of Li-salt doping. We begin by investigating the bulk electrolyte using quantum chemistry and ab initio molecular dynamics to elucidate the solvation structure of Li(+). MD simulations using the polarizable force field of Borodin and coworkers were then performed, from which we obtain an array of thermodynamic and transport properties. Excellent agreement is found with experiments for diffusion, ionic conductivity, and viscosity. Combining MD simulations with electronic structure computations, we computed the electrochemical window of the electrolytes across a range of Li(+)-doping levels and comment on the role of the liquid environment. Finally, we performed a suite of simulations of these Li-doped electrolytes at ideal electrified interfaces to evaluate the differential capacitance and the equilibrium Li(+) distribution in the double layer. The magnitude of differential capacitance is in good agreement with our experiments and exhibits the characteristic camel-shaped profile. In addition, the simulations reveal Li(+) to be highly localized to the second molecular layer of the double layer, which is supported by additional computations that find this layer to be a free energy minimum with respect to Li(+) translation.

molecular dynamics↗
DOE OSTI2024

ARJUNA: An Electrochemical Interface Mapping Probe for Solid-State Batteries

Solid-state batteries (SSBs) are promising candidates for next-generation energy storage, although their performance can be compromised by interfacial heterogeneity within the electrolyte. Furthermore, ensuring the quality of large form-factors electrolyte film is crucial for establishing a robust manufacturing platform for solid-state batteries. Herein, we report on the use of ARJUNA, an electrochemical interface mapping system, to characterize heterogeneities at solid electrolyte interfaces and to serve as a quality control system for SSB manufacturing. In addition to spatial mapping, the proposed system can also probe the interface behavior as a function of pressure and temperature. We present the operating principle, design, instrumentation, and evaluation of the system alongside a typical hybrid solid electrolyte produced using two common manufacturing processes. This report showcases the capability of ARJUNA to probe the heterogeneity and quality of processed solid electrolyte films.

25 ENERGY STORAGE↗
NASA NTRS1986

Measurements of electron attachment lineshapes and cross sections at ultralow electron energies for c-C6F10, c-C6F12, C8F16 and 1,1,2-C2Cl3F3

Electron attachment cross sections are reported in the electron energy range 0-160 meV and at energy resolutions of 4.5-7.5 meV (FWHM) for the molecules c-C6F10 (perfluorocyclohexene), c-C6F12 (perfluoro-1,2-dimethylcyclobutane), C8F16 (perfluoro-1,3-dimethylcyclohexane) and 1,1,2-C2Cl3F3 (1,1,2-trichlorotrifluoroethane). Use is made of the Kr photoionization technique, and measured attachment lineshapes are converted to cross sections by normalization through attachment rate constants. Comparisons are made with attachment cross sections derived from swarm-measured rate constants. Similar to previous results in eight other molecules, the present four molecules exhibit resolution-limited onsets at a threshold consistent with an s-wave attachment behavior and with a neutral-negative-ion curve crossing at zero energy.

Alajajian, S. H.↗
NASA NTRS2016

A Fragment-Cloud Model for Breakup of Asteroids with Varied Internal Structures

As an asteroid descends toward Earth, it deposits energy in the atmosphere through aerodynamic drag and ablation. Asteroid impact risk assessments rely on energy deposition estimates to predict blast overpressures and ground damage that may result from an airburst, such as the one that occurred over Chelyabinsk, Russia in 2013. The rates and altitudes at which energy is deposited along the entry trajectory depend upon how the bolide fragments, which in turn depends upon its internal structure and composition. In this work, we have developed an analytic asteroid fragmentation model to assess the atmospheric energy deposition of asteroids with a range of structures and compositions. The modeling approach combines successive fragmentation of larger independent pieces with aggregate debris clouds released with each fragmentation event. The model can vary the number and masses of fragments produced, the amount of mass released as debris clouds, the size-strength scaling used to increase the robustness of smaller fragments, and other parameters. The initial asteroid body can be seeded with a distribution of independent fragment sizes amid a remaining debris mass to represent loose rubble pile conglomerations, can be given an outer regolith later, or can be defined as a coherent or fractured monolith. This approach enables the model to represent a range of breakup behaviors and reproduce detailed energy deposition features such as multiple flares due to successive burst events, high-altitude regolith blow-off, or initial disruption of rubble piles followed by more energetic breakup of the constituent boulders. These capabilities provide a means to investigate sensitivities of ground damage to potential variations in asteroid structure.

Asteroid↗
DOE OSTI2023

Dynamic model-based feature extraction for fault detection and diagnosis of a supermarket refrigeration system

With the increasing concerns over climate change and carbon emissions, fault detection and diagnostics (FDD) of low–global warming potential (GWP) refrigerant supermarket refrigeration systems has gained great attention from academic and industrial sectors. Various FDD approaches have been developed to detect, identify, and diagnose faults to save energy, improve food quality, and protect the environment. Here, to mitigate the difficulty of collecting high-quality steady-state operational data in field operations faced by most model-based FDD methods, this study developed dynamic models of a low–GWP refrigerant (CO 2 ) supermarket refrigeration system. The model accuracy was validated using manufacturer data and experimental data. Simulations were conducted to predict the system dynamic response under two common operational faults—evaporator air path blockage fault and the display case door open fault—to identify fault patterns and define key dynamic behavior indexes for supporting FDD algorithm development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Part-scale evolution of fine-scale microstructural heterogeneity in solid-state additive manufacturing

Current solid-state additive manufacturing methods, refined through costly and time-consuming trial and error, have spurred interest in computational models that replicate material behavior under typical thermomechanical conditions (e.g., strain-rate ~ 102 s−1). These models, however, struggle to capture time-dependent microstructural evolution. In this work, Additive Friction-Stir Deposition (AFSD) is used as a representative case study for part-scale quantification of microstructure evolution at a fine spatial resolution (200 μm) by examining a liquid-nitrogen-cooled stop-action build via energy-dispersive X-ray diffraction coupled with a multi-channel detector. These results inform modeling efforts by linking process asymmetry to stored plastic strain, residual elastic strain, and texture development, and unlike current state-of-the-art characterization methods (e.g., EBSD or neutron diffraction), this approach provides both the spatial resolution and collection efficiency necessary to quantify fine-scale microstructural heterogeneity over large component volumes. As such, this technique provides essential validation data for computational models, e.g., crystal plasticity, enabling future prediction of heterogeneous behavior in AFSD and other additive manufacturing processes.

Franz, Cole [ORNL] (ORCID:0000000213465881)↗
DOE OSTI2020

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A miniaturized long-life acoustic transmitter for animal tracking in marine environments

Marine and hydrokinetic (MHK) energy devices are early-stage technologies and their impacts to ocean life and environment are largely unknown, which consequently limits the progress of deploying these devices at commercial scale. Optical and acoustic imaging have been the primary technologies for monitoring marine animals’ behavior at existing or potential MHK sites. However, these technologies lack the ability to identify individual animals, which makes them less useful in identifying protected species as well as studying the animals’ long-term behavior. Implantable acoustic transmitters with unique identification codes could be used as a complementary technology. To minimize the tag burden on the tagged animal (and hence the resulting bias) and to satisfy the long-term monitoring needs, it is desired that the implanted transmitter be as small and light as possible and can last for as long as possible. To date, the smallest acoustic transmitters suitable for marine animal tracking have relatively short service life and detection range (up to 200 m). In this report, we introduce a new miniaturized acoustic transmitter which operates at 208.3 kHz. It offers a smaller dry weight (0.45 g), a significantly longer service life (140 days at a 5-second ping rate) as well as a greater detection range than the existing commercial counterparts, which will help obtain more comprehensive movement and behavioral information of animals around MHK energy devices.

Li, Huidong↗
DOE OSTI2022

Complete characterization of a lithium battery electrolyte using a combination of electrophoretic NMR and electrochemical methods

Improving transport properties of electrolytes is important for developing lithium-ion batteries for future energy storage applications. In Newman's concentrated solution theory, electrolytes are characterized by three transport parameters, conductivity, diffusion coefficient, and transference number, in addition to the thermodynamic factor. In this work, these parameters are all determined for an exemplar liquid electrolyte, lithium bis(trifluoromethanesulfonyl)imide mixed in tetraethylene glycol dimethyl ether, using electrochemical methods. The intrinsic coupling between parameters obtained by electrochemical methods results in large error bars in the transference number that obscure the transport behavior of the electrolyte. Here, we use electrophoretic NMR (eNMR) to measure the electric-field-induced ion and solvent velocities to obtain the transference number directly, which enables determination of the thermodynamic factor with greater certainty. In conclusion, our work indicates that the combination of eNMR and electrochemical methods provides a robust approach for complete characterization of battery electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2023

The Role of SnO2 Processing on Ionic Distribution in Double-Cation-Double Halide Perovskites

Moving toward a future of efficient, accessible, and less carbon-reliant energy devices has been at the forefront of energy research innovations for the past 30 years. Metal-halide perovskite (MHP) thin films have gained significant attention due to their flexibility of device applications and tunable capabilities for improving power conversion efficiency. Serving as a gateway to optimize device performance, consideration must be given to chemical synthesis processing techniques. Therefore, how does common substrate processing techniques influence the behavior of MHP phenomena such as ion migration and strain? Here, we demonstrate how a hybrid approach of chemical bath deposition (CBD) and nanoparticle SnO2 substrate processing significantly improves the performance of (FAPbI3)0.97(MAPbBr3)0.03 by reducing micro-strain in the SnO2 lattice, allowing distribution of K+ from K-Cl treatment of substrates to passivate defects formed at the interface and produce higher current in light and dark environments. X-ray diffraction reveals differences in lattice strain behavior with respect to SnO2 substrate processing methods. Through use of conductive atomic force microscopy (c-AFM), conductivity is measured spatially with MHP morphology, showing higher generation of current in both light and dark conditions for films with hybrid processing. Additionally, time-of-flight secondary ionization mass spectrometry (ToF-SIMS) observed the distribution of K+ at the perovskite/SnO2 interface, indicating K+ passivation of defects to improve the power conversion efficiency (PCE) and device stability. We show how understanding the role of ion distribution at the SnO2 and perovskite interface can help reduce the creating of defects and promote a more efficient MHP device.

conductive atomic force microscopy↗
DOE OSTI2023

Effects of Aluminum and Molybdenum on the Phase Stability of Iron-Chromium Alloys: A First-Principles Study

The interaction between solute atoms is critical to the thermodynamic behavior of Fe-Cr alloys, but the effects of non-dilute Al and Mo on the Fe-Cr phase stability and vacancy formation energy are not clearly understood. In this study, density functional theory, cluster expansion, and Monte Carlo simulation are used to predict the effects of ternary solute elements on the thermodynamic properties in multicomponent Fe-Cr alloys. The machine learning regression approach is applied to train and construct energy models that accurately describe the complex chemical interactions. The computational outcomes include the prediction of the partial ternary phase diagram, mixing enthalpy, and vacancy formation energy for different compositions. The phase boundary calculation predicts a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from a’ precipitates. The mixing enthalpy calculation shows that Mo may also reduce the Cr solubility in bcc Fe. Additionally, the simulations show that non-dilute Cr decreases the vacancy formation energy in bcc Fe, while adding Al results in a less significant effect. Finally, the results demonstrate important applications of using machine learning energy models to study model or commercial alloys with multicomponent solute species and point defects.

36 MATERIALS SCIENCE↗