Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2023

Examining phase separation and crystallization in glasses with X-ray nano-computed tomography

X-ray nano-computed tomography (nano-CT) is a powerful technique to characterize and visualize 3 dimensional (3D) phenomena in complex glasses at the nanoscale. This technique can offer a unique opportunity to explore the intricate morphology of multicomponent glass and glass-ceramic samples, due to its low X-ray energy and high spatial resolution (down to 50 nm). In the current demonstration paper, nano-CT provided insight into crystallization and phase separation, as well as a variety of other phenomena, including corrosion, fracture, and porosity. At scales ranging from 100 to 20 μm, the microstructures of phase separation in borosilicate, complex silicate, and chalcogenide glass compositions were examined. In addition, crystallites present in simulant nuclear waste glasses, volcanic glasses, and other related samples were explored. Nano-CT analysis can add to the understanding of (1) the formation processes, (2) distributions, interactions, and compositions of phases, and (3) the mechanisms, structures, and growth of crystallization and phase separation. Nano-CT offers wide potential in the field of glass science, especially when (1) common characterization methods are insufficient, (2) simple sample preparation is required, (3) 3D rendering is needed, or (4) compelling images of small, complex features are desired.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Pilot-Scale Evaluation of Natural Gas-Based Foam at Elevated Pressure and Temperature Conditions

Abstract Nitrogen (N2) and Carbon Dioxide (CO2) foams have been used as hydraulic fracturing fluids for several decades to reduce water usage and minimize damage in water-sensitive reservoirs. These foam treatments require gases to be liquefied and transported to site. An alternative approach would be to use natural gas (NG) that is readily available from nearby wells, pipelines, and processing facilities as the internal, gaseous phase to create a NG-based foam. Hydraulic fracturing with NG foam is a relatively inexpensive option, makes use of an abundant and often wasted resource, and may even provide production benefits in certain reservoirs. As part of an ongoing development project sponsored by the Department of Energy (DOE), the surface process to create NG foam is being developed and the properties of NG foam are being explored. This paper presents recent results from a rigorous pilot-scale demonstration of NG foam over a range of operating scenarios relevant to surface and bottomhole conditions with a variety of base-fluid mixtures. The Pilot-scale Foam Test Facility (PFTF) used in these investigations is first described. The PFTF is capable of generating foamed fluids at pressures up to 7,500 psig and at temperatures in excess of 300°F. Then, results from several investigations aimed at proving NG foam at conditions relevant to the field are presented. NG foam was characterized using rheology measurements and flow visualization techniques. Experiments were performed to investigate the texture and stability of NG foam generated by two different mixing methods: one using a custom designed tee to match mixing velocities in the field where the gas phase is jetted into the aqueous stream, and another to ensure comprehensive mixing for laboratory analysis. Parametric studies were conducted to explore the effects of flow rate, foam quality, and temperature on the stability of NG foam. Moreover, different fluid preparations were used to investigate the effect of base fluid and additive concentrations on the stability of NG foams. Additional laboratory studies that investigated foam stability with produced water and multicomponent NG mixtures are also reported. The NG foams explored in these investigations exhibited typical, shear-thinning behavior observed in rheological studies of N2- and CO2-based foams. The measured viscosity and observed stability indicate that NG foams are well suited for fracturing applications. Like other foams, NG foam exhibits sensitivity to operating temperature characterized by a decrease in apparent viscosity as temperature increases. Rapid foam breakdown was observed at significantly elevated temperatures exceeding 290°F. In addition to fluid characterization, these investigations also yielded several key lessons that should be applied to future field demonstrations of NG foam.

Beck, Griffin↗
DOE OSTI2020

Recovery of a high-pressure phase formed under laser-driven compression

The recovery of metastable structures formed at high pressure has been a long-standing goal in the field of condensed matter physics. While laser-driven compression has been used as a method to generate novel structures at high pressure, to date no high-pressure phases have been quenched to ambient conditions. Here we demonstrate, using in situ x-ray diffraction and recovery methods, the successful quench of a high-pressure phase which was formed under laser-driven shock compression. We show that tailoring the pressure release path from a shock-compressed state to eliminate sample spall, and therefore excess heating, increases the recovery yield of the high-pressure ω phase of zirconium from 0% to 48%. Our results have important implications for the quenchability of novel phases of matter demonstrated to occur at extreme pressures using nanosecond laser-driven compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Phase diagram of the anisotropic triangular lattice Hubbard model

In a recent study [Phys. Rev. X 10, 021042 (2020)], we showed using large-scale density matrix renormalization group (DMRG) simulations on infinite cylinders that the triangular lattice Hubbard model has a chiral spin liquid phase. In this work, we introduce hopping anisotropy in the model, making one of the three distinct bonds on the lattice stronger or weaker compared with the other two. We implement the anisotropy in two inequivalent ways, one which respects the mirror symmetry of the cylinder and one which breaks this symmetry. In the full range of anisotropy, from the square lattice to weakly coupled one-dimensional chains, we find a variety of phases. Near the isotropic limit we find the three phases identified in our previous work: metal, chiral spin liquid, and 120$^\circ$ spiral order; we note that a recent paper suggests the apparently metallic phase may actually be a Luther-Emery liquid, which would also be in agreement with our results. Furthermore, when one bond is weakened by a relatively small amount, the ground state quickly becomes the square lattice N\'{e}el order. When one bond is strengthened, the story is much less clear, with the phases that we find depending on the orientation of the anisotropy and on the cylinder circumference. While our work is to our knowledge the first DMRG study of the anisotropic triangular lattice Hubbard model, the overall phase diagram we find is broadly consistent with that found previously using other methods, such as variational Monte Carlo and dynamical mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Synthesis method comparison of compositionally complex rare earth-based Ruddlesden–Popper n = 1 T'-type cuprates

The multicomponent approach has been successfully expanded to the Ruddlesden–Popper structure with the synthesis of two different high-entropy cuprate compositions: (La 0.2 Nd 0.2 Gd 0.2 Tb 0.2 Dy 0.2 ) 2 CuO 4 and (La 0.2 Pr 0.2 Nd 0.2 Sm 0.2 Eu 0.2 ) 2 CuO 4 . Additionally, the effect of synthesis method is explored using both solid-state reaction and polymeric steric entrapment (PSE) methods. It is found that PSE leads to more randomly distributed cation species, providing an advantageous method of synthesis for the growing field of high entropy oxides. In situ high-temperature x-ray diffraction tracks the amorphous to crystalline phase transformation in (La 0.2 Nd 0.2 Gd 0.2 Tb 0.2 Dy 0.2 ) 2 CuO 4 powder, synthesized using the PSE method. Using the High-Temperature XRD data, a method for gaining information on the kinetic behavior is also applied. Magnetometry of both compositions indicates ferrimagnetic behavior at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

King City Asbestos Corporation (KCAC) Mine Carbon Mineralization Field Test (Final Combined Report)

The overall objective of the KCAC carbon mineralization project is to determine the most successful approaches to CO 2 mineralization of asbestos-bearing serpentinite rocks and mine waste at the closed KCAC asbestos mine. The objectives of Budget Period 1 (BP1) are to characterize the site including, mineralogical and geochemical variability in site materials, baseline conditions including CO 2 flux and local climate, and to establish the designs for the second, experimental phase of the project. As described in Section 3 of this report, all tasks outlined in the Statement of Project Objective (SOPO) have been completed or are on schedule to be completed by the end of BP1. As described in Section 4, all success criteria and milestones have been met or will be met by the end of BP1. We conclude that the findings of BP1 strongly support moving to the second phase of the project to investigate methods to increase the rate of CO 2 mineralization. CO 2 flux monitoring during BP1 established an estimated baseline uptake of CO 2 over the 142,193 m 2 KCAC mine site of 147 and 179 metric tons CO 2 /yr using dynamic closed chamber and eddy covariance monitoring methods, respectively. Mineralogical analysis shows that the loose surficial waste material contains pyroaurite and minor hydromagnesite, indicating past carbon mineralization. Sediment samples at depths below 60 cm contained ~3 wt% brucite and localized high brucite contents of ~24 wt% were measured at depth, indicating significant unreacted material at relatively shallow depths. Together, these results demonstrate that despite prolonged exposure to the atmosphere, reactive material is still present and accessible and is naturally taking up atmospheric CO 2 . However, these natural rates are still well below rates that can be achieved by applying practices that can accelerate the CO 2 mineralization reaction. Designs for the Field-Testing Phase (BP2) of the project have been established, incorporating the findings of BP1. Based on results to date we have selected two methods to accelerate the rate of CO 2 mineralization, which we refer to as the Tilling Method and the Greenhouse Method. We hypothesize that these methods can increase the CO 2 mineralization rate by up to 5 times greater than the natural background rate measured in BP1. An Initial Lifecyle Assessment (LCA) and techno-economic analysis (TEA) were conducted for these two methods. The LCA concluded that the baseline version of the Greenhouse and Tilling methods at commercial scale are expected to achieve a CDR efficiency of 97% and 71%, respectively (cumulative, over 30 years of operation). Incorporating the net CO 2 removal rates calculated in the LCA, the TEA found the levelized, 30-year costs of the baseline Greenhouse and Tilling methods to be $\$$95/t-CO 2 and $\$$119/t-CO 2 , respectively, and that costs could be brought to well below $\$$100/t-CO 2 for either method by scaling up the volume of rock utilized. An initial EH&S Risk Assessment, evaluating the potential EH&S risks of the proposed carbon mineralization methods at commercial scale, was completed by LLNL EH&S professionals for these two methods.

58 GEOSCIENCES↗
DOE OSTI2021

Non-orthogonal computational grids for studying dislocation motion in phase field approaches

In this work, new non-orthogonal computational grids are implemented into a phase field model called Phase Field Dislocation Dynamics (PFDD). We demonstrate that the new non-orthogonal grid can accommodate multiple slip planes in either the face centered cubic (FCC) or body centered cubic (BCC) crystallographic systems. Here, we show that they avoid numerical errors induced when modeling glide on inclined slip planes in an orthogonal grid. The Gibbs effect that arises in the orthogonal or rotated orthogonal grids is substantially diminished when a non-orthogonal grid is employed. A few test cases demonstrate the effectiveness of using non-orthogonal grids in solving systems with multiple non-planar slip systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced polymorph metastability drives glycine nucleation in aqueous salt solutions

Crystal nucleation from aqueous solutions influences countless geological, biochemical, astrophysical, environmental, and materials science–related phenomena, including ice formation, the manufacturing of active pharmaceutical ingredients, development of diseases such as Alzheimer’s and the origin of life itself. Understanding and controlling nucleation is essential for designing materials with specific properties, developing strategies to inhibit or promote crystallization in various contexts and preventing pathological aggregation in neurodegenerative diseases. Similar to the protein structure prediction problem—where a single amino acid sequence can in theory adopt one most stable conformation but in practice may sample multiple competing conformations—crystal nucleation faces a parallel challenge: the same chemical species can form diverse polymorphs under different environmental conditions (e.g., temperature, pressure, solvent). Each polymorph presents its own set of physical and chemical properties, highlighting the importance of understanding and controlling polymorph selection in fields ranging from pharmaceuticals to materials design. Despite advances in experimental and computational methods for studying phase transitions and polymorph stability, nucleation remains challenging due to its nanoscale nature. Furthermore, in practical settings, salts and impurities can further influence crystal nucleation in diverse contexts, from scaling in pipelines and desalination plants to the durability of concrete and the efficiency of battery materials. This can lead to the formation of polymorphs that may differ from the most stable phase in pure solutions. Or, even though the final structure might appear same irrespective of whether the environment contained impurities or not, the mechanism through which it was formed might be completely different and not intuitive.

Wang, Ruiyu [University of Maryland, College Park,↗
DOE OSTI2021

Z 2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z 2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first order, in contrast to the conventional Z 2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic X X interaction while keeping the local gauge symmetry intact. We find that the Z 2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.

2-dimensional systems↗
DOE OSTI2022

On the S-matrix of Ising field theory in two dimensions

We explore the analytic structure of the non-perturbative S-matrix in arguably the simplest family of massive non-integrable quantum field theories: the Ising field theory (IFT) in two dimensions, which may be viewed as the Ising CFT deformed by its two relevant operators, or equivalently, the scaling limit of the Ising model in a magnetic field. Our strategy is that of collider physics: we employ Hamiltonian truncation method (TFFSA) to extract the scattering phase of the lightest particles in the elastic regime, and combine it with S-matrix bootstrap methods based on unitarity and analyticity assumptions to determine the analytic continuation of the 2 → 2 S-matrix element to the complex s-plane. Focusing primarily on the “high temperature” regime in which the IFT interpolates between that of a weakly coupled massive fermion and the E 8 affine Toda theory, we will numerically determine 3-particle amplitudes, follow the evolution of poles and certain resonances of the S-matrix, and exclude the possibility of unknown wide resonances up to reasonably high energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Covariant phase space with boundaries

The covariant phase space method of Iyer, Lee, Wald, and Zoupas gives an elegant way to understand the Hamiltonian dynamics of Lagrangian field theories without breaking covariance. The original literature however does not systematically treat total derivatives and boundary terms, which has led to some confusion about how exactly to apply the formalism in the presence of boundaries. In particular the original construction of the canonical Hamiltonian relies on the assumed existence of a certain boundary quantity “B”, whose physical interpretation has not been clear. We here give an algorithmic procedure for applying the covariant phase space formalism to field theories with spatial boundaries, from which the term in the Hamiltonian involving B emerges naturally. Our procedure also produces an additional boundary term, which was not present in the original literature and which so far has only appeared implicitly in specific examples, and which is already nonvanishing even in general relativity with sufficiently permissive boundary conditions. The only requirement we impose is that at solutions of the equations of motion the action is stationary modulo future/past boundary terms under arbitrary variations obeying the spatial boundary conditions; from this the symplectic structure and the Hamiltonian for any diffeomorphism that preserves the theory are unambiguously constructed. We show in examples that the Hamiltonian so constructed agrees with previous results. We also show that the Poisson bracket on covariant phase space directly coincides with the Peierls bracket, without any need for non-covariant intermediate steps, and we discuss possible implications for the entropy of dynamical black hole horizons.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A review of the metastable omega phase in beta titanium alloys: the phase transformation mechanisms and its effect on mechanical properties

Since its discovery in 1954, the omega (ω) phase in titanium and its alloys has attracted substantial attention from researchers. The β-to-ω and ω-to-α phase transformations are central to β-titanium alloy design, but the transformation mechanisms have been a subject of debate. With new generations of aberration-corrected transmission electron microscopy and atom probe tomography, both the spatial resolution and compositional sensitivity of phase transformation analysis have been rapidly improving. Furthermore, this review provides a detailed assessment of the new understanding gained and related debates in this field enabled by advanced characterization methods. Specifically, new insights into the possibility of a coupled diffusional-displacive component in the β-to-ω transformation and key nucleation driving forces for the ω-assisted α phase formation are discussed. Additionally, the influence of ω phase on the mechanical properties of β-titanium alloys is also reviewed. Lastly, a perspective on open questions and future direction for research is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Recent advances in describing and driving crystal nucleation using machine learning and artificial intelligence

With the advent of faster computer processors and especially graphics processing units (GPUs) over the last few decades, the use of data-intensive machine learning (ML) and artificial intelligence (AI) has increased greatly, and the study of crystal nucleation has been one of the beneficiaries. In this study, we outline how ML and AI have been applied to address four outstanding difficulties of crystal nucleation: how to discover better reaction coordinates (RCs) for describing accurately non-classical nucleation situations; the development of more accurate force fields for describing the nucleation of multiple polymorphs or phases for a single system; more robust identification methods for determining crystal phases and structures; and as a method to yield improved course-grained models for studying nucleation.

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

Sub-tesla On-Chip Nanomagnetic Metamaterial Platform for Angle-Resolved Photoemission Spectroscopy

Magnetically controlled states in quantum materials are central to their unique electronic and magnetic properties. However, direct momentum-resolved visualization of these states via angle-resolved photoemission spectroscopy (ARPES) has been hindered by the disruptive effect of magnetic fields on photoelectron trajectories. Here, we introduce an in situ method that is, in principle, capable of applying magnetic fields up to 1 T. This method uses substrates composed of nanomagnetic metamaterial arrays with alternating polarity. Such substrates can generate strong, homogeneous, and spatially confined fields applicable to samples with thicknesses up to the micron scale, enabling ARPES measurements under magnetic fields with minimal photoelectron trajectory distortion. Here, we demonstrate this minimal distortion with ARPES data taken on monolayer graphene. Our method paves the way for probing magnetic field-dependent electronic structures and studying field-tunable quantum phases with state-of-the-art energy-momentum resolutions.

Magnetic properties↗
DOE OSTI2023

Statistical Generation of Ocean Forcing With Spatiotemporal Variability for Ice Sheet Models

Melting of ice at the base of floating ice shelves that fringe the Antarctic ice sheet has been identified as a significant source of uncertainty in sea level rise projections. Part of this uncertainty derives from chaotic internal variability of the coupled ocean-atmosphere system. For numerical ice sheet model projections, this uncertainty has not previously been quantified because of the prohibitive computational expense of running large climate model ensembles. Here, we develop and demonstrate a technique that generates independent realizations of internal climate variability from a single climate model simulation. Building on previous developments in model emulation, this technique uses empirical orthogonal function decomposition and Fourier-phase randomization to generate statistically consistent realizations of spatiotemporal variability fields for the target climate variable. The method facilitates efficient sampling of a wide range of climate trajectories, which can also be incorporated within ice sheet or other physical models to represent feedback processes.

54 ENVIRONMENTAL SCIENCES↗