Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analysis of single-mode Richtmyer–Meshkov instability using high-order incompressible vorticity-streamfunction and shock-capturing simulations

In this study, two- and three-dimensional simulation results obtained using a new high-order incompressible, variable-density vorticity–streamfunction (VS) method and data from previous ninth-order weighted essentially nonoscillatory (WENO) shock-capturing simulations are used to investigate the nonlinear dynamics of single-mode Richtmyer–Meshkov instability using a model of a Mach 1.3 air(acetone)/SF6 shock tube experiment. A comparison of the density fields from both simulations with the experimental images demonstrates very good agreement in the large-scale structure with both methods but differences in the small-scale structure. The WENO method captures the small-scale disordered structure observed in the experiment, while the VS method partially captures such structure and yields a strong rotating core. The perturbation amplitude growth from the simulations generally agrees well with the experiment. The simulation bubble and spike amplitudes agree well at early times. At later times, the WENO bubble amplitude is smaller than the VS amplitude and vice versa for the spike amplitude. The predictions of nonlinear single-mode instability growth models are shown to agree with the simulation amplitudes at early-to-intermediate times but underpredict the amplitudes at later times in the nonlinear regime. Visualizations of the mass fraction and enstrophy isosurfaces, velocity and vorticity fields, and baroclinic vorticity production and vortex stretching terms from the three-dimensional simulations indicate that, with the exception of the small-scale structure within the rollups, the VS and WENO results are in good agreement.

42 ENGINEERING↗
DOE OSTI2021

An embedded variable step IMEX scheme for the incompressible Navier–Stokes equations

This report presents a series of implicit–explicit (IMEX) variable stepsize algorithms for the incompressible Navier–Stokes equations (NSE). Here, IMEX means the nonlinear term is treated fully explicitly while the remaining terms are treated implicitly. With the advent of new computer architectures there has been growing demand for low memory solvers of this type. The addition of time adaptivity improves the accuracy and greatly enhances the efficiency of the algorithm. We prove energy stability of an embedded first–second order IMEX pair. For the first order member of the pair, we prove stability for variable stepsizes, and analyze convergence. We believe this to be the first proof of this type for a variable stepsize IMEX scheme for the incompressible NSE. We then define and test a variable stepsize, variable order IMEX scheme using these methods. Our work contributes several firsts for IMEX NSE schemes, including an energy argument and error analysis of a two-step, variable stepsize method, and embedded error estimation for an IMEX multi-step method.

42 ENGINEERING↗
DOE OSTI2026

Coexistence of Synchronization and Stochasticity in Thermally Coupled Mott Oscillators

Synchronization is conventionally regarded as a mechanism for suppressing variability and enforcing order in coupled systems, from pendula and lasers to neurons and electronic oscillators. Here, we show that synchronization can also embed stochasticity at finer scales. We observe this phenomenon in thermally coupled VO 2 neuristors, where robust in-phase synchronization at the microsecond scale coexists with spike onset fluctuations at the nanosecond scale, with no fixed leader. The coexistence of order and disorder originates from stochastic domain-level physics of the insulator–metal and metal–insulator transitions, where local variations in transition temperature drive cycle-to-cycle randomness in nucleation, percolation, and relaxation. A stochastic domain model reproduces this effect by generating synchronized spike trains with random lead–lag jitter, and experimental interspike interval statistics confirm the persistence of fine-scale variability despite macroscopic phase locking. These findings establish that synchronization and stochasticity can coexist within the same physical platform, revealing hidden disorder within collective order. Furthermore, this insight reframes synchronization as not purely deterministic, but as a universal context where microscopic variability can persist, with implications for electronics, cryptography, and the fundamental physics of order–disorder coexistence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗
DOE OSTI2021

Modeling household online shopping demand in the U.S.: a machine learning approach and comparative investigation between 2009 and 2017

Despite the rapid growth of online shopping and research interest in the relationship between online and in-store shopping, national-level modeling and investigation of the demand for online shopping with a prediction focus remain limited in the literature. Here, this paper differs from prior work and leverages two recent releases of the U.S. National Household Travel Survey (NHTS) data for 2009 and 2017 to develop machine learning (ML) models, specifically gradient boosting machine (GBM), for predicting household-level online shopping purchases. The NHTS data allow for not only conducting nationwide investigation but also at the level of households, which is more appropriate than at the individual level given the connected consumption and shopping needs of members in a household. We follow a systematic procedure for model development including employing Recursive Feature Elimination algorithm to select input variables (features) in order to reduce the risk of model overfitting and increase model explainability. Among several ML models, GBM is found to yield the best prediction accuracy. Extensive post-modeling investigation is conducted in a comparative manner between 2009 and 2017, including quantifying the importance of each input variable in predicting online shopping demand, and characterizing value-dependent relationships between demand and the input variables. In doing so, two latest advances in machine learning techniques, namely Shapley value-based feature importance and Accumulated Local Effects plots, are adopted to overcome inherent drawbacks of the popular techniques in current ML modeling. The modeling and investigation are performed at the national level, with a number of findings obtained. The models developed and insights gained can be used for online shopping-related freight demand generation and may also be considered for evaluating the potential impact of relevant policies on online shopping demand.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗
DOE OSTI2024

Robust Combined Heat and Hybrid Power (CHHP) for High Electrical Efficiency Cogeneration

Georgia Tech (Prime Recipient), the University of Texas at El Paso (UTEP) and the National Energy Technology Laboratory (NETL) investigated a hybrid fuel cell/ gas turbine system concept as a combined heat and hybrid power (CHHP) system for both robust and high power-to-process heat ratio cogeneration. The novelty of the proposed system entailed the distinct, elevated electrical efficiencies it maintains while simultaneously supporting a broad span of heating needs (e.g., supply temperatures) demanded across variable heat loads. The scope included: 1) leveraging a pioneering national lab facility configured for dynamic system operability development of hybrid fuel cell/gas turbine cycles; 2) enabling technology development to adjust and modulate the quality and quantity of thermal supply to bottoming heat loads via novel extreme temperature gas bypass valves. Hybrid fuel cell/gas turbine systems have primarily been reduced-to-practice in a constrained (e.g., initial proof-of-concept) manner and have still demonstrated considerable electrical efficiencies. However, these pre-pilot systems have focused upon electrical efficiencies and electrical power generation as the exclusive energy demand. Such hybrid systems had not been extensively researched or developed for flexible and variable operation consisting of both power and heat demands; however, these variable combined power and heat demands are characteristic of many types of manufacturers such as animal/poultry processing, bakeries and milk/flour/pastry manufacturing, textile mills, and electrochemical processing. Commercially, developing the system into a working combined heat and power system benefits these types of manufacturers by allowing them to meet their power and heat demands at a lower cost, higher efficiency, and/or through onsite generation. Therefore, the technical scope of this project was largely to study and facilitate these hybrid systems as combined heat and hybrid power (CHHP) systems that include dynamic operability for variable heat and power loads and/or grid dynamics for various types of manufacturers. Simulation results were used to predict the performance of the CHHP system and conceptually develop it to achieve desired dynamic operability. Experimentally, the primary goal was to design, manufacture, and experiment upon a high-temperature bypass valve. Experimental data included air mass flow rates through the valve orifice when the valve was changed to variable extent between fully closed and fully open. The experimental data was then used to create a semi-empirical computational model of the bypass valve. Concluded simulation goals for the research included developing computational heat exchanger models for the hybrid system inclusive of the bottoming heat exchanger and the recuperative heat exchanger, and then combining the computational recuperator model with the computational valve model. Afterwards, the computational models were then integrated to predict the dynamic operation of hybrid fuel cell/gas turbine cycles throughout a design space and reporting such. The scope stated in the preceding paragraph was packaged into five specific goals: 1) enabling the simulation of dynamic combined heat and power through the creation of computational, modular heat exchanger models; 2) simulation and exploration of the CHHP system’s performance by integrating the heat exchanger models with the national lab’s pre-existing hybrid system (computational) simulation, but without the recuperator bypass valve concept in order to initially determine how the (baseline) system behaves and can be controlled in order to meet variable heat and power demands; 3) development and initial deployment of the high-temperature recuperator bypass valve technology in order to confirm and characterize the approach; 4) usage of the experimental data for the valve to create a semi-empirical computational model for the bypass valve which could then be combined with the heat exchanger computational models; 5) repeat of the second task of simulating and exploring the system’s performance, but this time including the bypass valve to resolve its efficacy. Tasks were successfully completed, and the general notion of flexibly operating, high electrical efficiency CHHP was further corroborated. Supportive details are provided in the report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

A Canonical Transformation to Eliminate Resonant Perturbations. I

We study dynamical systems that admit action-angle variables at leading order, which are subject to nearly resonant perturbations. If the frequencies characterizing the unperturbed system are not in resonance, the long-term dynamical evolution may be integrated by orbit-averaging over the high-frequency angles, thereby evolving the orbit-averaged effect of the perturbations. It is well known that such integrators may be constructed via a canonical transformation, which eliminates the high-frequency variables from the orbit-averaged quantities. An example of this algorithm in celestial mechanics is the von Zeipel transformation. However, if the perturbations are inside or close to a resonance, i.e., the frequencies of the unperturbed system are commensurate; these canonical transformations are subject to divergences. We introduce a canonical transformation that eliminates the high-frequency phase variables in the Hamiltonian without encountering divergences. This leads to a well-behaved symplectic integrator. We demonstrate the algorithm through two examples: a resonantly perturbed harmonic oscillator and the gravitational three-body problem in mean motion resonance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Practical gyrokinetics

We report thousands of gyrokinetic papers have been published since the introduction of the gyrokinetic change of variables. The intent here is not to review the field, but rather the goal is to present a tutorial providing insight into why gyrokinetics is an appropriate description of turbulent transport. The focus is on turbulent transport in axisymmetric tokamaks, but many of the ideas and techniques are applicable to stellarators and other magnetic fusion devices with nested surfaces of constant pressure. Besides the origins of gyrokinetics and recent insights, gyrokinetic orderings and gyrokinetic variables are summarized. Then a compact, but careful, derivation of the simplest electrostatic gyrokinetic equation for tokamaks is presented, along with a brief mention of gyrokinetics for stellarators. The advantages of assuming scale separation between the finer spatial scales and faster time variation of the turbulence and the global behaviour and slow temporal evolution of the background are stressed. Moreover, the procedure for removing the adiabatic or Maxwell–Boltzmann response is emphasized. Scale separation allows the near Maxwellian behaviour of the background and the rapid variation of the turbulence to be described by separate, local gyrokinetic equations. The turbulent fluctuations are found by solving the nonlinear gyrokinetic equation for the non-adiabatic portion of the fluctuating distribution function. Also, generalizations of the gyrokinetic equation so that it retains electromagnetic and flow modifications are considered in detail along with some symmetry properties. In addition, ambipolarity, particle transport and heat transport are addressed, along with a discussion of the complications associated with describing momentum transport and profile evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Practical Gyrokinetics

Thousands of gyrokinetic papers have been published since the introduction of the gyrokinetic change of variables. The intent here is not to review the field, but rather the goal is to present a tutorial providing insight into why gyrokinetics is an appropriate description of turbulent transport. The focus is on turbulent transport in axisymmetric tokamaks, but many of the ideas and techniques are applicable to stellarators and other magnetic fusion devices with nested surfaces of constant pressure. Besides the origins of gyrokinetics and recent insights, gyrokinetic orderings and gyrokinetic variables are summarized. Then a compact, but careful, derivation of the simplest electrostatic gyrokinetic equation for tokamaks is presented, along with a brief mention of gyrokinetics for stellarators. The advantages of assuming scale separation between the finer scales and faster time variation of the turbulence and the global behavior and slow temporal evolution of the background are stressed. Moreover, the procedure for removing the adiabatic or Maxwell-Boltzmann response is emphasized. Scale separation allows the near Maxwellian behavior of the background and the rapid variation of the turbulence to be described by separate local gyrokinetic equations. The turbulent fluctuations are found by solving the nonlinear gyrokinetic equation for the non-adiabatic portion of the fluctuating distribution function. Also, generalizations of the gyrokinetic equation so that it retains electromagnetic and flow modifications are considered in detail along with some symmetry properties. In addition, ambipolarity, particle transport, and heat transport are addressed, along with a discussion of the complications associated with describing momentum transport and profile evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Compact finite difference schemes for the time fractional diffusion equation with nonlocal boundary conditions

Fourth-order compact finite difference schemes for solving the time fractional diffusion equation with nonlocal boundary conditions are considered in this paper. The second-order derivative in one space variable is approximated by a fourth-order compact finite difference. Two different time derivatives, the Caputo derivative and Riemann–Liouville derivative are considered in the equations that we discuss. The Caputo time derivative is discretized by the L1-approximation and a second-order approximation is used for the Riemann–Liouville derivative. Stability analysis using the discrete energy method for Scheme I is given, and error estimation of this scheme is derived by conducting the local truncation error analysis, using the stability result. Numerical results are provided to show that both two fully discrete schemes have fourth-order accuracy in space, and 2−γ order of accuracy in time for Scheme I with γ being the order of the fractional derivative, while Scheme II has second-order accuracy in time when the theoretical solutions are smooth with zero initial values.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Coupled Multimodal Dynamics of Hydrogen-Containing Ion Networks in Water-Deficient, Sodium Hydroxide-Aluminate Solutions

The (meta)stability of low water activity sodium hydroxide/aluminate (Na + OH–/Al(OH) 4 –) electrolytes dictates kinetics in the Bayer process for aluminum refining and high-level nuclear waste processing. We utilized quasi-elastic neutron scattering (QENS) and proton nuclear magnetic resonance spectroscopy ( 1 H NMR) in extremely concentrated sodium aluminate solutions to investigate the picosecond (ps) to microsecond (ms) timescale motions of H-bearing species (Al(OH) 4 – monomers/clusters, OH– and H 2 O). In the QENS data, in contrast to typical liquids, no short-time translational diffusion was observed at 293 K, but two types of localized motions were found: (i) local backbone tumbling or a formation of large hydrated ion clusters on the order of 40–60 ps; and (ii) much slower, complex, and collective dynamics of the ensemble of H-bearing species on the order of 350–750 ps. Variable temperature, pulsed field gradient, diffusion-ordered 1H NMR was used to determine the ensemble translational motion along with relaxometry to calculate rotational correlation coefficients. The ensemble rotational correlation times were on the order of 184–300 ps from 1 H NMR, which is consistent with the timescale of the QENS components. Complementary molecular dynamics simulation of NaOH solutions exhibit extensive ion networks potentially responsible for the observed dynamical coupling of water with the motion of large hydrated ion clusters. Understanding these collective motions will aid in predicting the behavior of complex solutions during aluminum production and during nuclear waste processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Using Cosmic Ray Muons to Assess Geological Characteristics in the Subsurface

Cosmic rays are energetic nuclei and elementary particles that originate from stars and intergalactic events. The interaction of these particles with the upper atmosphere produces a wide range of secondary particles that reach the surface of the earth, of which muons are the most prominent. With enough energy, muons can travel up to a few kilometers beneath the surface of the earth before being stopped completely. The terrestrial muon flux profile and associated zenith angle can be utilized to determine geological characteristics of a location (e.g., rock overburden and density) without having to use conventional methods such as boreholes. This work uses a low-power plastic scintillator-based muon detection system as a prototype for this non-destructive geological assay methodology. Four custom designed 102 cm x 51 cm x 5 cm plastic scintillation panels are used to realize two orthogonal detection planes. Optical photons from each scintillation panel are read using OnSemi J-Series 4x4 silicon photomultiplier (SiPM) arrays in conjunction with preamplifiers. Simultaneous triggers between detectors from two planes indicate a coincidence event which is recorded using the QuarkNet data acquisition system (DAQ) from Fermi National Accelerator Laboratory. A custom detector holder was designed to securely mount the detection system and rotate the panels along the zenith to collect data at variable angles. In order to quantify the systematic uncertainties associated with the detector, such as energy depositions and angular resolution of the detector design, a Monte Carlo (MC) simulation using Geant4 is being developed. Cosmic ray flux prediction will be included in the project by adding the CORSIKA MC code to the simulation toolchain. Simulated and experimental data will drive the development and validation of a reconstruction algorithm that, upon completion, is expected to predict average overburden and rock density. Extended detector exposure to muons can be used as a means to understand changes in the surrounding environment like rock porosity. On the experimental front, muons will initially be measured at the surface, establishing the baseline flux. This is followed by recording the muon flux at variable depths and zenith angles, where the data will be used by the reconstruction algorithm to predict the overburden. The result will be benchmarked against geological surveys. The measured flux data will also be used to benchmark independent and established models. Successful proof-of-concept demonstration of this technology can open doors for long term non-invasive geological monitoring. The detector design, experimental methodology, and the benchmarking efforts are detailed in this work.

Gadey, Harish Reddy↗
DOE OSTI2019

Boundary-Value Problems with Impulsive Conditions for Parabolic Equations with Degenerations

We consider the first boundary-value problem and the unilateral boundary-value problem with impulsive conditions with respect to the time variable for a linear parabolic equation with power singularities of any order in the space variables. With the help of the maximum principle and a priori estimates, we establish the existence and uniqueness of the solutions of the posed problems in the Hölder spaces with power weights.

97 MATHEMATICS AND COMPUTING↗
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↗
DOE OSTI2018

Mass-preserving time second-order explicit–implicit domain decomposition schemes for solving parabolic equations with variable coefficients

In this paper, the new mass-preserving time second-order explicit–implicit domain decomposition (DD) schemes for solving parabolic equations with various coefficients are proposed. In the schemes, first, the interface values on the interfaces of sub-domains are explicitly calculated, in which two kinds of explicit difference schemes are used. Second, the interior values in sub-domains are determined by the implicit schemes. Third, the interface values are recalculated explicitly to preserve mass. The important feature is that the developed schemes are more than first order in time step and the correction operators at the third step ensure the schemes’ mass conservative. The stability and convergence are analyzed for the two schemes. The optimal error estimates are obtained in the H{sup 1} semi-norm under a stability condition. Numerical experiments confirm the theoretical results.

97 MATHEMATICS AND COMPUTING↗