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 793 records · Page 44

Anisotropic domains and antiferrodistortive-transition controlled magnetization in epitaxial manganite films on vicinal SrTiO 3 substrates

We studied the microstructural evolution and magnetism of ferroelastic La 0.9 Sr 0.1 MnO 3 (LSMO) epitaxial thin films grown on SrTiO 3 (001) substrates with different miscut angles. The substrate miscut angle plays a critical role in controlling the in-plane magnetic anisotropy. The microscopic origin of such magnetic anisotropy is attributed to the formation of anisotropic stripe domains along the surface step terraces. The magnetization in the LSMO films was found to be selectively modulated by the antiferrodistortive phase transition of the SrTiO 3 substrate. This phenomenon has been qualitatively explained by a strain modified Stoner–Wohlfarth model. We conclude that the magnetization modulation by the SrTiO 3 phase transition depends on h, the ratio of applied magnetic field to the saturation field. Such modulation is only visible with h < 1. The established domain microstructure–anisotropy–magnetism correlation in manganite films can be applied to a variety of complex oxide thin films on vicinal substrates.

36 MATERIALS SCIENCE↗

Constraints on early dark energy from the axion weak gravity conjecture

A popular proposal for resolving the Hubble tension involves an early phase of dark energy, driven by an axion field with a periodic potential. Here, in this paper, we argue that these models are tightly constrained by the axion weak gravity conjecture: for typical parameter values, the axion decay constant must satisfy f < 0.008,M Pl , which is smaller than the axion decay constants appearing in the vast majority of early dark energy models to date. We discuss possible ways to evade or loosen this constraint, arguing that its loopholes are small and difficult to thread. This suggests that it may prove challenging to realize early dark energy models in a UV complete theory of quantum gravity.

string theory and cosmology↗

Robustness of Vacancy-Bound Non-Abelian Anyons in the Kitaev Model in a Magnetic Field

Non-Abelian anyons in quantum spin liquids (QSLs) provide a promising route to fault-tolerant topological quantum computation. In the exactly solvable Kitaev honeycomb model, such anyons of the QSL state can be bound to nonmagnetic spin vacancies and endowed with non-Abelian statistics by an infinitesimal magnetic field. Here, we investigate how this approach for stabilizing non-Abelian anyons extends to a finite magnetic field represented by a proper Zeeman term. Through large-scale density-matrix renormalization group simulations, we compute the vacancy-anyon binding energy as a function of magnetic field for both the ferromagnetic and antiferromagnetic Kitaev models. Here, we find that anyon binding remains robust within the entire QSL phase for the ferromagnetic Kitaev model but breaks down already inside this phase for the antiferromagnetic Kitaev model. To compute a binding energy several orders of magnitude below the magnetic energy scale, we introduce both a refined definition and an extrapolation scheme based on carefully tailored perturbations.

Xiao, Bo [Oak Ridge National Laboratory (ORNL), Oa↗

Potential Vorticity Mixing in a Tangled Magnetic Field

A theory of potential vorticity (PV) mixing in a disordered (tangled) magnetic field is introduced. The analysis is in the context of β -plane MHD, with a special focus on the physics of momentum transport in the stably stratified, quasi-2D solar tachocline. A physical picture of mean PV evolution by vorticity advection and tilting of magnetic fields is proposed. In the case of weak field perturbations, quasi-linear theory predicts that the Reynolds and magnetic stresses balance as turbulence Alfvénizes for a larger mean magnetic field. Jet formation is explored quantitatively in the mean field–resistivity parameter space. However, since even a modest mean magnetic field leads to large magnetic perturbations for large magnetic Reynolds number, the physically relevant case is that of a strong but disordered field. We show that numerical calculations indicate that the Reynolds stress is modified well before Alfvénization—i.e., before fluid and magnetic energies balance. To understand these trends, a double-average model of PV mixing in a stochastic magnetic field is developed. Calculations indicate that mean-square fields strongly modify Reynolds stress phase coherence and also induce a magnetic drag on zonal flows. The physics of transport reduction by tangled fields is elucidated and linked to the related quench of turbulent resistivity. We propose here a physical picture of the system as a resisto-elastic medium threaded by a tangled magnetic network. Applications of the theory to momentum transport in the tachocline and other systems are discussed in detail.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of turbulence-induced refraction of lower hybrid waves using synthetic scrape-off layer filaments

Here, turbulence-induced refraction effects to lower hybrid (LH) wave propagation and current drive are studied using synthetic scrape-off layer (SOL) blob/filament fields. A synthetic 3D, field-following, blob turbulence model is implemented in the ray-tracing/Fokker-Planck (RTFP) codes GENRAY/CQL3D. In Alcator C-Mod, the blob field is shown to significantly affect LH ray-trajectories, leading to increased on-axis damping and smoother current profiles. This effect depends on the average blob size and amplitude. In addition, the diffusion of ray-trajectories in phase-space caused by turbulence increases the robustness of the RTFP model. A modified $N_{||}$ launch spectrum, acting as a proxy for parametric decay instability (PDI) effects, is included in simulations with the blob model. A synergy between the modified launch spectrum and turbulence-induced refraction results in synthetic hard x-ray profiles that agree with experiment. Lastly, the blob model is used to predict the effect of SOL turbulence on DIII-D high-field side (HFS) LH launch. Assuming low turbulence amplitude in the HFS SOL (~5%), turbulence-induced refraction is predicted to have little effect on current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Distinct critical behaviors from the same state in quantum spin and population dynamics perspectives

There is a deep connection between the ground states of transverse-field spin systems and the late-time distributions of evolving viral populations—within simple models, both are obtained from the principal eigenvector of the same matrix. However, that vector is the wave-function amplitude in the quantum spin model, whereas it is the probability itself in the population model. We show that this seemingly minor difference has significant consequences: Phase transitions that are discontinuous in the spin system become continuous when viewed through the population perspective, and transitions that are continuous become governed by new critical exponents. We introduce a more general class of models that encompasses both cases and that can be solved exactly in a mean-field limit. Numerical results are also presented for a number of one-dimensional chains with power-law interactions. We see that well-worn spin models of quantum statistical mechanics can contain unexpected new physics and insights when treated as population-dynamical models and beyond, motivating further studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Non-volatile electric-field control of inversion symmetry

Competition between ground states at phase boundaries can lead to significant changes in properties under stimuli, particularly when these ground states have different crystal symmetries. A key challenge is to stabilize and control the coexistence of symmetry-distinct phases. Using BiFeO 3 layers confined between layers of dielectric TbScO 3 as a model system, we stabilize the mixed-phase coexistence of centrosymmetric and non-centrosymmetric BiFeO 3 phases at room temperature with antipolar, insulating and polar semiconducting behaviour, respectively. Application of orthogonal in-plane electric (polar) fields results in reversible non-volatile interconversion between the two phases, hence removing and introducing centrosymmetry. Counterintuitively, we find that an electric field 'erases' polarization, resulting from the anisotropy in octahedral tilts introduced by the interweaving TbScO 3 layers. Consequently, this interconversion between centrosymmetric and non-centrosymmetric phases generates changes in the non-linear optical response of over three orders of magnitude, resistivity of over five orders of magnitude and control of microscopic polar order. Here our work establishes a platform for cross-functional devices that take advantage of changes in optical, electrical and ferroic responses, and demonstrates octahedral tilts as an important order parameter in materials interface design.

36 MATERIALS SCIENCE↗

Multidisciplinary design optimization in Architecture, Engineering, and Construction: a detailed review and call for collaboration

The design of buildings has become a complex and multidisciplinary problem involving multiple conflicting objectives as architects and designers address competing technical, economic, environmental, and societal concerns. This has been driving research in Architecture, Engineering, and Construction (AEC) toward rigorous multidisciplinary decision-making frameworks that generate and evaluate numerous design alternatives using multi-objective optimization in concert with simulation and analysis models of varying fidelity and computational expense. While such frameworks are well known and widely employed in the aerospace and systems engineering domains, efforts by design professionals and researchers in the AEC field are scattered at best. In this paper, we provide a detailed review of recent developments in optimization frameworks in the AEC field and subsequently highlight how such developments are largely compartmentalized into separate domains such as structural, energy, daylighting, and other performance factors. Herein, we further discuss the technical challenges involved in concurrent coupled multidisciplinary design optimization (MDO) in the AEC field such as interoperability issues between Building Information Modeling (BIM) environments, analysis/simulation tools, and optimization frameworks. We conclude by outlining research needed for more unified modeling and simulation-based optimization frameworks to aid in complex and multidisciplinary building design processes. We also highlight the need for the identification and development of multi-fidelity simulation tools for use across multiple design phases. As such, this paper contributes a novel roadmap to leverage aerospace and systems engineering research in MDO into the field of AEC, along with a call for researchers in the MDO community to seek collaborations in AEC field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Combined Effect of the Wegener–Bergeron–Findeisen Mechanism and Large Eddies on Microphysics of Mixed-Phase Stratiform Clouds

Abstract The process of glaciation in mixed-phase stratiform clouds was investigated by a novel Lagrangian–Eulerian model (LEM) in which thousands of adjoining Lagrangian parcels moved within a turbulent-like velocity field with statistical parameters typical of the Arctic boundary layer. We used detailed bin microphysics to describe the condensation/evaporation processes in each parcel, in which droplets, aerosols, and ice particles were described using size distributions of 500 mass bins. The model also calculated aerosol mass inside droplets and ice particles. Gravitational sedimentation of droplets and ice particles was also accounted for. Assuming that droplet freezing is the primary source of ice particles, the Arctic clouds observed in Indirect and Semi-Direct Aerosol Campaign (ISDAC) were successfully simulated. The model showed that at a low ice particle concentration typical of ISDAC, large vortices (eddies) led to a quasi-stationary regime, in which mixed-phase St existed for a long time. The large eddies controlled the water partitioning in the mixed-phase clouds. Droplets formed and grew in updrafts, typically reaching the cloud top, and evaporated in downdrafts. Ice particles grew in updrafts and downdrafts. The Wegener–Bergeron–Findeisen (WBF) mechanism was efficient in downdrafts and some parts of updrafts, depending on ice concentration and vertical velocity. At low ice concentrations, the effect of ice on the phase partitioning was negligible. In this regime, liquid droplets were found near the cloud top, whereas ice particles precipitated through the cloud base. When ice concentration exceeded about 10 L −1 , the WBF mechanism led to glaciation of almost the entire cloud, with the exception of narrow cloud regions associated with strong updrafts. At ice particle concentrations of a few tens per liter, the oscillatory regime took place due to the ice–liquid interaction. The microphysical structure of mixed-phase St forms as a combined effect of cloud dynamics (large eddies) and the WBF mechanism.

54 ENVIRONMENTAL SCIENCES↗

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

Enhancing superconductivity of Y 5 Rh 6 Sn 18 by atomic disorder

Here, we investigate the effect of enhancement of superconducting transition temperature $T_c$ by nonmagnetic atom disorder in the filled skutterudite-related tetragonal Y 5 Rh 6 Sn 18 compound doped with Ca. We documented experimentally that Y 5 Rh 6 Sn 18 and its Ca-doped alloys are electronically inhomogeneous at the nanoscale, when Ca content is smaller than ~1.5 in the Y 5– x Ca x Rh 6 Sn 18 system, while for $\textit{x}$ ≥ 1.5 much stronger chemical phase inhomogeneity is observed with an about 20% volume fraction of the second 3:4:13 cubic phase, which we interpret as a fluctuation of stoichiometry within the bulk sample. Then the enhancement of $T_c$ vs $\textit{x}$ could be modeled by a mechanism proposed in recent theoretical reports for increasing the mean-field transition temperature $T_c$ in the presence of nonmagnetic disorder. The increase in disorder with doping increases the sample inhomogeneity and causes a systematic increase in $T_c$, while for two-phase $\textit{x}$ ≥ 1.5 samples the critical temperature $T^{\star}_c ≈ 2 × T_c$ rapidly increases. Based on band structure calculations performed under pressure, we demonstrate how the change in density of states would affect $T_c$ of Y 5 Rh 6 Sn 18 . We obtained the Grüneisen parameter $γ_G$ larger for the inhomogeneous phase with respect to $γ_G$ of the bulk $T_c$ phase and attribute the enhancement of $T_c$ to larger stiffening of the high-temperature $T^{\star}_c$ phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ship-based Observations and Climate Model Simulations of Cloud Phase over the Southern Ocean

The Southern Ocean (SO) clouds exert a significant influence on the Earth’s radiation budget. Here we analyzed ship-based remote sensing observations of SO clouds over a five-month long DOE ARM Measurements of Aerosols, Radiation and Clouds over the Southern Oceans (MARCUS) field campaign to better understand cloud phase variability. We developed a method to classify eight categories of hydrometeors (ice, liquid, mixed phase, rain, drizzle, snow, aerosols, and clear sky) based on measurements of lidar, radar and radiosondes. Cloud thermodynamic phases (liquid, ice and mixed phase) at coarser scales were further derived to compare with the DOE Energy Exascale Earth System Model version 1 (E3SMv1) simulation and the Earth Model Column Collaboratory (EMC 2 ) instrument simulator. For a scale-aware comparison with climate model simulations, we found that spatially averaging the raw remote sensing data (e.g., backscatter, reflectivity) results in increased cloud cover and cloud liquid, whereas, directly averaging cloud phase from higher to lower resolution maintains clear air regions and is thus recommended for in-cloud frequency comparisons. For cloud thermodynamic phases in stratiform clouds, the E3SM underestimates cloud ice and overestimates cloud liquid at temperatures between -40 and 0°C. When latitudes increase, both observations and simulations show a transition of dominant phase from liquid to ice for cloud tops as well as for the entire cloud columns, but the model underestimation of ice phase is more severe at higher latitudes. Such model bias is unlikely caused by spatial scale differences or lack of heterogeneity in cloud vertical structure in the simulation.

54 ENVIRONMENTAL SCIENCES↗

Quarkyonic mean field theory

We discuss mean field theory of quarkyonic matter at zero temperature. We treat the nucleons with contact interactions in mean field approximation, discussing both vector and scalar mean field interactions. We treat the quarks without mean field vector interactions, but allow mass terms to be generated consistently from a scalar mean field consistent with the additive quark model for quark masses. Quarkyonic matter is composed of a shell of nucleons that under-occupy the total available phase space associated with the underlying quark degrees of freedom. Here, the fully occupied Fermi sphere beneath this shell of nucleons at high densities is thought of as quarks, but when this fully occupied distribution of states first appears, although the phase space is filled, the matter is at low density. For the transition between this low density and high density saturated matter, we advocate a dual description of the fully filled Fermi sea in terms of hadrons, and make a phenomenological hypothesis for the equation of state of this matter. We then proceed to an example where the mean field interactions are all vector and only associated with the nucleons, ignoring the effects of mass change associated with the scalar interactions. Except for the effects of Pauli blocking, the nucleons and quarks do not interact. To get a reasonable transition to quarkyonic matter the interaction of the quarks among themselves are assumed to be nonperturbative, and a simple phenomenological relation between quark Fermi energy and density is introduced.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗

A supersymmetric SYK model with a curious low energy behavior

We consider N = 2,4 supersymmetric SYK models that have a peculiar low energy behavior, with the entropy going like S = S 0 + (constant)T a , where a ≠ 1. The large N equations for these models are a generalization of equations that have been previously studied as an unjustified truncation of the planar diagrams describing the BFSS matrix quantum mechanics or other related matrix models. Here we reanalyze these equations in order to better understand the low energy physics of these models. We find that the scalar fields develop large expectation values which explore the low energy valleys in the potential. The low energy physics is dominated by quadratic fluctuations around these values. These models were previously conjectured to have a spin glass phase. We did not find any evidence for this phase by using the usual diagnostics, such as searching for replica symmetry breaking solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗