Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power coefficient”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2026

Applying Particle Swarm Optimization and Extended Kalman Filtering to Model Kaplan Generation Dynamics for Hydropower Systems

Variable renewable generation is increasing the need for hydropower plants to provide fast and flexible grid support, which places new demands on plant-level dynamic models used for monitoring, control, and operational decision-making. This need is especially important for hydroelectric systems, where turbine and generator dynamics are strongly coupled, nonlinear, and time-varying, making accurate real-time representation difficult. To address this problem, this paper develops a digital twin (DT) framework for a synchronous generator–Kaplan turbine system using an explicit separation of slow turbine dynamics and fast generator dynamics. The turbine subsystem is represented by a six-coefficient model, whose parameters are identified offline using particle swarm optimization, while the generator subsystem is updated online through an extended Kalman filter for real-time state and parameter estimation. These models are integrated within a closed-loop simulation that includes a proportional–integral–derivative–double-derivative governor and excitation system, allowing the DT to track plant behavior under realistic operating conditions. Unlike prior studies that treat turbine and generator modeling separately or rely mainly on simulated inputs, the proposed framework is validated using real operational data from a hydropower plant. Results show that the DT reproduces terminal voltage, active power, and reactive power with a normalized root mean square error of approximately 5%. This hybrid offline–online formulation constitutes the main contribution of the work, providing an adaptive and practically deployable DT for hydropower systems with direct relevance to control improvement, performance monitoring, and grid-support applications under high renewable penetration.

13 HYDRO ENERGY↗
DOE OSTI2026

Fast Growth of Ga2O3 on Highly Offcut (100)-Oriented Substrates

Beta-Ga2O3 has emerged as a leading candidate for next-generation power electronics, radio frequency (RF) switches, and extreme environment electronics due to a wide band gap (4.6 - 4.9 eV), high dopability (approximately 40 meV activation energy for an isolated silicon donor), and melt growth characteristics resulting in commercially available 4-inch substrates and commercial demonstrations of 6-inch substrates by multiple techniques. The (100) surface of Ga2O3 is highly desirable from a device and epitaxy standpoint - bulk growth of (100) material is more scalable than (010), the surface is nearly lattice-matched to p-type partner NiO, and Al2O3 incorporates at higher concentrations without phase separation. More importantly, the impact ionization coefficients along the [100] direction are low, leading to the highest possible critical fields. This is advantageous compared to the current state of the art, (001), due to reduced surface defects and increased possible breakdown voltage. However, the epitaxial growth rate on (100) surfaces is less than 10% of other faces due to weak bonding and favorable desorption, and on-axis (100) growth easily forms twin domains. Recent demonstrations have shown growth rate improvements from 0.4 nm/min to 1.5 nm/min by growing on (100) wafers that are offcut 6 degrees in the -c direction. These films show step-flow growth from (20-1) step-edges and high electron mobility due to suppressed twins. Despite these exciting results, offcuts greater than 6 degrees have not been explored due to the waste associated with grinding and polishing large offcuts. In this talk we will discuss the molecular beam epitaxy (MBE) growth and properties of Beta-Ga2O3 grown on (100) substrates offcut in the -c direction up to 13.4 degrees. These large offcuts are enabled by edge-fed film-defined growth (EFG) where the offcut is grown into the surface by pulling the crystal through the EFG die with the seed crystal rotated by the desired offcut angle. We will demonstrate that 13.4 degrees offcut substrates still exhibit a terraced (100) surface, and that a >10x increase (>5 nm/min) in growth rate is achieved. As previously reported on lower offcuts, we observe reversal of substrate twin domains around the (001) direction at the substrate-epilayer interface. We will discuss electrical properties including record-low (by MBE) unintentional doping densities of < 5E15 cm-3 and critical breakdown field in Schottky barrier diodes comparable with state-of-the-art (001) Ga2O3 without edge termination.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Performance Evaluation of Low-Cost Condensers Immersed in the Water Tank of a Heat Pump Water Heater

This paper reports the effect of low-cost immersed condensers on the performance of a heat pump water heater (HPWH) that does not require a water-circulating pump. The proposed design is based on the immersion of an L-style condenser coil through an opening in the top of the water tank. The novel design aims to eliminate the need for the water-circulating pumps, thereby substantially improving efficiency and reducing costs and maintenance of the HPWH system. Comprehensive computational fluid dynamics (CFD) simulations and experimental tests were performed, and the results showed that the CFD simulations were close to the experimental data. Further results confirm that the L-style condenser can introduce buoyancy-driven flow and eliminate temperature stratification in the studied HPWH water tank, thereby improving the heat transfer coefficient and coefficient of performance of the HPWH. The testing results further showed that the L-style condenser enabled more than 53% higher coefficient of performance and 400% higher heat transfer coefficient compared with the straight vertical condenser, and achieved much faster water heating. Furthermore, the effects of immersed condensers with different circuits on the HPWH performance were evaluated and compared. All results indicate that the current energy-saving and inexpensive HPWH technology is technically valid, and the benefits of performance improvement and low cost are attractive in future HPWH design.

Gao, Zhiming [ORNL] (ORCID:0000000271397995)↗
DOE OSTI2022

Wire-coil insert optimization for high-heat-load/flux synchrotron components

Many synchrotron components require high levels of internal-flow, forced-convection heat transfer to minimize surface temperatures, thermal gradients, and thermally induced stress on high-power beam-interacting surfaces. Wire-coil inserts, physically similar to a common spring, are mechanically fitted inside of component cooling passages to optimize heat transfer performance. They are routinely used in Advanced Photon Source (APS) front end and beamline high-heat-load/flux components to significantly enhance convection heat transfer-up to 400% compared to plain open passages. This has the additional benefit of greatly reducing coolant flow requirements for these components. Using several cooling passage sizes, five different wire sizes, and a range of pitch values, an experimental investigation conducted at the APS has determined the average heat transfer coefficient and resulting pressure loss as a function of water flow rate for 65 different wire-coil inserts. Data from this study have been non-dimensionalized and generalized to yield relationships that can be used to determine the heat transfer performance and resulting pressure loss for any given wire-coil insert that may be used at the APS. Through data reduction, the wire-coil insert characteristic dimensions have also been optimized to yield the highest heat transfer enhancement while minimizing the coolant flow requirements. These generalized expressions for wirecoil inserts will be presented, and they can be used by scientists and engineers during the component design process to evaluate achievable heat transfer performance and associated pressure loss, aiding in the establishment of optimized operating parameters and cooling passage flow distribution schemes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Sequential Doping of Ladder-Type Conjugated Polymers for Thermally Stable n-Type Organic Conductors

Doping of organic semiconductors is a powerful tool to optimize the performance of various organic (opto)electronic and bioelectronic devices. Despite recent advances, the low thermal stability of the electronic properties of doped polymers still represents a significant obstacle to implementing these materials into practical applications. Hence, the development of conducting doped polymers with excellent long-term stability at elevated temperatures is highly desirable. Here, we report on the sequential doping of the ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) with a benzimidazole-based dopant (i.e., N-DMBI). By combining electrical, UV–vis/infrared, X-ray diffraction, and electron paramagnetic resonance measurements, we quantitatively characterized the conductivity, Seebeck coefficient, spin density, and microstructure of the sequentially doped polymer films as a function of the thermal annealing temperature. Importantly, we observed that the electrical conductivity of N-DMBI-doped BBL remains unchanged even after 20 h of heating at 190 °C. This finding is remarkable and of particular interest for organic thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First lower hybrid current drive experiments on the WEST tokamak

The first lower hybrid current drive experiments in the full tungsten WEST tokamak are reported. Good wave coupling is found at rather low plasma current (q 95 ~4.3) and medium density (n e ~3×10 19 m -3 ). Reflection coefficient is in agreement with the expectation from the linear theory of coupling. With low reflection coefficients, 5MW was coupled for 2 seconds. High central electron temperature, up to 5keV, is achieved at ne = 3-4×10 19 m -3 . Flat and even hollow profiles of tungsten density are derived from the bolometry diagnostic. The stored thermal energy follows the H96-P scaling law with very slight degradation with density. The current drive efficiency has been assessed in low loop voltage (V L ~0.15V) discharges. Low plasma current operation and rather high effective charge (Z eff ~3) lead to modest current drive efficiency (η = 0.5-0.65 × 10 19 A.W -1 m -2 ). Long pulse operation (~30s) with high LHCD power (P LH =2.7MW) is achieved with stationary parameters, in particular electron and impurity densities, with the upper water-cooled tungsten divertor.

Goniche, M.↗
DOE OSTI2025

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗
DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗
DOE OSTI2024

Turbulent drifts of impurity ions as an explanation for anomalous radial transport in the far-SOL of DIII-D

Abstract Successful fusion reactor operation relies on minimal core contamination by impurities, otherwise too much power may be radiated and harm performance. This requires reliable predictions of impurity transport from the scrape-off layer (SOL) into the core, beyond the traditional ‘anomalous’ diffusion approach. We report a set of far-SOL tungsten transport simulations that demonstrate the role of turbulent drifts on radial impurity transport. A turbulent plasma background is simulated using the gyrokinetic SOL code Gkeyll. Tungsten ions are followed within the plasma background using only their drifts. We find that tungsten tends to travel radially outwards with velocities between v r = 300–1200 m s −1 primarily due to polarization drift. We also extract an anomalous radial diffusion coefficient that varies from D r anom = 5–20 m 2 s −1 . These results are compared to and agree with previous interpretive modeling results. We also show how the turbulent polarization drift can transport some tungsten ions from the wall inwards with effective pinch velocities up to 10 000 m s −1 . We conclude that turbulent drifts are a likely explanation for historically anomalous radial impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Development of an Integrated Multidiagnostic to Assess the High-Z Impurity Fluxes in the Metallic Environment of WEST Using IMAS

WEST is an actively cooled, long-pulse tokamak with nearly all plasma-facing components (PFC) made of tungsten. One of the aims of WEST is to study plasma operations with tungsten PFCs in preparation for long-pulse operations on high-Z divertor devices, such as ITER. For long-pulse operation, the high-Z impurity content and transport to the core plasma are critical concerns that require further measurement and interpretation to improve plasma performance and PFC durability. This work details the impurity influxes in WEST during a series of discharges in which the lower hybrid (LH) injected power was incrementally increased. An analysis has been performed of measurements collected from an array of edge diagnostics. Visible spectroscopy was utilized to measure the spectral radiances generated by fuel particles (D) and impurities (W, O, and C) at the divertor and at the antennas with a newly developed spectral peak-fitting tool used to analyze the data in WEST. The scrape-off layer (SOL) plasma conditions (n e and T e ) measured at the divertor target with flush Langmuir probes and near the outer mid-plane (OMP) using reciprocating Langmuir probes (RCPs) are used to evaluate the number of ionizations per photons (S/XB) coefficients required to estimate the impurity fluxes obtained with the collisional-radiative model ColRadPy. Here, the array of edge diagnostics discussed in this work, coupled with SOL plasma modeling tools, represents a multidiagnostic interpretative modeling workflow that will continue to be applied to upcoming experimental campaigns on the WEST experiment to assess high-Z impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A toolbox for $$q_{T}$$ and 0-jettiness subtractions at $$\hbox {N}^3\hbox {LO}$$

We derive the leading-power singular terms at three loops for both $q_T$ and 0-jettiness, $T_0$, for generic color-singlet processes. Our results provide the complete set of differential subtraction terms for $q_T$ and $T_0$ subtractions at N 3 LO, which are an important ingredient for matching N 3 LO calculations with parton showers. We obtain the full three-loop structure of the relevant beam and soft functions, which are necessary ingredients for the resummation of $q_T$ and $T_0$ at N 3 LL' and N 4 LL order, and which constitute important building blocks in other contexts as well. The nonlogarithmic boundary coefficients of the beam functions, which contribute to the integrated subtraction terms, are not yet fully known at three loops. By exploiting consistency relations between different factorization limits, we derive results for the $q_T$ and $T_0$ beam function coefficients at N 3 LO in the z→1 threshold limit, and we also estimate the size of the unknown terms beyond threshold.

Physics↗
DOE OSTI2021

Cost-Optimized Cold Climate Heat Pump Development and Field Test

Cold climate heat pumps (CCHPs) expand the heat pump market to climates where heating demand is dominant. They can achieve more than 70% energy savings compared with electric resistance heating and operate at lower cost than using tank-stored propane to fuel a furnace. A high-efficiency heat pump with a heating seasonal performance factor (HSPF)—as defined by the Air-Conditioning, Heating, and Refrigeration Institute (AHRI 2008)—greater than 10.0 would be more efficient than gas heating in terms of source energy. However, developing a cost-effective CCHP involves some challenges. A typical single-speed, air-source heat pump (ASHP) with an HSPF of 7.7 Btu/Wh does not work well under cold outdoor temperature conditions typical of cold climate locations for four major reasons: 1. Discharge temperature is too high—The low suction pressure and high compression pressure ratio at low ambient temperatures causes high compressor discharge temperatures in excess of the maximum limit for many of the compressors on the market. Furthermore, system charge of a heat pump is usually optimized in the cooling mode, which leads to overcharge conditions in the heating mode, further increasing the discharge temperature. 2. Heating capacity is insufficient if sized to meet the building design cooling load— Heating capacity of a single-speed heat pump decreases with ambient temperature. The heating capacity at -13°F (-25°C) typically decreases to 20%–40% of the rated heating capacity at 47°F (8.3°C) (~equivalent to the rated cooling capacity at 95°F [35°C]). Therefore, a single-speed heat pump, typically sized to match the building design cooling load, cannot provide adequate heating capacity to match the building heating load at low ambient temperatures. The capacity deficit is filled by inefficient resistance heat, thus lowering the system efficiency and significantly increasing power demand. 3. Cyclic loss is significant if sized to meet the building design heating load—If a single-speed heat pump is sized to meet the heating load, it will be significantly oversized relative to the cooling load in many cold climates. This will cause excessive on/off cyclic loss during the cooling and heating operations at moderately low ambient temperatures. Thus, capacity modulation capability (e.g., using a variable-speed or multi-stage compressor) is necessary for a CCHP, which uses its full capacity to meet the peak heating load and partial capacity to meet the cooling and part-load heating loads. 4. Coefficient of performance (COP) is low—Heating COP degrades significantly at low ambient temperatures owing to the large temperature difference between the heat source and sink. A target CCHP should be sized to meet the building design heating load while minimizing the cyclic loss for the cooling and heating operations at moderate ambient temperatures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Thermoelectric Performance of the 2D Bi 2 Si 2 Te 6 Semiconductor

Bi 2 Si 2 Te 6 , a 2D compound, is a direct band gap semiconductor with an optical band gap of 0.25 eV, and is a promising thermoelectric material. Single-phase Bi 2 Si 2 Te 6 is prepared by a scalable ball-milling and annealing process and the highly densified polycrystalline samples are prepared by spark plasma sintering. Bi 2 Si 2 Te 6 shows a p-type semiconductor transport behavior and exhibits an intrinsically low lattice thermal conductivity of ~0.48 Wm -1 K -1 (cross-plane) at 573 K. The first-principles density functional theory calculations indicate that such low lattice thermal conductivity is derived from the interactions between acoustic phonons and low-lying optical phonons, local vibrations of Bi, the low Debye temperature and strong anharmonicity result from the unique 2D crystal structure and metavalent bonding of Bi 2 Si 2 Te 6 . The Bi 2 Si 2 Te 6 exhibits an optimal figure of merit ZT of ~0.51 at 623 K, which can be further enhanced by the substitution of Bi with Pb. Pb doping leads to a large increase in power factor S 2 σ, from ~4.0 μWcm -1 K- 2 of Bi 2 Si 2 Te 6 to ~8.0 μWcm -1 K -2 of Bi 1.98 Pb 0.02 Si 2 Te 6 at 775 K, owing to the increase in carrier concentration. Moreover, Pb doping induces a further reduction in the lattice thermal conductivity to ~0.38 Wm -1 K -1 (cross-plane) at 623 K in Bi 1.98 Pb 0.02 Si 2 Te 6 ,. The simultaneous optimization of the power factor and lattice thermal conductivity achieves a peak ZT of ~0.90 at 723 K and a high average ZT of ~0.66 at 400–773 K in Bi 1.98 Pb 0.02 Si 2 Te 6 .

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2026

Locating the QCD critical point through contours of constant entropy density

We propose a new method to investigate the existence and location of the conjectured high-temperature critical point of strongly interacting matter via contours of constant entropy density. By approximating these lines as a power series in the baryon chemical potential 𝜇 𝐵 , one can extrapolate them from first-principle results at zero net-baryon density, and use them to locate the quantum chromodynamics (QCD) critical point, including the associated first-order and spinodal lines. As a proof of principle, we employ currently available continuum-extrapolated lattice data from the Wuppertal-Budapest collaboration to find a critical point at a temperature and a baryon chemical potential of 𝑇 𝑐 = 114.3 ± 6.9 MeV and 𝜇 𝐵,𝑐 = 602.1 ± 62.1 MeV, respectively, at expansion order 𝒪⁡(𝜇$^{2}_{𝐵}$). We advocate for a more precise determination of the required expansion coefficients via lattice QCD simulations as a means of pinpointing the location of the critical endpoint in the phase diagram of strongly interacting matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗