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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Validity of Machine Learning in the Quantitative Analysis of Complex Scanning Near-Field Optical Microscopy Signals Using Simulated Data

Scattering-type scanning near-field optical microscope (s-SNOM) is a modern technique for subdiffractional optical imaging and spectroscopy. Over the past two decades, tremendous efforts have been devoted to modeling complex tip-sample interactions in s-SNOM, aimed at understanding the electrodynamics of materials at the nanoscale. However, due to complexities in analytical methods and the limited computation power for fully numerical simulations, compromises must be made to facilitate the modeling of tip-sample interaction, such as using quasistatic approximation or unrealistic tip geometries. Here, we apply a variety of widely utilized machine-learning methods, including k nearest neighbor and feedforward neural network etc. to study the phase-resolved spectroscopic near-field response. With only a small set of training data, which is simulated using the finite-dipole model, we demonstrate that the relation between the experimental near-field signal and sample optical constant can be one to one mapped without the need for tip modeling: for a given material with a moderate dielectric function, its complex near-field spectrum can be accurately determined within the mid-IR spectral range, and vice versa. Our preliminary study sets the stage for future exploration using real experimental data. Our method is beneficial for processing the increasing amount of data accumulated across many research groups and especially useful for user facilities such as synchrotron-based national laboratories where a large amount of data is generated on a daily basis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Orbital-selective Mott phase and non-Fermi liquid in FePS 3

The layered metal phosphorous trisulfide FePS 3 is reported to be a Mott insulator at ambient conditions and to undergo structural and insulator-metal phase transitions under pressure. However, the character of the resulting metallic states has not been understood clearly so far. Here, we theoretically study the phase transitions of FePS 3 using first-principles methods based on density functional theory and embedded dynamical mean field theory. We find that the Mott transition in FePS 3 can be orbital selective, with t 2g states undergoing a correlation-induced insulator-to-metal transition while e g states remain gapped. We show that this orbital-selective Mott phase, which occurs only when nonhydrostatic pressure is used, is a bad metal (or non-Fermi liquid) with large fluctuating moments due to Hund's coupling. Further application of pressure increases the crystal-field splitting and converts the system to a conventional Fermi liquid with low-spin configurations dominant. Finally, our results show that FePS 3 is an example of a system that realizes an orbital-selective Mott phase, allowing tuning between correlated and uncorrelated metallic properties in an accessible pressure range (≤ 18 GPa).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Deformation heterogeneity and intragrain lattice misorientation in high strength contrast, dual-phase bridgmanite/periclase

A bridgmanite/periclase aggregate is investigated due to its prevalence in the Earth’s lower mantle and its importance for understanding geodynamic processes. In dual-phase polycrystalline aggregates, both strength contrast between phases and single crystal elastic and plastic anisotropy are known to influence the development of intragrain misorientation, and the evolution of crystallographic texture. In this study, full-field crystal plasticity simulations are performed with a finite element solution method, and applied to an aggregate of a mechanically strong orthorhombic phase with perovskite structure (bridgmanite, MgSiO 3 ), and a relatively weaker cubic phase (periclase, MgO). The relative strengths of the phases are parameterized to elucidate the effect that strength contrast has on texture evolution and single-phase simulations are performed for comparison. Altogether, results indicate that the relative strength between the two phases influences the development of plasticity, and the overall texture evolution. Results are discussed in light of trends related to the evolution of plasticity and misorientation in the aggregate, and their dependence on both the strength contrast between phases as well as the spatial distribution of grains and phases.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Insights into distorted lamellar phases with small-angle scattering and machine learning

Lamellar phases are essential in various soft matter systems, with topological defects significantly influencing their mechanical properties. In this report, we present a machine-learning approach for quantitatively analyzing the structure and dynamics of distorted lamellar phases using scattering techniques. By leveraging the mathematical framework of Kolmogorov–Arnold networks, we demonstrate that the conformations of these distorted phases – expressed as superpositions of complex waves – can be reconstructed from small-angle scattering intensities. Through the contour analysis of wave field phase singularities, we obtain the statistics of the spatial distribution of topological defects. Furthermore, we establish that the temporal evolution of these defects can be derived from the time-dependent traveling wave field, informed by the dispersion relation of spectral components. This method opens new avenues for investigating the dynamics of distorted lamellar phases using various dynamic scattering techniques such as neutron spin echo and X-ray photon correlation spectroscopy. These findings enhance our microscopic understanding of how defects influence the physical properties of lamellar materials, with implications for both equilibrium and non-equilibrium states in general lamellar systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Revealing the dynamics of ultrarelativistic non-equilibrium many-electron systems with phase space tomography

The description of physical processes with many-particle systems is a key approach to the modeling of numerous physical systems. For example in storage rings, where ultrarelativistic particles are agglomerated in dense bunches, the modeling and measurement of their phase-space distribution is of paramount importance: at any time the phase-space distribution not only determines the complete space-time evolution but also provides fundamental performance characteristics for storage ring operation. Here, we demonstrate a non-destructive tomographic imaging technique for the 2D longitudinal phase-space distribution of ultrarelativistic electron bunches. For this purpose, we utilize a unique setup, which streams turn-by-turn near-field measurements of bunch profiles at MHz repetition rates. To demonstrate the feasibility of our method, we induce a non-equilibrium state and show that the phase-space distribution microstructuring as well as the phase-space distribution dynamics can be observed in great detail. Our approach offers a pathway to control ultrashort bunches and supports, as one example, the development of compact accelerators with low energy footprints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

An eigenvalue-based method for computing the relaxed pressure in compressible multiphase flow with N phases

The modeling of compressible multiphase flows is a decades-old area of study with many applications across various fields. Many of these application areas use stiff pressure relaxation. This process involves the solution of a nonlinear system with N + 1 equations and N + 1 unknowns, where N is the number of phases. The resolution of this system with general equations of state (EOSs) is difficult. Furthermore, nonlinear systems can admit multiple solutions, and current solution methods do not address this possibility. Very recently, a thermodynamic relaxation method was introduced, which effectively maps a relatively simple predictor equation of state onto a more complex target equation of state. In this context, the target EOSs are the chosen EOSs for the thermodynamic model. Furthermore, this thermodynamic relaxation has the benefit of simplifying the stiff pressure relaxation system of equations. In this article, we show this system reduces to a polynomial of degree N, which can be recast as an eigenvalue problem through the use of the associated companion matrix. We show that although this eigenvalue method is generally less efficient than Newton–Raphson iteration, it does not suffer from convergence issues and finds all N roots of the polynomial. Hence, the method provides a fail-safe for root-finding iterative methods and a way to address the issue of multiple solutions to the nonlinear system of equations in stiff pressure relaxation.

Eigenvalue algorithm↗
DOE OSTI2026

Enhanced coercivity in Fe5C2/SiO2 core/shell nanocrystals

Rod-shaped Fe5C2 and core/shell Fe5C2/SiO2 nanocrystals were synthesized via a solution-based chemical method. Structural analysis confirmed the monoclinic phase of Fe5C2 with space group C2/c. Zero-field-cooling (ZFC) and field-cooling (FC) magnetization curves revealed distinct magnetic behaviors: uncoated Fe5C2 exhibited a low-temperature FC plateau indicative of strong dipolar interactions, while Fe5C2/SiO2 showed a monotonic increase in FC magnetization, suggesting reduced dipolar interactions due to SiO2 surface passivation. Isothermal remanent magnetization (IRM) and DC demagnetization (DCD) measurements supported this trend, with δM plots confirming weaker dipolar interactions in the coated sample. Bloch’s law fitting of temperature-dependent saturation magnetization showed a smaller Bloch’s constant for pure Fe5C2 and a larger value for Fe5C2/SiO2, reflecting enhanced surface disorder and reduced exchange coupling in the latter. Notably, Fe5C2/SiO2 demonstrated increased coercivity, attributed to decreased dipolar interaction and elevated surface anisotropy. Kneller’s law fitting yielded higher blocking temperatures for Fe5C2 (476 K) than Fe5C2/SiO2 (456 K), highlighting the impact of dipolar interactions on magnetic relaxation. These findings illustrate how SiO2 coatings effectively modulate dipolar interactions and enhance coercivity in Fe5C2 nanocrystals.

Joshi, Pramanand [Department of Physics, Universit↗
DOE OSTI2023

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Weak-Field Hall Resistivity and Spin-Valley Flavor Symmetry Breaking in Magic-Angle Twisted Bilayer Graphene

Near a magic twist angle, the lowest energy conduction and valence bands of bilayer graphene moiré superlattices become extremely narrow. The band dispersion that remains is sensitive to the moiré’s strain pattern, nonlocal tunneling between layers, and filling-factor-dependent Hartree and exchange band renormalizations. In this Letter, we analyze the influence of these band-structure details on the pattern of flavor symmetry breaking observed in this narrow band system and on the associated pattern of Fermi surface reconstructions revealed by weak-field Hall and Shubnikov–de Haas magnetotransport measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Ab Initio Study of Stability, Local Order, and Phase Diagram For a Series of bcc-based Transition Metal Alloys

In this work, a parameter-free electronic structure approach is applied to the study of stability and chemical order in the 15 substitutional body-centered cubic (bcc)-based alloys made of the six transition metals of groups 5 (V, Nb, Ta) and 6 (Cr, Mo, W) of the periodic table. The method is based on a Green’s function description of the electronic structure of the random alloys. Configurational order is treated within the generalized perturbation method, and temperature effects are examined with a generalized mean-field approach. In contrast to the results summarized in the assessed phase diagrams, stability and ordering trends are predicted in a broad range of alloy composition for at least seven alloys, and explanation is found in their electronic structure properties. Short-range order results, thermodynamics analysis, and bcc-based phase diagrams are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Shape evolution of neutron-rich 106,108,110 Mo isotopes in the triaxial degree of freedom

Background: Neutron-rich nuclei with mass number between 100 and 110 attract much attention, since several kinds of shapes, such as spherical, prolate, oblate, and triaxial shapes, are predicted. In particular, for neutron-rich Mo isotopes, different models predict different magnitudes and rigidity of triaxial deformation. Previous interpretations of experimental results based solely on low-lying $2^+_2$ states are insufficient to distinguish between the rigid triaxial shape, $\gamma$ vibration, or $\gamma$-soft rotor. Purpose: The shape evolution of 106 Mo, 108 Mo, and 110 Mo is investigated through their $2^+_1$-state lifetimes, decay-branching ratios of the $2^+_2$ state, and energies of the low-lying collective excited states with $K^π = 0^+, 2^+$, and $4^+$. Method: $\beta$-delayed $\gamma$-ray spectroscopy was employed for neutron-rich Nb and Zr isotopes produced at the RIKEN RI Beam Factory to populate excited states in 106 Mo, 108 Mo, and 110 Mo. The EUroball-RIKEN Cluster Array was used for high-resolution $\gamma$-ray detection and lifetimes of the $2^+_1$ states were determined using the UK fast-timing array of LaBr 3 (Ce) detectors. Results: New $\gamma$-ray transitions and levels are reported, including newly assigned $0^+_2$ states in 108,110 Mo. Quadrupole deformations were obtained for 106,108,110 Mo from their $2^+_1$ energies and lifetimes. The $\beta$-delayed neutron-emission probabilities of 108 Nb and 110 Nb were determined by examining the $\gamma$ rays of their respective daughter decays. Conclusions: In this work, the even-odd energy staggering in the $2^+_2$ band was compared with typical patterns of the $\gamma$-vibrational band, rigid triaxial rotor, and $\gamma$-soft rotor. The very small even-odd staggering of 106 Mo, 108 Mo, and 110 Mo favors a $\gamma$-vibrational band assignment. The kinematic moment of inertia for the $2^+_2$ band showed a trend similar to the ground-state band, which is as expected for the $\gamma$-vibrational band. Beyond-mean-field calculations employing the constrained Hartree-Fock-Bogoliubov and local quasiparticle-random-phase approximation method using the SLy5 + T interaction reproduced the ground and $2^+_2$ bands in 106 Mo and 108 Mo. The collective wave functions are consistent with the interpretation of the $2^+_2$ band as the $\gamma$-vibrational band of the prolate shape. However, the staggering pattern observed in 110 Mo differs from the one suggested in the calculations which predict a $\gamma$-soft rotor. There was no experimental indication of the oblate shape or the $\gamma$-soft rotor predicted in these Mo isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

AI-enabled high-resolution scanning coherent diffraction imaging

Ptychographic imaging is a powerful means of imaging beyond the resolution limits of typical x-ray optics. Capturing recovered images from raw ptychographic data, however, requires the solution of an inverse problem, namely, phase retrieval. Phase retrieval algorithms are computationally expensive, which precludes real-time imaging. In this work, we propose PtychoNN, an approach to solve the ptychography data inversion problem based on a deep convolutional neural network. We demonstrate how the proposed method can be used to predict real-space structure and phase at each scan point solely from the corresponding far-field diffraction data. Our results demonstrate the practical application of machine learning to recover high fidelity amplitude and phase contrast images of a real sample hundreds of times faster than current ptychography reconstruction packages. Furthermore, by overcoming the constraints of iterative model-based methods, we can significantly relax sampling constraints on data acquisition while still producing an excellent image of the sample. Besides drastically accelerating acquisition and analysis, this capability has profound implications for the imaging of dose sensitive, dynamic, and extremely voluminous samples.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

thornado-hydro: A Discontinuous Galerkin Method for Supernova Hydrodynamics with Nuclear Equations of State

This paper describes algorithms for non-relativistic hydrodynamics in the toolkit for high-order neutrino radiation hydrodynamics (thornado), which is being developed for multiphysics simulations of core-collapse supernovae (CCSNe) and related problems with Runge–Kutta discontinuous Galerkin (RKDG) methods. More specifically, thornado employs a spectral type nodal collocation approximation, and we have extended limiters — a slope limiter to prevent non-physical oscillations and a bound-enforcing limiter to prevent non-physical states — from the standard RKDG framework to be able to accommodate a tabulated nuclear equation of state (EoS). To demonstrate the efficacy of the algorithms with a nuclear EoS, we first present numerical results from basic test problems in idealized settings in one and two spatial dimensions, employing Cartesian, spherical-polar, and cylindrical coordinates. Then, we apply the RKDG method to the problem of adiabatic collapse, shock formation, and shock propagation in spherical symmetry, initiated with a 15 M ⊙ progenitor. Herein, we find that the extended limiters improve the fidelity and robustness of the RKDG method in idealized settings. The bound-enforcing limiter improves robustness of the RKDG method in the adiabatic collapse application, while we find that slope limiting in characteristic fields is vulnerable to structures in the EoS — more specifically, in the phase transition from nuclei and nucleons to bulk nuclear matter. The success of these applications marks an important step toward applying RKDG methods to more realistic CCSN simulations with thornado in the future.

Nuclear astrophysics (1129)↗
DOE OSTI2025

Large-𝑁 SU(4) Schwinger boson theory for coupled-dimer antiferromagnets

Here, we develop a systematic large-𝑁 expansion based on the Schwinger boson representation of SU(4) coherent states of dimers for the paradigmatic spin-1/2 bilayer square lattice Heisenberg antiferromagnet. This system exhibits a quantum phase transition between a quantum paramagnetic state and a Néel order state, driven by the coupling constant 𝑔 = 𝐽′/𝐽, which is defined as the ratio between the interdimer 𝐽′ and intradimer 𝐽 exchange interactions. We demonstrate that this approach accurately describes static and dynamic properties on both sides of the quantum phase transition. The critical coupling constant 𝑔 𝑐 ≈ 0.42 and the dynamic spin structure factor reproduce quantum Monte Carlo results with high precision. Notably, the 1/𝑁 corrections reveal the longitudinal mode of the magnetically ordered phase along with the overdamping caused by its decay into the two-magnon continuum. The present large-𝑁 SU(4) Schwinger boson theory can be extended to more general cases of quantum paramagnets that undergo a quantum phase transition into magnetically ordered states.

Schwinger boson method↗
DOE OSTI2026

Bootstrap embedding for interacting electrons in phonon coherent-state mean field

Here, we develop a Fermi–Bose bootstrap embedding framework for the ground state of interacting electrons coupled to a phonon mean field. The method combines bootstrap embedding for correlated electrons with a self-consistent coherent-state mean-field treatment for phonons. This method models the interacting electron–phonon problem as a system of correlated electrons traveling in a self-consistently specified potential landscape, allowing for efficient treatment of large lattice systems. Convergence of the methods for fragment size and total system size is demonstrated for the one-dimensional Hubbard–Holstein model for up to 350 sites. Finite-size scaling is performed to extrapolate to the infinite system size. Benchmarking against the density matrix renormalization group for a small 8-site system at half- and quarter-filling shows an orders-of-magnitude runtime advantage. The comparison further reveals that the method performs best in regimes dominated by localization, such as the Mott insulating phase and the strong-coupling tiny polaron regime, where the local embedding ansatz is still valid. However, due to the mean-field treatment for phonons, we find limitations of our methods in the weakly coupled delocalized region and at the Peierls transition, where quantum phonon fluctuations and long-range kinetic correlations become substantial.

Islam, Shariful [North Carolina State University, ↗