Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase field simulation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Emergence of high piezoelectricity from competing local polar order-disorder in relaxor ferroelectrics

Relaxor ferroelectrics are known for outstanding piezoelectric properties, finding a broad range of applications in advanced electromechanical devices. Decoding the origins of the enhanced properties, however, have long been complicated by the heterogeneous local structures. Here, we employ the advanced big-box refinement method by fitting neutron-, X-ray-based total scattering, and X-ray absorption spectrum simultaneously, to extract local atomic polar displacements and construct 3D polar configurations in the classical relaxor ferroelectric Pb(Mg 1/3 Nb 2/3 )O 3 –PbTiO 3 . Our results demonstrate that prevailing order-disorder character accompanied by the continuous rotation of local polar displacements commands the composition-driven global structure evolution. The omnidirectional local polar disordering appears as an indication of macroscopic relaxor characteristics. Combined with phase-field simulations, it demonstrates that the competing local polar order-disorder between different states with balanced local polar length and direction randomness leads to a flattening free-energy profile over a wide polar length, thus giving rise to high piezoelectricity. Our work clarifies that the critical structural feature required for high piezoelectricity is the competition states of local polar rather than relaxor.

42 ENGINEERING↗
DOE OSTI2023

Emergent chirality in a polar meron to skyrmion phase transition

Polar skyrmions are predicted to emerge from the interplay of elastic, electrostatic and gradient energies, in contrast to the key role of the anti-symmetric Dzyalozhinskii-Moriya interaction in magnetic skyrmions. Here, we explore the reversible transition from a skyrmion state (topological charge of –1) to a two-dimensional, tetratic lattice of merons (with topological charge of –1/2) upon varying the temperature and elastic boundary conditions in [(PbTiO 3 ) 16 /(SrTiO 3 ) 16 ] 8 membranes. This topological phase transition is accompanied by a change in chirality, from zero-net chirality (in meronic phase) to net-handedness (in skyrmionic phase). We show how scanning electron diffraction provides a robust measure of the local polarization simultaneously with the strain state at sub-nm resolution, while also directly mapping the chirality of each skyrmion. Using this, we demonstrate strain as a crucial order parameter to drive isotropic-to-anisotropic structural transitions of chiral polar skyrmions to non-chiral merons, validated with X-ray reciprocal space mapping and phase-field simulations.

42 ENGINEERING↗
DOE OSTI2023

Cell invasion during competitive growth of polycrystalline solidification patterns

Abstract Spatially extended cellular and dendritic array structures forming during solidification processes such as casting, welding, or additive manufacturing are generally polycrystalline. Both the array structure within each grain and the larger scale grain structure determine the performance of many structural alloys. How those two structures coevolve during solidification remains poorly understood. By in situ observations of microgravity alloy solidification experiments onboard the International Space Station, we have discovered that individual cells from one grain can unexpectedly invade a nearby grain of different misorientation, either as a solitary cell or as rows of cells. This invasion process causes grains to interpenetrate each other and hence grain boundaries to adopt highly convoluted shapes. Those observations are reproduced by phase-field simulations further demonstrating that invasion occurs for a wide range of misorientations. Those results fundamentally change the traditional conceptualization of grains as distinct regions embedded in three-dimensional space.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Giant electric field-induced second harmonic generation in polar skyrmions

Electric field-induced second harmonic generation allows electrically controlling nonlinear light-matter interactions crucial for emerging integrated photonics applications. Despite its wide presence in materials, the figures-of-merit of electric field-induced second harmonic generation are yet to be elevated to enable novel device functionalities. Here, we show that the polar skyrmions, a topological phase spontaneously formed in PbTiO3/SrTiO3 ferroelectric superlattices, exhibit a high comprehensive electric field-induced second harmonic generation performance. The second-order nonlinear susceptibility and modulation depth, measured under non-resonant 800 nm excitation, reach ~54.2 pm V −1 and ~664% V −1 , respectively, and high response bandwidth (higher than 10 MHz), wide operating temperature range (up to ~400 K) and good fatigue resistance (>10 10 cycles) are also demonstrated. Through combined in-situ experiments and phase-field simulations, we establish the microscopic links between the exotic polarization configuration and field-induced transition paths of the skyrmions and their electric field-induced second harmonic generation response. Our study not only presents a highly competitive thin-film material ready for constructing on-chip devices, but opens up new avenues of utilizing topological polar structures in the fields of photonics and optoelectronics.

Nonlinear optics↗
DOE OSTI2024

Manipulating the diffusion energy barrier at the lithium metal electrolyte interface for dendrite-free long-life batteries

Abstract Constructing an artificial solid electrolyte interphase (SEI) on lithium metal electrodes is a promising approach to address the rampant growth of dangerous lithium morphologies (dendritic and dead Li 0 ) and low Coulombic efficiency that plague development of lithium metal batteries, but how Li + transport behavior in the SEI is coupled with mechanical properties remains unknown. We demonstrate here a facile and scalable solution-processed approach to form a Li 3 N-rich SEI with a phase-pure crystalline structure that minimizes the diffusion energy barrier of Li + across the SEI. Compared with a polycrystalline Li 3 N SEI obtained from conventional practice, the phase-pure/single crystalline Li 3 N-rich SEI constitutes an interphase of high mechanical strength and low Li + diffusion barrier. We elucidate the correlation among Li + transference number, diffusion behavior, concentration gradient, and the stability of the lithium metal electrode by integrating phase field simulations with experiments. We demonstrate improved reversibility and charge/discharge cycling behaviors for both symmetric cells and full lithium-metal batteries constructed with this Li 3 N-rich SEI. These studies may cast new insight into the design and engineering of an ideal artificial SEI for stable and high-performance lithium metal batteries.

25 ENERGY STORAGE↗
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 OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructural impacts on ionic conductivity of oxide solid electrolytes from a combined atomistic-mesoscale approach

Abstract Although multiple oxide-based solid electrolyte materials with intrinsically high ionic conductivities have emerged, practical processing and synthesis routes introduce grain boundaries and other interfaces that can perturb primary conduction channels. To directly probe these effects, we demonstrate an efficient and general mesoscopic computational method capable of predicting effective ionic conductivity through a complex polycrystalline oxide-based solid electrolyte microstructure without relying on simplified equivalent circuit description. We parameterize the framework for Li 7- x La 3 Zr 2 O 12 (LLZO) garnet solid electrolyte by combining synthetic microstructures from phase-field simulations with diffusivities from molecular dynamics simulations of ordered and disordered systems. Systematically designed simulations reveal an interdependence between atomistic and mesoscopic microstructural impacts on the effective ionic conductivity of polycrystalline LLZO, quantified by newly defined metrics that characterize the complex ionic transport mechanism. Our results provide fundamental understanding of the physical origins of the reported variability in ionic conductivities based on an extensive analysis of literature data, while simultaneously outlining practical design guidance for achieving desired ionic transport properties based on conditions for which sensitivity to microstructural features is highest. Additional implications of our results are discussed, including a possible connection between ion conduction behavior and dendrite formation.

25 ENERGY STORAGE↗
DOE OSTI2022

Plateau–Rayleigh instability with a grain boundary twist

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

Hussein, Omar↗
DOE OSTI2021

Condensation of collective polar vortex modes

The dynamics of extended objects such as domain walls, domain bubbles, vortex structures, etc., can be described by their equations of motion associated with their effective mass and spring constant. Here we analytically derive the equations of motion for the polarization dynamics and elastodynamics for the structural responses of ferroelectric polar vortices, and theoretically extract their effective mass, spring constant, and mode frequencies. We demonstrate two subterahertz phonon modes and predicted their frequencies, both consistent with our recent experimental measurements and phase-field simulations. We show that elastic modulation of the energy function and spring constants leads to a condensation of a collective mode upon a second-order structural transition from symmetric to asymmetric vortices at a critical strain, analogous to the ferroelectric soft phonon mode at a ferroelectric transition. Furthermore, the present work offers a theoretical framework for predicting and manipulating the ultrafast collective dynamics of polar nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Anomalous strain-energy-driven macroscale translation of grains during nonisothermal annealing

We report a mode of grain growth, involving the macroscopic translation of grain centers during nonisothermal annealing. Through synchrotron high-energy x-ray diffraction microscopy, we find dissolution of semicoherent precipitates generates dislocations, thereby raising the stored strain energy within grains. Here, the subsequent evolution of grains shows unexpected grain translations over length scales of 10-100 mu m. Phase-field simulations reveal that such translations are not uncommon in strain-energy-driven grain growth, wherein different regions of a grain may grow and shrink simultaneously.

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

Phase transition enhanced superior elasticity in freestanding single-crystalline multiferroic BiFeO 3 membranes

The integration of ferroic oxide thin films into advanced flexible electronics will bring multifunctionality beyond organic and metallic materials. However, it is challenging to achieve high flexibility in single-crystalline ferroic oxides that is considerable to organic or metallic materials. Here, we demonstrate the superior flexibility of freestanding single-crystalline BiFeO 3 membranes, which are typical multiferroic materials with multifunctionality. They can endure cyclic 180° folding and have good recoverability, with the maximum bending strain up to 5.42% during in situ bending under scanning electron microscopy, far beyond their bulk counterparts. Such superior elasticity mainly originates from reversible rhombohedral-tetragonal phase transition, as revealed by phase-field simulations. This study suggests a general fundamental mechanism for a variety of ferroic oxides to achieve high flexibility and to work as smart materials in flexible electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low-voltage magnetoelectric coupling in membrane heterostructures

Strain-mediated magnetoelectric (ME) coupling in ferroelectric (FE)/ferromagnetic (FM) heterostructures offers a unique opportunity for both fundamental scientific research and low-power multifunctional devices. Relaxor-FEs, such as (1-x)Pb(Mg 1/3 Nb 2/3 )O 3 -(x)PbTiO 3 (PMN-xPT), are ideal FE layer candidates because of their giant piezoelectricity. However, thin films of PMN-PT suffer from substrate clamping, which substantially reduces piezoelectric in-plane strains. Here, we demonstrate low-voltage ME coupling in an all-thin-film heterostructure that uses the anisotropic strains induced by the (011) orientation of PMN-PT. We completely remove PMN-PT films from their substrate and couple with FM Ni overlayers to create membrane PMN-PT/Ni heterostructures showing 90° Ni magnetization rotation with 3 V PMN-PT bias, much less than the bulk PMN-PT ~100-V requirement. Scanning transmission electron microscopy and phase-field simulations clarify the membrane response. These results provide a crucial step toward understanding the microstructural behavior of PMN-PT thin films for use in piezo-driven ME heterostructures.

36 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 OSTI2022

An Investigation of Variables Affecting Plutonium Hydriding (20210946DI): Surface Phase Transitions and Reactions

Tusas was originally developed for phase field simulation of solidification (see [33, 6, 32]) and evolution of polycrystalline solid mechanics (see [16]). Tusas is particularly developed for PDEs compatible with structured or unstructured Lagrange (nodal) finite element discretizations and explicit (Euler) or implicit (Euler, Trapezoid, BDF2) temporal discretizations.

97 MATHEMATICS AND COMPUTING↗
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↗