SEARCH · Engineering Papers

Results for “elastic crystal”

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 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Composition dependence of atomic order in strain-relaxed, metastable GeSn alloys

Extended x-ray absorption fine structure (EXAFS) measurements of single-crystal Ge/GeSn radial heterostructure nanowires are used to examine the effects of composition on both short-range order (SRO) and longer-range disorder in GeSn alloys. GeSn has prompted significant interest because it can achieve a direct band gap for sufficient Sn concentrations beyond the equilibrium solid solubility limit in an all-group IV system. Short-range order in this material is particularly interesting as it has been predicted to affect the band gap independent of average composition or strain effects. By independently controlling the Sn composition and GeSn thickness during chemical vapor deposition of misfitting GeSn shells around ultrathin, elastically compliant, Ge core nanowires, the elastic misfit strain in the GeSn is minimized for Sn compositions over the studied range ≈Ge 0.96 Sn 0.04 to Ge 0.88 Sn 0.12 . The degree of SRO was found to decrease with increasing Sn composition. Additionally, damping of the EXAFS signal was observed as the Sn content increased, particularly for increasingly distant atomic shells about the absorbing atom, even for scattering paths not involving Sn atoms. This result is quantified as an increase in the mean-squared relative displacement parameters of the shells. These measurements reveal the accommodation of local strain due to the presence of the highly size-mismatched Sn atoms in the Ge diamond cubic lattice (≈14%), which may have effects on the band structure of the material in addition to the influence of short-range atomic order. Comparison among the nanowire samples allows for calculation of the topological rigidity parameter, a ∗∗ , for the first-neighbor bond lengths. Furthermore, these exhibit chemically distinct values for Ge-Ge, Ge-Sn, and Sn-Sn, and they are consistent with the value a ∗∗ = 0.75 ± 0.07 confirming the general applicability of the model to alloys with both large amounts of natural misfit strain and the potential for short-range order.

Crystal structure Semiconductors Transmission elec↗

DOE OSTI2021

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1989

Notched fatigue of single crystal PWA 1480 at turbine attachment temperatures

The focus is on the lower temperature, uncoated and notched features of gas turbine blades. Constitutive and fatigue life prediction models applicable to these regions are being developed. Fatigue results are presented which were obtained thus far. Fatigue tests are being conducted on PWA 1480 single crystal material using smooth strain controlled specimens and three different notched specimens. Isothermal fatigue tests were conducted at 1200, 1400, and 1600 F. The bulk of the tests were conducted at 1200 F. The strain controlled tests were conducted at 0.4 percent per second strain rate and the notched tests were cycled at 1.0 cycle per second. A clear orientation dependence is observed in the smooth strain controlled fatigue results. The fatigue lifes of the thin, mild notched specimens agree fairly well with this smooth data when elastic stress range is used as a correlating parameter. Finite element analyses were used to calculate notch stresses. Fatigue testing will continue to further explore the trends observed thus far. Constitutive and life prediction models are being developed.

Meyer, T. G.↗

DOE OSTI2020

Pressure-induced phase transition in barium hydride studied with neutron scattering

Barium hydride can undergo a structural phase transition from an orthorhombic phase to a hexagonal phase induced by high temperature or high pressure. This transition causes an immediate increase in the hydrogen diffusion rates by over an order of magnitude, and therefore, understanding the origin and details of such transition is of great interest not only for fundamental reasons but also for improving materials for future applications. In this work, the pressure evolution of the crystal structure was characterized using neutron powder diffraction up to a maximum pressure of 11.3 GPa. The pressure dependence of the unit cell volumes, lattice parameters, atomic sites, and compressibilities were determined for both phases. A structural phase transition occurred over a wide pressure range of P = 1.3 GPa–4.9 GPa. The transition to the higher density hexagonal phase reduced the volume per formula unit of BaD 2 by 13.6%, hence increasing the volumetric storage density. In addition, we investigated the hydrogen diffusion process using high pressure quasi-elastic neutron scattering up to 7.1 GPa. Our results show that the hydrogen mobility increases with pressure in the hexagonal phase. This work sheds light on the structural and dynamical aspects of barium hydride caused by the application of high pressure. Finally, the results may aid in the development of advanced metal hydride systems with increased hydrogen dynamics.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Phonon-mediated quantum efficiency measurement in semiconductors

Accurate quantum efficiency measurement not only provides crucial information for the photovoltaic cell industry but also supports experiments aimed at directly detecting dark matter and elastic neutrino interactions. The dark matter direct searches paradigm has recently expanded to include particles with masses below 1,MeV/c2, where the expected signal in an electron–recoil interaction is approximately in the eV range, just above the energy gap for silicon and germanium. A robust calibration method for ionization signals in this lower energy region is essential. This paper presents a method for measuring quantum efficiency and yield (q/E) in semiconductors using phonon-mediated calorimetry. The Neganov–Trofimov–Luke phonon amplification method in low-temperature semiconductor crystals has been employed to indirectly measure ionization down to single-electron accuracy. Specifically, at zero bias, the phonon readout directly quantifies the total energy deposited within the detector, independent of the ionization yield. This eliminates a significant source of systematic uncertainty in quantum efficiency estimates associated with total energy uncertainty. The paper includes results from an updated ionization efficiency measurement in a germanium detector.

Physics↗

NASA NTRS1985

Investigation of superlattice device structures

This report describes the investigation of growth properties, and the structure of epitaxial multilayer Si(Si(1x)Ge(x)) films grown on bulk Silicon Substrates. It also describes the fabrication and characterization of MOSFET and MESFET devices made on these epitaxial films. Films were grown in a CVD reactor using hydrides of Si and Ge with H2 and He as carrier gases. Growth temperatures were between 900 C and 1050 C with most films grown at 1000 C. Layer thickness was between 300A and 2000A and total film thickness was between 0.25 micro m and 7 micro m. The Ge content (X) in the alloy layers was between .05 and 0.2. N-type multilayer films grown on (100) p-type Si showed Hall mobility in the range 1000 to 1500 sq cm/v for an average carrier concentration of approx. 10 to the 16th power/cu cm. This is up to 50% higher than the Hall mobility observed in epitaxial Si films grown under the same conditions and with the same average carrier concentration. The mobility enhancement occurred in films with average carrier concentration (n) from 0.7 x 10 to the 16th power to 2 x 10 to the 17th power/cu cm, and total film thickness greater than 1.0 micro m. No mobility enhancement was seen in n-type multilayer films grown on (111) Si or in p-type multilayer films. The structure of the films was investigated was using SEM, TEM, AES, SIMS, and X-ray double crystal diffraction techniques. The film composition profile (AES, SIMS) showed that the transition region between layers is of the order of about 100A. The TEM examination revealed a well defined layered structure with fairly sharp interfaces and good crystalline quality. It also showed that the first few layers of the film (closest to the substrate) are uneven, most probably due to the initial growth pattern of the epitaxial film where growth occurs first in isolated islands that eventually growth and coalesce. The X-ray diffraction measurement determined the elastic strain and strain relief in the alloy layers of the film and the elastic strain in the intervening Si layers.

Gergis, I. S.↗

DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗

DOE OSTI2023

Crystal structure of calcium-ferrite type NaAlSiO4 up to 45 GPa

Abstract Alkali-rich aluminous high-pressure phases including calcium-ferrite (CF) type NaAlSiO4 are thought to constitute ~20% by volume of subducted mid-ocean ridge basalt (MORB) under lower mantle conditions. As a potentially significant host for incompatible elements in the deep mantle, knowledge of the crystal structure and physical properties of CF-type phases is therefore important to understanding the crystal chemistry of alkali storage and recycling in the Earth’s mantle. We determined the evolution of the crystal structure of pure CF-NaAlSiO4 and Fe-bearing CF-NaAlSiO4 at pressures up to ~45 GPa using synchrotron-based, single-crystal X-ray diffraction. Using the high-pressure lattice parameters, we also determined a third-order Birch-Murnaghan equation of state, with V0 = 241.6(1) Å3, KT0 = 220(4) GPa, and KT0′ = 2.6(3) for Fe-free CF, and V0 = 244.2(2) Å3, KT0 = 211(6) GPa, and KT0′ = 2.6(3) for Fe-bearing CF. The addition of Fe into CF-NaAlSiO4 resulted in a 10 ± 5% decrease in the stiffest direction of linear compressibility along the c-axis, leading to stronger elastic anisotropy compared with the Fe-free CF phase. The NaO8 polyhedra volume is 2.6 times larger and about 60% more compressible than the octahedral (Al,Si)O6 sites, with K0NaO8 = 127 GPa and K0(Al,Si)O6 ~304 GPa. Raman spectra of the pure CF-type NaAlSiO4 sample shows that the pressure coefficient of the mean vibrational mode, 1.60(7) cm–1/GPa, is slightly higher than 1.36(6) cm−1/GPa obtained for the Fe-bearing CF-NaAlSiO4 sample. The ability of CF-type phases to contain incompatible elements such as Na beyond the stability field of jadeite requires larger and less-compressible NaO8 polyhedra. Detailed high-pressure crystallographic information for the CF phases provides knowledge on how large alkali metals are hosted in alumina framework structures with stability well into the lowermost mantle.

Geochemistry & Geophysics↗

DOE OSTI2020

Beyond Graphene: Low-Symmetry and Anisotropic 2D Materials

Low-symmetry 2D materials—such as ReS 2 and ReSe 2 monolayers, black phosphorus monolayers, group-IV monochalcogenide monolayers, borophene, among others—have more complex atomistic structures than the honeycomb lattices of graphene, hexagonal boron nitride, and transition metal dichalcogenides. The reduced symmetries of these emerging materials give rise to inhomogeneous electron, optical, valley, and spin responses, as well as entirely new properties such as ferroelasticity, ferroelectricity, magnetism, spin-wave phenomena, large nonlinear optical properties, photogalvanic effects, and superconductivity. Novel electronic topological properties, nonlinear elastic properties, and structural phase transformations can also take place due to low symmetry. Lastly, the “Beyond Graphene: Low-Symmetry and Anisotropic 2D Materials” Special Topic was assembled to highlight recent experimental and theoretical research on these emerging materials.

2D materials↗

DOE OSTI2022

Lamellae of α phase in a metastable β -Ti alloy studied by small-angle x-ray scattering

The shapes and orientations of lamellae of the hexagonal close-packed α phase in single crystals of the body-centred cubic β phase in a metastable β -Ti alloy (Timetal LCB) were investigated by synchrotron small-angle x-ray scattering. Several steps of thermal treatment with various annealing times and temperatures were employed. The measured scattering data were fitted to a model of ellipsoidal lamellae of random sizes, and the mean lamella shapes and orientations of their axes with respect to the β -Ti matrix were determined from the fit. The number and relative volume of the lamellae were estimated from the integrated scattered intensity. The elastic energy of the lamellae as well as the energy of defects at the α / β interface (misfit dislocations) were calculated from the parameters of the lamellae. In conclusion, the results provide information on the kinetics of the β $ \rightarrow$ α transformation in metastable β -Ti alloys and allow assessing the degree of thermodynamic equilibrium of the final α / β composite.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Material Design Strategy for Halide Solid Electrolytes Li 3 MX 6 (X = Cl, Br, and I) for All-Solid-State High-Voltage Li-Ion Batteries

Although several solid electrolyte (SE) candidates have been explored, achieving the necessary combination of performance, stability, and processability has been challenging. Recently, several lithium ternary halides have attracted increasing attention for SEs because of their favorable combination of high ionic conductivity and wide electrochemical window. This study aims to provide a material design strategy for lithium halides Li 3 MX 6 (X = Cl, Br, and I) for high-voltage all-solid-state Li-ion batteries, achieved by the systematic investigation of crystal structures, phase and electrochemical stabilities, electronic and mechanical properties, and ionic conductivities. Calculation results reveal that the electronegativity difference between M and X affects structural properties and stabilities. Weak Coulomb interactions in Li 3 MX 6 result in the preference of the monoclinic phase, and the oxidation potential and chemical stability against the cathode materials of Li 3 MX 6 increase for relatively small X. Chlorides exhibit the highest oxidation potential (~4.3 V) among Li 3 MX 6 , suggesting that chlorides are appropriate SEs for high-voltage cathodes. The band gap and elastic moduli increase for relatively small X, suggesting the relatively low electronic conductivity and elastic deformability of chlorides. Chlorides with transition metals typically exhibit trigonal phases, a wider electrochemical stability window, a larger band gap, and higher elastic moduli compared to other types of halides. Additionally, chloride Li 3 MCl 6 is expected to have relatively high ionic conductivities with the aliovalent substitution of M 3+ to Zr 4+ and the anion mixing of Cl with Br. The findings of this study will provide fundamental guidelines for the development of lithium halide SEs for high-voltage all-solid-state Li-ion batteries.

25 ENERGY STORAGE↗

DOE OSTI2025

Elasticity of davemaoite as a primary contributor to lower-mantle heterogeneities

Geophysical detection of subducted mid–ocean ridge basalt (MORB) in the lower mantle is hindered by uncertainties in the elasticity of Fe,Al,Mg,Ti–bearing davemaoite, a key MORB component. Using Brillouin spectroscopy and x-ray diffraction, we determined the elasticity of a Ca 0.906(1) Fe 2+ 0.027(1) Fe 3+ 0.042(1) Mg 0.033(1) Al 0.072(1) Ti 0.020(1) Si 0.912(1) O 3 davemaoite up to 113 gigapascals and 2294 K. Here, we found that it exhibited a shear wave velocity 10 to 20% slower than end-member davemaoite, making it the slowest phase among major lower-mantle minerals. Our models show that MORB, containing 20 to 25 volume percent davemaoite, potentially contributes to large low-shear-velocity provinces (LLSVPs), whereas a cumulate layer enriched in davemaoite crystallized from basal magma ocean may comprise ultralow-velocity zones (ULVZs). Davemaoite’s ability to host incompatible and heat-producing elements possibly links LLSVPs and ULVZs to mantle plume initiation and geochemical signatures of ocean island basalts.

Zhou, Wen-Yi [Texas A & M Univ., College Station, ↗

NASA NTRS2004

NASA Tech Briefs, February 2004

Topics include: Simulation Testing of Embedded Flight Software; Improved Indentation Test for Measuring Nonlinear Elasticity; Ultraviolet-Absorption Spectroscopic Biofilm Monitor; Electronic Tongue for Quantitation of Contaminants in Water; Radar for Measuring Soil Moisture Under Vegetation; Modular Wireless Data-Acquisition and Control System; Microwave System for Detecting Ice on Aircraft; Routing Algorithm Exploits Spatial Relations; Two-Finger EKG Method of Detecting Evasive Responses; Updated System-Availability and Resource-Allocation Program; Routines for Computing Pressure Drops in Venturis; Software for Fault-Tolerant Matrix Multiplication; Reproducible Growth of High-Quality Cubic-SiC Layers; Nonlinear Thermoelastic Model for SMAs and SMA Hybrid Composites; Liquid-Crystal Thermosets, a New Generation of High-Performance Liquid-Crystal Polymers; Formulations for Stronger Solid Oxide Fuel-Cell Electrolytes; Simulation of Hazards and Poses for a Rocker-Bogie Rover; Autonomous Formation Flight; Expandable Purge Chambers Would Protect Cryogenic Fittings; Wavy-Planform Helicopter Blades Make Less Noise; Miniature Robotic Spacecraft for Inspecting Other Spacecraft; Miniature Ring-Shaped Peristaltic Pump; Compact Plasma Accelerator; Improved Electrohydraulic Linear Actuators; A Software Architecture for Semiautonomous Robot Control; Fabrication of Channels for Nanobiotechnological Devices; Improved Thin, Flexible Heat Pipes; Miniature Radioisotope Thermoelectric Power Cubes; Permanent Sequestration of Emitted Gases in the Form of Clathrate Hydrates; Electrochemical, H2O2-Boosted Catalytic Oxidation System; Electrokinetic In Situ Treatment of Metal-Contaminated Soil; Pumping Liquid Oxygen by Use of Pulsed Magnetic Fields; Magnetocaloric Pumping of Liquid Oxygen; Tailoring Ion-Thruster Grid Apertures for Greater Efficiency; and Lidar for Guidance of a Spacecraft or Exploratory Robot.

Source record↗

DOE OSTI2020

Imaging phonon dynamics with ultrafast electron microscopy: Kinematical and dynamical simulations

Ultrafast x-ray and electron scattering techniques have proven to be useful for probing the transient elastic lattice deformations associated with photoexcited coherent acoustic phonons. Indeed, femtosecond electron imaging using an ultrafast electron microscope (UEM) has been used to directly image the influence of nanoscale structural and morphological discontinuities on the emergence, propagation, dispersion, and decay behaviors in a variety of materials. Here, we describe our progress toward the development of methods ultimately aimed at quantifying acoustic-phonon properties from real-space UEM images via conventional image simulation methods extended to the associated strain-wave lattice deformation symmetries and extents. Using a model system consisting of pristine single-crystal Ge and a single, symmetric Lamb-type guided-wave mode, we calculate the transient strain profiles excited in a wedge specimen and then apply both kinematical- and dynamical-scattering methods to simulate the resulting UEM bright-field images. While measurable contrast strengths arising from the phonon wavetrains are found for optimally oriented specimens using both approaches, incorporation of dynamical scattering effects via a multi-slice method returns better qualitative agreement with experimental observations. Contrast strengths arising solely from phonon-induced local lattice deformations are increased by nearly an order of magnitude when incorporating multiple electron scattering effects. We also explicitly demonstrate the effects of changes in global specimen orientation on the observed contrast strength, and we discuss the implications for increasing the sophistication of the model with respect to quantification of phonon properties from UEM images.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Single-crystal X-ray diffraction of fluorapatite to 61 GPa

Apatite is a mineral of widespread importance in Earth and planetary science. Here we examine the behavior of a natural fluorapatite (FAp) crystal from Durango (Mexico) under compression to 61 GPa. Single-crystal X-ray diffraction experiments were carried out in a diamond-anvil cell using a synchrotron source. The apatite structure persists up to 32.4 GPa. Birch-Murnaghan equation of state parameters were fit to the pressure-volume data for fluorapatite for two cases: fixing V 0 at its measured ambient value resulted in a bulk modulus, K 0T , of 97.0(8) GPa and a pressure derivative of the bulk modulus, K' 0T , of 3.3(1), while fixing V 0 and K 0T at its ambient value 90.5 GPa (derived from ultrasonically measured elastic constants) resulted in a K' 0T value of 4.1(1). At 35.6 GPa, fluorapatite transforms to a triclinic phase (P$\bar{1}$, Z = 4), designated here as fluorapatite II (FAp-II). This phase persists up to at least 61 GPa. The major structural differences between FAp and FAp-II involve the buckling of the Ca polyhedra along the c-axis and changes in the number and coordination of the Ca sites. Furthermore, our study extends the pressure range over which fluorapatite has been examined by more than a factor of three, providing new insights into its structural response to high-pressure conditions.

58 GEOSCIENCES↗

DOE OSTI2022

Medial packing and elastic asymmetry stabilize the double-gyroid in block copolymers

Triply-periodic networks are among the most complex and functionally valuable self-assembled morphologies, yet they form in nearly every class of biological and synthetic soft matter building blocks. In contrast to simpler assembly motifs – spheres, cylinders, layers – networks require molecules to occupy variable local environments, confounding attempts to understand their formation. Here, we examine the double-gyroid network phase by using a geometric formulation of the strong stretching theory of block copolymer melts, a prototypical soft self-assembly system. The theory establishes the direct link between molecular packing, assembly thermodynamics and the medial map, a generic measure of the geometric center of complex shapes. We show that “medial packing” is essential for stability of double-gyroid in strongly-segregated melts, reconciling a long-standing contradiction between infinite- and finite-segregation theories. Additionally, we find a previously unrecognized non-monotonic dependence of network stability on the relative entropic elastic stiffness of matrix-forming to tubular-network forming blocks. The composition window of stable double-gyroid widens for both large and small elastic asymmetry, contradicting intuitive notions that packing frustration is localized to the tubular domains. This study demonstrates the utility of optimized medial tessellations for understanding soft-molecular assembly and packing frustration via an approach that is readily generalizable far beyond gyroids in neat block copolymers.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

DOE OSTI2022

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗