Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field”

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

Understanding coarsening of a post-corrosion microstructure in a molten salt by combining phase-field modeling and in situ tomography

Alloys corroding in molten salt have been observed to form bicontinuous, nanoporous microstructures via dealloying, which subsequently undergo coarsening due to facile transport in high-temperature conditions. In this work, we describe a methodology to elucidate the underlying transport mechanisms during coarsening of a bicontinuous microstructure via quantitative comparisons between phase-field simulations and four-dimensional in situ experiments, in this case X-ray nanotomography of the coarsening of a dealloyed 80 wt% Ni-20 wt% Cr microwire in molten KCl-MgCl 2 at 800°C. We conduct phase-field simulations initialized from experimental data to model coarsening via three different transport mechanisms: surface diffusion, solid bulk diffusion, and liquid bulk diffusion. These simulations reproduce key features of the experiment, such as the densification of the outer layer of the dealloyed wire and the reduction in radius over time. We quantitatively compare different microstructural characteristics between the simulations and experiment and extract temporal scaling factors that optimally match the time scales of the simulations to that of the experiment. This allows us to evaluate morphological similarity between the simulations and experiment and relate the experimental coarsening kinetics to fundamental material properties. We find that surface diffusion is most likely to be the dominant coarsening mechanism, and its kinetics imply a surface diffusivity of D S = 8.9 x 10 -20 m 3 /s, which is within the range of reported values for Ni-vacuum interfaces at 800°C. However, the difference between the experiment and the surface diffusion simulation increases substantially at late times, suggesting that other mechanisms, such as the dissolution of residual Cr, may be at play.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A phase field framework for capillary-induced fracture in unsaturated porous media: Drying-induced vs. hydraulic cracking

This work introduces a unified mathematical framework to replicate both desiccation-induced and hydraulic fracturing in low-permeable unsaturated porous materials observed in experiments. The unsaturated porous medium is considered as a three-phase solid–liquid–gas effective medium of which each constituent occupies a fraction of the representative elementary volume. As such, an energy-minimization-based phase-field model (PFM) is formulated along with the Biot’s poroelasticity theory to replicate the sub-critical crack growth in the brittle regime. Unlike hydraulic fracturing where the excess pore liquid pressure plays an important role at the onset and propagation of cracks, desiccation cracks are mainly driven by deformation induced by water retention. Therefore, the wettability of the solid skeleton may affect the evolution of the capillary pressure (suction) and change the path-dependent responses of the porous media. This air–water–solid interaction may either hinder or enhance the cracking occurrence. This difference of capillary effect on crack growth during wetting and drying is replicated by introducing retention-sensitive degradation mechanisms in our phase field fracture approach. To replicate the hydraulic behaviors of the pore space inside the host matrix and that of the cracks, the path-dependent changes of the intrinsic permeability due to crack growth and porosity changes are introduced to model the flow conduit in open and closed cracks. Numerical examples of drying-induced and hydraulic fracturing demonstrate the capability of the proposed model to capture different fracture patterns, which qualitatively agrees with the fracture mechanisms of related experiments documented in the literature.

42 ENGINEERING↗
DOE OSTI2022

The effects of dose, dose rate, and irradiation type and their equivalence on radiation-induced segregation in binary alloy systems via phase-field simulations

Radiation-induced segregation is a phenomenon commonly observed in many alloys which consists of the redistribution of elements (solute or interstitial impurities) under irradiation. The onset and development of radiation-induced segregation can only occur when a sufficient flux of defects is sustained and defect sinks are present. Irradiation dose, dose rate, and particle types all affect defect flux. In this work, we employ a phase-field model to examine the effects of dose, dose rate, and type of incident particles on radiation-induced segregation behavior in a model binary alloy. The phase-field model takes into account the formation and evolution of point defects as well as defect clusters, the diffusion and clustering of alloy species, the presence of additional extrinsic defect sinks in the form of dislocations, and two different methods of radiation-damage insertion, which are intended to simulate either light-ion/electron irradiation via Frenkel pairs or heavy-ion irradiation in the form of cascades. Our results show a dose-rate and particle-type dependence on the amount of solute segregation. We show that the material systems exposed to higher dose rates are less subjected to solute segregation at equivalent doses. We also show that such dose-rate-dependence behavior is due to a delay of the incubation dose at which radiation-induced segregation effectively starts. Particle type and the presence of dislocations can accentuate this behavior. Our model predictions correlate with many experimental observations made over the years on radiation-induced segregation providing credence to the simulation results. The methodology presented in this study allows for a first-order prediction of the dose rate at which proxy irradiation experiments could be performed to approximate radiation-induced segregation behaviors seen in targeted irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Recent progress on the mesoscale modeling of architected thin-films via phase-field formulations of physical vapor deposition

Thin-film coatings can be found everywhere in modern technological applications due to desirable electrical, mechanical, chemical, and optical properties. These properties directly depend upon the thin-film’s microstructural features, which are themselves influenced by the materials and vapor-deposition processing conditions used for fabrication. As such, understanding processing-microstructure relationships is essential to designing thin-films with optimized properties, and discovering new processing conditions that allow for novel thin-films with multifunctional microstructures. Here, a short review is presented on recent developments that utilize the phase-field method to simultaneously model the vapor-deposition process and corresponding microstructure formation at the mesoscale. Also phase-field-based vapor-deposition models that simulate thin-film growth of immiscible alloy and polycrystalline systems are highlighted in addition to machine-learning-based surrogate models that can facilitate accelerated high-fidelity simulations along with materials design and exploration studies.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A mixed-mode phase field fracture model in anisotropic rocks with consistent kinematics

Under a pure tensile loading, cracks in brittle, isotropic, and homogeneous materials often propagate such that pure mode I kinematics are maintained at the crack tip. However, experiments performed on geo-materials, such as sedimentary rock, shale, mudstone, concrete and gypsum, often lead to the conclusion that the mode I and mode II critical fracture energies/surface energy release rates are distinctive. This distinction has great influence on the formation and propagation of wing cracks and secondary cracks from pre-existing flaws under a combination of shear and tensile or shear and compressive loadings. To capture the mixed-mode fracture propagation, a mixed-mode I/II fracture model that employs multiple critical energy release rates based on Shen and Stephansson, IJRMMS, 1993 is reformulated in a regularized phase field fracture framework. Here we obtain the mixed-mode driving force of the damage phase field by balancing the microforce. Meanwhile, the crack propagation direction and the corresponding kinematics modes are determined via a local fracture dissipation maximization problem. Several numerical examples that demonstrate mode II and mixed-mode crack propagation in brittle materials are presented. Possible extensions of the model capturing degradation related to shear/compressive damage, as commonly observed in sub-surface applications and triaxial compression tests, are also discussed.

42 ENGINEERING↗
DOE OSTI2019

Phase-field model for diffusion-induced grain boundary migration: An application to battery electrodes

Diffusion-induced grain boundary migration is a phenomenon in which a grain boundary moves in response to the driving forces generated by diffusing solute species. For example, diffusing solute species change the atomic volume in the host material, either by filling a vacancy with a misfitting solute atom or by expanding host lattices through interstitial diffusion. These volume changes are inhomogeneous and are stored as elastic energy in the material that drives grain boundaries. In this paper, we introduce our previously developed Cahn-Hilliard–phase-field-crystal model (CH-PFC) as a computational tool to investigate diffusion-induced grain boundary migration in crystalline materials. This multiscale phase-field model couples the composition field of a diffusing species with the crystallographic texture of a host material. We apply the CH-PFC model to battery electrodes and investigate whether interstitial solute diffusion induces grain growth in FePO 4 /LiFePO 4 electrodes. To this end, we compute grain growth in 60 FePO 4 electrodes by conducting two parallel trials: In the first trial, we cycle the electrode and calculate diffusion-induced grain growth. In the second trial, we do not cycle the electrode and calculate curvature-driven grain growth. We find a statistically significant grain growth in the cycled electrodes and negligible grain growth in the noncycled electrodes. Altogether, we show that the CH-PFC model not only predicts electrode microstructures as a function of the Li composition, but also predicts the crystallographic features of an electrode during battery operation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the Effect of Electrochemical Properties and Microstructure on the Microgalvanic Corrosion of Mg Alloys via Phase-Field Simulations

In this work, we apply a phase-field model to gain insights into the corrosion behavior of Mg alloys. Specifically, we study the effect of electrochemical properties of the materials and their environment and the spatial distribution of second phases in a magnesium alloy. To this end, we perform sensitivity analyses in which we separately vary the exchange current density and corrosion potential of the second phase, as well as the electrolyte conductivity. We describe the effects of these parameters on the corrosion current and the resulting corrosion-front morphology. We also confirm that the ratio between the Wagner length, calculated for the anodic phase, and the initial width of the anodic domain is a good predictor of the morphological characteristic of the corroding surface. In addition, we examine the effect of the microstructure on the corrosion behavior in both 2D and 3D systems by varying the shape, size, and distribution of the second-phase regions. Our findings provide insights into the influence of the two-phase microstructure on the corrosion behavior in a magnesium alloy.

Electrochemistry↗
DOE OSTI2020

Structure and Dynamics of Domains in Ferroelectric Nanostructures – Phase-Field Modeling (Final Technical Report)

Domain pattern formation is one of the most common phenomena in nature and is a topic of immense interest in many fields ranging from materials science and physics to chemistry and biology. This DOE sponsored research project explored the basic science concerning the thermodynamic stability of mesoscale polarization domain patterns and their temporal evolution mechanisms during formation and subsequent switching in ferroelectric nanostructures and heterostructures. The project employed the computational phase-field method in combination of microelasticity and electrostatic theories. The research was carried out in close collaboration with a number of experimental groups who used High Resolution Transmission Electron Microscopy (HRTEM), in situ TEM with Scanning Probe Microscopy (SPM), or Piezoresponse Force Microscopy (PFM) to characterize the domain structures and dynamics of ferroelectric thin films and heterostructures and who grow high-quality ferroelectric and multiferroic thin films using advanced growth techniques such as Molecular Beam Epitaxy (MBE), Pulsed Laser Deposition (PLLD), and sputtering. The research efforts of the project help establish the phase-field method as the most powerful method for understanding and predicting domain structures in ferroelectric thin films and nanostructures. The findings of the project led to the basic understanding of stability and switching mechanisms of ferroelectric domains under different mechanical boundary conditions and under either homogeneous capacitor configurations or local fields using metallic probes as electrodes and the emergence of charged domain walls during domain switching. The project predicted the spatial length scales, temperature ranges, and electromechanical conditions for different polar states in heterostructures and guided the discovery of both transient and stable novel polarization states containing vortex lattices in oxide superlattices. The project resulted in 134 journal publications and 6 PhD theses with all the PhD graduates currently working in either academia or industry within the United States. The basic understanding on the stability of mesoscale polar states and pattern evolution achieved by the project improved our ability to control and engineer properties of ferroelectric thin films and heterostructures for potential applications in nanoscale electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electrochemical grand potential-based phase-field simulation of electric field-assisted sintering

Here, an electrochemical grand potential functional was proposed to describe the sintering of an ionic ceramic green body. The resultant phase-field description enables simulation of the consolidation of an arbitrary number of granular particles and their interactions with the surrounding void phase. The model includes the effects of charged vacancies and the associated interactions between internal and applied electric fields. Defect segregation to grain boundaries is also accounted for, as well as enhanced interfacial defect mobilities. The model was parameterized for Y 2 O 3 . Simulations of two-particle systems showed that the applied electric field had an increasingly important impact on neck growth as particle size increased. A sudden rapid increase in temperature occurred for larger field strengths, which has been reported to be correlated to the onset of a flash event in flash sintering. Simulations of many particles showed that internal heat generation by Joule heating was localized at particle–particle contacts (grain boundaries), even though their conductivities were lower than nearby internal particle-void interfaces. A percolative path for ionic charge across the green body and the ceramic sintered solid was thus defined, accelerating the Joule heating process as the porosity of the green body is removed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A combined kinetic Monte Carlo and phase field approach to model thermally activated dislocation motion

In this work, a kinetic Monte Carlo (KMC) scheme is integrated into a phase field dislocation dynamics (PFDD) model to account for thermal activation. While the implementation is general and could be applied to many different thermally activated processes, we apply the model to simulate a dislocation transmitting through an interface. Two cases are analyzed: a dislocation transmitting through a strong interface, and a dislocation transmitting through a more realistic grain boundary in tungsten. In conclusion, results and trends for both cases are discussed in detail, and ultimately show that the combined KMC-PFDD algorithm is effective in capturing thermally activated dislocation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2017

Phase-field based Multiscale Modeling of Heterogeneous Solid Electrolytes: Applications to Nanoporous Li 3 PS 4

Modeling the effective ion conductivities of heterogeneous solid electrolytes typically involves the use of a computer-generated microstructure consisting of randomly or uniformly oriented fillers in a matrix. But, the structural features of the filler/matrix interface, which critically determine the interface ion conductivity and the microstructure morphology, have not been considered during the microstructure generation. In using nanoporous β-Li 3 PS 4 electrolyte as an example, we develop a phase-field model that enables generating nanoporous microstructures of different porosities and connectivity patterns based on the depth and the energy of the surface (pore/electrolyte interface), both of which are predicted through density functional theory (DFT) calculations. Room-temperature effective ion conductivities of the generated microstructures are then calculated numerically, using DFT-estimated surface Li-ion conductivity (3.14×10 -3 S/cm) and experimentally measured bulk Li-ion conductivity (8.93×10 -7 S/cm) of β-Li 3 PS 4 as the inputs. We also use the generated microstructures to inform effective medium theories to rapidly predict the effective ion conductivity via analytical calculations. Furthemore, when porosity approaches the percolation threshold, both the numerical and analytical methods predict a significantly enhanced Li-ion conductivity (1.74×10 -4 S/cm) that is in good agreement with experimental data (1.64×10 -4 S/cm). The present phase-field based multiscale model is generally applicable to predict both the microstructure patterns and the effective properties of heterogeneous solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Machine learning materials physics: Surrogate optimization and multi-fidelity algorithms predict precipitate morphology in an alternative to phase field dynamics

Machine learning has been effective at detecting patterns and predicting the response of systems that behave free of natural laws. Examples include learning crowd dynamics, recommender systems and autonomous mobility. There also have been applications to the search for new materials that draw upon big-data classification problems. Though, when it comes to physical systems governed by conservation laws, the role of machine learning has been more limited. Here, we present our recent work in exploring the role of machine learning methods in discovering, or aiding, the search for physics. Specifically, we focus on using machine learning algorithms to represent high-dimensional free energy surfaces with the goal of identifying precipitate morphologies in alloy systems. Traditionally, this problem has been approached by combining phase field models, which impose first-order dynamics, with elasticity, to traverse a free energy landscape in search of minima. Equilibrium precipitate morphologies occur at these minima. Here, we exploit the machine learning methods to represent high-dimensional data, combined with surrogate optimization, sensitivity analysis and multifidelity modeling as an alternate framework to explore phenomena controlled by energy extremization. This combination of data-driven methods offers an alternative to the imposition of first-order dynamics via phase field methods, and represents one aspect of applying machine learning to materials physics.

42 ENGINEERING↗
DOE OSTI2022

Phase-field dislocation modeling of cross-slip

The phase-field dislocation dynamics (PFDD) model is extended to simulate cross-slip in body-centered cubic materials. The model makes use of a non-orthogonal, body-centered cubic numerical grid such that all calculable points lie on dislocation glide planes. We demonstrate that the advanced PFDD model predicts a low-energy, non-planar core for screw dislocations and a planar core for edge dislocations. These core differences result in screw-edge character dependence in the critical stress to initiate glide and dislocation drag coefficients. Finally, we show that this model can be used to study the propensity of screw dislocations to cross-slip around or transmit through a coherent, crystallographic obstacle, depending on obstacle strength.

36 MATERIALS SCIENCE↗