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

Anisotropic Cahn-Hilliard free energy and interfacial energies for binary alloys with pairwise interactions

The original Cahn-Hilliard derivation of the contribution of compositional inhomogeneity to the free energy of a binary alloy with pairwise interactions is extended to include higher-order inhomogeneity terms. For alloys on a cubic lattice, the coefficient of the first inhomogeneity is a second-rank tensor and reduces to a scalar, but it is shown that the second order and the third order inhomogeneity terms are weighted by fourth-rank and sixth-rank tensors, thus resulting in anisotropic contributions. Furthermore, each interaction shell generates a unique set of inhomogeneity coefficients that is determined by the set of vectors connecting an atom to its neighbors on that shell. These coefficients are calculated for fcc and bcc alloys with interactions up fourth nearest neighbors. Phase field simulations based on these extended Cahn-Hilliard free energies are performed to measure interface free energies along specific crystallographic directions as a function of temperature, and to obtain the equilibrium shape of precipitates. Interface free energies, and the resulting anisotropies, are compared to those obtained by discrete models and Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Flux-Closure Domain Structures in Ferroelectric K 0.5 Na 0.5 NbO 3 Thin Films

Topological domain structures in ferroelectric materials have garnered increasing attention due to their intriguing physical properties and promising applications. While most existing topological structures in ferroelectric perovskite oxides originate from tetragonal or rhombohedral bulk phases, much less is understood about their counterparts in orthorhombic ferroelectrics. Here, in this work, we employ ferroelectric K 0.5 Na 0.5 NbO 3 (KNN) thin films as a model system and leverage phase-field simulations to theoretically predict the static structures and dynamic behaviors of three types of flux-closure domain configurations: in-plane (Type-I), out-of-plane (Type-II), and superdomain (Type-III) flux-closure structures. We systematically investigate the effects of finite size, misfit strains, and electrical boundary conditions on the formation and switching of these topological structures. For the Type-I structure, size reduction or small misfit strain facilitates a transition of the flux-closure pattern to polar vortices. Type-II structures emerge under open-circuit electrical boundary conditions of the film, forming at the junctions of specific domain walls with the film surface or the film–substrate interface. The formation mechanisms of these two flux-closure structures are rationalized from an energy perspective. We demonstrated switching capabilities of Type-II and Type-III structures by obtaining polarization–electric field hysteresis loops. Our simulations also reveal a reversible electric-field-induced transition between the orthorhombic and rhombohedral ferroelectric phases with a checkerboard domain pattern, during which the integrity of the flux-closure structure is preserved. These findings provide theoretical insights and practical guidance for identifying and manipulating topological structures in low-symmetry ferroelectrics, paving the way for developing energy-efficient microelectronic devices based on topological structures.

P-E loop↗
DOE OSTI2021

C-FOG: Life of Coastal Fog

C-FOG is a comprehensive bi-national project dealing with the formation, persistence, and dissipation (life cycle) of fog in coastal areas (coastal fog) controlled by land, marine, and atmospheric processes. Given its inherent complexity, coastal-fog literature has mainly focused on case studies, and there is a continuing need for research that integrates across processes (e.g., air–sea–land interactions, environmental flow, aerosol transport, and chemistry), dynamics (two-phase flow and turbulence), microphysics (nucleation, droplet characterization), and thermodynamics (heat transfer and phase changes) through field observations and modeling. Central to C-FOG was a field campaign in eastern Canada from 1 September to 8 October 2018, covering four land sites in Newfoundland and Nova Scotia and an adjacent coastal strip transected by the Research Vessel Hugh R. Sharp . An array of in situ, path-integrating, and remote sensing instruments gathered data across a swath of space–time scales relevant to fog life cycle. Satellite and reanalysis products, routine meteorological observations, numerical weather prediction model (WRF and COAMPS) outputs, large-eddy simulations, and phenomenological modeling underpin the interpretation of field observations in a multiscale and multiplatform framework that helps identify and remedy numerical model deficiencies. An overview of the C-FOG field campaign and some preliminary analysis/findings are presented in this paper.

58 GEOSCIENCES↗
DOE OSTI2023

Phase Field Sintering Simulations of Tagged UO 2

Isotopic taggants are being considered as an additive for UO 2 manufacturing to improve nuclear security. Before taggants can be adopted, it must be determined whether taggants have any effect on the fuel microstructure or performance. As part of this effort, a grand potential formulation of the phase field method was used to model untagged, Cr-tagged, and Fe-tagged UO 2 . The simulations were compared to determine whether the taggants had any effect on fuel density or grain size. The results suggest that the taggants will not affect the density or grain size, but these findings may not be reliable because of high uncertainty. Causes of the uncertainty were discussed, and future work is proposed to reduce the model uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2022

Engineering infinite-range SU($\textit{n}$) interactions with spin-orbit-coupled fermions in an optical lattice

Here, we study multilevel fermions in an optical lattice described by the Hubbard model with on-site SU($\textit{n}$) -symmetric interactions. We show that in an appropriate parameter regime this system can be mapped onto a spin model with all-to-all SU($\textit{n}$)-symmetric couplings. Raman pulses that address internal spin states modify the atomic dispersion relation and induce spin-orbit coupling, which can act as a synthetic inhomogeneous magnetic field that competes with the SU($\textit{n}$) exchange interactions. We investigate the mean-field dynamical phase diagram of the resulting model as a function of $\textit{n}$ and different initial configurations that are accessible with Raman pulses. Consistent with previous studies for $\textit{n}$ = 2 , we find that for some initial states the spin model exhibits two distinct dynamical phases that obey simple scaling relations with $\textit{n}$. Moreover, for $\textit{n}$ > 2 we find that dynamical behavior can be highly sensitive to initial intraspin coherences. Our predictions are readily testable in current experiments with ultracold alkaline-earth-metal(-like) atoms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Primordial black holes and magnetic fields in conformal neutrino mass models

Sufficiently strong and long-lasting first-order phase transitions can produce primordial black holes (PBHs) that contribute substantially to the dark matter abundance of the Universe, and can produce large-scale primordial magnetic fields. Here, we study these mechanisms in a generic class of conformal U(1)' models that also explain active neutrino oscillation data via the type-I seesaw mechanism. We find that phase transitions that occur at seesaw scales between 10 4 GeV and 10 11 GeV produce gravitational wave signals (from the dynamics of the phase transition and from the decay of cosmic string loops) at LISA/ET that can be correlated with microlensing signals of PBHs at the Roman Space Telescope, while scales near 10 11 GeV can be correlated with Hawking evaporation signals at future gamma-ray telescopes. LISA can probe the entire range of PBH masses between 1 × 10 -16 M ⊙ and 8 × 10 -11 M ⊙ if PBHs fully account for the dark matter abundance. For Z' masses between 40 TeV and 10 4 TeV, and 10 TeV right-handed neutrinos, helical magnetic fields can be produced with magnitudes ≳ 0.5 pG and coherence lengths ≳ 0.008 Mpc, above current blazar lower bounds.

cosmological phase transitions↗
DOE OSTI2021

Topological Mott transition in a two band model of spinless fermions with on-site Coulomb repulsion

Highlights: • The effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling lead to the 2D topological semimetal. • The ground state phase diagram is calculated. • The topological semimetal phase is characterized by the zero energy Majorana states. In the framework of mean field approach, we study topological Mott transition in a two band model of spinless fermions on a square lattice at half filling. We consider the combined effect of the on-site Coulomb repulsion and the spin–orbit Rashba coupling. The ground state phase diagram is calculated as a function of the strength of the spin–orbit Rashba coupling and the Coulomb repulsion. The spin–orbit Rashba coupling leads to a distinct phase of matter, the topological semimetal. We study a new type of phase transition between the non-topological insulator and topological semimetal states. Topological phase state is characterized by the zero energy Majorana states, which there are in defined region of the wave vectors and are localized at the boundaries of the sample. The region of existence of the zero energy Majorana states tends to zero at the point of the Mott phase transition. The zero energy Majorana states are dispersionless (they can be considered as flat bands), the Chern number and Hall conductance are equal to zero (note in two dimensional model).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Theory of the Kitaev model in a [111] magnetic field

Recent numerical studies indicate that the antiferromagnetic Kitaev honeycomb lattice model undergoes a magnetic-field-induced quantum phase transition into a new spin-liquid phase. This intermediate-field phase has been previously characterized as a gapless spin liquid. By implementing a recently developed variational approach based on the exact fractionalized excitations of the zero-field model, we demonstrate that the field-induced spin liquid is gapped and belongs to Kitaev’s 16-fold way. Specifically, the low-field non-Abelian liquid with Chern number C = ±1 transitions into an Abelian liquid with C = ±4. The critical field and the field-dependent behaviors of key physical quantities are in good quantitative agreement with published numerical results. Furthermore, we derive an effective field theory for the field-induced critical point which readily explains the ostensibly gapless nature of the intermediate-field spin liquid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Grand potential sintering simulations of doped UO 2 accident-tolerant fuel concepts

Chromium-doped UO is a widely-studied near-term deployable accident-tolerant fuel concept because it results in a dense, large-grain structure that increases the fuel resistance to densification, swelling, and fission gas release. A new charged-interstitial mechanism was recently proposed to describe the behavior of dopants like chromium in sintered UO 2 . Based on that mechanism, manganese was suggested as an even stronger dopant than chromium. In the current work we use mesoscale sintering simulations in an effort to validate the new mechanism. We compare the relative behavior of Cr-doped and undoped UO 2 against experimental data in the literature. We also make predictions of the relative behavior of Mn-doped UO 2 . Simulations are done using the phase field-based grand potential sintering model. Dopants have two effects on sintered UO 2 . They marginally increase the densification rate and greatly increase the average grain size. Both of these effects are individually tested in small-scale simulations. Then large 3D sintering simulations are used to test the combined behavior of both effects. The results for the densification rate simulations are consistent with experiments for Cr-doped fuels. However, we note the grain growth rates are lower than what is found in the literature for Cr-doped fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Modeling impacts of ice-nucleating particles from marine aerosols on mixed-phase orographic clouds during 2015 ACAPEX field campaign

Abstract. A large fraction of annual precipitation over the western United States comes from wintertime orographic clouds associated with atmospheric rivers (ARs). Transported African and Asian dust and marine aerosols from the Pacific Ocean may act as ice-nucleating particles (INPs) to affect cloud and precipitation properties over the region. Here we explored the effects of INPs from marine aerosols on orographic mixed-phase clouds and precipitation at different AR stages for an AR event observed during the 2015 ACAPEX field campaign under low dust (<0.02 cm−3) conditions. Simulations were conducted using the chemistry version of the Weather Research and Forecasting Model coupled with the spectral-bin microphysics at 1 km grid spacing, with ice nucleation connected with dust and marine aerosols. By comparing against airborne and ground-based observations, accounting for marine INP effects improves the simulation of AR-precipitation. The marine INPs enhance the formation of ice and snow, leading to less shallow warm clouds but more mixed-phase and deep clouds, as well as to a large spillover effect of precipitation after AR landfall. The responses of cloud and precipitation to marine INPs vary with the AR stages, with more significant effects before AR landfall and post-AR than after AR landfall, mainly because the moisture and temperature conditions change with the AR evolution. This work suggests weather and climate models need to consider the impacts of marine INPs since their contribution is notable under low dust conditions despite the much lower relative ice nucleation efficiency of marine INPs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Coupled magnetostructural continuum model for multiferroic BiFeO 3

Here, a continuum approach to study magnetoelectric multiferroic BiFeO 3 (BFO) is proposed. Our modeling effort marries the ferroelectric (FE) phase field method and micromagnetic simulations to describe the entire multiferroic order parameter sector (polarization, oxygen antiphase tilts, strain, and magnetism) self-consistently on the same time- and length scale. In this paper, we discuss our choice of FE and magnetic energy terms and demonstrate benchmarks against known behavior. We parametrize the lowest order couplings of the structural distortions against previous predictions from density functional theory calculations giving access to simulations of the FE domain wall (DW) topology. This allows us to estimate the energetic hierarchy and thicknesses of the numerous structural DWs. We then extend the model to the canted antiferromagnetic order and demonstrate how the FE domain boundaries influence the resulting magnetic DWs. We also highlight some capabilities of this model by providing two examples relevant for applications. We demonstrate spin-wave transmission through the multiferroic domain boundaries which identify rectification in qualitative agreement with recent experimental observations. As a second example of application, we model fully dynamical magnetoelectric switching, where we find a sensitivity on the Gilbert damping with respect to switching pathways. We envision that this modeling effort will set the basis for further work on properties of arbitrary 3D nanostructures of BFO (and related multiferroics) at the mesoscale.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

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↗