SEARCH · Engineering Papers

Results for “Single-phase approach”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Boundary Element Model for Assessing Large‐Scale Pressurization in Faulted Geological Storage Systems

Assessing large-scale pressurization at the regional scale—a possible outcome of large subsurface storage applications such as wastewater injection and geological carbon sequestration—presents significant computational challenges. These challenges are particularly pronounced when accounting for complex geologic structures with multiple reservoir and caprock layers, fault zones, and wells. This study introduces a computationally efficient model that integrates single-phase semi-analytical solutions with a boundary element (BE) approach. The model simulates pressure propagation in multilayered 3D systems, including vertical faults, caprock, basement, and confining units. We apply this new model to a representative scenario involving CO 2 injection near a partially sealing fault with verification against an independent two-phase flow model. Results demonstrate that our model accurately captures far-field pressure responses and that, outside the CO 2 plume zone, pressure predictions from single-phase and two-phase models are nearly identical. This supports the use of single-phase models like ours for efficient estimation of far-field pressure changes. Additionally, we demonstrate its effectiveness at a large scale, incorporating multiple wells and faults. With its ability to represent multiple wells, fault zones, and geological heterogeneity, our model is well suited for assessments of basin-scale pressurization. Its computational efficiency also makes it a promising tool for integration with optimization frameworks aimed at designing and managing injection strategies in faulted storage systems.

Cihan, A. [Lawrence Berkeley National Laboratory (↗

DOE OSTI2025

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Holistic energy analysis method for thermal management architectures of data centers

Modern high-performance computing (HPC) data centers (DCs), particularly those supporting energy-intensive artificial intelligence (AI) workloads, face escalating thermal management challenges that degrade performance through thermal throttling and drive up cooling power consumption and operational costs. To address this challenge, many have developed a wide variety of thermal management solutions (single-phase, two-phase, direct, indirect, hybrid, and more) which attempt to cool HPC DCs effectively while attempting to minimize overall system power consumption. However, the analysis of these solutions and methods to effectively compare one with another is lacking. Overall power usage effectiveness (PUE) and total-power usage effectiveness (TUE) provide a metric to quantify power consumption but fail to identify components in the system which require further optimization. To address this, we propose a holistic analytical framework – the waterfall diagram (WFD) – which leverages a waterfall chart methodology, offering a comprehensive visualization of both the thermal management system loop and heat flow pathways from individual server components to the outdoor ambient. Use of the WFD enables graphical estimations of power efficiency and cooling performance across each component of a DC cooling system and complements Sankey-style energy flow visualizations by additionally resolving stage-wise temperature changes and incremental TUE contributions. The framework is used in conjunction with simulation-based approaches, to conduct a detailed pressure drop and flow distribution analysis aimed at identifying the optimal coolant distribution architecture for a single-phase direct-to-chip water-cooled DC, which serves as the baseline for subsequent WFD analysis. Among the evaluated architectures, the 3 U modular coolant distribution architecture is found to demonstrate the best performance, considering minimal pressure drop and uniform flow distribution. In addition, TUE is calculated for each cooling loop component based on its associated pressure drop and corresponding pumping power, which are integrated into the WFD. This correlation between TUE and local temperature offers immediate insight into the power efficiency and thermal performance contributions of individual components, facilitating further development and optimization. Examples of WFD applications are presented under varying thermal loads and ambient conditions, demonstrating reasonable cooling strategies. Notably, the 3 U modular architecture maintains a consistent chip case temperature of 85°C, achieving a TUE of 1.016 at ambient temperature of 47°C, and a TUE of 1.026 at ambient temperature of 52°C. The WFD methodology provides an efficient, holistic, and streamlined framework for DC thermal management architecture assessment and enables design optimization which is important for addressing the thermal-fluidic energy challenges of current and next-generation DCs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Achieving high tensile strength and ductility in refractory alloys by tuning electronic structure

The energy efficiency of heat engines (gas and steam turbines) for electricity production and propulsion is determined by the Carnot cycle and scales with operating temperature. Commercial nickel- and cobalt-based superalloys melt near 1,500 °C and rapidly lose mechanical strength beyond 1,000 °C. Refractory metals melt well above 2,000 °C but have inherent manufacturability challenges that are barriers to adoption, such as high ductile-to-brittle transition temperatures. Using density functional theory-guided design, we demonstrate tailored local lattice distortions that promote phase-stable, non-equiatomic refractory concentrated solid solutions with both high ductility and strength. Here, we exemplify this for single-phase, body-centred cubic Nb 4 Ta 4 V 3 Ti that exhibits castability, excellent room-temperature tensile yield strength (∼1 GPa) and ductility (approaching 20% uniform strain), and exceptional high-temperature tensile strength (500 MPa at 1,000 °C). These findings illustrate a path for designing materials that hold great potential for advancing next-generation technologies such as Generation IV fission reactors, first-generation fusion-plasma reactors, and more efficient gas turbines for electricity generation and propulsion.

DFT↗

DOE OSTI2025

Microstructurally Strained Pyrochlore–Perovskite Biphasic Electrocatalysts for the Oxygen Evolution Reaction

Efficiency of water splitting for hydrogen production is often limited by the sluggish kinetics of multiple electronic transfers required in the heterogeneous oxygen evolution reaction (OER). Catalyst design for reducing the high OER overpotential remains a major scientific challenge. Lattice-strain engineering, a method for tuning the electronic structure and surface geometric configuration of active sites, may greatly affect the interaction between adsorbates and catalytic surfaces for high activity and stability. Here, in this study, we present the synthesis of biphasic oxides of YPrSrRuMnO x , which consists of distinct phases of Y 2 Ru 2 O 7 pyrochlore and (Pr 0.7 Sr 0.3 )MnO 3 perovskite, and the development of a suitable analytical approach to study the strain–catalytic property relationship. Linear sweep voltammetry results reveal that the biphasic oxide exhibits approximately 3.1 times greater mass activity and 2.4 times larger turnover frequency (TOF) than single-phase Y 2 Ru 2 O 7 in the 0.1 M HClO 4 electrolyte. The biphasic catalyst is also about 3 times more stable than the single-phase oxide under acidic conditions. X-ray photoelectron spectroscopy, nitrogen isotherm, and electrochemical surface area analyses indicate that the oxidation state, specific surface area, and electrochemical surface area do not cause enough difference in the observed enhancement of OER performance. We examined the effects of microstrain on electrocatalysis, originating from lattice mismatch between different phases, using three different structural models. Specifically, we compared the Williamson–Hall method, standard stress–strain analysis, and Rietveld refinement in analyzing the structure–property relationship. Strain mapping using geometric phase analysis (GPA) further revealed significant microstrain and lattice dislocations localized near phase boundaries in the biphasic oxide, in contrast to the uniform strain in single-phase materials. The results reveal that the increased microstrain correlates well with the improved OER performance, as the biphasic oxide catalyst exhibits 2–3 times greater microstrain than Y 2 Ru 2 O 7 pyrochlore.

electrocatalysts↗

DOE OSTI2025

A General Strategy for Bandgap Engineering Via Anion‐Lattice Doping in High‐Entropy Oxides

Bandgap engineering is a critical tool for tailoring the electronic properties of functional materials, traditionally achieved by modifying the cation sublattice. Here, a generalizable strategy is introduced that leverages facile anion-lattice doping in high entropy materials to modulate the bandgap in high-entropy metal oxides (HEMOs). By incorporating nitrogen into a single-phase high-entropy metal oxide/nitride (HEMO:HEMN) solid solution, a substantial bandgap reduction is achieved from 3.55 eV (HEMO) to ≈2.46 eV (HEMO:HEMN), significantly enhancing electronic conductivity. Unlike conventional bandgap tuning approaches that rely on cation substitution or heterojunction formation, this method exploits anion-mediated entropy stabilization, enabling uniform bandgap narrowing across the entire solid solution. This anion-lattice engineering strategy is broadly applicable to high-entropy systems, providing a new pathway for designing energy materials with tailored electronic properties. The resulting HEMO:HEMN solid solution exhibits a tenfold increase in capacitance and capacity compared to HEMO in supercapacitor and lithium-ion battery tests, demonstrating the transformative potential of anion-driven bandgap modulation for next-generation energy storage and conversion technologies.

energy storage↗

DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗

DOE OSTI2025

Pseudo-viscous modeling of transport in dense granular flows for thermal energy storage applications

Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.

14 SOLAR ENERGY↗

DOE OSTI2025

Discovery, Design, Synthesis and Testing of High Performance Structural Alloys (Final Technical Report)

The overarching goal of this project is to understand the phase stability and mechanical behavior of non-stoichiometric multi-principal element alloy (MPEA) materials. In order to identify suitable alloys, we plan to use a combinatorial thin film screening approach, in collaboration with scientists at Lawrence Berkeley National Laboratory who are performing computational work as well as complementary experimental work. Specific tasks within the scope of this project include the fabrication, using thin film deposition from six sputtering targets, of combinatorial samples with multi-dimensional gradients in composition and microstructure. These samples are studied to screen MPEA systems for promising candidate alloys with specific composition(s), based on characterization of composition, structure and mechanical behavior across the thin film. We want to produce single-phase MPEAs with chemical homogeneity in a given thin film region, simple grain structures, and no intermetallic phases present. Gradient films facilitate first-pass screening for desirable characteristics and inform the next stage of work that involves fabrication of bulk MPEA specimens for (tensile) mechanical testing and characterization. To make the bulk alloys, metal (elemental) pieces are melted to form MPEAs, followed by heat treatment to homogenize the composition and microstructure. A subset of alloys is also cast, using vacuum arc melting, to yield larger samples (diameter ~1 cm and length ~5-10 cm) and these allow us to assess viability of scale-up for the alloys in structural applications. Further processing plans include rolling and heat treatment to recrystallize selected bulk MPEAs and grow grains to different extents, in order to investigate size effects in the mechanical behavior of MPEAs. Microspecimen testing will be performed (primarily in tension) to assess the mechanical behavior over a range of temperatures. The deformation microstructure of mechanically tested alloys will be characterized using transmission electron microscopy (TEM) to provide a scientific basis for understanding the structure-property relationships in MPEA mechanical behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Multiscale Thermal-hydraulic analysis of the MARVEL microreactor using a coupled SAM and SCM simulation. (PPTX)

This is .pptx document presenting a summary of the paper of the same name that has been submitted and accepted to NURETH-21: MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor's primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL?s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

DOE OSTI2025

Facile One-Pot Block Copolymer-Mediated Solvothermal Approach for Synthesis of High-Entropy Alloy with Enhanced OER Activity

Herein, we introduce a straightforward synthesis approach for highly active dendritic multimetallic high-entropy alloy (DMHEA@PtIrPdAgRu) nanoparticles with sufficient entropic mixing, featuring uniform distribution of five noble group metals (Pt, Ir, Pd, Ag, and Ru) via a block copolymer-mediated one-pot solvothermal reduction method for oxygen evolution reaction (OER). In this synthesis, N,N -dimethylformamide (DMF) is used as a reductant as well as solvent and core–shell-corona-type (poly(styrene)- block -poly(vinylpyridine)- block -poly(ethylene oxide)) (PS-PVP-PEO) block copolymer as a structure directing agent. The cooperative effect between the copolymer architecture and the reducing environment of DMF promoted a confined nucleation mechanism for forming a single-phase dendritic structure HEA with high compositional uniformity, thereby mitigating phase segregation, a common challenge in the synthesis of multimetallic nanoparticles. This prepared DMHEA@PtIrPdAgRu catalyst exhibits a low overpotential of 490 mV to attain a high current density of 100 mA cm –2 with a Tafel slope of 442 mV dec –1 for oxygen evolution. The superior OER performance is attributed to the synergistic cooperation among its active and coordinated metal centers as well as the incorporation of corrosion-resistant metal like platinum.

block copolymer↗

DOE OSTI2025

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗

DOE OSTI2026

Phase change material-to-refrigerant heat exchangers: Experimental validation and uneven melting analysis

Thermal energy storage (TES) using phase change materials (PCMs) enables load shifting and reduces mismatches between the building's thermal demand and the heat pump (HP) system's thermal capacity. While most TES studies have focused on single-phase heat transfer fluids, PCM-to-refrigerant heat exchangers (PRHX), in which the refrigerant undergoes liquid-vapor phase change, remain largely underexplored. Here, to address this gap, this study systematically investigates a shell-and-tube PRHX operating as the condenser in an HP through experiments and simulations. A two-dimensional enthalpy-based finite-volume PRHX model was developed and validated against experimental data, including refrigerant and PCM temperature profiles, PRHX capacity, and system-level performance. Results show three characteristic stages: Stage I, the condenser temperature rapidly increases with PCM in the solid state absorbing sensible heat; Stage II, stable operation during PCM phase change; and Stage III, performance decline once part of the PCM becomes fully melted. Heat transfer analysis revealed that uneven PCM melting along the condenser length was driven primarily by variations in the refrigerant-side heat transfer coefficient, rather than local approach temperature differences. Further study demonstrated that minimizing refrigerant outlet subcooling led to more uniform PCM melting and extended Stages I and II by up to 89%. This paper provides a validated PRHX model, clarifies the mechanisms of uneven PCM melting, and highlights subcooling control as an effective strategy to improve TES performance in HP systems without requiring secondary loops.

Graphite matrix↗

DOE OSTI2026

Deterministic Fabrication of Large-Area, High-Crystallinity Oxide Moiré Superlattices

Oxide twistronics extends moiré engineering beyond van der Waals materials, offering a promising platform for accessing emergent interfacial phenomena arising from the strong coupling of lattice, charge, and orbital degrees of freedom in complex oxides. However, deterministic fabrication of high-crystallinity oxide moiré superlattices over large lateral dimensions remains challenging due to the three-dimensional bonding network of oxides. Here, we demonstrate a scalable, generalized fabrication strategy that enables the formation of high-crystallinity oxide moiré superlattices with clean, chemically bonded interfaces and precisely controlled twist angles down to nominal values of 0.1°, achieving subdegree twist-angle accuracy across large contiguous lateral dimensions approaching the millimeter scale. Using NaNbO3 as a model system, we show that the resulting interlayer coupling drives pronounced structural reconstruction that modifies both the phase structure and ferroelectric domain configuration. Synchrotron-based X-ray 3D reciprocal space mapping reveals the emergence of a single-phase state in twisted bilayers, in contrast to the mixed-phase structure observed in single-layer membranes prior to twist assembly. The structural signatures are further consistent with gradual lattice rotation distributed along the thickness direction that may accommodate interfacial shear strain, distinct from reconstruction observed in van der Waals moiré systems which primarily occurs through in-plane stacking rearrangement. This collective lattice response is correlated with twist-dependent nanoscale electromechanical modulations observed by piezoresponse force microscopy. These results establish a scalable materials platform for oxide twistronics and support the implementation of twist-engineered functionalities in practical, macroscale device architectures.

Ghanbari, Reza [North Carolina State University (N↗