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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Transitions in the morphology and critical stresses of gliding dislocations in multiprincipal element alloys

Refractory multiprincipal element alloys (MPEAs) are promising material candidates for high-temperature, high-strength applications. However, their deformation mechanisms, particularly at the dislocation scale, are unusual compared to those of conventional alloys, making it challenging to understand the origin of their high strength. Here, using an atomistically informed phase-field dislocation dynamics model, we study transitions in the morphology and critical stresses of long gliding screw dislocations over extended distances in a refractory MPEA. The model MPEA crystal accounts for the atomic-scale fluctuations in chemical composition across the glide planes via spatially correlated lattice energies and for the differences in glide resistances between screw and edge dislocations of unit length. Here, we show that the dislocation moves in a stop-start motion, alternating between wavy morphology in free flight and nearly recovered straight screw orientation in full arrest. The periods of wavy glide are due to variable kink-pair formation and migration rates along the length, where portions with higher rates glide more quickly. The critical stress to initiate motion corresponds to the stress required to form and migrate a kink pair at the weakest region along the length of the dislocation. Heterogeneity in lattice energy leads to variability in the local stress-strain response and to a strain hardening-like response, in which the critical stress to reactivate glide increases with glide distance. Statistical assessment of hundreds of realizations of dislocations indicates that the amount of hardening directly scales with the dispersion in underlying lattice energy.

3-dimensional systems↗
DOE OSTI2014

Periodic ordering of clusters and stripes in a two-dimensional lattice model. I. Ground state, mean-field phase diagram and structure of the disordered phases

The short-range attraction and long-range repulsion between nanoparticles or macromolecules can lead to spontaneous pattern formation on solid surfaces, fluid interfaces, or membranes. In order to study the self-assembly in such systems we consider a triangular lattice model with nearest-neighbor attraction and third-neighbor repulsion. At the ground state of the model (T = 0) the lattice is empty for small values of the chemical potential μ, and fully occupied for large μ. For intermediate values of μ periodically distributed clusters, bubbles, or stripes appear if the repulsion is sufficiently strong. At the phase coexistences between the vacuum and the ordered cluster phases and between the cluster and the lamellar (stripe) phases the entropy per site does not vanish. As a consequence of this ground state degeneracy, disordered fluid phases consisting of clusters or stripes are stable, and the surface tension vanishes. For T > 0 we construct the phase diagram in the mean-field approximation and calculate the correlation function in the self-consistent Brazovskii-type field theory.

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

A phase-field formulation for dynamic cohesive fracture

We extend a phase-field/gradient damage formulation for cohesive fracture to the dynamic case. The model is characterized by a regularized fracture energy that is linear in the damage field, as well as non-polynomial degradation functions. Two categories of degradation functions are examined, and a process to derive a given degradation function based on a local stress-strain response in the cohesive zone is presented. The resulting model is characterized by a linear elastic regime prior to the onset of damage, and controlled strain-softening thereafter. The governing equations are derived according to macro- and microforce balance theories, naturally accounting for the irreversible nature of the fracture process by introducing suitable constraints for the kinetics of the underlying microstructural changes. The model is complemented by an efficient staggered solution scheme based on an augmented Lagrangian method. Numerical examples demonstrate that the proposed model is a robust and effective method for simulating cohesive crack propagation, with particular emphasis on dynamic fracture.

42 ENGINEERING↗
DOE OSTI2020

Quantum phase crossover in the spherical mean-field plus quadrupole–quadrupole and pairing model with two j-orbits

The shape phase crossover in the mean-field plus the geometric quadrupole–quadrupole and pairing model within two [Formula: see text]-orbits is analyzed, for which a simple description of [Formula: see text]Sn confined in the lowest [Formula: see text] and [Formula: see text] orbits above the [Formula: see text]Sn core is demonstrated to reveal the crossover behavior of the model. It is shown that [Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text] in this case may serve as effective order parameters for even–even and odd-A systems, respectively, from which the shape phase crossover from the rotation-like phase to the superconducting-like phase can be observed with variation of the number of valence neutrons.

Physics↗
DOE OSTI2015

Multi-scale modeling of microstructure dependent intergranular brittle fracture using a quantitative phase-field based method

In this study, the fracture behavior of brittle materials is strongly influenced by their underlying microstructure that needs explicit consideration for accurate prediction of fracture properties and the associated scatter. In this work, a hierarchical multi-scale approach is pursued to model microstructure sensitive brittle fracture. A quantitative phase-field based fracture model is utilized to capture the complex crack growth behavior in the microstructure and the related parameters are calibrated from lower length scale atomistic simulations instead of engineering scale experimental data. The workability of this approach is demonstrated by performing porosity dependent intergranular fracture simulations in UO 2 and comparing the predictions with experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field simulations of fission gas bubbles in high burnup UO2 to inform engineering-scale fuel performance modeling

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). The cause of pulverization has been hypothesized to be 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. Simulation 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 OSTI2023

A multiscale phase field fracture approach based on the non-affine microsphere model for rubber-like materials

Rubber-like materials have a broad scope of applications due to their unique properties like high stretchability and increased toughness. Hence, computational models for simulating their fracture behavior are paramount for designing them against failures. In this study, the phase field fracture approach is integrated with a multiscale polymer model for predicting the fracture behavior in elastomers. At the microscale, damaged polymer chains are modeled to be made up of a number of elastic chain segments pinned together. Using the phase field approach, the damage in the chains is represented using a continuous variable. Both the bond stretch internal energy and the entropic free energy of the chain are assumed to drive the damage, and the advantages of this assumption are expounded. A framework for utilizing the non-affine microsphere model for damaged systems is proposed here by considering the minimization of a hypothetical undamaged free energy, ultimately connecting the chain stretch to the macroscale deformation gradient. At the macroscale, a thermodynamically consistent formulation is derived in which the total dissipation is assumed to be mainly due to the rupture of molecular bonds. Using a monolithic scheme, the proposed model is numerically implemented and the resulting three-dimensional simulation predictions are compared with existing experimental data. The capability of the model to qualitatively predict the propagation of complex crack paths and quantitatively estimate the overall fracture behavior is verified. Additionally, the effect of the length scale parameter on the predicted fracture behavior is studied for an inhomogeneous system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Development of a grain growth model for U 3 Si 2 using experimental data, phase field simulation and molecular dynamics

The purpose of this work is to develop a model for normal grain growth in U 3 Si 2 . The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. Here, the mobility was determined by comparing phase field grain growth simulations to the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase-field simulations to inform nuclear fuel performance modeling

Software tools to simulate nuclear fuel performance at the engineering scale, such as Idaho National Laboratory (INL)’s BISON code, are increasingly relied upon in regulatory and economic decision-making. However, accurate results from these tools depends on the availability of materials parameters that are used as input. In recent years, atomistic and mesoscale simulation methods have emerged as a cost-effective, expedient means to obtain such input parameters. Phase-field simulations using INL’s Marmot application have been used to obtain microstructure-level parameters and to improve material models for fuel performance modeling using BISON. In this talk, recent examples of this process are given, including applications in UO2, U3Si2, and UZr fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Investigation of Stress Field and Fracture Development During Shale Maturation Using Analog Rock Systems

The emergence of hydrocarbons within shale as a major recoverable resource has sparked interest in fluid transport through these tight mudstones. Recent studies suggest the importance to recovery of microfracture networks that connect localized zones with large organic content to the inorganic matrix. This paper presents a joint modeling and experimental study to examine the onset, formation, and evolution of microfracture networks as shale matures. Both the stress field and fractures are simulated and imaged. A novel laboratory-scale, phase-field fracture propagation model was developed to characterize the material failure mechanisms that play a significant role during the shale maturation process. The numerical model developed consists of coupled solid deformation, pore pressure, and fracture propagation mechanisms. Benchmark tests were conducted to validate model accuracy. Laboratory-grade gelatins with varying Young’s modulus were used as scaled-rock analogs in two dimensional Hele-Shaw cell apparatus. Yeast within the gelatin generates gas in a fashion analogous to hydrocarbon formation as shale matures. These setups allow study and visualization of host rock elastic-brittle fracture and fracture network propagation mechanisms. The experimental setup was fitted to utilize photoelasticity principles coupled with birefringence properties of gelatin to explore visually the stress field of the gelatin as the fracture network developed. Stress optics image analysis and Linear Elastic Fracture Mechanics (LEFM) principles for crack propagation were used to monitor fracture growth for each gelatin type. Observed and simulated responses suggest gas diffusion within and deformation of the gelatin matrix as predominant mechanisms for energy dissipation depending on gelatin strength. LEFM, an experimental estimation of principal stress development with fracture growth, at different stages was determined for each gelatin rheology. Here, the interplay of gas diffusion and material deformation determines the resulting frequency and pattern of fractures. Results correlate with Young’s modulus. Experimental and computed stress fields reveal that fractures resulting from internal gas generation are similar to, but not identical to, type 1 opening mode.

58 GEOSCIENCES↗
DOE OSTI2025

Interface morphology and dislocation-mediated processes during rapid solidification of thin films

Rapid solidification experiments have, in recent years, revealed a wealth of new microstructural phenomena that suggest a strong connection between the kinetics of solidification and the crystalline structures that emerge as a result. In this work, we investigate the interplay between interface morphology and defect-mediated processes during rapid solidification conditions using a Phase Field Crystal (PFC) model, enabling us to simultaneously and efficiently explore the physics of solidification and elasto-plasticity in the formalism of a single-field theory. We predict that there are two mechanisms by which dislocations emitted directly from the solid–liquid interface induce orientation gradients as well as the formation of subgrain boundaries within a single solidifying cell. We relate these mechanisms to the morphology of the moving solid–liquid interface and identify a suitable control parameter in the PFC model with which we can go between said morphologies by effectively changing the relative strength of the capillary length and kinetic coefficients of the solid–liquid interface. Thus, we are able to provide mechanistic explanations for several microstructural features (with an emphasis on orientation gradients and subgrain boundaries) observed during the rapid solidification of pure materials. We also provide a simple explanation for the formation of “jagged” subgrain boundaries, which is consistent with our experimental observations in rapidly solidified samples of Aluminum, whose mechanisms have thus far been unknown.

Interface morphology↗