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

Magnetic field-temperature phase diagram of the spin-$\frac{1}{2}$ triangular lattice antiferromagnet KYbSe 2

A quantum spin liquid (QSL) is a state of matter characterized by fractionalized quasiparticle excitations, quantum entanglement, and a lack of long-range magnetic order. However, QSLs have evaded definitive experimental observation. Several Yb 3+ -based triangular lattice antiferromagnets with effective 𝑆 = $\frac{1}{2}$ have been suggested to stabilize the QSL state as the ground state. Here, in this work, we build a comprehensive magnetic temperature phase diagram of a high-quality single crystalline KYbSe 2 via heat capacity and magnetocaloric effect down to 30 mK with magnetic field applied along the 𝑎 axis. At zero magnetic field, we observe the magnetic long-range order at 𝑇 N =0.29 K entering 120 degrees ordered state in heat capacity, consistent with neutron scattering studies. Analysis of the low-temperature (𝑇) specific heat (𝐶) at zero magnetic field indicates linear 𝑇 dependence of 𝐶/𝑇 and a broad hump of 𝐶/𝑇 in the proximate QSL region above 𝑇 N . By applying magnetic field, we observe the up-up-down phase with 1/3 magnetization plateau and oblique phases, in addition to two new phases. These observations strongly indicate that while KYbSe 2 closely exhibits characteristics resembling an ideal triangular lattice, deviations may exist, such as the effect of the next-nearest-neighbor exchange interaction, calling for careful consideration for spin Hamiltonian modeling. Further investigations into tuning parameters, such as chemical pressure, could potentially induce an intriguing QSL phase in the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Reduced-order model for microstructure evolution prediction in the electrodes of solid oxide fuel cell with dynamic discrepancy reduced modeling

Microstructure evolution in the electrodes of solid oxide fuel cell is an important degradation mechanism which reduces active sites for redox reaction and the electric conductivity. Phase field models for microstructure evolution simulation are usually expensive for large scale simulations. In this paper, a reduced-order coarsening model is developed using dynamic discrepancy reduced modeling, which reduces the model order by inserting Gaussian process stochastic functions into the dynamic equations of Ostwald ripening. The reduced order model has been calibrated on a dataset generated by a phase field model that has been well validated to experiments. A validating dataset has also been generated with which the model prediction show good agreement. This model is further applied to predict long term microstructure evolution in different SOFC electrodes. Finally, this work is the first attempt of building a degradation model of SOFC using data science techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Coupled Lattice Polarization and Ferromagnetism in Multiferroic NiTiO 3 Thin Films

Polarization-induced weak ferromagnetism (WFM) was demonstrated a few years back in LiNbO3-type compounds, MTiO3 (M = Fe, Mn, Ni). Although the coexistence of ferroelectric polarization and ferromagnetism has been demonstrated in this rare multiferroic family before, first in bulk FeTiO3, then in thin-film NiTiO3, the coupling of the two order parameters has not been confirmed Here, we report the stabilization of polar, ferromagnetic NiTiO3 by oxide epitaxy on a LiNbO3 substrate utilizing tensile strain and demonstrate the theoretically predicted coupling between its polarization and ferromagnetism by X-ray magnetic circular dichroism under applied fields. The experimentally observed direction of ferroic ordering in the film is supported by simulations using the phase-field approach. Our work validates symmetry-based criteria and first-principles calculations of the coexistence of ferroelectricity and WFM in MTiO3 transition metal titanates crystallizing in the LiNbO3 structure. It also demonstrates the applicability of epitaxial strain as a viable alternative to high-pressure crystal growth to stabilize metastable materials and a valuable tuning parameter to simultaneously control two ferroic order parameters to create a multiferroic. Multiferroic NiTiO3 has potential applications in spintronics where ferroic switching is used, such as new four-stage memories and electromagnetic switches.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Mesoscale model for fission-induced recrystallization in U-7Mo alloy

A mesoscale model is developed by integrating the rate theory and phase-field models and is used to study the fission-induced recrystallization in U-7Mo alloy. The rate theory model is used to predict the dislocation density and the recrystallization nuclei density due to irradiation. The predicted fission rate and temperature dependences of the dislocation density are in good agreement with experimental measurements. This information is used as input for the multiphase phase-field model to investigate the fission-induced recrystallization kinetics. The simulated recrystallization volume fraction and bubble induced swelling agree well with experimental data. The effects of the fission rate, initial grain size, and grain morphology on the recrystallization kinetics are discussed based on an analysis of recrystallization growth rate using the modified Avrami equation. Here, we conclude that the initial microstructure of the U-Mo fuels, especially the grain size, can be used to effectively control the rate of fission-induced recrystallization and therefore swelling.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Effect of patterned inclusions on the fracture behavior of ceramic composites

Developments in additive manufacturing have enabled the reinforcement of ceramic composites across the scales. Understanding the fracture behavior of such composites is critical for structural design and engineering. However, systematic investigations on the relationship between the patterning characteristics and the effective fracture properties have been scarce. In this work, we propose a mixed-mode phase field model to investigate the effects of elastic modulus and pattern characteristics of circular inclusions on the fracture behavior of a ceramic composite beam. The modeling results have shown that the stiff and compliant inclusions can either repel or attract the crack. The extent of repelling and attraction can be controlled by altering the pattern characteristics. In addition, utilizing the relationship between the pattern characteristics and resulting fracture path, fracture growth can be controlled with appropriate inclusion properties. The targeted enhancement in fracture toughness can also be predicted using an effective fracture model linked with the underlying fracture path computed by the proposed mixed-mode phase field model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unraveling the Determinant Mechanisms in Flow-Mediated Crystal Growth and Phase Behaviors

To uncover the critical mechanisms responsible for mesoscopic level development during flow-mediated crystal growth, we develop a semi-two-way hydrodynamic coupled structural phase-field crystal formalism (HXPFC-s2). The new formalism, inspired by previous attempts at coupling hydrodynamic and phase-field crystal (PFC) equations, allows for studying mesoscopic flow-mediated crystallization at diffusive timescales pertinent to industrial applications. Unlike previous efforts, the devised coupling to the structural PFC (XPFC) equations allows generalization to more complex crystal structures through explicit parameterization of the direct correlation function (DCF). Utilizing the HXPFC-s2 formalism, we seek to uncover the determinant physical mechanisms in crystallization under simple shear flows by comparing temperature-driven crystallization to flow-mediated crystallization under varying flow-strengths. Parallels and deviations of under-cooling and flow-strength effects on crystal growth are drawn using the crystal cluster-size and system ordering time evolutions. In doing so, we identify scaling behaviors with a Peclet-like number, Pe∼, a critical Peclet-like number, Pe∼*, and flow-field-crystal plane-dependent interactions. Our findings may be relevant for controlling crystal growth and phase behaviors in flow applications.

Willis, L. Connor (ORCID:0009000961321848)↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2024

Magnetic field induced phases and spin Hamiltonian in Cs 2 CoBr 4

Magnetic structures and spin excitations are studied across the phase diagram of the geometrically frustrated S = 3/2 quantum antiferromagnet Cs 2 CoBr 4 in magnetic fields applied along the magnetic easy axis, using neutron diffraction, inelastic neutron scattering, and terahertz absorption spectroscopy. Here, the data are analyzed, where appropriate, using extended SU(4) linear spin wave theory. A minimal magnetic Hamiltonian is proposed based on measurements in the high-field paramagnetic state. It deviates considerably from the previously considered models. Additional dilatometry experiments highlight the importance of magnetoelastic coupling in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Microstructure Experiments-Enabled MARMOT Simulations of SiC/SiC-based Accident Tolerant Nuclear Fuel System

We have undertaken an experimental-computational project that addresses a few key technology gaps associated with the use of SiC/SiC composites for light water reactor fuel cladding. Two principal endeavors in this project are to assess the irradiation-induced microstructural changes and swelling in SiC/SiC composites, and to characterize and model the porous oxide surface layer (when exposed to steam) along SiC recession that is detrimental to the clad integrity under accident conditions. The project tasks are: (i) ion irradiation and characterization (University of Tennessee, Knoxville), (ii) mechanical tests and analysis (University of South Carolina), (iii) steam exposure tests (Oak Ridge National Laboratory), (iv) electron microscopy and spectroscopic characterization (NC State University), (v) x-ray microscopy and reconstruction (University of South Carolina and NC State University), (vi) phase field modeling and simulations (NC State University and Idaho National Laboratory). The effects of 10 MeV Au ion irradiation at 350°C on the microstructure evolution in SiC/SiC composites are investigated at doses up to 400 displacements per atom (dpa) at the University of Tennessee, Knoxville. Atomic force microscopy and optical profilometry reveal irradiation induced axial and radial shrinkage of the fibers for doses greater than 10 dpa. Based on detailed electron microscopy characterization, the primary cause of the fiber shrinkage is attributed to irradiation-induced loss of carbon packets. Additionally, the multilayer PyC interface is observed to portray high resistance to irradiation damage. The mechanical response and failure mechanisms of un-irradiated samples is also assessed through loading tests and X-ray imaging at the University of South Carolina. Steam exposure tests are performed at the Oak Ridge National Laboratory in a facility that represents a reactor pressure vessel under a loss-of-coolant-accident scenario. The samples that are analyzed methodically in this report are tested for 32/31 hours at 1200°C with a velocity of 0.25 cm/s for the pressures: 0.1 MPa, 0.45 MPa, 0.92 MPa and 1.38 MPa. Scanning/transmission electron microscopy analysis conducted at the NC State University (NCSU) shows that the oxide layer thickness increases with the steam pressure. While the oxide layer is crystalline (α- cristobalite) for the pressures 0.45 MPa, 0.92 MPa and 1.38 MPa, the layer is amorphous at 0.1 MPa. Results from Raman spectroscopy have also confirmed the formation of α-cristobalite phase of SiO₂. The abrupt increase in the integrated Raman intensity ratio between steam pressures 0.1 and 0.45 MPa suggests the onset of accelerated crystallization. Non-destructive three dimensional X-ray microscopy/tomography (XCT) and computational image processing techniques are employed by the University of South Carolina to probe the porous oxidation features on SiC samples at varying pressures. Interestingly, most pores are observed to be located away from the surface as well as the oxide-SiC interface. The average oxide layer thickness assessed from the XCT analysis is seen to be in excellent agreement with the values determined through scanning electron microscopy analysis for the highest pressures where the oxide layer is relatively more uniform. A phase-field model developed by the NC State University and Idaho National Laboratory for simulating oxidation of SiC by steam captures the paralinear kinetics of SiC oxidation with a high degree of fidelity. Results from quasi-one-dimensional and two dimensional simulations show that the pores with oxidizing species lead to a higher volatilization rate. These simulations indicate that the enhanced apparent volatilization rates at higher pressures observed in experiments can be rationalized by the increased volatilization from the pores. The project team has also successfully developed the capability to import images from experiments for realistic evolution of the oxidizing microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

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.

17 WIND ENERGY↗
DOE OSTI2015

LESSONS FROM NUCLEAR ENERGY: BUILDING SNOW AND AVALANCHE SIMULATION TOOLS USING A COMMON OPEN-SOURCE FRAMEWORK

The Department of Energy (DOE) is currently investing millions of dollars annually into various modeling and simulation tools for all aspects of nuclear energy. An important part of this effort includes developing applications based on the open-source Multiphysics Object Oriented Simulation Environment (MOOSE; mooseframework.org) from Idaho National Laboratory (INL). Through the DOE and outside collaborators, MOOSE includes a large set of physics modules, including phase field, level set, heat conduction, tensor mechanics, Navier-Stokes, and porous media, among others. The phase-field module, in particular, is well-suited for micro-structure evolution simulations, including solidification, as shown in the figure. It was also used, in conjunction with the heat conduction module, to build a multiscale snow model (Pika) as presented at the 2014 International Snow Science Workshop (Slaughter et al., 2014). The snow science community should learn from the nuclear industry, and harness the enormous effort underway to build simulation tools that are open, modular, and share a common framework. In particular, MOOSE-based multi-physics solvers are inherently parallel, dimension agnostic, adaptive in time and space, fully-coupled, and capable of interacting with other applications. The snow science community must build on existing tools to enable collaboration between researchers and practitioners throughout the world, and advance the state-of-the-art in line with other scientific research efforts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Field induced phase transition in Sm{sub 0.5}(Ca{sub 1−x}Sr{sub x}){sub 0.5}MnO{sub 3}

We demonstrate a fascinating consequence of A-site disorder at x = 0.5 in a new series Sm{sub 0.5}(Ca{sub 1−x}Sr{sub x}){sub 0.5}MnO{sub 3}. An unusual magnetic field-induced ultra-sharp transition to ferromagnetic metallic (FMM) state and collapse of charge ordering are observed both in magnetization and magnetoresistance curves at much lower field compared to end compounds. When magnetic field is applied above a critical field, FMM state appears and it retains the high-field state, although magnetic field decreases to zero. With decreasing field sweep rate the transition shifts to higher field indicating a meta-magnetic transition. The atomic-scale local inhomogeneity/distortion arising from the difference in ionic radii and/or the random Coulomb potential from ion mixture is supposed to be critical behind such striking phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Visualization of dielectric constant-electric field-temperature phase maps for imprinted relaxor ferroelectric thin films

The dielectric phase transition behavior of imprinted lead magnesium niobate–lead titanate relaxor ferroelectric thin films was mapped as a function of temperature and dc bias. To compensate for the presence of internal fields, an external electric bias was applied while measuring dielectric responses. The constructed three-dimensional dielectric maps provide insight into the dielectric behaviors of relaxor ferroelectric films as well as the temperature stability of the imprint. The transition temperature and diffuseness of the dielectric response correlate with crystallographic disorder resulting from strain and defects in the films grown on strontium titanate and silicon substrates; the latter was shown to induce a greater degree of disorder in the film as well as a dielectric response lower in magnitude and more diffuse in nature over the same temperature region. Strong and stable imprint was exhibited in both films and can be utilized to enhance the operational stability of piezoelectric devices through domain self-poling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The Contributions of Polar Nanoregions to the Dielectric and Piezoelectric Responses in Domain-Engineered Relaxor-PbTiO 3 Crystals

The existence of polar nanoregions is the most important characteristic of relaxor-based ferroelectric materials. Recently, the contributions of polar nanoregions to the shear piezoelectric property of relaxor-PbTiO 3 (PT) crystals are confirmed in a single domain state, accounting for 50%–80% of room temperature values. For electromechanical applications, however, the outstanding longitudinal piezoelectricity in domain-engineered relaxor-PT crystals is of the most significance. In this paper, the contributions of polar nanoregions to the longitudinal properties in [001]-poled Pb(Mg 1/3 Nb 2/3 )O 3 -0.30PbTiO 3 and [110]-poled Pb(Zn 1/3 Nb 2/3 )O3-0.15PbTiO 3 (PZN-0.15PT) domain-engineered crystals are studied. Taking the [110]-poled tetragonal PZN-0.15PT crystal as an example, phase-field simulations of the domain structures and the longitudinal dielectric/piezoelectric responses are performed. According to the experimental results and phase-field simulations, the contributions of polar nanoregions (PNRs) to the longitudinal properties of relaxor-PT crystals are successfully explained on the mesoscale, where the PNRs behave as “seeds” to facilitate macroscopic polarization rotation and enhance electric-field-induced strain. The results reveal the importance of local structures to the macroscopic properties, where a modest structural variation on the nanoscale greatly impacts the macroscopic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Unique Helical Magnetic Order and Field-Induced Phase in Trillium Lattice Antiferromagnet EuPtSi

Magnetic transition phenomena in cubic chiral antiferromagnet EuPtSi with T N = 4.0 K were investigated by means of single crystal neutron diffraction. At 0.3 K in the ground state, magnetic peaks emerge at positions represented by an ordering vector q 1 =(0.2,0.3,0) and its cyclic permutation. Upon heating, an additional magnetic peak splitting with hysteresis was uncovered at around T$*\atop{N}$~2.5 K, indicating the presence of a first-order commensurate–incommensurate transition with q$*\atop{1}$=(0.2,0.3,δ) (δ max ≃ 0.04) at T$*\atop{N}$. A half-polarized neutron scattering experiment for polarization parallel to the scattering vector revealed that polarization antiparallel to the scattering vector has stronger intensity in both magnetic phases. This feature clarifies the single chiral character of the helical structure with moments lying perpendicular to the ordering vector in both ordered states. Under a vertical magnetic field of 1.2 T for B ∥[1,1,1] at 1.9 K entering into the so-called A phase, magnetic peaks form characteristic hexagonal patterns in the equatorial scattering plane around nuclear peaks. An ordering vector q A ≃(±0.09,±0.20,∓0.28) of the A-phase has similar periodic length as q 1 , and could be the hallmark of a formation of skyrmion lattice in EuPtSi.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Phase Behavior of Grafted Polymer Nanocomposites from Field-Based Simulations

There are limited theoretically predicted phase diagrams for polymer nanocomposites (PNCs) because conventional modeling techniques are largely unable to predict the macroscale phase behavior of PNCs. Here, we show that recent field-based methods, including PNC field theory (PNC-FT) and theoretically informed Langevin dynamics, can be used to calculate phase diagrams for polymer-grafted nanoparticles (gNPs) incorporated into a polymer matrix. We calculate binodals for the transition from the miscible, dispersed phase to the macrophase separated state as functions of important experimental parameters, including the ratio of the matrix chain degree of polymerization ( P ) to the grafted chain degree of polymerization ( N ), the enthalpic repulsion between the matrix and grafted chains, the grafting density (σ), the size of the NPs, and the NP volume fraction. We demonstrate that thermal and polymer conformational fluctuations enhance the degree of phase separation in gNP-PNCs, a result of depletion interaction effects. We support this conclusion by experimentally investigating the phase separation of poly(methyl methacrylate)-grafted silica NPs in a polystyrene matrix as a function of P/N . The simulations only agree with experiments when fluctuations are included because fluctuations are needed to properly capture the depletion interactions between the gNPs. Here, we clarify the role of conformational entropy in driving depletion interactions in PNCs and suggest that inconsistencies in the literature may be due to polymer chain length effects since conformational entropy increases with increasing chain length.

36 MATERIALS SCIENCE↗
DOE OSTI2018

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

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

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

36 MATERIALS SCIENCE↗