Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Field Model”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Morphology selection in dealloying: A phase field study of the coupling among kinetic mechanisms

A phase field model is used to investigate how the corrosion rate and morphology of a dealloying metal is controlled by the kinetic rates of bulk solid diffusion, interfacial diffusivity, and interface dissolution. A specific focus is investigating the manner by which variations in these kinetic parameters synergistically influence the onset of a morphological instability at the solid-liquid interface, which is known to precede rapid dealloying corrosion towards highly topologically complex structures such as bicontinuous porous solids. This work demonstrates that this instability can be suppressed by either an increase of interdiffusion in the solid alloy, or by a decreasing rate of dissolution at the solid-liquid interface. When the instability is suppressed, the result is a relatively planar, passivated dealloying front. Importantly, however, coupled changes in kinetic parameters can also promote the morphological instability and sustain fine-scale ligament and pore formation. Finally, the relevance of the solid-liquid interfacial diffusion for determining the rate of dealloying and ligament formation is highlighted. By elucidating the mechanistic influence of bulk and interfacial kinetic factors on the dealloyed morphology, this work helps to rationalize the contrasting morphologies observed in various dealloying systems.

Diffusion↗
DOE OSTI2020

Transformation-induced fracture toughening in CuAlBe shape memory alloys: A phase-field study

Stress-induced transformation ahead of a crack tip can relax the stress concentration and produce fracture toughness in shape memory alloys (SMAs). In this manuscript, we utilize a nonisothermal phase-field model (PFM) to study the martensitic transformation induced crack tip toughening in CuAlBe SMA. The force-displacement curve, transformation zone and high stress zone in single crystalline samples show high dependency on the grain orientation with respect to the crack alignment. Comparison between isothermal and non-isothermal simulations reveals that the transformation-induced self-heating decreases the toughening capability by increasing the critical transformation stress. Investigating the high stress zones in CuAlBe and the nontransforming CuAl alloy shows that the toughening is obtained by redistributing the stress concentration far from the crack tip, as the high stress zone follows the tip of the transformation zone. Increasing the acuity of the crack tip is found to generate more symmetric martensite wings on both side of the crack axis and a more localized high stress zone. The polycrystal specimen displays higher toughening due to the internal constraints related to the presence of various grains with difference orientations. Depending on the orientation of the grain inclosing the crack tip, the toughening can be lower or higher. Finally, coincidence of the crack tip with a triple junction is found to improve the toughening behavior.

42 ENGINEERING↗
DOE OSTI2020

Domain Dynamics under Ultrafast Electric-Field Pulses

Exploring the dynamic responses of a material is of importance to both understanding its fundamental physics at high frequencies and potential device applications. Here we develop a phase-field model for predicting the dynamics of ferroelectric materials and study the dynamic responses of ferroelectric domains and domain walls subjected to an ultrafast electric field pulse. We discover a transition of domain evolution mechanisms from pure domain growth at a relatively low field to combined nucleation and growth of domains at a high field. We derive analytical models for the two regimes which allow us to extract the effective mass and damping coefficient of ferroelectric domain walls. The exhibition of two regimes for the ferroelectric domain dynamics at low and high electric fields is expected to be a general phenomenon that would appear for ferroic domains under other ultrafast stimuli. Finally, the present work also offers a general framework for studying domain dynamics and obtaining fundamental properties of domain walls and thus for manipulating the dynamic functionalities of ferroelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Intermittent Defect Fluctuations in Oxide Heterostructures

Abstract The heterogeneous nature, local presence, and dynamic evolution of defects typically govern the ionic and electronic properties of a wide variety of functional materials. While the last 50 years have seen considerable efforts into development of new methods to identify the nature of defects in complex materials, such as the perovskite oxides, very little is known about defect dynamics and their influence on the functionality of a material. Here, the discovery of the intermittent behavior of point defects (oxygen vacancies) in oxide heterostructures employing X‐ray photon correlation spectroscopy is reported. Local fluctuations between two ordered phases in strained SrCoO x with different degrees of stability of the oxygen vacancies are observed. Ab‐initio‐informed phase‐field modeling reveals that fluctuations between the competing ordered phases are modulated by the oxygen ion/vacancy interaction energy and epitaxial strain. The results demonstrate how defect dynamics, evidenced by measurement and modeling of their temporal fluctuations, give rise to stochastic properties that now can be fully characterized using coherent X‐rays, coupled for the first time to multiscale modeling in functional complex oxide heterostructures. The study and its findings open new avenues for engineering the dynamical response of functional materials used in neuromorphic and electrochemical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2017

An Atomistic Introduction to Orientation Relations Between Phases in the Face-centered Cubic to Body-centered Cubic Phase Transition in Iron and Steel

We establish an atomistic view of the high- and low-temperature phases of iron/steel as well as some elements of the phase transition between these phases on cooling. In particular we examine the 4 most common orientation relationships between the high temperature austenite and low-temperature ferrite phases seen in experiment. With a thorough understanding of these relationships we are prepared to set up various atomistic simulations, using techniques such as Density Functional Theory and Molecular Dynamics, to further study the phase transition, in particular, quantities needed for Phase Field Modeling, such as the free energies of bulk phases and the phase transition front propagation velocity.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Alpha – omega and omega – alpha phase transformations in zirconium under hydrostatic pressure: A 3D mesoscale study

Here, a three dimensional (3D) elastoplastic phase-field model is developed for modeling the hydrostatic pressure-induced alpha – omega phase transformation and the reverse phase transformation, i.e. omega – alpha, in zirconium (Zr). Plastic deformation and strain hardening of the material are also considered in the model. The microstructure evolution during both phase transformations is studied. The transformation start pressures at different temperatures are predicted and are plotted as a phase diagram. The effect of phase transformations on the mechanical properties of the material is also studied. The input data corresponding to pure Zr are acquired from experimental studies as well as by using the CALPHAD method. Our simulations show that three different omega variants form as laths. On release of pressure, reverse phase transformation initiates at lath boundaries. We observe that both phase transformations are martensitic in nature and also occur at the same pressure, i.e. little hysteresis. The transformation start pressures and the kinetics of the transformation predicted by our model are in good agreement with experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Multi-scale modeling of inter-granular fracture in UO2

A hierarchical multi-scale approach is pursued in this work to investigate the influence of porosity, pore and grain size on the intergranular brittle fracture in UO2. In this approach, molecular dynamics simulations are performed to obtain the fracture properties for different grain boundary types. A phase-field model is then utilized to perform intergranular fracture simulations of representative microstructures with different porosities, pore and grain sizes. In these simulations the grain boundary fracture properties obtained from molecular dynamics simulations are used. The responses from the phase-field fracture simulations are then fitted with a stress-based brittle fracture model usable at the engineering scale. This approach encapsulates three different length and time scales, and allows the development of microstructurally informed engineering scale model from properties evaluated at the atomistic scale.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Crystal plasticity including a phase-field deformation twinning model for the high-rate deformation of cyclotetramethylene tetranitramine

A new finite strain thermomechanical model for the high-rate deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) has been developed and applied to simulations of plate impact experiments. The crystal plasticity model is based on a model developed previously for RDX (Luscher et al. 2017), which is extended to incorporate deformation twinning. Twinning during normal plate impacts is simulated with a phase-field twin model. First, material parameters governing the kinetics of dislocation slip are calibrated on the subset of simulations which had negative Schmid factors for the twin system. Second, a parametric study of the twin material parameters was performed to find suitable values. Overall, the results of the simulations with the phase-field twinning model are reported for impacts on several crystal orientations. We find that the twin growth decreases with increasing distance from the impact surface because of dissipation of the shock front via dislocation-mediated plasticity, and that the simulated interface velocity with and without phase-field twinning do not show appreciable differences. These modeling results suggest that the significance of twinning in β-HMX cannot be determined with traditional loading configurations and diagnostics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mesoscale Modeling of Kinetic Phase Behaviors in Mg-B-H (Subcontract Report)

Storage of hydrogen on board vehicles is one of the critical enabling technologies for creating hydrogenfueled transportation systems that can reduce oil dependency and mitigate the long-term effects of fossil fuels on climate change. Stakeholders in developing hydrogen infrastructure are currently focused on highpressure storage at 350 bar and 700 bar, in part because no viable solid-phase storage material has emerged. Nevertheless, solid-state materials, including high-density hydrides, remain of interest because of their unique potential to meet all DOE targets and deliver hydrogen at lower pressures and higher on-board densities. A successful solution would significantly reduce costs and ensure the economic viability of a U.S. hydrogen infrastructure. The Mg(BH 4 ) 2 -MgB 2 system represents a highly promising solution because of its reasonable reaction enthalpy, high intrinsic capacity, and demonstrated reversibility, yet suffers from poor reaction kinetics. This subcontract aims to deliver a phase-field model for the kinetics of the evolution of the relevant phases within the Mg-B-H system during hydrogenation and dehydrogenation. This model will be used within a broader theory, synthesis, and characterization framework to study the properties of geometry-selected nanoparticles of pristine and doped MgB 2 /Mg(BH 4 ) 2 with two aims: (1) understand the intrinsic limitations in (de)hydrogenation; (2) devise strategies for improving thermodynamics and kinetics through nanostructuring.

08 HYDROGEN↗
DOE OSTI2018

Influence of surface nano-patterning on the placement of InAs quantum dots

We have examined the influence of spontaneous nano-patterning on the placement of InAs quantum dots (QDs) on (Al)GaAs surfaces using an experimental-computational approach. Both atomically flat and mounded surfaces, generated via a surface instability induced by the Ehrlich- Schwoebel barrier, are employed as templates for the subsequent deposition of InAs QDs. Using height profiles from atomic-force micrographs, we simulate QD deposition with a 2D phase field model, which describes the time evolution of the InAs layer driven by a chemical potential gradient. For flat surfaces, phase-field simulations result in QD densities comparable to experimental observations. For mounded surfaces, the simulations reveal QDs preferentially positioned in regions of positive curvature (substrate valleys), e.g., at the edge of surface mounds, consistent with the anisotropic QD placement observed experimentally. We discuss the role of curvaturedriven diffusion in the spontaneous ordering of QDs, demonstrating the applicability of this mechanism to AlGaAs mounds.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Hybrid models for the simulation of microstructural evolution influenced by coupled, multiple physical processes

Most materials microstructural evolution processes progress with multiple processes occurring simultaneously. In this work, we have concentrated on the processes that are active in nuclear materials, in particular, nuclear fuels. These processes are coarsening, nucleation, differential diffusion, phase transformation, radiation-induced defect formation and swelling, often with temperature gradients present. All these couple and contribute to evolution that is unique to nuclear fuels and materials. Hybrid model that combines elements from the Potts Monte Carlo, phase-field models and others have been developed to address these multiple physical processes. These models are described and applied to several processes in this report. An important feature of the models developed are that they are coded as applications within SPPARKS, a Sandiadeveloped framework for simulation at the mesoscale of microstructural evolution processes by kinetic Monte Carlo methods. This makes these codes readily accessible and adaptable for future applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Diffuse interface approach to modeling crystal plasticity with accommodation of grain boundary sliding

In diffuse-interface (phase field) models of polycrystalline solids, the grain boundaries naturally possess a narrow region of finite thickness. A Ginzburg-Landau type kinetic equation for plastic deformation of a solid is derived that combines crystal plasticity and J 2 plasticity. Within such a framework for a polycrystal the grain interior can be defined by crystal plasticity, and the grain boundary region can be assigned J 2 type plasticity. The two are connected smoothly. This enables the modeling framework to accommodate grain boundary sliding (GBS), an important deformation mechanism for creep. The relaxation of elastic modulus of a polycrystalline solid and the power-law creep (both accommodated by GBS) are reported on by 2D and 3D simulations, respectively, and the results are compared to the 2D finite-element simulations of Ghahremani and Crossman-Ashby. A strong grain shape dependence and orientation dependence (for non-equiaxed grains) of the GBS effect are predicted by this model.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Modeling the Restraint of Liquid Jets by Surface Tension in Microgravity

An axisymmetric phase field model is developed and used to model surface tension forces on liquid jets in microgravity. The previous work in this area is reviewed and a baseline drop tower experiment selected 'for model comparison. A mathematical model is developed which includes a free surface. a symmetric centerline and wall boundaries with given contact angles. The model is solved numerically with a compact fourth order stencil on a equally spaced axisymmetric grid. After grid convergence studies, a grid is selected and all drop tower tests modeled. Agreement was assessed by comparing predicted and measured free surface rise. Trend wise agreement is good but agreement in magnitude is only fair. Suspected sources of disagreement are suspected to be lack of a turbulence model and the existence of slosh baffles in the experiment which were not included in the model.

Chato, David J.↗
DOE OSTI2019

Structural deformation and transformation of θ'-Al 2 Cu precipitate in Al matrix via interfacial diffusion

Metastable θ'-Al 2 Cu precipitates enhance mechanical properties of Al alloy, but gradually transform to θ phase due to diffusion related coarsening at higher temperatures (>200°C). To improve understanding of the phase instability and transformation of θ'-Al 2 Cu, we investigated the interfacial atomic mobility of θ'-Al2Cu precipitate in Al matrix via molecular dynamics (MD). To characterize the interfacial atomic mobility, diffusion properties (activation energy, diffusivity, and jump attempt frequency) with respect to atomic species, interfacial structure, and temperature are calculated using atomic trajectories from MD. For the enhanced accuracy of this analysis, especially at low temperatures, temporal scale of molecular dynamics is significantly extended by employing the parallel replica dynamics. This study, as the first MD investigation on atomic diffusion of θ'-Al 2 Cu in Al matrix, suggests and discusses (i) structural deformation activated at the semi-coherent θ'-Al 2 Cu/Al interfaces, (ii) transformation initiated at the edge of the coherent θ'-Al 2 Cu/Al interface, and (iii) directional dependence of diffusion properties. Overall, these findings are expected to be employed in a larger scale phase-field modeling of precipitate growth, and to contribute to the understanding of phase transformation process and the development of alloys with θ'-phase retained to higher temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2016

Exploring Polarization Rotation Phase Instabilities in Super-Tetragonal BiFeO3 Epitaxial Thin Films and Their Technological Implications

Many functional properties of ferroelectrics are underlain by structural instabilities, which render these materials very susceptible to small alternating applied fields (electric, mechanical, etc.) through certain constitutive coupling relations, e.g., elastic compliance and piezoelectric response, and often such instabilities can be shifted by static applied fields thus meaning tunable dynamic properties. Structural instabilities are naturally accommodated on the brink of morphotropic phase boundaries (MPB s) where multiple phases of small energy difference coexist in different crystallographic forms. Canonical MPB is realized through compositional mixture, as is typically exemplified by Pb(Zr1-xTix)O3 solid solutions and relaxor ferroelectrics of (1-x)PbMg1/3Nb2/3O3-xPbTiO3. More recently, a strain-driven MPB has been discovered in BiFeO3 (BFO) thin films epitaxially grown on LaAlO3 (LAO) crystal substrates (which imposes about -4.5% in-plane strains). Such an MPB is in between a rhombohedral (R) phase that bulk BFO exhibits and a so-called super-tetragonal (T) phase, which name hints at its giant lattice axial ratio (c/a ~ 1.23) and accordingly high electric polarization (~1.5 C m-2). The discovery of an MPB in BFO has revealed another facet of this multiferroic system, further adding opportunities to its many exotic functionalities such as domain wall conduction, magnetoelectric and photovoltaic effects As with other MPB s, large electric-field induced strains as well as more underlying lattice softening effects are observed near this MPB promising piezoelectric-based applications. In addition, T-phase BFO itself shows distinct properties, e.g., electronic band gap and optical absorption, from the R-phase and the resultant switching effects between them may also be exploitable. However, unlike conventional ferroelectric oxides where the phases across an MPB usually have subtle difference caused by rotations of an ion off-centering polarization, the BFO system bears multiple structural degrees of freedom, in particular antiferrodistortive modes of oxygen octahedral tilt, and a multitude of structural transition paths are thereby facilitated. Moreover, since the MPB of BFO is driven by epitaxial strain, it is sensitive to the strain relaxation related to film thickness and growth conditions, and above some critical thickness the films appear in the form of a hierarchical mixed-phase microstructure involving several coexistent polymorphs with distinct lattice distortion (monoclinic phases) and tilts. Elastic and electrical heterogeneities are necessary consequences of such complex microstructure, which couple to the intrinsic order parameters and expectedly have profound influence on the structural dynamics and material properties. All these make it a demanding task to obtain a deep understanding of this MPB system on par with its application prospect. The thin-film material form also brings in experimental restrictions to the MPB phase transition studies of BFO since many pertinent techniques fail to operate at the nanoscale or suffer from formidably weak signals. Therefore, static structural characterizations using e.g. X-ray diffraction, electron microscopy and Raman scattering have prevailed thus far. Recently, we adapted band-excitation piezoresponse spectroscopy (BEPS) to probe the bias-induced R/T phase transition dynamics of BFO/SrTiO3 (STO) thin films, and revealed the soft-mode elastic behavior of the transition. The efficacy of our method, along with standard piezoresponse force microscopy (PFM), can be significantly leveraged by integrating a variety of local and/or global excitations, e.g., tip pressure (10 s GPa attainable), heating, photo-irradiation, available to modern scanning probe platforms. This can thus afford unique opportunities to survey the structural dynamics of ferroelectric materials coupled to those field variables, enabling rapid discovery or validation of their functional properties. In this study, we focus on the T-phase BFO/LAO system; we have examined its phase transition behavior not only due to local coaction of tip bias and loading force, but also under device-level global electric fields in a coplanar capacitor. The intrinsic elastic softening phenomena near the structural transitions are comprehensively revealed by BEPS and corroborated by phase-field modeling. Our findings may open a new pathway for technological utilization of the MPB phase instabilities of BFO.

BiFeO3↗