Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “order”

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 1,639 records · Page 91

Untangling charge-order dependent bulk states from surface effects in a topological kagome metal ScV 6 Sn 6

Kagome metals with charge density wave (CDW) order exhibit a broad spectrum of intriguing quantum phenomena. The recent discovery of the novel kagome CDW compound ScV 6 Sn 6 has spurred significant interest. However, understanding the interplay between CDW and the bulk electronic structure has been obscured by a profusion of surface states and terminations in this quantum material. Here, in this study, we employ photoemission spectroscopy and potassium dosing to elucidate the complete bulk band structure of ScV 6 Sn 6 , revealing multiple van Hove singularities near the Fermi level. We surprisingly discover a robust spin-polarized topological Dirac surface resonance state at the M point within the twofold van Hove singularities. Assisted by first-principles calculations, the temperature dependence of the k z -resolved angle-resolved photoemission spectroscopy spectrum provides unequivocal evidence for the proposed $\sqrt{3}$×$\sqrt{3}$×3 charge order over other candidates. Our work not only enhances the understanding of the CDW-dependent bulk and surface states in ScV 6 Sn 6 , but also establishes an essential foundation for potential manipulation of the CDW order in kagome materials.

36 MATERIALS SCIENCE↗

Coupling-dependent antiferromagnetic-ferromagnetic ordering in a pinwheel artificial spin ice

Nanopatterned magnetic thin films offer a platform for exploration of tailored magnetic properties such as emergent long-range order. A prominent example is artificial spin ice (ASI), where an arrangement of nanoscale magnetic elements, acting as macrospins, interact via their dipolar fields. In this study, we discuss the transition from antiferromagnetic (AFM) to ferromagnetic (FM) long-range order in a square lattice ASI as the magnetic elements are gradually rotated through 45⁢° to a “pinwheel” configuration. The AFM-FM transition is observed experimentally using synchrotron radiation x-ray spectromicroscopy and occurs for a certain rotation angle of the nanomagnets, dependent on the dipolar coupling strength determined by the separation of the magnets in the lattice. Large-scale magnetic dipole simulations show that the point-dipole approximation fails to capture the correct AFM-FM transition angle. However, excellent agreement with experimental data is obtained using a dumbbell-dipole model, which better reflects the actual dipolar fields of the magnets. This model also explains the coupling dependence of the transition angle, another feature not captured by the point-dipole model. Our findings resolve a discrepancy between measurement and theory in previous work on “pinwheel” ASIs and establish the coupling dependence of the AFM-FM transition. The revised dipole model, with a more accurate representation of the stray field, offers more precise control of magnetic order in artificial spin systems.

Artificial spin ice↗

Structural and metal-insulator transitions in rhenium-based double perovskites via orbital ordering

Re-based double perovskites (DPs) have garnered substantial attention due to their high Curie temperatures (T C ) and display of complex interplay of structural and metal-insulator transitions (MIT). In this report we systematically study the ground-state electronic and structural properties for a family of Re-based DPs A 2 BReO 6 (A=Sr, Ca and B=Cr, Fe), which are related by a common low-energy Hamiltonian, using density functional theory +U calculations. We show that the on-site interaction U of Re induces orbital ordering (denoted C-OO), with each Re site having an occupied d xy orbital and a C-type alternation among d xz /d yz , resulting in an insulating state consistent with experimentally determined insulators Sr 2 CrReO 6 , Ca 2 CrReO 6 , and Ca 2 FeReO 6 . The threshold value of U Re for orbital ordering is reduced by inducing E g octahedral distortions of the same C-type wavelength (denoted C-OD), which serves as a structural signature of the orbital ordering; octahedral tilting also reduces the threshold. The C-OO and the concomitant C-OD are a spontaneously broken symmetry for the Sr-based materials (i.e., a 0 a 0 c - tilt pattern), while not for the Ca-based systems (i.e., a - a - b + tilt pattern). Spin-orbit coupling does not qualitatively change the physics of the C-OO/C-OD, but can induce relevant quantitative changes. We prove that a single set of U Cr , U Fe , U Re capture the experimentally observed metallic state in Sr 2 FeReO 6 and insulating states in other three systems. We predict that the C-OO is the origin of the insulating state in Sr 2 CrReO 6 , and that the concomitant C-OD may be experimentally observed at sufficiently low temperatures (i.e., space group P4 2 /m) in pure samples. Additionally, given our prescribed values of U, we show that the C-OO induced insulating state in Ca 2 CrReO 6 will survive even if the C-OD amplitude is suppressed (e.g., due to thermal fluctuations). The role of the C-OO/C-OD in the discontinuous, temperature driven MIT in Ca 2 FeReO 6 is discussed.

36 MATERIALS SCIENCE↗

Local-order metric for condensed-phase environments

Here we introduce a local order metric (LOM) that measures the degree of order in the neighborhood of an atomic or molecular site in a condensed medium. The LOM maximizes the overlap between the spatial distribution of sites belonging to that neighborhood and the corresponding distribution in a suitable reference system. The LOM takes a value tending to zero for completely disordered environments and tending to one for environments that perfectly match the reference. The site-averaged LOM and its standard deviation define two scalar order parameters, S and δS, that characterize with excellent resolution crystals, liquids, and amorphous materials. We show with molecular dynamics simulations that S, δS, and the LOM provide very insightful information in the study of structural transformations, such as those occurring when ice spontaneously nucleates from supercooled water or when a supercooled water sample becomes amorphous upon progressive cooling.

36 MATERIALS SCIENCE↗

Correlations and electronic order in a two-orbital honeycomb lattice model for twisted bilayer graphene

The recent observation of superconductivity in proximity to an insulating phase in twisted bilayer graphene (TBG) at small “magic” twist angles has been linked to the existence of nearly flat bands, which make TBG a fresh playground to investigate the interplay between correlations and superconductivity. The low-energy narrow bands were shown to be well described by an effective tight-binding model on the honeycomb lattice (the dual of the triangular Moiré superlattice) with a local orbital degree of freedom. In this paper, we perform a strong-coupling analysis of the proposed (p x ,p y ) two-orbital extended Hubbard model on the honeycomb lattice. By decomposing the interacting terms in the particle-particle and particle-hole channels, we classify the different possible superconducting, magnetic, and charge instabilities of the system. In the pairing case, we pay particular attention to the two-component (d–wave) pairing channels, which admit vestigial phases with nematic or chiral orders, and study their phenomenology. Furthermore, we explore the strong-coupling regime by obtaining a simplified spin-orbital exchange model which may describe a putative Mott-type insulating state at quarter-filling. Our mean-field solution reveals a rich intertwinement between ferromagnetic and antiferromagnetic orders with different types of nematic and magnetic orbital orders. Altogether, our work provides a solid framework for further investigations of the phase diagram of the two-orbital extended Hubbard model in both strong- and weak-coupling regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Charge-Order-Induced Ferroelectricity in LaVO 3 / Sr VO 3 Superlattices

In this paper, the structure and properties of the $1{\mathbin:}1$ superlattice of ${\mathrm{LaVO}}_{3}$ and ${\mathrm{SrVO}}_{3}$ are investigated with a first-principles density-functional-theory-plus-$U$ ($\mathrm{DFT}+U$) method. The lowest energy states are antiferromagnetic charge-ordered Mott-insulating phases. In one of these insulating phases, layered charge ordering combines with the layered La/Sr cation ordering to produce a polar structure with a large nonzero spontaneous polarization normal to the interfaces. This polarization, comparable to that of conventional ferroelectrics, is produced by electron transfer between the ${\mathrm{V}}^{3+}$ and ${\mathrm{V}}^{4+}$ layers. The energy of this normal-polarization state relative to the ground state is only 3 meV per vanadium. Under tensile strain, this energy difference can be further reduced, suggesting that the normal-polarization state can be induced by an electric field applied normal to the superlattice layers, yielding an antiferroelectric double-hysteresis loop. Finally, if the system does not switch back to the ground state on removal of the field, a ferroelectric-type hysteresis loop could be observed.

36 MATERIALS SCIENCE↗

Suppressed Fluctuations as the Origin of the Static Magnetic Order in Strained Sr 2 RuO 4

Combining first-principles density-functional calculations and Moriya's self-consistent renormalization theory, we explain the recently reported counterintuitive appearance of an ordered magnetic state in uniaxially strained Sr 2 RuO 4 beyond the Lifshitz transition. We show that strain weakens the quantum spin fluctuations, which destroy the static order, more strongly than the tendency to magnetism. Furthermore, a different rate of decrease of the spin fluctuations vs magnetic stabilization energy promotes the onset of a static magnetic order beyond a critical strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Antiferromagnetic order in Co-doped Fe 5 GeTe 2 probed by resonant magnetic x-ray scattering

The quasi-two-dimensional van der Waals magnet Fe 5-δ GeTe 2 has emerged as a promising platform for electronic and spintronic functionalities at room temperature, owing to its large ferromagnetic ordering temperature T C ≈ 315 K. Interestingly, by cobalt (Co) substitution of iron in F5GT, i.e., (Fe 1-x Co x ) 5-δ GeTe 2 (Co-F5GT), not only can its magnetic transition temperature be further enhanced, but the magnetic and structural ground states can also be tuned. Specifically, an antiferromagnetic (AFM) order is induced beyond the Co doping level x ≥ 0.4. Here, in this work, we investigate the magnetic properties of a Co-F5GT single crystal at x = 0.45(1), by utilizing the element-specific, resonant magnetic x-ray scattering technique. Our study reveals an A-type, Ising-like AFM ground state, with a transition temperature T N ≈ 340 K. In addition, our work unveils an important contribution from Co magnetic moments to the magnetic order. The application of the in-plane magnetic fields gradually polarizes the spin moments along the field direction, but without inducing incommensurate spin texture(s).

36 MATERIALS SCIENCE↗

Efficient Scalable Parallel Higher Order Direct MoM-SIE Method With Hierarchically Semiseparable Structures for 3-D Scattering

A novel fast scalable parallel algorithm is proposed for the solution of large 3-D scattering problems based on: 1) the double (geometrical and current-approximation) higher order (DHO) method of moments (MoM) in the surface integral equation (SIE) formulation and 2) a direct solver for dense linear systems utilizing hierarchically semiseparable (HSS) structures. Namely, an HSS matrix representation is used for compression, factorization, and solution of the system matrix. In addition, a rank-revealing QR decomposition for memory compression is used, with a stopping criterion in terms of the relative rank tolerance value. One method for geometrical preprocessing of the scatterers based on the cobblestone distance sorting technique is employed in order to enhance the HSS algorithm accuracy and parallelization. Numerical examples show how the accuracy of the DHO HSS-MoM-SIE method is easily controllable by using the relative tolerance for the matrix compression. Furthermore, the examples demonstrate low memory consumption, as well as much faster simulation time, when compared to the direct LU decomposition. The method enables dramatically faster monostatic scattering computations than iterative solvers and reduced number of unknowns when compared to low-order discretizations. Finally, great scalability of the algorithm is demonstrated on more than one thousand processes.

97 MATHEMATICS AND COMPUTING↗

Network-Level Traffic Signal Cooperation: A Higher-Order Conflict Graph Approach

Traffic signal control and cooperation are extremely important to alleviate traffic congestion in a large traffic network. This study develops a higher-order conflict graph approach for network-wide traffic signal control and cooperation. A conflict graph is applied to model the traffic signal configurations, which identifies the conflict and unconflicted movements for each intersection. In conflict graph, the node represents each movement. The weight of each node can be defined as traffic volume, queue length, fuel consumption, or any weighted combinations of these measurements. The calculation of the optimal green light duration and green light sequence (for different movements) is equivalent to sequentially finding the maximum weight independent set (MWIS) in the conflict graph. The conflict graph also provides a uniform and efficient way to connect traffic signal operations among nearby intersections spatially. Then, we introduced the concept of the k -th order neighborhood to model the degree of connectivity between each movement to the movements at upstream or downstream intersections. The weight of each node in the higher-order conflict graph not only represents its own congestion level, but also relates to the traffic conditions of nearby intersections. Through this approach, the cooperation of multiple intersections can be realized by incorporating their spatial connectivity into conflict graph and solving the MWIS problem. A simulation network is built in SUMO to test the effectiveness of the proposed method. Results suggested that the proposed model outperformed other state-of-the-art signal control methods. Also, the scheme maintains good performance under varying traffic demands.

42 ENGINEERING↗

Reduced Order Model of Transactive Bidding Loads

Transactive energy (TE) has been identified to provide better grid efficiency and reliability by market-based transactive exchanges between energy producers and energy consumers. Simulations of TE systems are crucial to evaluate the benefits and impacts of different transactive mechanisms. However, such simulations can be time consuming due to the information exchange between various participants and complex co-simulation environments. In this paper, we develop a reduced order model to speed up the simulation of transactive systems in TE simulation platform (TESP) while achieving very low error between the reduced order and full model results. Specifically, the developed reduced order model consists of an aggregate responsive load agent which utilizes two Recurrent Neural Networks (RNNs) with Long Short-Term Memory units (LSTMs) to enable transactive elements to collectively participate in the TE system. The proposed aggregate responsive load (ARL) agent is able to produce similar transactive behaviors to the full simulation model while achieving significant simulation time reduction. Finally, we also show that the developed model enables generalization of simulation results across different dates and across different number of loads included in the simulations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

High-Order Multi-Material ALE Hydrodynamics

Here, we present a new approach for multi-material arbitrary Lagrangian--Eulerian (ALE) hydrodynamics simulations based on high-order finite elements posed on high-order curvilinear meshes. The method builds on and extends our previous work in the Lagrangian and remap phases of ALE, and depends critically on a functional perspective that enables subzonal physics and material modeling. Curvilinear mesh relaxation is based on node movement, which is determined through the solution of an elliptic equation. The remap phase is posed in terms of advecting state variables between two meshes over a fictitious time interval. The resulting advection equation is solved by a discontinuous Galerkin (DG) formulation, combined with a customized Flux Corrected Transport (FCT) type algorithm. Because conservative fields are remapped, additional synchronization steps are introduced to preserve bounds with respect to primal fields. These steps include modification of the low-order FCT solutions, definition of conservative FCT fluxes based on primal field bounds, and monotone transitions between primal and conservative fields.

97 MATHEMATICS AND COMPUTING↗

SNS: A Solution-Based Nonlinear Subspace Method for Time-Dependent Model Order Reduction

Several reduced order models have been successfully developed for nonlinear dynamical systems. To achieve a considerable speed-up, a hyper-reduction step is needed to reduce the computational complexity due to nonlinear terms. Many hyper-reduction techniques require the construction of nonlinear term basis, which introduces a computationally expensive offline phase. A novel way of constructing nonlinear term basis within the hyper-reduction process is introduced. In contrast to the traditional hyper-reduction techniques where the collection of nonlinear term snapshots is required, the SNS method avoids collecting the nonlinear term snapshots. Instead, it uses the solution snapshots that are used for building a solution basis, which enables avoiding an extra data compression of nonlinear term snapshots. As a result, the SNS method provides a more efficient offline strategy than the traditional model order reduction techniques, such as the DEIM, GNAT, and ST-GNAT methods. The SNS method is theoretically justified by the conforming subspace condition and the subspace inclusion relation. It is useful for model order reduction of large-scale nonlinear dynamical problems to reduce the offline cost. It is especially useful for ST-GNAT that has shown promising results, such as a good accuracy with a considerable online speed-up for hyperbolic problems in a recent paper by Choi and Carlberg [SIAM J. Sci. Comput., 41 (2019), pp. A26--A58], because ST-GNAT involves an expensive offline cost related to collecting nonlinear term snapshots. Error analysis for the SNS method is presented. Numerical results support that the accuracy of the solution from the SNS method is comparable to the traditional methods and a considerable speed-up (i.e., a factor of two to a hundred) is achieved in the offline phase.

97 MATHEMATICS AND COMPUTING↗

Conservative High-Order Time Integration for Lagrangian Hydrodynamics

In tis study, we develop novel time integration methods for the compressible Euler equations in the Lagrangian frame that are of arbitrary high order and exactly preserve the mass, momentum, and total energy of the system. The equations are considered in nonconservative form, that is, common for staggered grid hydrodynamics (SGH) methods; namely, the evolved quantities are mass, momentum, and internal energy. A general family of time integration schemes is formulated, and practical pairs for orders three and four are derived. Numerical results on standard hydrodynamics benchmarks confirm the high-order convergence on smooth problems and the exact numerical preservation of all physically conserved quantities.

97 MATHEMATICS AND COMPUTING↗

An adaptive multiblock high-order finite-volume method for solving the shallow-water equations on the sphere

We present a high-order finite-volume approach for solving the shallow-water equations on the sphere, using multiblock grids on the cubed-sphere. This approach combines a Runge--Kutta time discretization with a fourth-order accurate spatial discretization, and includes adaptive mesh refinement and refinement in time. Results of tests show fourth-order convergence for the shallow-water equations as well as for advection in a highly deformational flow. Hierarchical adaptive mesh refinement allows solution error to be achieved that is comparable to that obtained with uniform resolution of the most refined level of the hierarchy, but with many fewer operations.

97 MATHEMATICS AND COMPUTING↗

A fourth-order Cartesian grid embeddedboundary method for Poisson’s equation

In this paper, we present a fourth-order algorithm to solve Poisson's equation in two and three dimensions. We use a Cartesian grid, embedded boundary method to resolve complex boundaries. We use a weighted least squares algorithm to solve for our stencils. We use convergence tests to demonstrate accuracy and we show the eigenvalues of the operator to demonstrate stability. We compare accuracy and performance with an established second-order algorithm. We also discuss in depth strategies for retaining higher-order accuracy in the presence of nonsmooth geometries.

97 MATHEMATICS AND COMPUTING↗

Emissions minimization on road networks via Generic Second Order Models

In this paper we consider the problem of estimating emissions due to vehicular traffic on complex networks, and minimizing their effect by regulating traffic at junctions. For the traffic evolution, we consider a Generic Second Order Model, which encompasses the majority of two-equations (i.e., second-order) models available in the literature, and extend it to road networks with merge and diverge junctions. The dynamics on the whole network is determined by selecting a solution to the Riemann Problems at junctions, i.e., the Cauchy problems with constant initial data on each incident road. The latter are solved by assuming the maximization of the flow and assigning a traffic distribution coefficient for outgoing roads of diverges, and a priority rule for incoming roads of merges. A general emission model is considered and its parameters are tuned to the $ {\mathrm{NO_{x}}} $ emission rate. The minimization of emissions is then formulated in terms of the traffic distribution and priority parameters, taking into account travel times. A comparison is provided between roundabouts with optimized parameters and traffic lights, which correspond to time-varying traffic priorities. Our approach can be adapted to manage traffic in complex networks in order to reduce emissions while keeping travel time at acceptable levels.

33 ADVANCED PROPULSION SYSTEMS↗