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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Xenon Hydrate as an Analog of Methane Hydrate in Geologic Systems Out of Thermodynamic Equilibrium

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas-liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas-facing side to the liquid-facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling Ni Coarsening and Migration in the Hydrogen Electrode of Solid oxide Cells under Operating Conditions

Ni coarsening and migration is the most important degradation in the hydrogen electrode of solid oxide cells. This presentation reviews our works on Ni coarsening and migration in NETL over the past years. We used simulation techniques including phase-field modeling and density function theory to investigate possible mechanisms of Ni coarsening and migration, including self-diffusion of Ni, diffusion of gaseous Ni(OH)2, surface diffusion of Ni(OH)x and Ni-YSZ wettability change, and examined the effect of operating conditions on these mechanisms. So far, none of the mechanisms can fully explain the experiments. Remaining questions and possible paths forward are summarized and discussed.

Lei, Yinkai↗
DOE OSTI2022

Electrochemically-Induced Phase Transformations in Battery Storage Compounds (Final Technical Report)

Compounds of interest for ion storage in advanced batteries frequently exhibit phase transformations as the working ion concentration varies. Under large electrochemical driving forces inherent to practical use, systems are often driven far from equilibrium. This program combines experiments and theory to understand the phase transition behavior of ion insertion compounds when electrochemically driven far from equilibrium. As model systems, we focus on alkaline metal phosphates AMPO4 (A = alkali; M = first row transition metal) of olivine structure, which are both technologically interesting and ideally suited for fundamental study due to the ability to systematically tune transformation strain, and along with it, the phase transformation pathway. Behavior in compositions having large transformation strains (~15 vol%) requiring plasticity for strain accommodation is emphasized. Experimental techniques include operando characterization of structure while simultaneously varying electrokinetic parameters, and high resolution microscopy of nanoscale and interfacial phenomena. Phase-field modeling is used to model the thermodynamics and kinetics of competing transformation pathways, extended to include the effects of plasticity, and integrated with porous electrode kinetic theory to treat multi-particle effects. Success in this project will lead to an ability to design ion storage compounds with predictable transformation pathways, electrochemical kinetics, capacity utilization, and durability. New technologically important compounds may also be discovered.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗
DOE OSTI2022

Optimization of WAAM Process to Produce AUSC Components with Increased Service Life

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. Wire Arc Additive manufacturing (WAAM) process is a class of directed energy deposition process enabling higher build rate and allowing custom wire feedstock allowing spatial variation of microstructure. The larger spot size and lower speed creates a larger melt pool, reducing residual stress and often time creating directional/columnar microstructure for Ni-superalloy. The use of flexible platform provides freedom in deposition strategy, which can accommodate complex substrates, including feature addition onto existing structures, non-flat layers, and repair methods. However, the certification of the final component needs to match the strength requirement. Hence, the quality of the component produced is stringently monitored to avoid buildup of residual stress, cracks, porosity and to reduce detrimental segregated phases commonly observed during alloy solidification. Simulation plays a huge role in predicting the melt pool dimension and can be used to optimize the process parameter. Similar development is also required to perform physics-based modeling of microstructural development in WAAM that can predict the microstructural features during solidification and can be used to optimize the process more effectively. To move toward this goal, Raytheon Technologies Research Center together with Siemens worked to create a set of computational tools to control the process parameters, enable on-line measurements and acquisition with feedback to the optimized process parameters, and eventually track material evolution through each step of the additive process. Computational fluid dynamics is used for accurate prediction and calibration of the thermal field during WAAM process. and phase field models for microstructure evolution as a function of processing parameters to establish a connection between additive parameters and the final microstructure. Here we report cellular automata (CA) model development to predict the dendritic microstructure evolution with surface and bulk nuclei for single track and multiple layers. The CA model was developed to account for secondary element addition and predict segregation, local melting, and latent heat release as well as prediction of Euler angles from orientation information and validated against experiments. This framework was utilized to tailor spatially-varying composition in a part by appropriately controlling the microstructure evolution during the additive process. Functionally graded Haynes 282 alloy with high Cr content at the surface was tested for oxidation and mechanical properties. An advanced physics-based reaction-diffusion model predicting the simultaneous creation of chromium oxide and alumina is developed and validated to extend life expectancy of the WAAM manufactured high temperature part. A machine-learning data-driven framework establishing the process-structure relationship from a dataset of real microstructure images and corresponding process history data has been developed and implemented in the NX Siemens design system. The digital twin configuration along with the tool path generation enabled prediction of WAAM component buildup time and the techno-economic analysis provided a favorable option for all 4 cases with 15-40% cost reduction.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Modeling Ni Coarsening under Humid Atmosphere in the Electrode of Solid Oxide Cells

Ni coarsening is an important degradation mechanism in solid oxide cells. It is reported that Ni coarsening is faster under a humid atmosphere experimentally, but the underline mechanism is not well understood. In this work, Ni coarsening through Ni(OH)x surface diffusion was investigated by a combination of density-functional theory (DFT) and phase-field modeling (PFM). DFT is used to evaluate the surface coverage and diffusivity of Ni(OH)x on Ni (111) surface and the results are used as an input to the PFM to simulate Ni coarsening under humid atmospheres on both reconstructed and synthetic microstructures. It is found that steam partial pressure and microstructure effect alone cannot explain the faster Ni coarsening under a humid atmosphere observed experimentally. Large local overpotential in fuel cell mode is needed to explain the fast Ni coarsening.

Lei, Yinkai↗
DOE OSTI2017

TA [B] Predicting Microstructure-Creep Resistance Correlation in High Temperature Alloys over Multiple Time Scales

DoE-NETL partnered with Purdue University to predict the creep and associated microstructure evolution of tungsten-based refractory alloys. Researchers use grain boundary (GB) diagrams, a new concept, to establish time-dependent creep resistance and associated microstructure evolution of grain boundaries/intergranular films GB/IGF controlled creep as a function of load, environment, and temperature. The goal was to conduct a systematic study that includes the development of a theoretical framework, multiscale modeling, and experimental validation using W-based body-centered-cubic alloys, doped/alloyed with one or two of the following elements: nickel, palladium, cobalt, iron, and copper—typical refractory alloys. Prior work has already established and validated a basic theory for W-based binary and ternary alloys; the study conducted under this project extended this proven work. Based on interface diagrams phase field models were developed to predict long term microstructural evolution. In order to validate the models nanoindentation creep data was used to elucidate the role played by the interface properties in predicting long term creep strength and microstructure evolution.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Simulating microgalvanic corrosion in alloys using the PRISMS phase-field framework

In this prospective paper, we first review the existing simulation tools to simulate microgalvanic corrosion during free immersion. Then, we describe a recently developed application that employs PRISMS-PF, an open-source, high-performance phase-field modeling framework. The model employed in the application accounts for the electrochemical reaction at the metal/electrolyte interface and ionic migration in the electrolyte to determine the evolution of the corrosion front. We present the implementation details for the application and discuss its features such as super-linear parallel scaling performance for a sufficiently large system. Finally, we demonstrate the capability of the application by simulating corrosion of the matrix phase of an alloy near a secondary phase particle in two and three dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2016

A phase-field approach to model multi-axial and microstructure dependent fracture in nuclear grade graphite

The fracture behavior of nuclear grade graphites is strongly influenced by underlying microstructural features such as the character of filler particles, and the distribution of pores and voids. These microstructural features influence the crack nucleation and propagation behavior, resulting in quasi-brittle fracture with a tortuous crack path and significant scatter in measured bulk strength. This paper uses a phase-field method to model the microstructural and multi-axial fracture in H-451, a historic variant of nuclear graphite that provides the basis for an idealized study on a legacy grade. The representative volume elements are constructed from randomly located pores with random size obtained from experimentally determined log-normal distribution. The representative volume elements are then subjected to simulated multi-axial loading, and a reasonable agreement of the resulting fracture stress with experiments is obtained. Finally, quasi-brittle stress-strain evolution with a tortuous crack path is also observed from the simulations and is consistent with experimental results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A variational phase-field framework for thermal softening and dynamic ductile fracture

A variational phase field model for dynamic ductile fracture is presented. Here, the model is designed for elasto-viscoplastic materials subjected to rapid deformations in which the effects of heat generation and material softening are dominant. The variational framework allows for the consistent inclusion of plastic dissipation in the heat equation as well as thermal softening. It employs a coalescence function to degrade fracture energy during regimes of high plastic flow. A variationally consistent form of the Johnson–Cook model is developed for use with the framework. Results from various benchmark problems in dynamic ductile fracture are presented to demonstrate capabilities. In particular, the ability of the model to regularize shear band formation and subsequent damage evolution in two- and three-dimensional problems is demonstrated. Importantly, these phenomena are naturally captured through the underlying physics without the need for phenomenological criteria such as stability thresholds for the onset of shear band formation.

42 ENGINEERING↗
DOE OSTI2017

A New Superalloy Enabling Heavy Duty Gas Turbine Wheels for Improved Combined Cycle Efficiency

The drive to increase combined cycle turbine efficiency from 62% to 65% for the next-generation advanced cycle requires a new heavy duty gas turbine wheel material capable of operating at 1200°F and above. Current wheel materials are limited by the stability of their major strengthening phase (gamma double prime), which coarsens at temperatures approaching 1200°F, resulting in a substantial reduction in strength. More advanced gamma prime superalloys, such as those used in jet engine turbine disks, are also not suitable due to size constraints; the gamma prime phase overages during the slow cooling rates inherent in processing thick-section turbine wheels. The current program addresses this need by screening two new alloy design concepts. The first concept exploits a gamma prime/gamma double prime coprecipitation reaction. Through manipulation of alloy chemistry, coprecipitation is controlled such that gamma double prime is used only to slow the growth of gamma prime during slow cooling, preventing over-aging, and allowing for subsequent heat treatment to maximize strength. In parallel, phase field modeling provides fundamental understanding of the coprecipitation reaction. The second concept uses oxide dispersion strengthening to improve on two existing alloys that exhibit excellent hold time fatigue crack growth resistance, but have insufficient strength to be considered for gas turbine wheels. Mechanical milling forces the dissolution of starting oxide powders into a metal matrix allowing for solid state precipitation of new, nanometer scale oxides that are effective at dispersion strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Real-space imaging of periodic nanotextures in thin films via phasing of diffraction data

New properties and exotic quantum phenomena can form due to periodic nanotextures, including Moire patterns, ferroic domains, and topologically protected magnetization and polarization textures. Despite the availability of powerful tools to characterize the atomic crystal structure, the visualization of nanoscale strain-modulated structural motifs remains challenging. Here, we develop nondestructive real-space imaging of periodic lattice distortions in thin epitaxial films and report an emergent periodic nanotexture in a Mott insulator. Specifically, we combine iterative phase retrieval with unsupervised machine learning to invert the diffuse scattering pattern from conventional X-ray reciprocal-space maps into real-space images of crystalline displacements. Our imaging in PbTiO 3 /SrTiO 3 superlattices exhibiting checkerboard strain modulation substantiates published phase-field model calculations. Furthermore, the imaging of biaxially strained Mott insulator Ca 2 RuO 4 reveals a strain-induced nanotexture comprised of nanometer-thin metallic-structure wires separated by nanometer-thin Mott-insulating-structure walls, as confirmed by cryogenic scanning transmission electron microscopy (cryo-STEM). The nanotexture in Ca 2 RuO 4 film is induced by the metal-to-insulator transition and has not been reported in bulk crystals. We expect the phasing of diffuse X-ray scattering from thin crystalline films in combination with cryo-STEM to open a powerful avenue for discovering, visualizing, and quantifying the periodic strain-modulated structures in quantum materials.

Science & Technology - Other Topics↗
DOE OSTI2021

Development of High-Performance Ni-base Alloys for Gas Turbine Wheels Using a Coprecipitation Approach (Final Technical Report)

In this project, a synergistic approach involving experiments informed by modelling was used to develop γ'/γ'' coprecipitation-strengthened Ni-base superalloys for land-based gas turbine wheels. In thick sections of such turbine wheels, the slow cooling rate normally observed, as well as variations in cooling rate could potentially result in drastically inhomogeneous and coarsened microstructure. Thus, ten alloy compositions were selected based on coprecipitation and sluggish γ’ growth and coarsening during three iterations of alloy development, assisted by extensive CALPHAD modeling. Phase equilibrium calculations, informed by the characterization of early alloys, were successful in predicting suitable aging treatments to obtain bimodal coprecipitate distribution in several alloys. These aging treatments were applied to the subsequent round of alloys and the resulting microstructures were characterized in detail. Two different cooling rates were used during each of these heat treatments to simulate the range of cooling rates observed in a thick section forging. The microstructures were used to calibrate a multi-phase field model (MPF), which was used to study the formation mechanisms of various coprecipitate morphologies observed in the experiments. Following this, the most promising alloys in each round were selected and studied using high temperature tensile and creep tests. Considering that a slower cooling rate would lead to a strength debit, while a faster cooling rate might lead to a decline in creep rupture life of a full-scale wheel, tensile tests were carried out for slower cooled samples, while creep tests were carried out for faster cooled samples. The deformation microstructures were then characterized using advanced characterization techniques such as Diffraction Contrast-Scanning Transmission Electron Microscopy (DC-STEM) and Electron Channeling Contrast Imaging (ECCI). Further, a mean-field creep model was developed and calibrated using these experimental findings to predict long-term creep behavior of these coprecipitation-strengthened alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Simulation studies of nucleation of ferroelectric polarization reversal.

Electric field-induced reversal of spontaneous polarization is the defining characteristic of a ferroelectric material, but the process(es) and mechanism(s) associated with the initial nucleation of reverse-polarity domains are poorly understood. This report describes studies carried out using phase field modeling of LiTaO 3 , a relatively simple prototype ferroelectric material, in order to explore the effects of either mechanical deformation or optically-induced free charges on nucleation and resulting domain configuration during field-induced polarization reversal. Conditions were selected to approximate as closely as feasible those of accompanying experimental work in order to provide not only support for the experimental work but also ensure that additional experimental validation of the simulations could be carried out in the future. Phase field simulations strongly support surface mechanical damage/deformation as effective for dramatically reducing the overall coercive field (Ec) via local field enhancements. Further, optically-nucleated polarization reversal appears to occur via stabilization of latent nuclei via the charge screening effects of free charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Development of a Gibbs Energy Minimiser for the MOOSE-based Corrosion Modelling App Yellowjacket and Validation of MSTDB

Nuclear materials are highly complex multiscale, multiphysics systems,and an effective prediction of nuclear reactor performance and safety requires simulation capabilities that tightly couple different physical phenomena. The Idaho National Laboratory’s Multiphysics Object Oriented Simulation Environment (MOOSE) provides the computational foundation for performing such simulations. With the move towards advanced reactors, such as the Molten Salt Reactor (MSR), that employ high temperature fluids compared to conventional reactors, corrosion has become a problem of great interest. A new application called Yellowjacket is currently under development to directly couple thermodynamic equilibrium and kinetics with phase field models in order to model corrosion in MSRs. As part of Yellowjacket, a Gibbs energy minimiser is being developed to perform thermochemical equilibrium calculations for a range of different materials, which is currently in its infancy. This report describes the further progress towards the development of Yellowjacket Gibbs energy minimiser. Ontario Tech University is developing a new Gibbs energy minimiser for Yellowjacket which is the primary contribution of this work. The aim to develop a thermochemistry solver for the MOOSE framework following the same development philosophy and using the same tools and libraries. A special focus is on performance, documentation and SQA. Furthermore through a scope extension partway through the fiscal year a thorough assessment of the MSTDB-TC v1.3 was performed at Ontario Tech University in the context of continuous improvement and quality assurance. The objective of the work was to have an arm’s length review of the database to give confidence that the database is performing as it was intended while assessing its current state to give recommendations to future developments. This assessment involved two parts: A) a quantitative assessment, and B) a qualitative assessment. Part A involved developing an automated test-suite that would compute values from the database using Thermochimica with comparisons to experimental measurements for validation purposes, which gives confidence to the database’s stakeholders that its functioning properly. Part B involved reviewing all binary systems in the database and making a qualitative assessment with two performance indicators: comprehensiveness and overall confidence. It is important to note that the models in the database are empirical, which is to say that the quality of any model is highly dependent on the experimental data used to inform its development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

In-situ Study of Dynamic Phenomena at Metal Nanosolder Interfaces Using Aberration Corrected Scanning Transmission Electron Microcopy.

Controlling metallic nanoparticle (NP) interactions plays a vital role in the development of new joining techniques (nanosolder) that bond at lower processing temperatures but remain viable at higher temperatures. The pr imary objective of this project is t o develop a fundamental understanding of the actual reaction processes, associated atomic mechanisms, and the resulting microstructure that occur during thermally - driven bond formation concerning metal - metal nano - scale (%3C50nm) interfaces. In this LDRD pr oject, we have studied metallic NPs interaction at the elevated temperatures by combining in - situ transmission electron microscopy (TEM ) using an aberration - corrected scanning transmission electron microscope (AC - STEM) and atomic - scale modeling such as m olecular dynamic (MD) simulations. Various metallic NPs such as Ag, Cu and Au are synthesized by chemical routines. Numerous in - situ e xperiments were carried out with focus of the research on study of Ag - Cu system. For the first time, using in - situ STEM he ating experiments , we directly observed t he formation of a 3 - dimensional (3 - D) epitaxial Cu - Ag core - shell nanoparticle during the thermal interaction of Cu and Ag NPs at elevated temperatures (150 - 300 o C). The reaction takes place at temperatures as low as 150 o C and was only observed when care was taken to circumvent the effects of electron beam irradiation during STEM imaging. Atomic - scale modeling verified that the Cu - Ag core - shell structure is energetically favored, and indicated that this phenomenon is a nano - scale effect related to the large surface - to - volume ratio of the NPs. The observation potentially can be used for developing new nanosolder technology that uses Ag shell as the "glue" that stic ks the particles of Cu together. The LDRD has led to several journal publications and numerous conference presentations, and a TA. In addition, we have developed new TEM characterization techniques and phase - field modeling tools that can be used for future materials research at Sandia. Acknowledgeme nts This work was supported by the Laboratory Directed Research and Development (LDRD) program of Sandia National Laboratories. Sandia National Laboratories is a multi - program laboratory managed and operated by Sandia Corporation, a wholly owned subsidia ry of Lockheed Martin Corporation, for the US Department of Energy's National Nuclear Security Administration under contract DE - AC04 - 94AL85000.

Lu, Ping↗
DOE OSTI2016

Tunable thermal conductivity via domain structure engineering in ferroelectric thin films: A phase-field simulation

Effective thermal conductivity as a function of domain structure is studied by solving the heat conduction equation using a spectral iterative perturbation algorithm in materials with inhomogeneous thermal conductivity distribution. Using this proposed algorithm, the experimentally measured effective thermal conductivities of domain-engineered {001} p- BiFeO 3 thin films are quantitatively reproduced. In conjunction with two other testing examples, this proposed algorithm is proven to be an efficient tool for interpreting the relationship between the effective thermal conductivity and micro-/domain-structures. By combining this algorithm with the phase-field model of ferroelectric thin films, the effective thermal conductivity for PbZr 1-x Ti x O 3 films under different composition, thickness, strain, and working conditions is predicted. It is shown that the chemical composition, misfit strain, film thickness, film orientation, and a Piezoresponse Force Microscopy tip can be used to engineer the domain structures and tune the effective thermal conductivity. Furthermore, we expect our findings will stimulate future theoretical, experimental and engineering efforts on developing devices based on the tunable effective thermal conductivity in ferroelectric nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗