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

Inhomogeneous mean-field approach to collective excitations near the superfluid–Mott glass transition

Highlights: • Unconventional collective excitations across the superfluid–Mott glass transition. • Higgs mode localized in superfluid and insulating phases. • Broad uncritical amplitude response function at quantum phase transition. • Goldstone mode delocalizes in the superfluid phase. We develop an inhomogeneous quantum mean-field approach to the behavior of collective excitations across the superfluid–Mott glass quantum phase transition in two dimensions, complementing recent quantum Monte Carlo simulations (Puschmann et al. 2020). In quadratic (Gaussian) approximation, the Goldstone (phase) and Higgs (amplitude) modes completely decouple. Each is described by a disordered Bogoliubov Hamiltonian which can be solved by an inhomogeneous multi-mode Bogoliubov transformation. We find that the Higgs mode is spatially localized in both phases. The corresponding scalar spectral function shows a broad peak that is noncritical in the sense that its peak frequency does not soften but remains nonzero across the quantum phase transition. In contrast, the lowest-energy Goldstone mode delocalizes in the superfluid phase, leading to a zero-frequency spectral peak. We compare these findings to the results of the quantum Monte Carlo simulations. We also relate them to general results on the localization of bosonic excitations, and we discuss the limits and generality of our approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Optical spectroscopy and modeling of uranium gas-phase oxidation: Progress and perspectives

We report that studies related to U gas-phase oxidation through plasma- and thermo-chemistry are important for many fields, including environmental monitoring, forensic analysis, debris analysis in a weapon detonation event, and nucleation physics. Recently, significant efforts have been made to understand the chemical pathways involved in the progression from U atoms to diatoms (UO) and polyatomic molecules (UxOy), employing optical spectroscopy tools and computational modeling. In many studies, laser ablation of U or a U-containing flow reactor are used as a highly resource-efficient, repeatable, tunable, and lab-scale testbed for studying gas-phase oxidation in U plasmas. The spectroscopic analysis of high-temperature gas-phase oxidation of U is challenging due to the congested U spectra, resolution limitations of instrumentation, and the numerous chemical reaction pathways possible. This article focuses on the current understanding and challenges related to studying U plasma chemistry, specifically U gas-phase oxidation and molecular formation, via optical spectroscopy of plasmas and associated computational and spectral modeling. The physical and chemical processes involved in the evolution from U atoms to U oxide molecules to nanoparticles and agglomerates (i.e., debris) are discussed in the context of optical spectroscopic studies. The article concludes by highlighting opportunities for future research efforts based on existing knowledge published in the literature.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Microwave Heating-Induced Temperature Gradients in Liquid–Liquid Biphasic Systems

Microwaves (MWs) can enable the electrification and intensification of chemical manufacturing. They have been applied to various unit separations, such as drying, distillation, and extraction, entailing gas–liquid and solid–liquid systems. However, a limited quantitative understanding of MW-heated liquid–liquid biphasic systems related to extraction exists. This work measures the temporal and spatial temperature difference between an aqueous and an organic phase in batch and continuous microfluidic modes. We demonstrate permanent temperature differences between phases over 35 °C and spatiotemporal periodic and quasiperiodic oscillations modulated by the flow patterns. The temperature differences are primarily driven by the faster absorption rate of MW irradiation by the aqueous phase versus the slower heat transfer from the aqueous phase to the organic phase. These are amplified by low specific interfacial area and modifications of the electromagnetic field. We employ a multiphysics model to predict the temperature difference in a batch system. The model is in good agreement with the experiments. We demonstrate a strong effect of input power, dielectric properties of organic solvents, the volume of solvents, and the volume ratio between phases on the temperature difference. A simple analytical model describes the temperature difference and provides design principles. Furthermore, the combined approach offers new insights into the design and optimization of the MW-heated biphasic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2026

Benchmarking of massively parallel phase-field codes for directional solidification

We present a detailed benchmark comparing two state-of-the-art phase-field implementations for simulating alloy solidification under experimentally relevant conditions. The study investigates the directional solidification of Al-3wt%Cu under high-velocity solidification conditions and SCN-0.46wt% camphor under microgravity conditions from National Aeronautics and Space Administration (NASA) DECLIC-DSI-R experiments. Both codes, one employing finite-difference discretization with uniform mesh and GPU-acceleration (GPU-PF) and the other one employing finite-element discretization with adaptive-mesh and CPU-parallelization (PRISMS-PF), solve the same quantitative phase-field formulation that incorporates an anti-trapping current for the solidification of dilute alloys. We evaluate the predictions of each code for dendritic morphology, primary spacing, and tip dynamics in both 2D and 3D, as well as their numerical convergence and computational performance. While existing benchmark problems have primarily focused on simplified or small-scale simulations, they do not reflect the computational and modeling challenges posed by employing experimentally relevant time and length scales. Our results provide a practical framework for assessing phase-field code performance as well as validating and facilitating their application in integrated computational materials engineering (ICME) workflows that require integration with realistic experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Magnetic phase transitions in quantum spin-orbital liquids

Here, we investigate the spin and orbital correlations of a superexchange model with spin $\textit{S}$ = 1 and orbital $\textit{L}$ = 1 relevant for $5d^4$ transition-metal Mott insulators, using exact diagonalization and density matrix renormalization group (DMRG). For spin-orbit coupling λ = 0, the orbitals are in an entangled state that is decoupled from the spins. We find two phases with increasing λ: (I) the $\textit{S}$2 phase with two peaks in the structure factor for λ ≤ λ$_{c1}$ ≈ 0.34$\textit{J}$ where $\textit{J}$ is the ferromagnetic exchange; and, (II) the $\textit{S}$1 phase for λ$_{c1}$ < λ ≤ λ$_{c2}$ ≈ 1.2$\textit{J}$ with emergent antiferromagnetic correlations. Both $\textit{S}$1 and $\textit{S}$2 phases are shown to exhibit power-law correlations, indicative of a gapless spectrum. Upon increasing λ > λ$_{c2}$ leads to a product state of local spin-orbital singlets that exhibit exponential decay of correlations, indicative of a gapped phase. We obtain insights into the phases from the well-known Uimin-Lai-Sutherland model in an external field that provides an approximate description of our model within mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Cosmic Coincidences of Primordial-Black-Hole Dark Matter

If primordial black holes (PBHs) contribute more than 10% of the dark matter (DM) density, their energy density today is of the same order as that of the baryons. Such a cosmic coincidence might hint at a mutual origin for the formation scenario of PBHs and the baryon asymmetry of the Universe. Baryogenesis can be triggered by a sharp transition of the rolling rate of inflaton from slow-roll to (nearly) ultraslow-roll phases that produce large curvature perturbations for PBH formation in single-field inflationary models. We show that the baryogenesis requirement drives the PBH contribution to DM, along with the inferred PBH mass range, the resulting stochastic gravitational wave background frequency window, and the associated cosmic microwave background tensor-to-scalar ratio amplitude, into potentially observable regimes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Validation of Maui PSCAD Model: Motivation, Methodology, and Lessons Learned: Preprint

As power systems become increasingly dependent on inverter-based resources, the dominant physical and control interactions that determine system stability may no longer be accurately captured in positive sequence dynamic simulation tools; electromagnetic transient (EMT) simulations may be needed. This paper describes the process of creating and validating an EMT model of the entire Maui island transmission system in PSCAD software. The model is validated against the utility's PSSE model and against field data for an event that consisted of a single-phase fault followed by a generation trip. The methodology and lessons learned are summarized here.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2024

Particle-Hole Asymmetric Ferromagnetism and Spin Textures in the Triangular Hubbard-Hofstadter Model

In a lattice model subject to a perpendicular magnetic field, when the lattice constant is comparable to the magnetic length, one enters the “Hofstadter regime,” where continuum Landau levels become fractal magnetic Bloch bands. Strong mixing between bands alters the nature of the resulting quantum phases compared to the continuum limit; lattice potential, magnetic field, and Coulomb interaction must be treated on equal footing. Using determinant quantum Monte Carlo and density matrix renormalization group techniques, we study this regime numerically in the context of the Hubbard-Hofstadter model on a triangular lattice. In the field-filling phase diagram, we find a broad wedge-shaped region of ferromagnetic ground states for filling factor ν ≤ 1 , bounded below by filling factor ν = 1 and bounded above by half filling the lowest Hofstadter subband. We observe signatures of SU(2) quantum Hall ferromagnetism at filling factors ν = 1 and ν = 3 . The phases near ν = 1 are particle-hole asymmetric, and we observe a rapid decrease in ground-state spin polarization consistent with the formation of skyrmions only on the electron doped side. At large fields, above the ferromagnetic wedge, we observe a low-spin metallic region with spin correlations peaked at small momenta. We argue that the phenomenology of this region likely results from exchange interaction mixing fractal Hofstadter subbands. The phase diagram derived beyond the continuum limit points to a rich landscape to explore interaction effects in magnetic Bloch bands. Published by the American Physical Society 2024

Ding, Jixun K.↗
DOE OSTI2024

On CCN Effects upon Convective Cold Pool Timing and Features

Abstract Cold pools produced by deep convection can initiate new convection, and their representation in larger-scale weather and climate models could improve prediction of the extent and timing of upscale growth. Cold pools originate from latent cooling from precipitation changing phase, but little attention has been paid to microphysical influences on cold pool characteristics, particularly CCN effects. Datasets obtained from the CACTI and RELAMPAGO field campaigns, along with idealized numerical modeling, are utilized to investigate the hypothesis that convective storms forming in higher-CCN environments generate their first surface rainfall later, delaying cold pool initiation. Aircraft observations of CCN and shallow convection on 9 days do suggest a CCN effect. Those ingesting more CCN contained fewer drizzle drops, although a decreased cloud depth with increasing CCN was also likely a limiting factor. In three of those cases that later developed into deep convection, the timing of cold pool onset was not ubiquitously delayed in environments with more CCN. Idealized numerical simulations suggest that an ordinary thunderstorm can experience small delays in cold pool onset with increasing CCN due to changes in graupel production, but CCN effects on the cold pool from a supercell thunderstorm can be easily overpowered by its unique dynamics. A strong inverse relationship between cold pool strength, expansion rate, and depth with increasing CCN is suggested by the results of the ordinary thunderstorm simulation. Further consideration of CCN appears warranted for future cold pool parameterization development, but other environmental factors affecting storm morphology and precipitation cannot be ignored.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Gapless to gapless phase transitions in quantum spin chains

We investigate spin chains with bilinear-biquadratic (BLBQ) spin interactions as a function of an applied magnetic field h. At the Uimin-Lai-Sutherland (ULS) critical point we find a gapless to gapless transition revealed by the dynamical structure factor S(q,ω) as a function of h. At h=0, the envelope of the lowest-energy excitations goes soft at two points, q 1 =2π/3 and q 2 =4π/3, dubbed the phase A. With increasing field, the spectral peaks at each of the gapless points bifurcate, making in total four soft modes, and combine to form a new set of excitations that soften at a single point q=π at h c1 ≈0.94. Beyond h c1 the system enters another gapless B phase until the transition at h c2 =4 to the fully polarized phase. We compare the ULS model results with those for the Affleck-Kennedy-Lieb-Tasaki model as a representative of the gapped Haldane phase. We explain the mechanism of the gapless to gapless transition in the ULS model using its conserved charges and a spinon band picture. We also discuss the universality of central charges of the BLBQ family of models subjected to a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Many-body reactive force field development for carbon condensation in C/O systems under extreme conditions

In this paper, we describe the development of a reactive force field for C/O systems under extreme temperatures and pressures, based on the many-body Chebyshev Interaction Model for Efficient Simulation (ChIMES). The resulting model, which targets carbon condensation under thermodynamic conditions of 6500 K and 2.5 g cm –3 , affords a balance between model accuracy, complexity, and training set generation expense. We show that the model recovers much of the accuracy of density functional theory for the prediction of structure, dynamics, and chemistry when applied to dissociative condensed phase systems at 1:1 and 1:2 C:O ratios, as well as molten carbon. Our C/O modeling approach exhibits a 104 increase in efficiency for the same system size (i.e., 128 atoms) and a linear system size scalability over standard quantum molecular dynamics methods, allowing the simulation of significantly larger systems than previously possible. We find that the model captures the condensed-phase reaction-coupled formation of carbon clusters implied by recent experiments, and that this process is susceptible to strong finite size effects. Overall, we find the present ChIMES model to be well suited for studying chemical processes and cluster formation at pressures and temperatures typical of shock waves. We expect that the present C/O modeling paradigm can serve as a template for the development of a broader high pressure–high temperature force-field for condensed phase chemistry in organic materials.

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

Field-induced transitions in the Kitaev material α - RuCl 3 probed by thermal expansion and magnetostriction

High-resolution thermal expansion and magnetostriction measurements were performed on single crystals of α - RuCl 3 in magnetic fields applied parallel to the Ru-Ru bonds. The length changes were measured in the direction perpendicular to the honeycomb planes. Here, our data show clear thermodynamic characteristics for the field-induced phase transition at the critical field μ 0 H c 1 = 7.8 ( 2 ) T where the antiferromagnetic zigzag order is suppressed. At higher fields, a kink in the magnetostriction coefficient signals an additional transition or crossover around μ 0 H c 2 ≈ 11 T . The extracted Grüneisen ratio shows typical hallmarks for quantum criticality near H c 1 but also displays anomalous behavior above H c 1 . We compare our experimental data with spin-wave calculations employing a minimal Kitaev-Heisenberg model in the semiclassical limit. Most of the salient features are in agreement with each other, however, the peculiar features in the region above H c 1 cannot be accounted for in our semiclassical modeling and hence suggest a genuine quantum nature. We construct a phase diagram for α - RuCl 3 in a magnetic field along the Ru-Ru bonds, displaying a zigzag ordered state below H c 1 , a quantum paramagnetic regime between H c 1 and H c 2 , and a semiclassical partially polarized state above H c 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Tensor network study of the m = 1/2 magnetization plateau in the Shastry-Sutherland model at finite temperature

The two-dimensional infinite projected entangled pair state tensor network is evolved in imaginary time with the full update (FU) algorithm to simulate the Shastry-Sutherland model in a magnetic field at finite temperature directly in the thermodynamic limit. We focus on the phase transition into the m = 1/2 magnetization plateau, which was observed in experiments on SrCu 2 (BO 3 ) 2 . For the largest simulated bond dimension, the early evolution in the high-temperature regime is simulated with the simple update (SU) scheme and then, as the correlation length increases, continued with the FU scheme towards the critical regime. We apply a small symmetry-breaking bias field and then extrapolate towards zero bias using a simple scaling theory in the bias field. The combined SU + FU scheme provides an accurate estimate of the critical temperature, even though the results could not be fully converged in the bond dimension in the vicinity of the transition. Here, the critical temperature estimate is improved with a generalized scaling theory that combines two divergent length scales: One due to the bias, and the other due to the finite bond dimension. The obtained results are consistent with the transition being in the universality class of the two-dimensional classical Ising model. The estimated critical temperature is 3.5(2) K, which is well above the temperature 2.1 K used in the experiments.

2-dimensional systems↗
DOE OSTI2023

Network model for periodically strained graphene

The long-wavelength physics of monolayer graphene in the presence of periodic strain fields has a natural chiral scattering network description. When the strain field varies slowly compared to the graphene lattice and the effective magnetic length of the induced valley pseudomagnetic field, the low-energy physics can be understood in terms of valley-polarized percolating domain-wall modes. Inspired by a recent experiment, we consider a strain field with threefold rotation and mirror symmetries but without twofold rotation symmetry, resulting in a system with the connectivity of the oriented kagome network. Scattering processes in this network are captured by a symmetry-constrained phenomenological S matrix. Here, we analyze the phase diagram of the kagome network and show that the bulk physics of the strained graphene can be qualitatively captured by the network when we account for a percolation transition at charge neutrality. We also discuss the limitations of this approach to properly account for boundary physics.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2020

Confinement, reduced entanglement, and spin-glass order in a random quantum spin-ice model

In this work, we study an effective spin model derived perturbatively from a random transverse-field Ising model on the pyrochlore lattice. The model consists of spin configurations on the pyrochlore lattice, restricted to the spin-ice subspace, with spins interacting with random Ising exchange couplings as well as ring exchanges along the hexagons of the lattice. This model is studied by exact diagonalization up to N=64 site systems. Furthermore, we calculate spin-glass correlation functions and local entanglement entropy S T between spins in a single tetrahedron and the rest of the system. We find that the model undergoes two phase transitions. At weak randomness the model is in a quantum spin-ice phase where S T =ln 6. Increasing randomness, at low transverse fields, first leads to a frozen phase, with long-range spin-glass order and S T =ln 2 corresponding to the Cat states associated with Ising order. Further increase in randomness leads to a random resonating-hexagon phase with a frozen backbone of spins and a broad distribution of entanglement entropies. The implications of these studies for non-Kramers rare-earth pyrochlores are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Industrial PLC Network Modeling and Parameter Identification Using Sensitivity Analysis and Mean Field Variational Inference

A multiple input multiple output (MIMO) power line communication (PLC) model for industrial facilities was developed that uses the physics of a bottom-up model but can be calibrated like top-down models. The PLC model considers 4-conductor cables (three-phase conductors and a ground conductor) and has several load types, including motor loads. The model is calibrated to data using mean field variational inference with a sensitivity analysis to reduce the parameter space. The results show that the inference method can accurately identify many of the model parameters, and the model is accurate even when the network is modified.

47 OTHER INSTRUMENTATION↗