Search NASA⌕ Search

SEARCH · Search NASA

Results for “multi-component expansion”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

First Principles Prediction of Structure, Structure Selectivity, and Thermodynamic Stability under Realistic Conditions

Novel materials are often the enabler for new energy technologies. In ab-initio computational materials science, method are developed to predict the behavior of materials starting from the laws of physics, so that properties can be predicted before compounds have to be synthesized and tested. As such, a virtual materials laboratory can be constructed, saving time and money. The objectives of this program were to develop first-principles theory to predict the structure and thermodynamic stability of materials. Since its inception the program focused on the development of the cluster expansion to deal with the increased complexity of complex oxides. This research led to the incorporation of vibrational degrees of freedom in ab-initio thermodynamics, developed methods for multi-component cluster expansions, included the explicit configurational degrees of freedom of localized electrons, developed the formalism for stability in aqueous environments, and culminated in the first ever approach to produce exact ground state predictions of the cluster expansion. Many of these methods have been disseminated to the larger theory community through the Materials Project, pymatgen software, or individual codes. We summarize three of the main accomplishments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A General Framework for Progressive Data Compression and Retrieval

In scientific simulations, observations, and experiments, the transfer of data to and from disk and across networks has become a major bottleneck for data analysis and visualization. Compression techniques have been employed to tackle this challenge, but traditional lossy methods often demand conservative error tolerances to meet the numerical accuracy requirements of both anticipated and unknown data analysis tasks. Progressive data compression and retrieval has emerged as a promising solution, where each analysis task dictates its own accuracy needs. However, few analysis algorithms inherently support progressive data processing, and adapting compression techniques, file formats, client/server frameworks, and APIs to support progressivity can be challenging. Here, this paper presents a framework that enables progressive-precision data queries for any data compressor or numerical representation. Our strategy hinges on a multi-component representation that successively reduces the error between the original and compressed field, allowing each field in the progressive sequence to be expressed as a partial sum of components. We have implemented this approach with four established scientific data compressors and assessed its effectiveness using real-world data sets from the SDRBench collection. The results show that our framework competes in accuracy with the standalone compressors it is based upon. Additionally, (de)compression time is proportional to the number of components requested by the user. Finally, our framework allows for fully lossless compression using lossy compressors when a sufficient number of components are employed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

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

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Error controlling of the combined Cluster-Expansion and Wang–Landau Monte-Carlo method and its application to FeCo

Simulating order–disorder phase transitions in magnetic materials requires the accurate treatment of both the atomic and magnetic interactions, which span a vast configuration space. Taking FeCo as a prototype system, we demonstrate that this can be accurately and efficiently addressed by the combined Cluster-Expansion method and Wang–Landau algorithm. All involved key parameters in this method are carefully tested. Based on these tests, error controlling approaches are also developed. Furthermore our method is readily applicable to the study of similar second-order phase transitions in other binary and multi-component systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. II. UV, Optical, and Near-infrared Light Curves and Comparison to Kilonova Models

We present UV, optical, and near-infrared (NIR) photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at 0.47–18.5 days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the Hubble Space Telescope (HST). The spectral energy distribution (SED) inferred from this photometry at 0.6 days is well described by a blackbody model with T≈8300 K, a radius of R≈4.5×10{sup 14} cm (corresponding to an expansion velocity of v≈0.3c), and a bolometric luminosity of L{sub bol}≈5×10{sup 41} erg s{sup −1}. At 1.5 days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set, we find that models with heating from radioactive decay of {sup 56}Ni, or those with only a single component of opacity from r-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data; the resulting “blue” component has M{sub ej}{sup blue}≈0.01 M{sub ⊙} and v{sub ej}{sup blue}≈0.3 c, and the “red” component has M{sub ej}{sup red}≈0.04 M{sub ⊙} and v{sub ej}{sup red}≈0.1 c. These ejecta masses are broadly consistent with the estimated r-process production rate required to explain the Milky Way r-process abundances, providing the first evidence that binary neutron star (BNS) mergers can be a dominant site of r-process enrichment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗
NASA NTRS2006

Hubble Space Telescope Observations of Oxygen-rich Supernova Remnants in the Magellanic Clouds. III. WFPC2 Imaging of the Young, Crab-like Supernova Remnant SNRO540-69.3

Hubble Space Telescope images with the Wide Field Planetary Camera 2 of the young, oxygen-rich, Crab-like supernova remnant SNR0540-69.3 in the Large Magellanic Cloud (LMC) reveal details of the emission distribution and the relationship between the expanding ejecta and synchrotron nebula. The emission distributions appear very similar to those seen in the Crab nebula, with the ejecta located in a thin envelope surrounding the synchrotron nebula. The [O III] emission is more extended than other tracers, forming a faint "skin" around the denser filaments and synchrotron nebula, as also observed in the Crab. The [O III] exhibits somewhat different kinematic structure in long-slit spectra, including a more extended high-velocity emission halo not seen in images. Yet even the fastest expansion speeds in SNR 0540 s halo are slow when compared to most other young supernova remnants, though the Crab nebula has similar slow expansion speeds. We show a striking correspondence between the morphology of the synchrotron nebula observed in an optical continuum filter with that recently resolved in Chandra X-ray images. We argue that the multi-component kinematics and filamentary morphology of the optical emission-line features likely result from magnetic Rayleigh-Taylor instabilities that form as the synchrotron nebula expands and sweeps up ejecta, as seen in the Crab nebula. Our images and spectra help to refine our understanding of SNR 0540 in several more detailed respects: they confirm the identification of H(alpha)+[N II] in the red spectrum, they show that the systemic velocity of SNR 0540 is not significantly different from that of the LMC, and they hint at a lower Ne abundance than the Crab (potentially indicating a more massive progenitor star).

Morse, Jon A.↗
DOE OSTI2014

Intrinsic fluctuations of dust grain charge in multi-component plasmas

A master equation is formulated to model the states of the grain charge in a general multi-component plasma, where there are electrons and various kinds of positive or negative ions that are singly or multiply charged. A Fokker-Planck equation is developed from the master equation through the system-size expansion method. The Fokker-Planck equation has a Gaussian solution with a mean and variance governed by two initial-value differential equations involving the rates of the attachment of ions and electrons to the dust grain. Also, a Langevin equation and a discrete stochastic method are developed to model the time variation of the grain charge. Grain charging in a plasma containing electrons, protons, and alpha particles with Maxwellian distributions is considered as an example problem. The Gaussian solution is in very good agreement with the master equation solution numerically obtained for this problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Containerless Measurement of Thermophysical Properties of Ti-Zr-Ni Alloys

The surface tension, viscosity, density, and thermal expansion of Ti-Zr-Ni alloys were measured for a number of compositions by electrostatic levitation methods. Containerless methods greatly reduce heterogeneous nucleation, increasing access to the undercooled liquid regime at finite cooling rates. The density and thermal expansion are measured optically, while the surface tension and viscosity are measured by the oscillating drop method. The measured alloys include compositions which form a metastable quasicrystal phase from the undercooled liquid, and alloys close to the composition of several multi-component bulk metallic glass-forming alloys. Measurements of surface tension show behavior typical of transition metals at high temperature, but a sudden decrease in the deeply undercooled liquid for alloys near the quasicrystal-forming composition range, but not for compositions which form the solid-solution phase first.

Hyers, Robert↗
DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Thin film combinatorial sputtering of Al-Ce alloys: Investigating the phase separation of as-deposited solid solutions and determining the coefficient of thermal expansion

Al x Ce 100–x thin films with a composition range of ~75.0 < x < 99.5 at% (36.5 < x < 97.5 wt%) were synthesized via combinatorial co-sputtering from an Al and an Al 50 Ce 50 target. The crystal structure, phase fraction, film morphology, electrical resistivity, and temperature-dependent coefficients of thermal expansion (CTE) are all correlated to the Al x Ce 100–x composition. The as-deposited films form a metastable solid-solution, and annealing leads to the formation of the thermodynamically stable two-phase system of Al and the α-Al 11 Ce 3 intermetallic. Temperature dependent x-ray diffraction (XRD) reveals that the two phases expand independently of one another, and the thin film Al temperature-dependent CTE is similar to bulk Al. The thin film Al 11 Ce 3 intermetallic phase has a nearly constant CTE of ~1.5 × 10 –5 /°C within the temperature range studied (25–550 °C). To confirm the thin film Al 11 Ce 3 results, bulk stoichiometric Al 11 Ce 3 and +/- 1 wt% Ce samples were prepared and the CTE of each was measured with the same conditions. A Rietveld analysis of the bulk data enabled an estimation of the CTE in each of the 3 orthorhombic lattice parameters, which displayed anisotropic behavior. The thin film and bulk CTE measurements were in very good agreement. Estimations of the temperature dependent CTE of the two-phase alloys are made via the Reuss and Voigt models. Finally, by demonstrating the efficacy of the approach, more complex multi-component rapid materials discovery of low CTE Al-alloys can be pursued via the combinatorial thin film synthesis and XRD measurement.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study

The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.

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

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2019

Pronghorn

Pronghorn is designed as a multi-dimensional coarse-mesh conjugate heat transfer simulator (tightly coupled fluid and solid phase heat transfer) for nuclear reactor core applications. Pronghorn’s design philosophy of is that: 1. The tool will be agnostic of reactor concept, including VHTRs, SFRs, FHRs, and LWRs. 2. Must be higher fidelity than traditional nuclear systems/safety analysis (0D-1D). 3. Avoid lower-length scale phenomena, such as boundary layer theory, i.e., rely on closure relations. 4. Apply mixture theory (homogenization) wherever possible. 5. Conserve mass and energy (solve all physics on single mesh). 6. Must be easily coupled to both MOOSE-based and external physics applications. 7. Does not require a super computer. Thus, Pronghorn’s algorithmic philosophy is to resolve large scale flow features and model lower-length scale phenomena, such as boundary layers. The coarse mesh approach for lower-length scale modeling is to develop a multi-component multi-phase homogenization of solid and liquid phases into a mixture formulation inside a finite element. For example, a pebble bed core is represented as a two-phase mixture of coolant and fuel pebbles, where the fuel phase is stationary. Thus, the solid phase exerts a drag force on the fluid phase, composed of form and viscous drag, and both phases exchange internal energy through conduction and convective heat transfer. For Pronghorn, the volume fraction of fluid can vary in space, but not in time. Pronghorn has been validated against the VHTR PBMR400 benchmark problems, used for rod ejection simulations have been conducted for thermal-fluids/neutronics for both pebble-bed and prismatic gas-cooled reactors and simple thermal-neutronics coupled LWR benchmark problems. Pronghorn flow models cover a wide range of transient and steady flow conditions, including, porous fluid flow, compressible flow, and nearly incompressible, thermally expansive flow. The physics can be solved in three-dimensional Cartesian coordinates and two-dimensional cylindrical (r-z) coordinates.

Martineau, RichardC.↗