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

Toward a Theory of the Evolution of Drop Morphology and Splintering by Freezing

Abstract The drop freezing process is described by a phase-field model. Two cases are considered: when the freezing is triggered by central nucleation and when nucleation occurs on the drop surface. Depending on the environmental temperature and drop size, different morphological structures develop. Detailed dendritic growth was simulated at the first stage of drop freezing. Independent of the nucleation location, a decrease in temperature within the range from ~−5° to −25°C led to an increase in the number of dendrites and a decrease in their width and the interdendritic space. At temperatures lower than about −25°C, a planar front developed following surface nucleation, while dendrites formed a granular-like structure with small interdendritic distances following bulk nucleation. An ice shell grew in at the same time (but slower) as dendrites following surface nucleation, while it started forming once the dendrites have reached the drop surface in the case of central nucleation. The formed ice morphology at the first freezing stage predefined the splintering probability. We assume that stresses needed to break the ice shell arose from freezing of the water in the interdendritic spaces. Under this assumption, the number of possible splinters/fragments was proportional to the number of dendrites, and the maximum rate of splintering/fragmentation occurred within a temperature range of about −10° to −20°C, is in agreement with available laboratory and in situ measurements. At temperatures < −25°C, freezing did not lead to the formation of significant stresses, making splintering unlikely. The number of dendrites increased with drop size, causing a corresponding increase in the number of splinters. Examples of morphology that favors drop cracking are presented, and the duration of the freezing stages is evaluated. Sensitivity of the freezing process to the surface fluxes is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

VTO FY20Q4 quarterly report

Project Objective.This project has three objectives: (1) tuning microstructures of 3D-printed SSE separators,(2) determining material and processing compatibilities with cathode printing,and (3)3D printing of sintering-free SSE separators.Project Impact.All-solid-state lithium batteries are difficult to process due to the brittleness of ceramic materials, poor solid-solid contact, and electrolyte-electrode stability issues. As a result, the energy and power density and also cycling stability are far from satisfying. This project will address fabrication difficulties by using state-of-the-art 3D-printing techniques that can introduce 3D interfaces and architectures to enhance solid-solid contact and reduce charge transfer resistance. Success will benefit the DOE by establishing the best manufacturing methods for ASSBsto achieve VTO goals on the performance of beyond Li-ion batteries. Approach.The project employs 3D-printing techniques to manufacture SSEs and related components for ASSBs. The team starts with direct ink writing (DIW) to develop ink recipes with desired rheological properties and explore post-sintering approaches to achieve high densification. In parallel with DIW3D printing, the team also explores other 3D printing options, such as projection microstereolithography (PuSL). Both sintering and sintering-free approaches can gain benefit from 3D printing and therefore will be investigated.Out-Year Goals.The team will determine particle sizes and morphologies as well as post-processing conditions that deliver good ionic conductivity and charge transfer properties of LLZTO separators. They will examine co-sintering of LLZTO/cathode and determine the effects of material, particle size, surface coating, and sintering conditions.Collaborations.Microstructures, ionic conductivities, and mechanical properties will be provided to the LLNL simulation team, led by PI B. Wood, for establishing and validating phase-field modeling methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational Tools for Additive Manufacture of Tailored Microstructure and Properties

Additive manufacturing has the potential to revolutionize industrial hardware and unlock efficiency gains through the fabrication of geometries and architectures not possible by conventional processing. Currently most additive builds use a single set of process parameters (e.g. laser power and scan speed) which results in a part with a homogenous microstructure that provides a singular performance level. To move beyond this state, Raytheon Technologies Research Center worked to create a set of computational tools to track material evolution through each step of the additive process. Computational fluid dynamics and phase field models for microstructure evolution as a function of processing parameters, and crystal plasticity models fully coupling microstructure and mechanical properties for performance predictions were leveraged to establish a connection between additive parameters and the final microstructure. This framework was utilized to tailor spatially-varying mechanical properties in a part by appropriately controlling the microstructure evolution during the additive process. Specifically, a turbine blade was 3D printed from nickel superalloy IN718 using laser powder bed fusion with coarse grains in the airfoil section which experiences the highest temperatures and is creep limited while finer grains were printed in the root of the blade which experience higher stresses but at lower temperatures and is therefore fatigue limited. The benefit of being able to intentionally insert coarse grains in the high temperature region of the blade was showcased with a microstructure sensitive creep model that indicates longer creep life for coarser grains.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Effect of Processing on Precipitation Kinetics in Nano Polycrystalline 7075 Aluminum Alloy

Fairmount Technologies (FT) seeks to optimize the thermomechanical processing of Aluminum Alloy 7075. The target yield strength is 20% higher than commercially available 7075-T6 temper sheet which will enable structural light weighting by a similar magnitude and reduce the carbon footprint of the transportation sector. To reduce processing cost and make this commercially viable, it is important to understand physical mechanisms driving microstructural changes. A mesoscale phase-field model seeded with experimentally measured microstructures and assessed thermodynamic and kinetic properties will be developed and exercised using HPC at PNNL. It will account for the effect of deformation, temperature, and chemistry on precipitation kinetics in nanocrystalline microstructures. Our objective is to understand this “stress-aging” coupling by performing parametric simulations with mesoscale models for a better understanding of GP zone and η' formation in nanocrystalline 7075 Al. This will help FT optimize processing schedules to achieve desired microstructures that increase both strength and ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray movies 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 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, 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 7. 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.

Chueh, William↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2024

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
DOE OSTI2021

Accelerating phase-field-based microstructure evolution predictions via surrogate models trained by machine learning methods

Abstract The phase-field method is a powerful and versatile computational approach for modeling the evolution of microstructures and associated properties for a wide variety of physical, chemical, and biological systems. However, existing high-fidelity phase-field models are inherently computationally expensive, requiring high-performance computing resources and sophisticated numerical integration schemes to achieve a useful degree of accuracy. In this paper, we present a computationally inexpensive, accurate, data-driven surrogate model that directly learns the microstructural evolution of targeted systems by combining phase-field and history-dependent machine-learning techniques. We integrate a statistically representative, low-dimensional description of the microstructure, obtained directly from phase-field simulations, with either a time-series multivariate adaptive regression splines autoregressive algorithm or a long short-term memory neural network. The neural-network-trained surrogate model shows the best performance and accurately predicts the nonlinear microstructure evolution of a two-phase mixture during spinodal decomposition in seconds, without the need for “on-the-fly” solutions of the phase-field equations of motion. We also show that the predictions from our machine-learned surrogate model can be fed directly as an input into a classical high-fidelity phase-field model in order to accelerate the high-fidelity phase-field simulations by leaping in time. Such machine-learned phase-field framework opens a promising path forward to use accelerated phase-field simulations for discovering, understanding, and predicting processing–microstructure–performance relationships.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A brief review on strain engineering of ferroelectric K x Na 1- x NbO 3 epitaxial thin films: Insights from phase-field simulations

Strains play a pivotal role in determining the phase equilibrium, domain configuration, and functional properties of the low-dimensional ferroelectrics. There is growing interest in the strain engineering of ferroelectric K x Na 1- x NbO 3 (KNN) epitaxial thin films, which exhibit excellent physical properties and promise as eco-friendly alternatives to lead-based ferroelectrics for microdevice applications. Further, advances have been made in understanding the phase equilibria and transitions, domains and domain walls, and their relations to the physical properties of KNN epitaxial thin films using a combination of experiments and theoretical modeling, particularly phase-field simulations. Here, we review recent progress in these aspects and showcase the phase-field method for establishing strain phase diagrams, elucidating the domain and domain wall structures at equilibrium, and predicting the structure–property relationships in ferroelectric KNN thin films. We also discuss challenges and opportunities to further advance our understanding of KNN thin films and potentially unlock new functionalities by leveraging phase-field simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2021

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Utah FORGE 2-2446: Closing the Loop Between In-situ Stress Complexity and EGS Fracture Complexity - Workshop Presentation

This is a presentation on the Closing the Loop Between In-situ Stress Complexity and EGS Fracture Complexity project by Lawrence Livermore National Laboratory, presented by Dr. Matteo Cusini. The project's objective was to employ a combination of high-fidelity simulations and true-triaxial block fracturing tests at high temperature to explore the intricate relationship between in-situ stress and hydraulic fracture patterns and better characterize the in-situ stress at Utah FORGE. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗