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

The contribution of black carbon to global ice nucleating particle concentrations relevant to mixed-phase clouds

Significance The Earth’s climate is affected by the radiative properties of clouds. In clouds containing both supercooled water droplets and ice particles (mixed-phase clouds), cloud lifetime and albedo are affected by the presence of ice nucleating particles (INP). After decades of research, the importance of black carbon (BC) as an INP remains unresolved. In this work, we determined the contribution of BC to INP concentrations ([INP]) from real-world wildfires and prescribed burns. We found that BC does not contribute significantly to [INP] in near-field biomass burning plumes. These results were applied as parameterizations to global climate model output. Combined with previous work on diesel-engine exhaust, we find that BC does not contribute significantly to global [INP] relevant to mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Structural modeling of high-entropy oxides battery anodes using x-ray absorption spectroscopy

High-entropy oxides (HEOs) are single phase solid solutions where five or more metals share the same sublattice, giving rise to unexpected features in various fields of applications. Recently, HEOs have emerged as an alternative conversion electrode anode material for next-generation Li-ion batteries, where the combination of several different elements in a single solid solution can synergistically act to overcome some of its main drawbacks, improving performance. Due to their chemical complexity, x-ray absorption spectroscopy (XAS) emerges as an appropriate technique to study the electronic (x-ray absorption near edge structure, XANES) and local structure (extended x-ray absorption fine structure, EXAFS) of these compounds as a function of cycling. This work aims to highlight the capabilities of XAS as an element-specific probe to understand a material’s structure at the atomistic level through EXAFS modeling of (MgFeCoNiCuZn)O high-entropy system and how to extract valuable information about the bond distance, number of near neighbors, and local disorder, which are crucial to a full understanding of the electrochemical reaction mechanisms of such battery electrodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Anomalous Hall effect at the Lifshitz transition in ZrTe 5

Zirconium pentatelluride ZrTe 5 is a topological semimetal. The presence of a temperature-induced Lifshitz transition, in which the Fermi level goes from the conduction band to the valence band with increasing temperature, provides unique opportunities to study the interplay between Fermi-surface topology, dynamics of Dirac fermions, and Berry curvature in one system. Here we present a combined experimental and theoretical study and show that a low-energy model can be used to understand the complicated Hall response and large anomalous Hall effect observed in ZrTe 5 over a wide range of temperature and magnetic field. We found that the anomalous Hall contribution dominates the Hall response in a narrow temperature window around the Lifshitz transition, away from which the orbital contribution dominates. Moreover, our results indicate that a topological phase transition coexists with the Lifshitz transition. To conclude, our model provides a unifying framework to understand the Hall effect in semimetals with large Zeeman splitting and nontrivial topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Modeling brittle fracture due to anisotropic thermal expansion in polycrystalline materials

Here, this work investigated brittle fracture of polycrystalline materials due to thermal stresses arising from anisotropic thermal expansion. We used phase-field fracture simulations with the properties of alpha-uranium (α-U) and assumed a linear elastic mechanical response. Three-dimensional simulations were used to predict fracture for various conditions and crystallographic textures. We found that fracture was more pronounced during cooling than during heating because the anisotropy increased with temperature. We also found that the total crack surface area increased with increasing average misorientation, while the net shape change of the material decreased with increasing misorientation. Two-dimensional simulations in which one crystallographic coefficient of thermal expansion (CTE) was set to zero indicated that the expansion behavior in the crystallographic direction with the smallest CTE was the primary cause of fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

SCExAO/CHARIS Near-infrared Integral Field Spectroscopy of the HD 15115 Debris Disk

We present new, near-infrared (1.1–2.4 μm) high-contrast imaging of the debris disk around HD 15115 with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). The SCExAO/CHARIS resolves the disk down to ρ ∼ 0.″2 (r {sub proj} ∼ 10 au), a factor of ∼3–5 smaller than previous recent studies. We derive a disk position angle of PA ∼ 279.°4–280.°5 and an inclination of i ∼ 85.°3–86.2.°. While recent SPHERE/IRDIS imagery of the system could suggest a significantly misaligned two-ring disk geometry, CHARIS imagery does not reveal conclusive evidence for this hypothesis. Moreover, optimizing models of both one- and two-ring geometries using differential evolution, we find that a single ring having a Hong-like scattering phase function matches the data equally well within the CHARIS field of view (ρ ≲ 1″). The disk’s asymmetry, well evidenced at larger separations, is also recovered; the west side of the disk appears, on average, around 0.4 mag brighter across the CHARIS bandpass between 0.″25 and 1″. Comparing Space Telescope Imaging Spectrograph (STIS) 50CCD optical photometry (2000–10500 Å) with CHARIS near-infrared photometry, we find a red (STIS/50CCD−CHARIS broadband) color for both sides of the disk throughout the 0.″4–1″ region of overlap, in contrast to the blue color reported at similar wavelengths for regions exterior to ∼2″. Further, this color may suggest a smaller minimum grain size than previously estimated at larger separations. Finally, we provide constraints on planetary companions and discuss possible mechanisms for the observed inner disk flux asymmetry and color.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Topological order in Mott insulators

Here in this commemorative issue for Phil Anderson, I reproduce parts of a joint unpublished paper that first got me thinking about topological order in Mott insulators. With Phil, we explored aspects of the momentum distribution function and the non-classical rotational inertia in the Bose Hubbard model close to the quantum phase transition. My interest on this question continued, and later in the article I take up the question of topological order in the Kitaev Mott insulator on the honeycomb lattice with bond-dependent interactions in a finite magnetic field. Kitaev obtained the exact solution for zero field that demonstrated an exotic Z 2 quantum spin liquid phase with topological order. In a finite magnetic field, we discovered a new intermediate U(1) gapless quantum spin liquid phase, sandwiched between the known gapped Kitaev spin liquid and a polarized phase, that harbors fractionalized neutral fermionic excitations, dubbed spinons, that remarkably form a Fermi surface in a Mott insulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ultrafast nonlinear absorption of Haldane model quantum dots

We study theoretically the nonlinear absorbance of Haldane model quantum dots (QDs) placed in the field of an ultrashort and strong optical pulse. The absorbance strongly depends on the frequency of the pulse. When the frequency of the pulse is much less than the QD bandgap, the absorbance shows strong dependence on the pulse amplitude and, as a function of an internal phase of the Haldane model, the absorbance has maxima at intermediate values of the phase. When the frequency of the pulse becomes closer to, but still less than, the bandgap, the absorbance has a weak dependence on the pulse amplitude and, as a function of the internal phase, it has a maximum at the phase of 900 when the QD bandgap also has the smallest value. Furthermore, nonlinear electron dynamics in such QD systems changes from almost reversible one at small pulse frequencies to highly irreversible dynamics at large frequencies of the pulse.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modeling Sequences of Earthquakes and Aseismic Slip (SEAS) in Elasto-Plastic Fault Zones With a Hybrid Finite Element Spectral Boundary Integral Scheme

We present a coupled finite element spectral boundary integral framework for modeling sequences of earthquakes and aseismic slip on a 2-D planar rate-and-state fault with off-fault visco-plastic response in the plane strain approximation. The model resolves both slow aseismic deformation and inertia effects during rapid slip. As an application, we perform two sets of simulations with different choices of cohesion to explore the co-evolution of fault slip, bulk plasticity and local stress fields. The first set implements a relatively large value of the cohesion parameter, which results in limiting inelastic strain accumulation to dynamic rupture phases. The second set implements a smaller cohesion, allowing for plastic strain to accumulate in both seismic and aseismic phases. For the first model, our results indicate that the extent and distribution of plastic strain depend on the angle of maximum compressive principal stress. At larger angles, inelastic strain accumulates on the extensional side of a dynamically propagating rupture. At smaller angles, the extent of plasticity is limited to the compressional side of the domain. At smaller cohesion values, off-fault plasticity may occur during aseismic slip, which alters the nucleation characteristics and earthquake sequence pattern. Furthermore, our results at lower cohesion values indicate that plastic strain accumulation may occur in both the extensional and compressional sides of the off-fault bulk even at higher angles of maximum compression. In conclusion, this produces damage patterns that deviate from the traditional off-fault fan-like distribution observed in dynamic rupture simulations and emphasizes the significance of long-term deformation in interpreting observations.

58 GEOSCIENCES↗
DOE OSTI2026

Onset of High-Harmonic Fast Wave Propagation in Annulus Geometry

Previous experiments on NSTX have reported high-harmonic fast wave (HHFW) power losses in the scrape-off layer (SOL), with dependencies on the wave parallel wavenumber, plasma density, and magnetic field. This study presents an analysis of HHFW propagation in the SOL plasma, focusing on parasitic coupling to the waveguide mode bounded by the steep gradient pedestal and the outer wall boundary. The analytic and finite element wave models are discussed using a 2D annulus geometry, which serves as a simplified representation of the NSTX-U plasma. The model shows that a minimum radial gap width is required to support wave propagation, analogous to the higher-order modes in coaxial geometries. The impact of several factors, including gap width, wave parallel phasing, and magnetic field strength, on parasitic wave coupling is examined, highlighting the role of the geometric dimension in relation to the wave perpendicular wavelength. The findings suggest that controlling the edge/SOL density profile (gap control) is crucial for mitigating parasitic coupling when the density in front of the antenna exceeds the density cutoff.

Baek, Seung-Gyou [Massachusetts Inst. of Technolog↗
DOE OSTI2020

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magnetic and singlet phases in the three-dimensional periodic Anderson Model

Heavy-fermion materials are compounds in which localized f orbitals hybridize with delocalized d ones, leading to quasiparticles with large renormalized masses. The presence of strongly correlated f electrons at the Fermi level may also lead to long-range order, such as magnetism, or unconventional superconductivity. From a theoretical point of view, the “standard model” for heavy-fermion compounds is the periodic Anderson model (PAM). Despite being extensively scrutinized, its thermodynamic properties in three-dimensional (3D) lattices have not been carefully addressed by unbiased methodologies. Here we investigate the 3D PAM employing state-of-the-art finite-temperature auxiliary-field quantum Monte Carlo simulations. Furthermore, we present the behavior of the kinetic energy, the entropy, the specific heat, and the double occupancy as functions of the temperature and the hybridization strength. From these quantities, and by the analysis of the spin-spin correlation functions, we investigate the occurrence of magnetic phase transitions at finite temperatures, and we determine the phase diagram of the model, including the behavior of the Néel temperature as a function of the external parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

pyTDGL: Time-dependent Ginzburg-Landau in Python

Time-dependent Ginzburg-Landau (TDGL) theory is a phenomenological model for the dynamics of superconducting systems. Due to its simplicity in comparison to microscopic theories and its effectiveness in describing the observed properties of the superconducting state, TDGL is widely used to interpret or explain measurements of superconducting devices. Here, we introduce pyTDGL, a Python package that solves a generalized TDGL model for superconducting thin films of arbitrary geometry, enabling simulations of vortex and phase dynamics in mesoscopic superconducting devices. pyTDGL can model the nonlinear magnetic response and dynamics of multiply connected films, films with multiple current bias terminals, and films with a spatially inhomogeneous critical temperature. We demonstrate these capabilities by modeling quasi-equilibrium vortex distributions in irregularly shaped films, and the dynamics and current-voltage-field characteristics of nanoscale superconducting quantum interference devices (nanoSQUIDs).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electromagnetic turbulence in increased β plasmas in the Large Plasma Device

The variation of pressure-gradient-driven turbulence with plasma $\beta$ (up to $\beta \approx 15\,\%$ ) is investigated in linear, magnetized plasma. The magnitude of magnetic fluctuations is observed to increase substantially with increasing $\beta$ . More importantly, parallel magnetic fluctuations are observed to dominate at higher $\beta$ values, with $\delta B_\parallel / \delta B_\perp \approx 2$ and $\delta B / B_0 \approx 1\,\%$ . Parallel magnetic fluctuations are strongly correlated with density fluctuations and the two are observed to be out of phase. The relative magnitude of and cross-phase between density and parallel magnetic field fluctuations are consistent with the dynamic pressure balance ( $P+{B_{0}^2}/{2\mu _0} = \textrm {constant}$ ). A local slab model theory for electromagnetic, modified drift Alfvén waves, including parallel magnetic fluctuations, shows partial agreement with experimental observations.

Physics↗
DOE OSTI2023

Positional ordering induced by dynamic steric interactions in superparamagnetic rods

The dynamic motion produced by precessing magnetic fields can drive matter into far-from-equilibrium states. We predict 1D periodic ordering in systems of precessing rods when magnetic interactions between rods remain insignificant. The precession angle of the rods is completely determined by the field's precession angle and the ratio of the field's precession frequency and the characteristic response frequency of the rods. Here, we develop a molecular dynamics model that explicitly calculates magnetic interactions between particles, and we also simulate rods in the limit of a strong and fast precessing magnetic field where inter-rod magnetic interactions are negligible, using a purely steric model. Our simulations show how steric interactions drive the rods from a positionally disordered phase (nematic) to a layered (smectic) phase. As the rod precession angle increases, the nematic–smectic transition density significantly decreases. The minimization of unfavorable steric interactions also induces phase separation in binary mixtures of rods of different lengths. This effect is general to any force that produces precession in elongated particles. This work will advance the understanding and control of out-of-equilibrium soft matter systems.

42 ENGINEERING↗
DOE OSTI2021

Automatic Differentiation in MetaPhysicL and Its Applications in MOOSE

Efficient solution via Newton’s method of nonlinear systems of equations requires an accurate representation of the Jacobian, corresponding to the derivatives of the component residual equations with respect to the degrees of freedom. In practice these systems of equations often arise from spatial discretization of partial differential equations used to model physical phenomena. These equations may involve domain motion or material equations that are complex functions of the systems’ degrees of freedom. Computing the Jacobian by hand in these situations is arduous and prone to error. Finite difference approximations of the Jacobian or its action are prone to truncation error, especially in multiphysics settings. Symbolic differentiation packages may be used, but often result in an excessive number of terms in realistic model scenarios. An alternative to symbolic and numerical differentiation is automatic differentiation (AD), which propagates derivatives with every elementary operation of a computer program, corresponding to continual application of the chain rule. Automatic differentiation offers the guarantee of an exact Jacobian at a relatively small overhead cost. In this work, we outline the adoption of AD in the Multiphysics Object Oriented Simulation Environment (MOOSE) via the MetaPhysicL package. We describe the application of MOOSE’s AD capability to several sets of physics that were previously infeasible to model via hand-coded or Jacobian-free simulation techniques, including arbitrary Lagrangian-Eulerian and level-set simulations of laser melt pools, phase-field simulations with free energies provided through neural networks, and metallic nuclear fuel simulations that require inner Newton loop calculation of nonlinear material properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Fractional vortices, Z 2 gauge theory, and the confinement-deconfinement transition

Here, in this Letter we discuss the classical three-dimensional XY model whose nearest-neighbor interaction is a mixture of cos (θ i -θ j ) (ferromagnetic) and cos 2 (θ i -θ j ) (nematic). This model is dual to a theory with integer and half-integer vortices. While both types of vortices interact with a noncompact U(1) gauge field, the half-integer vortices interact with an extra interaction mediated by a Z 2 gauge field. We shall discuss the confinement-deconfinement transition of the half-integer vortices, the Wilson and the 't Hooft loops, and their mutual statistics in path integral language. In addition, we shall present a quantum version of the classical model which exhibits these physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Reentrant Correlated Insulators in Twisted Bilayer Graphene at 25 T ( 2 π Flux)

Twisted bilayer graphene (TBG) is remarkable for its topological flat bands, which drive strongly interacting physics at integer fillings, and its simple theoretical description facilitated by the Bistritzer-MacDonald Hamiltonian, a continuum model coupling two Dirac fermions. Because of the large moiré unit cell, TBG offers the unprecedented opportunity to observe reentrant Hofstadter phases in laboratory-strength magnetic fields near 25 T. This Letter is devoted to magic angle TBG at 2π flux where the magnetic translation group commutes. Here we use a newly developed gauge-invariant formalism to determine the exact single-particle band structure and topology. We find that the characteristic TBG flat bands reemerge at 2π flux, but, due to the magnetic field breaking C 2z T, they split and acquire Chern number ±1. We show that reentrant correlated insulating states appear at 2π flux driven by the Coulomb interaction at integer fillings, and we predict the characteristic Landau fans from their excitation spectrum.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthetic tensor gauge fields

Synthetic gauge fields have provided physicists with a unique tool to explore a wide range of fundamentally important phenomena. However, most experiments have focused on synthetic vector gauge fields. The very rich physics brought about by coupling tensor gauge fields to fracton phases of matter remains unexplored in laboratories. Here, we propose schemes to realize synthetic tensor gauge fields that address dipoles instead of single particles. A lattice tilted by a strong linear potential and a weak quadratic potential yields a rank-2 electric field for a dipole formed by a particle-hole pair. Such a rank-2 electric field leads to a new type of Bloch oscillations, which modulates the quadrupole moment and preserves the dipole moment of the system. In higher dimensions, the interplay between interactions and vector gauge potentials imprints a phase to the ring-exchange interaction and thus generates synthetic tensor gauge fields. Such tensor gauge fields make it possible to realize a dipolar Harper-Hofstadter model in laboratories. The resultant dipolar Chern insulators feature chiral edge currents of dipoles in the absence of net charge currents. Published by the American Physical Society 2025

Zhang, Shaoliang (ORCID:000000016635044X)↗