Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field simulations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reactive Carbide‐Based Synthesis and Microstructure of NASICON Sodium Metal All Solid‐State Electrolyte

Reactive carbide precursor-based synthesis of NASICON-type NZSP (Na 1+x Zr 2 Si x P 3-x O 12 ) solid-state electrolyte (SSE) is demonstrated, in contrast to the established oxide-based approach. Exothermic decomposition of ZrC and SiC in air homogenizes microstructure, yielding 98% compact density after conventional sintering at 1200 °C. Quantitative stereology demonstrates that significant microstructural differences are present. Compacts of carbide-derived Carb-NZSP are 98% dense with a secondary zirconium oxide (ZrO 2 ) volume fraction of 0.2% ± 0.3%, versus 93% dense and 3% ± 1% for oxide-derived baseline. For Carb-NZSP, the secondary glassy phosphate phase is agglomerated, while for baseline, it is dispersed and percolated. Electrochemical testing combined with post-mortem analysis demonstrates how microstructural control of secondary phases is critical for dendrite suppression: Carb-NZSP critical current density (CCD) is 3.1 ± 0.8 mA cm −2 at 0.1 mAh cm −2 , versus 1.0 ± 0.7 mA cm −2 at 0.1 mAh cm −2 . Cryogenic focused ion beam (cryo-FIB) analysis demonstrates that in both materials, the porous 2D sheet-like sodium metal dendrites propagate around and subsume NZSP grains, likely following a path enriched with glassy phase and with porosity. Dendrites also flow around isolated zirconia particles. Phase field simulation reveals deflection of dendrites by mechanically tough zirconia, while brittle glassy phase accelerates dendrite growth, especially when finely distributed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Strain phase equilibria and phase‐field method of ferroelectric polydomain: A case study of monoclinic K x Na 1 − x NbO 3 thin films

Abstract Knowledge of the thermodynamic equilibria and domain structures of ferroelectrics is critical to establishing their structure–property relationships that underpin their applications from piezoelectric devices to nonlinear optics. Here, we establish the strain condition for strain phase separation and polydomain formation and analytically predict the corresponding domain volume fractions and wall orientations of, relatively low symmetry and theoretically more challenging, monoclinic ferroelectric thin films by integrating thermodynamics of ferroelectrics, strain phase equilibria theory, microelasticity, and phase‐field method. Using monoclinic K x Na 1 − x NbO 3 (0.5 < x < 1.0) thin films as a model system, we establish the polydomain strain–strain phase diagrams, from which we identify two types of monoclinic polydomain structures. The analytically predicted strain conditions of formation, domain volume fractions, and domain wall orientations for the two polydomain structures are consistent with phase‐field simulations and in good agreement with experimental results in the literature. The present study demonstrates a general, powerful analytical theoretical framework to predict the strain phase equilibria and domain wall orientations of polydomain structures applicable to both high‐ and low‐symmetry ferroelectrics and provide fundamental insights into the equilibrium domain structures of ferroelectric K x Na 1 − x NbO 3 thin films that are of technology relevance for lead‐free dielectric and piezoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Asymmetry of the Ferroelectric Phase Transition in BaTiO 3

Symmetry changes during phase transformations fundamentally determine the behavior of thermodynamic systems, governing phenomena as diverse as water evaporation, fermion condensation, and epidemic spreading. Phase transitions are conventionally divided into two classes, and this classification is typically assumed to remain invariant upon reversing the transition. Here, an asymmetric phase transformation is uncovered in the ferroelectric–paraelectric transition of single-crystal BaTiO 3 , a model system for first-order transitions. Under slow temperature variation (≤0.1 °C min −1 ), thermodynamic, dielectric, and domain-structure measurements reveal that the ferroelectric-to-paraelectric transition exhibits latent heat, phase coexistence, and a discontinuous order parameter, while these signatures are absent upon cooling. Complementary phase-field simulations demonstrate similar behavior, attributing it to distinct elastic strain energy accumulation and release during heating and cooling. These findings reveal a first-order character upon heating but a second-order-like behavior upon cooling, challenging the conventional paradigm of symmetric phase-transition classification and suggesting new possibilities for ferroelectric-based energy storage.

Hershkovitz, Asaf [Technion-Israel Institute of Te↗
DOE OSTI2026

Nonequilibrium defect-phase nanostructures stabilized by irradiation in undersaturated Ni-Si nanocrystalline alloy

Nanocrystalline thin films of the undersaturated alloy Ni-8.5 at% Si were subjected to 2 MeV Ti irradiation at temperatures ranging from 450˚C to 550˚C. Correlative microscopy combining transmission electron microscopy (TEM), scanning-TEM and atom probe tomography (APT revealed that large dose irradiation at 450˚C of samples with initial grain sizes below 100 nm stabilized a novel nanostructure which surprisingly contained three co-existing phases, the γ face-centered-cubic (FCC) matrix, γ' L12 ordered precipitates on intragranular dislocation loops and Ni 31 Si 12 precipitates at triple junctions (TJs). In contrast, irradiation at 550˚C and irradiation of larger grain-size samples at 450˚C only produced a γ-γ' two-phase coexistence. Analysis of the three-phase nanostructure and phase field simulations indicates that radiation-induced segregation is most pronounced at TJs, thus triggering the formation of Ni 31 Si 12 precipitates. These incoherent precipitates, in turn, are expected to stabilize the grain size under irradiation. The results are generalized using the concept of driven defect-phases. It is suggested that the stabilization of driven defect-phases may impart radiation resilience by providing localized relaxation modes to the microstructure evolution during and after temporary perturbations in irradiation conditions.

36 MATERIALS SCIENCE↗
NASA NTRS2019

In Situ Observation of Growth Dynamics in DECLIC Directional Solidification Insert Onboard ISS: DSI-R Flight Campaign

The study of solidification microstructure formation is of utmost importance for materials design and processing, as solid-liquid interface patterns largely govern mechanical and physical properties. Pattern selection occurs under dynamic conditions of growth in which the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The dynamic and nonlinear nature of this instability renders in situ observation of the interface an invaluable tool to gain knowledge on the time-evolution of the interface pattern. Transparent organic analogs, which solidify like metallic alloys, allow direct visualization of interface dynamics. Extensive ground-based studies of both metallic and organic bulk samples have established the presence of significant convection during solidification processes that alters the formation of interfacial microstructures. A reduced-gravity environment is therefore mandatory for fluid flow elimination in bulk samples. In the framework of the CNES project MISOL3D (MIcrostrutures de SOLidification 3D) and the NASA projects DSIP (Dynamical Selection of 3D Interface Patterns), SPADES (SPAtiotemporal Evolution of three-dimensional DEndritic array Structures) and CAMUS (ComputAtional Studies of MicrostrUcture Formation During Alloy Solidification in Microgravity), we participated in the development of the Directional Solidification Insert (DSI) of the DEvice for the study of Critical Liquids and Crystallization (DECLIC). The DECLIC-DSI is dedicated to in situ and real-time characterization of solid-liquid interface patterns during directional solidification of transparent alloys in diffusive transport regime. Between April 2010 and March 2011, the first ISS campaign (DSI) explored the entire range of microstructures resulting in unprecedented observations. A second campaign (DSI-R), performed between October 2017 and December 2018, in which the insert contained an alloy of higher solute concentration, allowed to complete the benchmark database. The increase of solute concentration resulted in well-developed dendritic patterns at lower velocities (lower interface curvature and larger tip radius). The microstructure resulting from dendritic growth is dominant in metallurgy so that it is fundamental to understand the mechanisms of its formation. The main aims of this experimental campaign are to understand: the mechanisms of the cell to dendrite transition, the fundamental mechanisms of sidebranching formation, the dependence of dendrite tip shapes on growth conditions, the interaction of primary array and secondary sidebranches, and the influence of subgrain boundaries on the spatiotemporal organization of the array structure. Preparation, analysis and interpretation of the experiments performed onboard ISS are considerably enhanced by experiments performed on ground using thin-samples (Pr. Trivedi’s group) and phase-field simulations of microstructure formation in a diffuse growth regime (Pr. Karma’s group). In this summary, we will present an initial assessment of the results obtained during the DSI-R campaign. Nomenclature Vp: pulling velocity G: thermal gradient L: solidified length rtip: tip curvature radius Acronyms/Abbreviations CADMOS: Centre d'Aide au Développement des activités en Micro-pesanteur et des Opérations Spatiales CNES: Centre National d’Études Spatiales DECLIC: Device for the study of Critical LIquids and Crystallization DSI: Directional Solidification Insert DSI-R: Directional Solidification Insert - Refurbish ISS: International Space Station NASA: National Aeronautics and Space Administration PF: phase-field SCN: Succinonitrile 3D: three-dimensional

Mota, F.L.↗
DOE OSTI2025

Morphologies of dealloying corrosion attack at grain boundaries

Dealloying corrosion at grain boundaries severely compromises the performance of polycrystalline materials across a wide variety of technological applications. The impact of this phenomenon depends upon the morphology and rate of intergranular dealloying, which can range from planar to wormhole-like patterns that rapidly advance into the alloy. Using 2D and 3D multi-phase field simulations, we reveal how diverse microstructures result from a fundamental interplay between alloy composition and a grain boundary migration mechanism that alters diffusional pathways of dealloying. Inside intergranular dealloying channels, corrosion product buildup can spawn new channels that branch into grain interiors, such that alloys can be degraded from the inside out. These processes further lead to an atypical coarsening mode assisted by diffusion in the dealloying agent. We summarize a unifying explanation for dealloying morphology selection in polycrystalline alloys, which provides an important step towards their optimization for dealloying corrosion environments.

Corrosion↗
DOE OSTI2025

Holistic Microstructure Control Strategies in Photopolymerization‐Induced Phase Separation of Acrylate Systems

Open porous materials, known for their large surface area and interconnected structures, are essential in various applications, including batteries, ion exchange, catalysis, filtration, and electronic waste recycling. A critical aspect of the functionality of porous membranes is the precise control of pore size and morphology. Photopolymerization-induced phase separation (photo-PIPS) offers a convenient and versatile methods for creating porous structures. However, controlling the porous morphology remains challenging due to the complex interplay between thermodynamics, polymerization kinetics, and monomer structures, which makes it difficult to establish the relationship between processing conditions and resulting morphology in photo-PIPS. Herein, a physics-based phase-field model capable of generating and characterizing the microstructures of porous materials based on both average and localized features is developed. Using the phase-field simulations as test bed, the effects of polarity, light intensity, and curing temperature, as well as the previously unexplored roles of chain transfer agents and substrates, on the morphology of the resulting porous microstructure are investigated. Experiments are performed to verify the results predicted by the simulations. This work lays out a comprehensive guide for designing PIPS-derived porous microstructures and offers practical engineering strategies for tailoring microstructure-level topology and size of pores for application-specific needs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2024

Optical Control of Adaptive Nanoscale Domain Networks

Adaptive networks can sense and adjust to dynamic environments to optimize their performance. Understanding their nanoscale responses to external stimuli is essential for applications in nanodevices and neuromorphic computing. However, it is challenging to image such responses on the nanoscale with crystallographic sensitivity. Here, the evolution of nanodomain networks in (PbTiO 3 ) n /(SrTiO 3 ) n superlattices (SLs) is directly visualized in real space as the system adapts to ultrafast repetitive optical excitations that emulate controlled neural inputs. The adaptive response allows the system to explore a wealth of metastable states that are previously inaccessible. Their reconfiguration and competition are quantitatively measured by scanning x-ray nanodiffraction as a function of the number of applied pulses, in which crystallographic characteristics are quantitatively assessed by assorted diffraction patterns using unsupervised machine-learning methods. The corresponding domain boundaries and their connectivity are drastically altered by light, holding promise for light-programable nanocircuits in analogy to neuroplasticity. Phase-field simulations elucidate that the reconfiguration of the domain networks is a result of the interplay between photocarriers and transient lattice temperature. The demonstrated optical control scheme and the uncovered nanoscopic insights open opportunities for the remote control of adaptive nanoscale domain networks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unlocking Electrostrain in Plastically Deformed Barium Titanate

Achieving substantial electrostrain alongside a large effective piezoelectric strain coefficient (d 33 *) in piezoelectric materials remains a formidable challenge for advanced actuator applications. Here, in this work, a straightforward approach to enhance these properties by strategically designing the domain structure and controlling the domain switching through the introduction of arrays of ordered {100}<100> dislocations is proposed. This dislocation engineering yields an intrinsic lock-in steady–state electrostrain of 0.69% at a low field of 10 kV cm -1 without external stress and an output strain energy density of 5.24 J cm -3 in single-crystal BaTiO 3 , outperforming the benchmark piezoceramics and relaxor ferroelectric single-crystals. Additionally, applying a compression stress of 6 MPa fully unlocks electrostrains exceeding 1%, yielding a remarkable d 33 * value over 10 000 pm V -1 and achieving a record-high strain energy density of 11.67 J cm -3 . Optical and transmission electron microscopy, paired with laboratory and synchrotron X-ray diffraction, is employed to rationalize the observed electrostrain. Phase-field simulations further elucidate the impact of charged dislocations on domain nucleation and domain switching. These findings present an effective and sustainable strategy for developing high-performance, lead-free piezoelectric materials without the need for additional chemical elements, offering immense potential for actuator technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Transient Terahertz Oscillations During Photoinduced Polarization Topology Reconfiguration in Ferroelectric Superlattices

Terahertz resonances embedded in crystalline heterostructures could close a spectral gap between conventional electronics and photonics while opening new windows on non-equilibrium lattice dynamics. We show that femtosecond optical screening of the depolarization field in epitaxial PbTiO3/SrTiO3 superlattices launches a collective polar mode that oscillates near 1 THz and coherently spans the entire mini-Brillouin zone. Wave-vector-resolved pump–probe X-ray diffraction resolves a nearly dispersion-less oscillation at 0.87 THz and 0.94 THz at the zone boundary and zone center, respectively, persisting for ~2.5 ps, corresponding to a weakly damped resonance. Dynamical phase-field simulations reveal the origin of the mode to mesoscopic rotation of closure-domain textures during the photo-excited transition from an unscreened to a screened electrostatic state. Varying the PbTiO3 and SrTiO3 ratio tunes the mode frequency continuously from 0.9 to 1.4 THz, providing a quantitative design rule for frequency-selectable THz oscillators in ferroelectric heterostructures. By coupling nanoscale polarization reconfiguration to long-wavelength coherent dynamics, this work establishes depolarization-field engineering to topology-driven THz functionality and expanding the landscape of collective lattice dynamics.

Sri Gyan, Deepankar [Univ. of Wisconsin, Madison, ↗
DOE OSTI2024

Electrical de-poling and re-poling of relaxor-PbTiO 3 piezoelectric single crystals without heat treatment

Re-poling of unexpected partially depoled piezoelectric materials conventionally needs to be first fully depoled through annealing above their Curie temperature to revive piezoelectric performances. Here, we investigated de-poling and re-poling of Pb(In 1/2 Nb 1/2 )O 3 -Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 single crystals under electric fields at room temperature. We found that alternating current electric fields with amplitudes near the coercive field at low frequencies (<10 Hz) can be employed to successfully depolarize poled crystals at room temperature. We also demonstrated a reversible polarization switching process with a relaxor-PbTiO 3 single crystal ultrasound transducer without device performance degradations. This experimental observation is supported by phase-field simulation, showing that alternating current electric fields can readily induce de-poling at room temperature, while direct current electric fields induce a transient depoled state only within an uncontrollable short period of time. The findings suggest new strategies for unprecedented in-device tailoring of the polarization states of ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reversible long-range domain wall motion in an improper ferroelectric

Reversible ferroelectric domain wall movements beyond the 10 nm range associated with Rayleigh behavior are usually restricted to specific defect-engineered systems. Here, we demonstrate that such long-range movements naturally occur in the improper ferroelectric ErMnO 3 during electric-field-cycling. We study the electric-field-driven motion of domain walls, showing that they readily return to their initial position after having traveled distances exceeding 250 nm. By applying switching spectroscopy band-excitation piezoresponse force microscopy, we track the domain wall movement with nanometric spatial precision and analyze the local switching behavior. Phase field simulations show that the reversible long-range motion is intrinsic to the hexagonal manganites, linking it to their improper ferroelectricity and topologically protected structural vortex lines, which serve as anchor point for the ferroelectric domain walls. Our results give new insight into the local dynamics of domain walls in improper ferroelectrics and demonstrate the possibility to reversibly displace domain walls over much larger distances than commonly expected for ferroelectric systems in their pristine state, ensuring predictable device behavior for applications such as tunable capacitors or sensors.

Zahn, Manuel [Norwegian Univ. of Science and Techn↗
DOE OSTI2026

Parity breaking at faceted crystal growth fronts during ice templating

Directional solidification of water-based solutions has emerged as a versatile technique to template hierarchical porous materials, but this nonequilibrium process remains incompletely understood. Here, in this study, we use phase-field simulations to shed light on the mechanism that selects the growth direction of the lamellar ice structure that templates those materials. Our results show that this selection can be understood within the general framework of spontaneous parity breaking, yielding quantitative predictions for the tilt angle of lamellae with respect to the thermal axis. The results provide a theoretical basis to interpret a wide range of experimental observations.

Materials science↗
DOE OSTI2024

Elucidating the complex interplay between thermodynamics, kinetics, and electrochemistry in battery electrodes through phase-field modeling

Abstract This article highlights applications of phase-field modeling to electrochemical systems, with a focus on battery electrodes. We first provide an overview on the physical processes involved in electrochemical systems and applications of the phase-field approach to understand the thermodynamic and kinetic mechanisms underlying these processes. We employ two examples to highlight how realistic thermodynamics and kinetics can naturally be incorporated into phase-field modeling of electrochemical processes. One is a composite battery cathode with an intercalation compound (Li x FePO 4 ) as the electrochemically active material, and the other is a displacement reaction compound (Li–Cu–TiS 2 ). With the input parameters mostly from atomistic calculations and experimental measurements, phase-field simulations allowed us to untangle the interactions among transport, reaction, electricity, chemistry, and thermodynamics that lead to highly complex evolution of the materials within battery electrodes. The implications of these observations for battery performance and degradation are discussed. Graphical abstract

Andrews, W. Beck (ORCID:0000000287824621)↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗