Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase field”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

Magnetic structure and magnetoelectric properties of the spin-flop phase in LiFePO 4

We investigate the magnetic structure and magnetoelectric(ME) effect in the high-field phase of the antiferromagnet LiFePO 4 above the critical field of 31 T. A neutron diffraction study in pulsed magnetic fields reveals the propagation vector to be q = 0 for the high-field magnetic structure. Pulsed-field electric polarization measurements show that, at the critical field, the low-field off-diagonal ME coupling α ab is partially suppressed, and the diagonal element α bb emerges. These results are consistent with a spin-flop transition where the spin direction changes from primarily being along the easy b axis below the transition to being along a above. The persistence of off-diagonal ME tensor elements above the critical field suggests a lowering of the magnetic point-group symmetry and hence a more complex magnetic structure in the high-field phase. In addition, neutron diffraction measurements in low magnetic fields show no observable field-induced spin canting, which indicates a negligible Dzyaloshinskii-Moriya interaction. The observed spin-flop field supports the Hamiltonian recently deduced from inelastic neutron studies and indicates that the system is less frustrated and with a larger single-ion anisotropy than originally thought. Our results demonstrate the effectiveness of combining pulsed-field neutron diffraction and electric polarization measurements to elucidate the magnetic structures and symmetries at the highest attainable field strengths.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structural phase transitions in the Ag{sub 2}Nb{sub 4}O{sub 11}-Na{sub 2}Nb{sub 4}O{sub 11} solid solution

The phase transitions between various structural modifications of the natrotantite-structured system xAg{sub 2}Nb{sub 4}O{sub 11}-(1-x)Na{sub 2}Nb{sub 4}O{sub 11} have been investigated and a phase diagram constructed as a function of temperature and composition. This shows three separate phase transition types: (1) paraelectric-ferroelectric, (2) rhombohedral-monoclinic and (3) a phase transition within the ferroelectric rhombohedral zone between space groups R3c and R3. The parent structure for the entire series has space group R3{sup Macron }c. Compositions with x>0.75 are rhombohedral at all temperatures whereas compositions with x<0.75 are all monoclinic at room temperature and below. At x=0.75, rhombohedral and monoclinic phases coexist with the phase boundary below room temperature being virtually temperature-independent. The ferroelectric phase boundary extends into the monoclinic phase field. No evidence was found for the R3-R3c phase boundary extending into the monoclinic phase field and it is concluded that a triple point is formed. - Graphical abstract: Phase diagram for xAg{sub 2}Nb{sub 4}O{sub 11}-(1-x)Na{sub 2}Nb{sub 4}O{sub 11} solid solution showing changes in crystal symmetry as a function of temperature and composition. The crystal structure is depicted. Highlights: Black-Right-Pointing-Triangle Ferroelectric, rhombohedral Ag{sub 2}Nb{sub 4}O{sub 11} in solid solution with monoclinic Na{sub 2}Nb{sub 4}O{sub 11}. Black-Right-Pointing-Triangle Three phase boundaries were studied as a function of composition and temperature. Black-Right-Pointing-Triangle Both rhombohedral and monoclinic variants exhibit ferroelectricity. The parent phase of the series has space group R3{sup Macron }c.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of High Magnetic Fields on Phase Transformations in Amorphous Nd2Fe14B

We briefly summarize the results from a set of experiments designed to demonstrate the effects of high magnetic fields applied during thermal annealing of amorphous Nd2Fe14B produced through melt-spinning. A custom-built differential scanning calorimeter was used to determine the crystallization temperatures in zero-field and in applied fields of 20 kOe and 90 kOe, which guided subsequent heat treatments to evaluate phase evolution. X-ray diffraction was used for phase identification and transmission electron microscopy was employed for observation of the crystallite size and morphology. Magnetization measurements were also used to evaluate the resulting magnetic phases after thermomagnetic processing. While the applied magnetic fields do not appear to affect the crystallization temperature, significant effects on the kinetics of phase evolution are observed and correlated strongly to the magnetic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

HPC for Optimizing Process Parameters to Control Material Evolution in Seamless Induction Hardening of Wind Turbine Main Shaft Bearings

Work proposed in this project focused on understanding the effect of martensitic transformation in the steel on the potential for cracking during seamless induction hardening (SIH) as a function of process conditions to allow the process to optimally scale up. Large-scale, three-dimensional phase-field simulations of martensitic transformation were performed using MEUMAPPS-SS (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid State) code developed at Oak Ridge National Laboratory. The simulations were guided by location-specific thermal history generated by experimental measurements of time-temperature history generated at The Timken Company. The simulations were able to capture the morphological evolution of the martensite variants in an Fe-1.0C-1.5Cr steel based on the Nishiyama-Wasserman (NW) orientation relationship. The simulations were also able to quantify the stress-state at the interface between impinging martensite variants. The simulations indicated that the magnitude of the various stress and strain components were dependent on the sizes of the impinging plates with a reduction in these quantities with reduced plate size in agreement with experimental findings. The results obtained from the simulations will be used to guide the optimization of the alloy thermal conditions to eliminate quench cracking during SIH of bearing steels.

99 GENERAL AND MISCELLANEOUS↗

External-Field-Induced Phase Transformation and Associated Properties in a Ni 50 Mn 34 Fe 3 In 13 Metamagnetic Shape Memory Wire

Metamagnetic shape memory alloys exhibit a series of intriguing multifunctional properties and have great potential for applications in magnetic actuation, sensing and magnetic refrigeration. However, the poor mechanical properties of these alloys with hardly any tensile deformability seriously limit their practical application. In the present work, we developed a Ni-Fe-Mn-In microwire that exhibits both a giant, tensile superelasticity and a magnetic-field-induced first-order phase transformation. The recoverable strain of superelasticity is more than 20% in the temperature range of 233–283 K, which is the highest recoverable strain reported heretofore in Ni-Mn-based shape memory alloys (SMAs). Moreover, the present microwire exhibits a large shape memory effect with a recoverable strain of up to 13.9% under the constant tensile stress of 225 MPa. As a result of the magnetic-field-induced first-order phase transformation, a large reversible magnetocaloric effect with an isothermal entropy change ΔS m of 15.1 J kg –1 K –1 for a field change from 0.2 T to 5 T was achieved in this microwire. The realization of both magnetic-field and tensile-stress-induced transformations confers on this microwire great potential for application in miniature multi-functional devices and provides an opportunity for multi-functional property optimization under coupled multiple fields.

36 MATERIALS SCIENCE↗

Holistic Microstructure Control Strategies in Photopolymerization‐Induced Phase Separation of Acrylate Systems

Open porous materials, known for their large surface area and interconnected structures, are essential in various applications, including batteries, ion exchange, catalysis, filtration, and electronic waste recycling. A critical aspect of the functionality of porous membranes is the precise control of pore size and morphology. Photopolymerization-induced phase separation (photo-PIPS) offers a convenient and versatile methods for creating porous structures. However, controlling the porous morphology remains challenging due to the complex interplay between thermodynamics, polymerization kinetics, and monomer structures, which makes it difficult to establish the relationship between processing conditions and resulting morphology in photo-PIPS. Herein, a physics-based phase-field model capable of generating and characterizing the microstructures of porous materials based on both average and localized features is developed. Using the phase-field simulations as test bed, the effects of polarity, light intensity, and curing temperature, as well as the previously unexplored roles of chain transfer agents and substrates, on the morphology of the resulting porous microstructure are investigated. Experiments are performed to verify the results predicted by the simulations. This work lays out a comprehensive guide for designing PIPS-derived porous microstructures and offers practical engineering strategies for tailoring microstructure-level topology and size of pores for application-specific needs.

36 MATERIALS SCIENCE↗

Preferential Paths of Air-water Two-phase Flow in Porous Structures with Special Consideration of Channel Thickness Effects

Accurate understanding and predicting the flow paths of immiscible two-phase flow in rocky porous structures are of critical importance for the evaluation of oil or gas recovery and prediction of rock slides caused by gas-liquid flow. A 2D phase field model was established for compressible air-water two-phase flow in heterogenous porous structures. The dynamic characteristics of air-water two-phase interface and preferential paths in porous structures were simulated. The factors affecting the path selection of two-phase flow in porous structures were analyzed. Transparent physical models of complex porous structures were prepared using 3D printing technology. Tracer dye was used to visually observe the flow characteristics and path selection in air-water two-phase displacement experiments. The experimental observations agree with the numerical results used to validate the accuracy of phase field model. The effects of channel thickness on the air-water two-phase flow behavior and paths in porous structures were also analyzed. The results indicate that thick channels can induce secondary air flow paths due to the increase in flow resistance; consequently, the flow distribution is different from that in narrow channels. This study provides a new reference for quantitatively analyzing multi-phase flow and predicting the preferential paths of immiscible fluids in porous structures.

42 ENGINEERING↗

Multi-scale modeling of the electric field assisted sintering process

The electric field assisted sintering (EFAS) process involves tightly coupled physics that influence microstructural evolution in the particles being compacted. It is also an inherently multi-scale phenomenon, with the microstructure of the compact influencing the subsequent engineering-scale response of the sintering system. To improve understanding of how processing parameters influence microstructural evolution, we have developed a multi-scale modeling approach that couples a continuum-level model of the sintering system with a phase-field model for microstructural evolution of particles within the compact. The phase-field model couples the effect of chemical and electrical driving forces on microstructural evolution and includes the effect of charged defect segregation to surfaces and grain boundaries; this segregation leads to enhanced defect transport and heat generation at these interfaces in response to applied electric field. The effect of enhanced heat generation on particle neck growth and the influence of microstructural evolution on the engineering-scale model are demonstrated.

36 - MATERIALS SCIENCE↗

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Development of Macroscale Models of UO 2 Fuel Sintering and Densification using Multiscale Modeling and Simulation

Light water reactor fuel pellets are fabricated using sintering to final densities of 95% or greater. During reactor operation, the porosity remaining in the fuel after fabrication decreases further due to irradiation-assisted densification. While empirical models have been developed to describe this densification process, a mechanistic model is needed as part of the ongoing work by the NEAMS program to develop a more predictive fuel performance code. In this work we will develop a phase field model of sintering of UO 2 in the MARMOT code, and validate it by comparing to published sintering data. We will then add the capability to capture irradiation effects into the model, and use it to develop a mechanistic model of densification that will go into the BISON code and add another essential piece to the microstructure-based materials models. The final step will be to add the effects of applied fields, to model field-assisted sintering of UO 2 . The results of the phase field model will be validated by comparing to data from field-assisted sintering. Tasks over three years: 1. Develop a sintering model for UO 2 in MARMOT 2. Expand model to account for irradiation effects 3. Develop a mechanistic macroscale model of densification for BISON

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

High-field magnetic phase diagrams of the 𝑅⁢Mn 6 ⁢Sn 6 (𝑅=Gd–Tm) kagome metals

𝑅⁢Mn 6 ⁢Sn 6 (𝑅=Y, Gd–Lu) kagome metals are promising materials hosting flat electronic bands and Dirac points that interact with magnetism. The coupling between the two magnetic 𝑅 and Mn sublattices can drive complex magnetic states with potential consequences for spin and charge transport and other topological properties. Here, in this work, we use a detailed magnetic Hamiltonian to calculate and predict the magnetic phase diagrams for 𝑅⁢Mn 6 ⁢Sn 6 kagome metals within the mean-field approximation. These calculations reveal a variety of collinear, noncollinear, and noncoplanar phases that arise from competition between various interlayer magnetic exchange interactions and magnetic anisotropies of the 𝑅 and Mn ions. We enumerate these phases and their magnetic space groups for future analysis of their impact on topological and trivial bands near the Fermi surface.

kagome metal↗

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗

Quantification of Dead Lithium on Graphite Anode under Fast Charging Conditions

A series of computational and experimental studies were conducted to understand the onset of lithium plating and subsequent quantification of dead lithium on graphite electrodes in the design of fast charging batteries. The experiments include titration and relaxation studies for detecting initiation of lithium metal plating for various SOC and C-rates, which are compared against the thermodynamically consistent phase field computational results. The collaborative study on "model graphite electrode" with 2.18 mAh cm-2 nominal capacity at 25 degrees C demonstrates: (1) the macroscopic voltage response during relaxation studies indicate the reintercalation of plated lithium into the graphite anode; (2) for SOC below 60% and low C-Rates, there is no dead lithium; (3) for SOC between 60% to 80%, and C-Rates in the range of 4C-6C show dead lithium both in experiments and simulations; (4) at 100% SOC and 4C-6C rates, large amounts of dead lithium are observed. The study presented here allows us to evaluate the effects of the physical properties of the electrochemical system on plating and stripping kinetics and the amount of dead lithium on graphite electrodes, which determines the cell capacity loss under fast charge.

batteries↗

Magnetic field-temperature phase diagram of multiferroic (NH4)2FeCl5·H2O

Abstract Owing to their overall low energy scales, flexible molecular architectures, and ease of chemical substitution, molecule-based multiferroics are extraordinarily responsive to external stimuli and exhibit remarkably rich phase diagrams. Even so, the stability and microscopic properties of various magnetic states in close proximity to quantum critical points are highly under-explored in these materials. Inspired by these opportunities, we combined pulsed-field magnetization, first-principles calculations, and numerical simulations to reveal the magnetic field–temperature ( B – T ) phase diagram of multiferroic (NH 4 ) 2 FeCl 5 ⋅H 2 O. In this system, a network of intermolecular hydrogen and halogen bonds creates a competing set of exchange interactions that generates additional structure in the phase diagram—both in the vicinity of the spin flop and near the 30 T transition to the fully saturated state. Consequently, the phase diagrams of (NH 4 ) 2 FeCl 5 ⋅H 2 O and its deuterated analog are much more complex than those of other molecule-based multiferroics. The entire series of coupled electric and magnetic transitions can be accessed with a powered magnet, opening the door to exploration and control of properties in this and related materials.

36 MATERIALS SCIENCE↗

Multi-scale modeling of the electric field assisted sintering process

The electric field assisted sintering (EFAS) process involves tightly coupled physics that influence microstructural evolution in the particles being compacted. It is also an inherently multi-scale phenomenon, with the microstructure of the compact influencing the subsequent engineering-scale response of the sintering system. To improve understanding of how processing parameters influence microstructural evolution, we have developed a multi-scale modeling approach that couples a continuum-level model of the sintering system with a phase-field model for microstructural evolution of particles within the compact. The phase-field model couples the effect of chemical and electrical driving forces on microstructural evolution and includes the effect of charged defect segregation to surfaces and grain boundaries; this segregation leads to enhanced defect transport and heat generation at these interfaces in response to applied electric field. The effect of enhanced heat generation on particle neck growth and the influence of microstructural evolution on the engineering-scale model are demonstrated.

36 MATERIALS SCIENCE↗

Use of Microgravity to Control the Microstructure of Eutectics

This grant began in June of 1996. Its long term goal is to be able to control the microstructure of directionally solidified eutectic alloys, through an improved understanding of the influence of convection. The primary objective of the present projects is to test hypotheses for the reported influence of microgravity on the microstructure of three fibrous eutectics (MnBi-Bi, InSb-NiSb, Al3Ni-Al). A secondary objective is to determine the influence of convection on the microstructure of other eutectic alloys. Two doctoral students and a masters student supported as a teaching assistant were recruited for this research. Techniques were developed for directional solidification of MnBi-Bi eutectics with periodic application of current pulses to produce an oscillatory freezing rate. Image analysis techniques were developed to obtain the variation in MnBi fiber spacing, which was found to be normally distributed. The mean and standard deviation of fiber spacing were obtained for several freezing conditions. Eighteen ampoules were prepared for use in the gradient freeze furnace QUELD developed at Queen's University for use in microgravity. Nine of these ampoules will be solidified soon at Queen's in a ground-based model. We hope to solidify the other nine in the QUELD that is mounted on the Canadian Microgravity Isolation Mount on MIR. Techniques are being developed for directional solidification of the Al-Si eutectic at different freezing rates, with and without application of accelerated crucible rotation to induce convection. For the first time, theoretical methods are being developed to analyze eutectic solidification with an oscillatory freezing rate. In a classical sharp-interface model, we found that an oscillatory freezing rate increases the deviation of the average interfacial composition from the eutectic, and increases the undercooling of the two phases by different amounts. This would be expected to change the volume fraction solidifying and the fiber spacing. Because of difficulties in tracking the freezing interfaces of the two solid phases, a phase-field model is also being developed. A paper demonstrating application of phase field methods to periodic structures has been submitted for publication.

Wilcox, William R.↗