Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field modeling”

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

Optimal design of chemoepitaxial guideposts for the directed self-assembly of block copolymer systems using an inexact Newton algorithm

Directed self-assembly (DSA) of block copolymers (BCPs) is one of the most promising developments in the cost-effective production of nanoscale devices. The process makes use of the natural tendency for BCP melts to form nanoscale structures upon phase separation. The phase separation can be directed through the use of chemically patterned substrates to promote the formation of morphologies that are essential to the production of semiconductor devices. Moreover, the design of substrate pattern can be formulated as an optimization problem for which we seek optimal substrate designs that effectively produce given target morphologies. In this paper, we adopt a phase field model given by a nonlocal Cahn–Hilliard partial differential equation (PDE) based on the minimization of the Ohta–Kawasaki free energy, and present an efficient PDE-constrained optimization framework for the optimal design problem. The design variables are the locations of circular- or strip-shaped guiding posts that are used to model the substrate chemical pattern. To solve the ensuing optimization problem, we propose a variant of an inexact Newton conjugate gradient algorithm tailored to this problem. Additionally, we demonstrate the effectiveness of our computational strategy on numerical examples that span a range of target morphologies. Owing to our second-order optimizer and fast state solver, the numerical results demonstrate five orders of magnitude reduction in computational cost over previous work. The efficiency of our framework and the fast convergence of our optimization algorithm enable us to rapidly solve the optimal design problem in not only two, but also three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Disordering of helium gas bubble superlattices in molybdenum under ion irradiation and thermal annealing

Self-organization of gas bubbles causes the formation of an ordered array of nanoscale gas bubbles (a gas bubble superlattice), a highly efficient mechanism for gas storage under irradiation. Here, the stability of helium (He) gas bubble superlattices in molybdenum (Mo) under krypton (Kr) ion irradiation and thermal annealing has been investigated. The He gas bubble superlattices gradually become disordered under Kr ion irradiation at 300° C, and the order-disorder transformation process completes at 2.5 dpa. Both transmission electron microscopy (TEM) and synchrotron-based small-angle X-ray scattering (SAXS) reveal that the order-disorder transformation of He gas bubble superlattices is associated with a slight increase in the average bubble size. Phase-field modeling indicates that the inhomogeneous growth and coarsening of bubbles/voids cause the disordering of imperfect superlattices under irradiation and implies that highly ordered superlattices could potentially exhibit much stronger resistance to irradiation damage. Under thermal annealing, the He gas bubble superlattices in Mo become unstable and disordered at 1000° C with the bubble size increasing from ~1.1 to ~1.6 nm. The finding in this research provides insights into the disordering mechanisms of defect superlattices as well as guidance for designing stable defect superlattices in harsh environments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Kinetics and mechanism of light-induced phase separation in a mixed-halide perovskite

Halide ion phase separation is a barrier to the application of mixed-halide perovskites, whereby the presence of large populations of photo-generated or injected carriers causes undesirable changes in the local band gap. Here, we investigate the mechanism of phase separation in CsPbI x Br 3-x perovskite single crystals driven by light. The phase separation process and its dynamics are visualized at the nanometer scale at cryogenic temperatures using in-situ scanning transmission electron microscopy and cathodoluminescence. Combined with phase field modeling that accounts for the coupling between electronic and chemical driving forces, our observations point to a spinodal decomposition mechanism in which both the amplitude of composition fluctuation and the characteristic length scale grow non-linearly with time. Our findings provide microscopic insights that can assist in further engineering mixed-halide perovskites either for stability or for intentional programming of the local halide ion composition, opening pathways to a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fragmentation of the edge of a terminated Cu nanolayer within a Nb matrix upon annealing

Here, we investigate the evolution of the tip of a terminated Cu nano-layer embedded within a Nb matrix upon annealing at 800°C. The tip of the terminated Cu layer fragments progressively into a series of single-crystal Cu particles upon continued annealing. This finding shows that layer fragmentation is the second step, after layer pinchoff, in the thermal coarsening of nanolayered metal composites. Comparison to prior phase-field modeling permits us to infer that the mobility at 800°C of Cu-Nb, as modeled by the Cahn-Hillard equation, is M(c)=(M 0 = 7.3 ± 3.6 nm 5 / s • eV) |1 – c 2 |, where the order parameter, c, is 1 for Cu and -1 for Nb.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thin-Film Paradigm to Probe Interfacial Diffusion during Solid-State Metathesis Reactions

We highlight a paradigm for studying complex reaction mechanisms that guide the synthesis of materials. Thin films of FeCl 2 and Na 2 S 2 were deposited to study the metathesis reaction to form FeS 2 and NaCl. In situ X-ray reflectivity was used to monitor the interface between materials, which revealed a slow, impeded reaction at high temperatures as compared to previous studies using powder samples. AC impedance and X-ray photoelectron spectroscopy provided insight into distribution of elements and conductivity of the phases present during the reaction, and phase-field modeling was used to elucidate the diffusion of ions throughout the thin-film bilayers. Here, the use of thin-film bilayers provides a simplified system to study solid-state metathesis reactions and highlights the complexity of diffusion at solid-state interfaces.

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

Probing Metastable Domain Dynamics via Automated Experimentation in Piezoresponse Force Microscopy

The dynamics of complex topological defects in ferroelectric materials is explored using automated experimentation in piezoresponse force microscopy. Specifically, a complex trigger system (i.e., “FerroBot”) is employed to study metastable domain-wall dynamics in Pb 0.6 Sr 0.4 TiO 3 thin films. Several regimes of superdomain wall dynamics have been identified, including smooth domain-wall motion and significant reconfiguration of the domain structures. We have further demonstrated that microscopic mechanisms of the domain-wall dynamics can be identified; i.e., domain-wall bending can be separated from irreversible domain reconfiguration regimes. In conjunction, phase-field modeling was used to corroborate the observed mechanisms. As such, the observed superdomain dynamics can provide a model system for classical ferroelectric dynamics, much like how colloidal crystals provide a model system for atomic and molecular systems.

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

Terahertz-field activation of polar skyrons

Unraveling collective modes arising from coupled degrees of freedom is crucial for understanding complex interactions in solids and developing new functionalities. Unique collective behaviors emerge when two degrees of freedom, ordered on distinct length scales, interact. Polar skyrmions, three-dimensional electric polarization textures in ferroelectric superlattices, disrupt the lattice continuity at the nanometer scale with nontrivial topology, leading to previously unexplored collective modes. Here, using terahertz-field excitation and femtosecond x-ray diffraction, we discover subterahertz collective modes, dubbed “skyrons”, which appear as swirling patterns of atomic displacements functioning as atomic-scale gearsets. The key to activating skyrons is the use of the THz field that couples primarily to skyrmion domain walls. Momentum-resolved time-domain measurements of diffuse scattering reveal an avoided crossing in the dispersion relation of skyrons. Atomistic simulations and dynamical phase-field modeling provide microscopic insights into the three-dimensional crystallographic and polarization dynamics. The amplitude and dispersion of skyrons are demonstrated to be controlled by sample temperature and electric-field bias. The discovery of skyrons and their coupling with terahertz fields opens avenues for ultrafast control of topological polar structures.

ferroelectrics↗
DOE OSTI2020

High-throughput phase-field simulations and machine learning of resistive switching in resistive random-access memory

Metal oxide-based Resistive Random-Access Memory (RRAM) exhibits multiple resistance states, arising from the activation/deactivation of a conductive filament (CF) inside a switching layer. Understanding CF formation kinetics is critical to achieving optimal functionality of RRAM. Here a phase-field model is developed, based on materials properties determined by ab initio calculations, to investigate the role of electrical bias, heat transport and defect-induced Vegard strain in the resistive switching behavior, using MO 2- x systems such as HfO 2- x as a prototypical model system. It successfully captures the CF formation and resultant bipolar resistive switching characteristics. High-throughput simulations are performed for RRAMs with different material parameters to establish a dataset, based on which a compressed-sensing machine learning is conducted to derive interpretable analytical models for device performance (current on/off ratio and switching time) metrics in terms of key material parameters (electrical and thermal conductivities, Vegard strain coefficients). These analytical models reveal that optimal performance (i.e., high current on/off ratio and low switching time) can be achieved in materials with a low Lorenz number, a fundamental material constant. This work provides a fundamental understanding to the resistive switching in RRAM and demonstrates a computational data-driven methodology of materials selection for improved RRAM performance, which can also be applied to other electro-thermo-mechanical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Grain boundary effects in high-temperature liquid-metal dealloying: a multi-phase field study

Abstract A multi-phase field model is employed to study the microstructural evolution of an alloy undergoing liquid dealloying, specifically considering the role of grain boundaries. A semi-implicit time-stepping algorithm using spectral methods is implemented, which enables simulating large 2D and 3D domains over long time scales while still maintaining a realistic interfacial thickness. Simulations reveal a mechanism of coupled grain–boundary migration to maintain equilibrium contact angles with the topologically complex solid–liquid interface, which locally accelerates diffusion-coupled growth of a liquid channel into the precursor. This mechanism asymmetrically disrupts the ligament connectivity of the dealloyed structure in qualitative agreement with published experimental observations. The grain boundary migration-assisted corrosion channels form even for precursors with small amounts of the dissolving alloy species, below the parting limit . The activation of this grain boundary dealloying mechanism depends strongly on grain boundary mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-equilibrium pathways to emergent polar supertextures

Ultrafast stimuli can stabilize metastable states of matter inaccessible by equilibrium means. Establishing the spatiotemporal link between ultrafast excitation and metastability is crucial to understand these phenomena. Here, in this work, we utilize single-shot optical pump–X-ray probe measurements to capture snapshots of the emergence of a persistent polar vortex supercrystal in a heterostructure that hosts a fine balance between built-in electrostatic and elastic frustrations by design. By perturbing this balance with photoinduced charges, an initially heterogeneous mixture of polar phase disorders within a few picoseconds, leading to a state composed of disordered ferroelectric and suppressed vortex orders. On the picosecond–nanosecond timescales, transient labyrinthine fluctuations develop, accompanied by the recovery of the vortex order. On longer timescales, these fluctuations are progressively quenched by dynamical strain modulations, which drive the collective emergence of a single vortex supercrystal phase. Our results, corroborated by dynamical phase-field modelling, reveal non-equilibrium pathways following the ultrafast excitation of designer systems to persistent metastability.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Subterahertz collective dynamics of polar vortices

The collective dynamics of topological structures are of interest from both fundamental and applied perspectives. For example, studies of dynamical properties of magnetic vortices and skyrmions have not only deepened our understanding of many-body physics but also offered potential applications in data processing and storage. Topological structures constructed from electrical polarization, rather than electron spin, have recently been realized in ferroelectric superlattices, and these are promising for ultrafast electric-field control of topological orders. However, little is known about the dynamics underlying the functionality of such complex extended nanostructures. In this study, using terahertz-field excitation and femtosecond X-ray diffraction measurements, we observe ultrafast collective polarization dynamics that are unique to polar vortices, with orders-of-magnitude higher frequencies and smaller lateral size than those of experimentally realized magnetic vortices. A previously unseen tunable mode, hereafter referred to as a vortexon, emerges in the form of transient arrays of nanoscale circular patterns of atomic displacements, which reverse their vorticity on picosecond timescales. Its frequency is considerably reduced (softened) at a critical strain, indicating a condensation (freezing) of structural dynamics. We use first-principles-based atomistic calculations and phase-field modelling to reveal the microscopic atomic arrangements and corroborate the frequencies of the vortex modes. The discovery of subterahertz collective dynamics in polar vortices opens opportunities for electric-field-driven data processing in topological structures with ultrahigh speed and density.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗