SEARCH · Engineering Papers

Results for “Two-dimensional phase field equation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2024

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Wave-packet manipulation of He Rydberg states by a seeded free-electron laser

We report a two-dimensional pump-control-probe spectroscopic study of the dynamics of singly excited He Rydberg-state wave packets with a seeded extreme ultraviolet (XUV) free-electron laser (FEL) source. A pair of coherent XUV pulses, defined by their coarse time separation and relative phase, created and manipulated the wave packets. The He atoms were postionized by infrared (IR) pulses, and the ion yield was measured as a function of XUV phase and IR arrival time. We tagged and sorted the relative phase of the XUV pulse pair on a single-shot basis by fitting each FEL spectrum with a suitable function that accounts for nonidealities of the XUV pulse pairs, associated with the seeding process; more generally, the fit returns the time-dependent electric field of the FEL spectra. The experimental two-dimensional maps of ion yields, measured as a function of IR (probe) delay and of XUV (pump-control) phase, were compared with the solution of the first-order time-dependent Schrödinger equation for this field. Despite the fact that the experimental conditions imply strong excitation, beyond the approximations of first-order perturbation theory, the simulated map satisfactorily reproduces the experimental one for temporally well-separated pulses. We show that by selecting data at appropriate values of pump-control phase, we enhance or suppress the amplitude of chosen wave-packet components consisting of two or more Rydberg states. When the temporal overlap of the pulse pair cannot be neglected, the phase reconstruction is underdetermined, and we provide a simplified comparison between data and simulations. Published by the American Physical Society 2024

Dumergue, Mathieu (ORCID:0000000236984406)↗

DOE OSTI2021

Dynamical vertex approximation for many-electron systems with spontaneously broken SU(2) symmetry

We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.

2-dimensional systems↗

DOE OSTI2021

A phase-field method for boiling heat transfer

Here we present a phase field method for heat transfer in two-phase flow with boiling. The vapor/liquid interface evolution is modeled by the Cahn-Hilliard equation. The phase change rate is determined by accounting for the heat conduction balance on either vapor or liquid side of the interfacial area, depending on which side the temperature is assumed to be maintained at the saturation temperature during boiling. The velocity correction scheme proposed by Dong & Shen [27] is extended to solve the Navier-Stokes equations for a non-solenoidal velocity field, and the entropy viscosity method is employed for stabilization. The phase change model is verified by two-dimensional simulations of a vapor bubble growing in super-heated liquid and in film boiling. In both cases, mesh independence of the results is systematically performed. Subsequently, the method is applied to predict the growth of three-dimensional vapor bubble in a rectangular microchannel with boiling flow, achieving good agreement with experimental measurements and available simulation results using the level-set method. The numerical experiments demonstrate that the required mesh resolution for the phase field method is comparable with that of volume of fluid (VoF) and level-set methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

DOE OSTI2021

Disentangling observable dependence in SCET I and SCET II anomalous dimensions: angularities at two loops

The resummation of radiative corrections to collider jet observables using soft collinear effective theory is encoded in differential renormalization group equations (RGEs), with anomalous dimensions depending on the observable under consideration. This observable dependence arises from the ultraviolet (UV) singular structure of real phase space integrals in the effective field theory. We show that the observable dependence of anomalous dimensions in SCET I problems can be disentangled by introducing a suitable UV regulator in real radiation integrals. Resummation in the presence of the new regulator can be performed by solving a two-dimensional system of RGEs in the collinear and soft sectors, and resembles many features of resummation in SCET II theories by means of the rapidity renormalization group. We study the properties of SCET I with the additional regulator and explore the connection with the system of RGEs in SCET II theories, highlighting some universal patterns that can be exploited in perturbative calculations. As an application, we compute the two-loop soft and jet anomalous dimensions for a family of recoil-free angularities and give new analytic results. This allows us to study the relations between the SCET I and SCET II limits for these observables. We also discuss how the extra UV regulator can be exploited to calculate anomalous dimensions numerically, and the prospects for numerical resummation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Room-temperature polariton quantum fluids in halide perovskites

Abstract Quantum fluids exhibit quantum mechanical effects at the macroscopic level, which contrast strongly with classical fluids. Gain-dissipative solid-state exciton-polaritons systems are promising emulation platforms for complex quantum fluid studies at elevated temperatures. Recently, halide perovskite polariton systems have emerged as materials with distinctive advantages over other room-temperature systems for future studies of topological physics, non-Abelian gauge fields, and spin-orbit interactions. However, the demonstration of nonlinear quantum hydrodynamics, such as superfluidity and Čerenkov flow, which is a consequence of the renormalized elementary excitation spectrum, remains elusive in halide perovskites. Here, using homogenous halide perovskites single crystals, we report, in both one- and two-dimensional cases, the complete set of quantum fluid phase transitions from normal classical fluids to scatterless polariton superfluids and supersonic fluids—all at room temperature, clear consequences of the Landau criterion. Specifically, the supersonic Čerenkov wave pattern was observed at room temperature. The experimental results are also in quantitative agreement with theoretical predictions from the dissipative Gross-Pitaevskii equation. Our results set the stage for exploring the rich non-equilibrium quantum fluid many-body physics at room temperature and also pave the way for important polaritonic device applications.

42 ENGINEERING↗

DOE OSTI2023

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗

DOE OSTI2024

Multidimensional coherent spectroscopy of light-driven states and their collective modes in multiband superconductors

We present a comprehensive theory of light-controlled multiband superconductivity, and apply it to predict distinctive signatures of light-driven superconducting (SC) states in terahertz multidimensional coherent spectroscopy (THz-MDCS) experiments. We first derive gauge-invariant Maxwell-Bloch equations for multiband BCS superconductors with spatial fluctuations. We consider driving electromagnetic fields determined self-consistently by Maxwell's equations. By calculating the THz-MDCS spectra measured experimentally in the clean SC limit, we identify unique signatures of finite-momentum Cooper-pairing states that live longer than the laser pulse. They are controlled by a pair of THz laser pulses with well-defined relative phase (pulse pair). The pseudospin oscillators that describe the properties of these SC states are parametrically driven by both finite-momentum Cooper pairing and by time oscillations of the order-parameter relative phase. We show that such strong parametric driving leads to drastic changes in the THz-MDCS spectral shape from the predictions of third-order nonlinear susceptibility calculations. These spectral changes strongly depend on the interband-to-intraband interaction ratio and on the collective modes of the light-driven state. For negligible interband interaction, the spectra show a transition with increasing field, from traditional pump-probe, four-wave-mixing, and third-harmonic generation peaks determined by the laser frequency to sidebands determined by the excitations of the driven system. These sidebands emerge from difference-frequency Raman processes in the nonequilibrium SC state. For interband couplings weaker than the intraband pairing, we show that the Leggett phase collective mode leads to harmonic sidebands around the traditional pump-probe peaks. Additional Higgs collective mode peaks result from light-induced inversion-symmetry breaking in a thin-film geometry. For strong interband coupling, we find a transition from a nonequilibrium finite Cooper-pair momentum state characterized by hybrid-Higgs amplitude mode peaks in THz-MDCS spectra to a driven state identified experimentally by the emergence of Floquet-type sidebands at bi-Higgs frequencies. Those dominant bi-Higgs-frequency satellites are manifestations of a new order parameter relative phase collective mode that characterizes the nonequilibrium SC state. The predicted interaction- and field-dependent transitions in the spectral profile allow us to propose THz-MDCS experiments for quantum tomography of light-driven superconductivity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Squeezing quantum states in three-dimensional twisted crystals

Bloch's theorem provides a conventional starting point for describing wave propagation in periodic media, but in ordered materials where competing spatial periods coexist it is rendered ineffective, often with dramatic consequences. Here we develop an alternate approach that uses coherent free-particle vortex states to study quantum states in supertwisted crystals: three-dimensional stacks of atomically thin two-dimensional layers. Here, this formalism leads naturally to the representation of the spectrum using squeezed coherent states, and it reveals the crucial role of a Coriolis coupling in the equations of motion. This identifies an underlying noncommutative geometry and novel edge state structure in a family of complex ordered structures.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Interlayer coherence in superconductor bilayers

We investigate the possibility and implications of interlayer coherence in electrically isolated superconducting bilayers. We find that in a mean-field approximation bilayers can have superconducting and excitonic order simultaneously if repulsive interactions between layers are sufficiently strong. The excitonic order implies interlayer phase coherence and can be studied in a representation of symmetric and antisymmetric bilayer states. When both orders are present we find several solutions of the mean-field equations with different values of the the symmetric and antisymmetric state pair amplitudes. The mixed state necessarily has nonzero pair amplitudes for electrons in different layers in spite of the repulsive interlayer interactions, and these are responsible for spatially indirect Andreev reflection processes in which an incoming electron in one layer can be reflected as a hole in the opposite layer. Furthermore, we evaluate layer diagonal and off-diagonal current-voltage relationships that can be used to identify this state experimentally.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Multiband superconductivity in strongly hybridized 1 T ' - WTe 2 / NbSe 2 heterostructures

The interplay of topology and superconductivity has become a subject of intense research in condensed-matter physics for the pursuit of topologically nontrivial forms of superconducting pairing. An intrinsically normal-conducting material can inherit superconductivity via electrical contact to a parent superconductor via the proximity effect, usually understood as Andreev reflection at the interface between the distinct electronic structures of two separate conductors. However, at high interface transparency, strong coupling inevitably leads to changes in the band structure, locally, owing to hybridization of electronic states. Here, we investigate such strongly proximity-coupled heterostructures of monolayer $1T$' - $\mathrm{WTe_2}$ grown on NbSe 2 by van der Waals epitaxy. The superconducting local density of states, resolved in scanning tunneling spectroscopy down to 500 mK, reflects a hybrid electronic structure well described by a multiband framework based on the McMillan equations which captures the multiband superconductivity inherent to the NbSe 2 substrate and that is induced by proximity to $\mathrm{WTe_2}$ self-consistently. Our material-specific tight-binding model captures the hybridized heterostructure quantitatively and confirms that strong interlayer hopping gives rise to a semimetallic density of states in the two-dimensional $\mathrm{WTe_2}$ bulk, even for nominally band-insulating crystals. The model further accurately predicts the measured order parameter Δ ≅ 0.6 meV induced in the $\mathrm{WTe_2}$ monolayer bulk, stable beyond a 2 T magnetic field. So we believe that our detailed multiband analysis of the hybrid electronic structure provides a useful tool for sensitive spatial mapping of induced order parameters in proximitized atomically thin topological materials.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Crystalline solutions of the Kohn-Sham equations in the fractional quantum Hall regime

A Kohn-Sham density functional approach has recently been developed for the fractional quantum Hall effect, which maps the strongly interacting electrons into a system of weakly interacting composite fermions subject to an exchange correlation potential as well as a density dependent gauge field that mimics the "flux quanta" bound to composite fermions. To get a feel for the role of various terms, we study the behavior of the self-consistent solution as a function of the strength of the exchange correlation potential, which is varied through an ad hoc multiplicative factor. Here, we find that a crystal phase is stabilized when the exchange correlation interaction is sufficiently strong relative to the composite-fermion cyclotron energy. Various properties of this crystal are examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗