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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2023

Complete lowest order radiative corrections in semi-inclusive scattering of polarized particles

For this work, the lowest order radiative corrections to the cross section and asymmetries measured in experiments on semi-inclusive deep inelastic scattering of polarized particles were calculated. Both exact and leading log expressions were presented and discussed for the total correction that include the contributions from the processes of i) real photon emission with semi-inclusive processes, ii) loop diagrams, and iii) real photon emission with exclusive processes. Radiative corrections to the Sivers and Collins asymmetries in $\pi^+$ electroproduction were studied numerically within the kinematical conditions of modern experimental environments at Jefferson Laboratory (JLab). The Wandzura-Wilczek approximation for the semi-inclusive structure functions and MAID2007 parameterization for the six amplitudes of exclusive processes were used in numeric analyses. The results show that i) radiative effects can generate a correction comparable to the size of Sivers and Collins asymmetries at the Born level, ii) there is good agreement between the exact and leading-order corrections, iii) external functions (that is, other than the Sivers and Collins functions in the respective asymmetries) can generate a contribution to the radiative correction up to 20%, and iv) there exists a strong dependence of the radiative correction on the models for semi-inclusive and exclusive structure functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hypernetwork science via high-order hypergraph walks

We propose high-order hypergraph walks as a framework to generalize graph-based network science techniques to hypergraphs. Edge incidence in hypergraphs is quantitative, yielding hypergraph walks with both length and width. Graph methods which then generalize to hypergraphs include connected component analyses, graph distance-based metrics such as closeness centrality, and motif-based measures such as clustering coefficients. We apply high-order analogs of these methods to real world hypernetworks, and show they reveal nuanced and interpretable structure that cannot be detected by graph-based methods. Lastly, we apply three generative models to the data and find that basic hypergraph properties, such as density and degree distributions, do not necessarily control these new structural measurements. Our work demonstrates how analyses of hypergraph-structured data are richer when utilizing tools tailored to capture hypergraph-native phenomena, and suggests one possible avenue towards that end.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗
DOE OSTI2022

ROMNet: Neural-Network-Based Surrogates for Reduced-Order Model Dynamics

SAND2023-05829O ROMNet is a toolbox built on the Tensorflow library for training reduced-order, neural operator surrogate models. The software is used for training reduced order neural operator surrogate models for dynamical systems applications. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗
DOE OSTI2024

Nonlinear manifold reduced order model

Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a lowdimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physics-informed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LS-ROMs. Our software takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers' equations. A speedup of up to 2.6 for 1D Burgers' and a speedup of 11.7 for 2D Burgers' equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique.

Choi, Youngsoo↗
DOE OSTI2021

Simultaneous type-I and type-II phase matching for second-order nonlinearity in integrated lithium niobate waveguide

Second-order optical nonlinearity is widely used for both classical and quantum photonic applications. Due to material dispersion and phase matching requirements, the polarization of optical fields is pre-defined during the fabrication. Only one type of phase matching condition is normally satisfied, and this limits the device flexibility. Here, we demonstrate that phase matching for both type-I and type-II second-order optical nonlinearity can be realized simultaneously in the same waveguide fabricated from thin-film lithium niobate. This is achieved by engineering the geometry dispersion to compensate for the material dispersion and birefringence. The simultaneous realization of both phase matching conditions is verified by the polarization dependence of second-harmonic generation. Correlated photons are also generated through parametric down conversion from the same device. This work provides a novel approach to realize versatile photonic functions with flexible devices.

Briggs, Ian↗
DOE OSTI2024

Enhancement of harmonic generation by an intense driving laser with high-order waveguide modes in a high-pressure gas-filled hollow waveguide

High-order harmonics have been widely used as reliable tabletop coherent radiation sources recently, but their applications have often been limited by the available pulse energy. Here, we report that by using an overdriven intense laser in a long waveguide with high-pressure gas, phase matching can be achieved in three distinct “regimes”. In the third regime, favorable phase matching is achieved at near-axis positions to enhance harmonic yields. Our results are supported by a full theoretical analysis, and we demonstrate that coupling of the driving laser with the high-order waveguide modes (instead of the fundamental mode used in most prior experiments) is responsible for achieving phase matching. Furthermore, we establish that this phase matching (and harmonic enhancement) is robust, and a scaling relation is derived for the necessary waveguide and gas parameters, allowing our predictions to be tested immediately in any laboratory today.

Li, Baochang [Nanjing University of Science and Te↗
DOE OSTI2021

Multifidelity computing for coupling full and reduced order models

Hybrid physics-machine learning models are increasingly being used in simulations of transport processes. Many complex multiphysics systems relevant to scientific and engineering applications include multiple spatiotemporal scales and comprise a multifidelity problem sharing an interface between various formulations or heterogeneous computational entities. To this end, we present a robust hybrid analysis and modeling approach combining a physics-based full order model (FOM) and a data-driven reduced order model (ROM) to form the building blocks of an integrated approach among mixed fidelity descriptions toward predictive digital twin technologies. At the interface, we introduce a long short-term memory network to bridge these high and low-fidelity models in various forms of interfacial error correction or prolongation. The proposed interface learning approaches are tested as a new way to address ROM-FOM coupling problems solving nonlinear advection-diffusion flow situations with a bifidelity setup that captures the essence of a broad class of transport processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Emergent pair density wave order across a Lifshitz transition

We numerically investigate the telltale signs of pair-density-wave order (PDW) in the Kondo-Heisenberg chain by focusing on the momentum resolved spectrum in different parameter regimes. Density matrix renormalization group calculations reveal that this phase is characterized by a dispersion with two minima and four Fermi points, indicating the emergence of an effective next-nearest-neighbor hopping that arises as a second-order effect to avoid magnetic frustration. The pairs appear in the spectrum as in-gap bound states with weight concentrated in the hole pockets. The low-energy physics can be understood by means of a generalized t−J model with next-nearest-neighbor hopping. Our results offer a guide for searching for experimental signatures, and for other models that can realize PDW physics.

Yang, Luhang [University of Tennessee Knoxville]↗
DOE OSTI2022

PMC$_\infty$: Infinite-Order Scale-Setting method using the Principle of Maximum Conformality and preserving the Intrinsic Conformality

We show results for Thrust and C-parameter in e^+ e^- e + e − annihilation to 3 jets obtained using the recently developed new method for eliminating the scale ambiguity and the scheme dependence in pQCD namely the Infinite-Order Scale-Setting method using the Principle of Maximum Conformality (PMC _\infty ∞ ). This method preserves an important underlying property of gauge theories: intrinsic Conformality (iCF). It leads to a remarkably efficient method to eliminate the conventional renormalization scale ambiguity at any order in pQCD. A comparison with Conventional Scale Setting method (CSS) is also shown.

Di Giustino, Leonardo↗
DOE OSTI2020

Nonlinear manifold-based reduced order model

Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations, in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a low-dimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physicsinformed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LSROMs. Our method takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers’ equations. A speedup of up to 2.6 for 1D Burgers’ and a speedup of 11.7 for 2D Burgers’ equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique. Finally, a posteriori error bounds for the NM-ROMs are derived that take account of the hyper-reduced operators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multiscale modeling high-order methods and data-driven modeling

Projection-based reduced-order models (ROMs) comprise a promising set of data-driven approaches for accelerating the simulation of high-fidelity numerical simulations. Standard projection-based ROM approaches, however, suffer from several drawbacks when applied to the complex nonlinear dynamical systems commonly encountered in science and engineering. These limitations include a lack of stability, accuracy, and sharp a posteriori error estimators. This work addresses these limitations by leveraging multiscale modeling, least-squares principles, and machine learning to develop novel reduced-order modeling approaches, along with data-driven a posteriori error estimators, for dynamical systems. Theoretical and numerical results demonstrate that the two ROM approaches developed in this work - namely the windowed least-squares method and the Adjoint Petrov - Galerkin method - yield substantial improvements over state-of-the-art approaches. Additionally, numerical results demonstrate the capability of the a posteriori error models developed in this work.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Primer on FERC Order No. 2222: Insights for International Power Systems

The U.S. Federal Energy Regulatory Commission (FERC) Order No. 2222 is heralded as a landmark rule on distributed energy resources (DER). Issued in September 2020, it is specifically designed for organized wholesale power markets. The rule is nonbinding for utilities that have not chosen to join an organized wholesale power market. This primer, which analyzes the order, focuses on potential insights for similar changes in the international power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Preconditioning for Hyper-reduction in Reduced Order Models

Many projection-based reduced order models (pROM) that utilize the governing equation and data to accelerate physical simulations can be applied to nonlinear dynamical systems. To avoid full order model (FOM) scale update for each time step, hyperreduction techniques are developed to sample high dimensional nonlinear terms. Our study aims to investigate if preconditioning the least-squares problem used for the nonlinear approximation can improve the robustness of the condition number of the problem while achieving high accuracy. In our study, we use the row-normalization matrix motivated by the Christoffel function as the preconditioner and solve the corresponding weighted least-squares problem. Numerical results for Lagrangian hydrodynamics examples are analyzed to explore how the preconditioner works compared to existing hyper-reduction techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Impact of FERC Order 2222 on DER Participation Rules in US Electricity Markets

Electricity markets in the bulk grid are beginning to implement market mechanisms that support the procurement of flexible capabilities from wide range of technologies, including distributed energy resources (DERs). The flexibility of these resources will help counterbalance supply uncertainties from large-scale integration of variable renewable generation. To encourage development of distributed and aggregated market participants, FERC Order 2222 was issued in September 2020 to require each Independent System Operator (ISO) in the US to implement rules that enable broader participation from aggregations of DERs in the bulk market. The following paper first describes the generic design of ISO markets before introducing the new market participation rules that ISOs have proposed for compliance with Order 2222. The paper then describes how software performance issues may continue to affect the eligibility requirements and offer structures for DER aggregations participating in ISOs, noting that continued research on computational methods may help reduce burdens for DER integration. The prospects for transmission and distribution system coordination is second major issue discussed, which will require minor changes to existing processes in the short term. In the longer term, there is more opportunity for more wide-ranging reforms, such as the development of a Distribution System Operator (DSO) framework. Newly proposed market rules may affect how Transactive Energy Systems (TES) will help facilitate efficient formation of DER aggregations and operation of the individual DERs within an aggregation. Within the TES context, the challenge is to fully understand how resource eligibility and operational and planning coordination methods will affect the design and implementation of TES.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

The mapping & avoidance of high order cavity modes

In this paper a technique developed in 2016 is presented for mapping the multidimensional, multi-mode high order mode landscape of a beam driven super conducting RF cavity. The initial use case for this technique was intended to allow longitudinal gradient control of a beam driven super conducting cavity while simultaneously avoiding many dangerous high order modes and excessive component stresses. This paper will focus on the stability and component stress considerations of the 56MHz superconducting cavity, located in the 4 o’clock sector interaction region of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory. The 56MHz cavity is intended to work in conjunction with the Stochastic Cooling system to improve the longitudinal confinement in the central 197MHz RF bucket and thus increase the integrated luminosity during heavy ion beam runs at RHIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Short-Depth QAOA circuits and Quantum Annealing on Higher-Order Ising Models (Rev.2)

The Quantum Alternating Operator Ansatz (QAOA) and Quantum Annealing (QA) are quantum algorithms that are both based on the adiabatic theorem and both have the goal of sampling the optimal solution(s) of combinatorial optimization problems. Quantum annealing has been physically instantiated on D-Wave devices using superconducting flux qubits, and QAOA can be programmed on digital gate-model quantum computers such as the programmable superconducting transmon qubits devices of the IBMQ series, for instance ibm washington. QAOA and QA address the same types of problems, but it is unclear how they will scale to large problem sizes and to larger and higher-fidelity quantum computers. In this article, we present a direct comparison between QAOA, one and two rounds, run on all 127 qubits of ibm washington and QA run on D-Wave Advantage system4.1 and Advantage system6.1. The problems which allow for this comparison are random Ising model problems whose connectivity matches the heavy hexagonal lattice topology of ibm washington and the Pegasus graph connectivity of the two D-Wave devices. We create two classes of problem instances for this comparison: one with higher order terms (ZZZ variable interactions), linear terms, and quadratic terms, and a separate problem type with only linear and quadratic terms. Our QAOA circuits are novel and extremely short depth, with a CNOT depth of 6 per round, which allows whole chip usage of ibm washington’s heavy hexagonal lattice and can be applied to future heavy-hex chips. We also test the effectiveness of the error suppression technique digital dynamical decoupling on the QAOA circuits. The QAOA circuits compiled to ibm washington are composed of several thousand circuit instructions, approximately 3, 000 depending on the details of the circuit, making these some the largest quantum circuits ever executed on a digital quantum processor. QAOA and QA are compared against the classical heuristic algorithm of simulated annealing and all problem instances are exactly solved using CPLEX in order to evaluate which samplers, if any, correctly found the ground state solution(s) of the problem instances. We find that (i) QA outperforms QAOA on all problem instances, (ii) QAOA samples the problems better than random sampling, and (iii) QAOA angle computation exhibits clear parameter concentration across the ensemble of Ising models.

127 Qubits↗