Search NASA⌕ Search

SEARCH · Search NASA

Results for “backward”

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 271 records · Page 15

A Hybrid Lagrangian/Eulerian Collocated Velocity Advection and Projection Method for Fluid Simulation

Abstract We present a hybrid particle/grid approach for simulating incompressible fluids on collocated velocity grids. Our approach supports both particle‐based Lagrangian advection in very detailed regions of the flow and efficient Eulerian grid‐based advection in other regions of the flow. A novel Backward Semi‐Lagrangian method is derived to improve accuracy of grid based advection. Our approach utilizes the implicit formula associated with solutions of the inviscid Burgers’ equation. We solve this equation using Newton's method enabled by C 1 continuous grid interpolation. We enforce incompressibility over collocated, rather than staggered grids. Our projection technique is variational and designed for B‐spline interpolation over regular grids where multiquadratic interpolation is used for velocity and multilinear interpolation for pressure. Despite our use of regular grids, we extend the variational technique to allow for cut‐cell definition of irregular flow domains for both Dirichlet and free surface boundary conditions.

Gagniere, S.↗

Attributes of Bi-Directional Turbomachinery for Pumped Thermal Energy Storage

Abstract In this paper, we (i) present a methodology for determining the aerodynamic performance of bi-directional turbomachines for pumped thermal energy storage, i.e., turbomachines designed to operate as a compressor in one direction, and then as a turbine in the opposite direction, (ii) carry out performance computations for such turbomachines, and (iii) propose principles for conceptual design of these devices. Focus is placed on using the energy storage cycle not only to identify the novel requirements placed on bi-directional turbomachines but also to estimate the effect of these requirements on the efficiency of the energy storage process. In particular, the difference between aerodynamic loading in forward and in backward operation causes the blading to work at incidences leading to the performance below maximum efficiency, resulting in a lower round-trip efficiency. The description of the design principles includes determination of the number of stages, definition of nondimensional parameters for blading selection, and optimization of two-dimensional blading for bi-directional operation. The assessment of stage count shows the relationship between relative Mach number, pressure ratio, and round-trip efficiency. The nondimensional parameters are assessed through a bi-directional analogue to existing “Smith charts,” for the efficiency of single-direction turbomachines, as a function of camber and stagger. The blade shape evaluation and optimization show how the blade profile can be modified to address the requirements of a bi-directional turbomachine, enabling an increase in round-trip efficiency of two percentage points compared to a baseline double circular arc configuration.

Engineering↗

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

Layer-Parallel Training of Deep Residual Neural Networks

Residual neural networks (ResNets) are a promising class of deep neural networks that have shown excellent performance for a number of learning tasks, e.g., image classification and recognition. Mathematically, ResNet architectures can be interpreted as forward Euler discretizations of a nonlinear initial value problem whose time-dependent control variables represent the weights of the neural network. Hence, training a ResNet can be cast as an optimal control problem of the associated dynamical system. For similar time-dependent optimal control problems arising in engineering applications, parallel-in-time methods have shown notable improvements in scalability. This paper demonstrates the use of those techniques for efficient and effective training of ResNets. The proposed algorithms replace the classical (sequential) forward and backward propagation through the network layers with a parallel nonlinear multigrid iteration applied to the layer domain. This adds a new dimension of parallelism across layers that is attractive when training very deep networks. From this basic idea, we derive multiple layer-parallel methods. The most efficient version employs a simultaneous optimization approach where updates to the network parameters are based on inexact gradient information in order to speed up the training process. Finally, using numerical examples from supervised classification, we demonstrate that the new approach achieves a training performance similar to that of traditional methods, but enables layer-parallelism and thus provides speedup over layer-serial methods through greater concurrency.

97 MATHEMATICS AND COMPUTING↗

Stability Analysis of Inline ZFP Compression for Floating-Point Data in Iterative Methods

Currently, the dominating constraint in many high performance computing applications is data capacity and bandwidth, in both internode communications and even moreso in intranode data motion. A new approach to address this limitation is to make use of data compression in the form of a compressed data array. Storing data in a compressed data array and converting to standard IEEE-754 types as needed during a computation can reduce the pressure on bandwidth and storage. However, repeated conversions (lossy compression and decompression) introduce additional approximation errors, which need to be shown to not significantly affect the simulation results. Here, we extend recent work that analyzed the error of a single use of compression and decompression of the ZFP compressed data array representation to the case of time-stepping and iterative schemes, where an advancement operator is repeatedly applied in addition to the conversions. We show that the accumulated error for iterative methods involving fixed-point and time evolving iterations is bounded under standard constraints. An upper bound is established on the number of additional iterations required for the convergence of stationary fixed-point iterations. An additional analysis of traditional forward and backward error of stationary iterative methods using ZFP compressed arrays is also presented. The results of several 1D, 2D, and 3D test problems are provided to demonstrate the correctness of the theoretical bounds.

97 MATHEMATICS AND COMPUTING↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

An Algebraic Quantum Circuit Compression Algorithm for Hamiltonian Simulation

Quantum computing is a promising technology that harnesses the peculiarities of quantum mechanics to deliver computational speedups for some problems that are intractable to solve on a classical computer. Current generation noisy intermediate-scale quantum (NISQ) computers are severely limited in terms of chip size and error rates. Shallow quantum circuits with uncomplicated topologies are essential for successful applications in the NISQ era. In this work, based on matrix analysis, we derive localized circuit transformations to efficiently compress quantum circuits for simulation of certain spin Hamiltonians known as free fermions. The depth of the compressed circuits is independent of simulation time and grows linearly with the number of spins. The proposed numerical circuit compression algorithm behaves backward stable and scales cubically in the number of spins enabling circuit synthesis beyond O(10 3 ) spins. The resulting quantum circuits have a simple nearest-neighbor topology, which makes them ideally suited for NISQ devices.

Hamiltonian simulation↗

Two-Level Sketching Alternating Anderson Acceleration for Complex Physics Applications

We present a novel two-level sketching extension of the Alternating Anderson–Picard (AAP) method for accelerating fixed-point iterations in challenging single- and multiphysics simulations governed by discretized PDEs. Our approach combines a static, physics-based projection that reduces the least-squares (LS) problem to the most informative field (e.g., via Schur-complement insight) with a dynamic, algebraic sketching stage driven by a backward stability analysis under Lipschitz continuity. We introduce inexpensive estimators for stability thresholds and cache-aware randomized selection strategies to balance computational cost against memory access overhead. The resulting algorithm solves reduced LS systems in place, minimizes memory footprints, and seamlessly alternates between low-cost Picard updates and Anderson mixing. Implemented in Julia, our two-level sketching AAP achieves up to 50% time-to-solution reductions compared to standard Anderson acceleration—without degrading convergence rates—on benchmark problems including Stokes, 𝑝-Laplacian, bidomain, and Navier–Stokes formulations at varying problem sizes. These results demonstrate the method’s robustness, scalability, and potential for integration into high-performance scientific computing frameworks. Our implementation is available open source in the AAP.jl library.

Barnafi, Nicolas [University of Chile, Santiago]↗

Transverse energy–energy correlations of jets in the electron–proton deep inelastic scattering at HERA

We study the event shape variables, transverse energy–energy correlation TEEC (cosφ) and its asymmetry ATEEC (cosφ) in deep inelastic scattering (DIS) at the electron–proton collider HERA, where φ is the angle between two jets defined using a transverse-momentum (k T ) jet algorithm. At HERA, jets are defined in the Breit frame, and the leading nontrivial transverse energy–energy correlations arise from the 3-jet configurations. With the help of the NLOJET++, these functions are calculated in the leading order (LO) and the next-to-leading order (NLO) approximations in QCD at the electron–proton center-of-mass energy $\sqrt{s}$ = 314 GeV. We restrict the angular region to -0.8 ≤ cosφ ≤ 0.8, as the forward- and backward-angular regions require resummed logarithmic corrections, which we have neglected in this work. Following experimental jet-analysis at HERA, we restrict the DIS-variables x, y = Q 2 /(xs), where Q 2 = -q 2 is the negative of the momentum transfer squared q 2 , to 0 ≤ x ≤ 1, 0.2 ≤ y ≤ 0.6, and the pseudo-rapidity variable in the laboratory frame (η lab ) to the range -1 ≤ η lab ≤ 2.5. The TEEC and ATEEC functions are worked out for two ranges in Q 2 , defined by 5.5GeV 2 ≤ Q 2 ≤ 80GeV 2 , called the low-Q 2 -range, and 150GeV 2 ≤ Q 2 ≤ 1000GeV 2 , called the high-Q 2 -range. We show the sensitivity of these functions on the parton distribution functions (PDFs), the factorization (μ F ) and renormalization (μ R ) scales, and on α s (M$^2_Z$). Of these the correlations are stable against varying the scale μ F and the PDFs, but they do depend on μ R . For the choice of the scale μ R =$\sqrt{\langle ET \rangle^2+Q^2}$, advocated in earlier jet analysis at HERA, the shape variables TEEC and ATEEC are found perturbatively robust. These studies are useful in the analysis of the HERA data, including the determination of α s (M$^2_Z$) from the shape variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radiative three-body D -meson decays in and beyond the standard model

We study radiative charm decays D → P 1 P 2γ , P 1,2 = π, K in QCD factorization at leading order and within heavy hadron chiral perturbation theory. Branching ratios including resonance contributions are around ~10 –3 for the Cabibbo-favored modes into Kπγ and ~10 –5 for the singly Cabibbo-suppressed modes into π + π – γ,K + K – γ, and thus in reach of the flavor factories BES III and Belle II. Dalitz plots and forward–backward asymmetries reveal significant differences between the two QCD frameworks; such observables are therefore ideally suited for a data-driven identification of relevant decay mechanisms in the standard-model dominated D → Kπγ decays. This increases the potential to probe new physics with the D → π + π – γ and D → K + K – γ decays, which are sensitive to enhanced dipole operators. CP asymmetries are useful to test the SM and look for new physics in neutral |ΔC| = 1 transitions. Cuts in the Dalitz plot enhance the sensitivity to new physics due to the presence of both s- and t, u-channel intermediate resonances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An experiment for electron-hadron scattering at the LHC

Abstract Novel considerations are presented on the physics, apparatus and accelerator designs for a future, luminous, energy frontier electron-hadron ( eh ) scattering experiment at the LHC in the thirties for which key physics topics and their relation to the hadron-hadron HL-LHC physics programme are discussed. Demands are derived set by these physics topics on the design of the LHeC detector, a corresponding update of which is described. Optimisations on the accelerator design, especially the interaction region (IR), are presented. Initial accelerator considerations indicate that a common IR is possible to be built which alternately could serve eh and hh collisions while other experiments would stay on hh in either condition. A forward-backward symmetrised option of the LHeC detector is sketched which would permit extending the LHeC physics programme to also include aspects of hadron-hadron physics. The vision of a joint eh and hh physics experiment is shown to open new prospects for solving fundamental problems of high energy heavy-ion physics including the partonic structure of nuclei and the emergence of hydrodynamics in quantum field theory while the genuine TeV scale DIS physics is of unprecedented rank.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hint of a new scalar interaction in LHCb data?

Abstract We explain recent LHCb measurements of the lepton universality ratios, $$R_{D^{(*)}}^{\tau /\ell }\equiv \frac{\mathcal {B}(\bar{B} \rightarrow D^{(*)+} \tau ^- \bar{\nu }_\tau )}{\mathcal {B}(\bar{B} \rightarrow D^{(*)+}\ell ^- \bar{\nu }_\ell )}$$ R D ( ∗ ) τ / ℓ ≡ B ( B ¯ → D ( ∗ ) + τ - ν ¯ τ ) B ( B ¯ → D ( ∗ ) + ℓ - ν ¯ ℓ ) and $${R(\Lambda _c^+)}^{\tau /\ell } \equiv \frac{\mathcal {B}(\Lambda _b \rightarrow \Lambda _c^+ \tau ^- \bar{\nu }_{\tau })}{\mathcal {B}(\Lambda _b \rightarrow \Lambda _c^+ \ell ^- \bar{\nu }_{\ell })}$$ R ( Λ c + ) τ / ℓ ≡ B ( Λ b → Λ c + τ - ν ¯ τ ) B ( Λ b → Λ c + ℓ - ν ¯ ℓ ) with $$\ell =\mu $$ ℓ = μ , via new physics that affects $$R_D^{\tau /\ell }$$ R D τ / ℓ and $$R(\Lambda _c^+)^{\tau /\ell }$$ R ( Λ c + ) τ / ℓ but not $$R_{D^*}^{\tau /\ell }$$ R D ∗ τ / ℓ . The scalar operator in the effective theory for new physics is indicated. We find that the forward-backward asymmetry and $$\tau $$ τ polarization in $$\bar{B} \rightarrow D^+ \tau ^{-} \bar{\nu }_{\tau }$$ B ¯ → D + τ - ν ¯ τ and $$\Lambda _b \rightarrow \Lambda _c^+ \tau ^- \bar{\nu }_{\tau }$$ Λ b → Λ c + τ - ν ¯ τ decays are significantly affected by the scalar interaction. We construct a simple two Higgs doublet model as a realization of our scenario and consider lepton universality in semileptonic charm and top decays, radiative B decay, B -mixing, and $$Z \rightarrow b \bar{b}$$ Z → b b ¯ .

Physics↗

Measurement of the production cross section of prompt $\Xi ^0_{\textrm{c}}$ baryons in p–Pb collisions at $\sqrt{s_{{\textrm{NN}}}}=5.02$ TeV

The transverse momentum (p T ) differential production cross section of the promptly produced charm-strange baryon $\Xi ^0_{\textrm{c}}$ (and its charge conjugate $\Xi ^0_{\textrm{c}}$) is measured at midrapidity via its hadronic decay into π + $\Xi$ – in p–Pb collisions at a centre-of-mass energy per nucleon–nucleon collision $\sqrt{s_{{\textrm{NN}}}}=5.02$ TeV with the ALICE detector at the LHC. The $\Xi ^0_{\textrm{c}}$ nuclear modification factor (R pPb ), calculated from the cross sections in pp and p–Pb collisions, is presented and compared with the R pPb of Λ$^{+}_{c}$ baryons. The ratios between the p T -differential production cross section of $\Xi ^0_{\textrm{c}}$ baryons and those of D 0 mesons and Λ$^{+}_{c}$ baryons are also reported and compared with results at forward and backward rapidity from the LHCb Collaboration. The measurements of the production cross section of prompt $\Xi ^0_{\textrm{c}}$ baryons are compared with a model based on perturbative QCD calculations of charm-quark production cross sections, which includes only cold nuclear matter effects in p–Pb collisions, and underestimates the measurement by a factor of about 50. This discrepancy is reduced when the data is compared with a model that includes string formation beyond leading-colour approximation or in which hadronisation is implemented via quark coalescence. The p T -integrated cross section of prompt $\Xi ^0_{\textrm{c}}$-baryon production at midrapidity extrapolated down to p T = 0 is also reported. These measurements offer insights and constraints for theoretical calculations of the hadronisation process. Additionally, they provide inputs for the calculation of the charm production cross section in p–Pb collisions at midrapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring Hadron Structure with GPDs at EIC: New Topics in Theory, Experiment, Interpretation

Exploring hadron structure with generalized parton distributions (GPDs) and related concepts is an essential part of the physics program of the planned Electron-Ion Collider (EIC). Here, we discuss new topics in GPD physics (theory, experiment, interpretation) that have emerged in the last years and could be included in the EIC program. This includes: (a) D-term and QCD forces in the nucleon; (b) Transversity GPDs in pseudoscalar meson production and high-mass pair production; (c) Pion and kaon form factors and GPDs; (d) Transition distribution amplitudes and backward meson production; (e) Neutron GPDs and nuclear modifications; (f) Nuclear GPDs and quark/gluon imaging of nuclei; (g) Twist-3 GPDs and quark-gluon correlations; (h) Wigner functions and potential observables. We also comment on future theoretical developments needed for the GPD program at EIC and emerging connections with other areas of nuclear physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

TorchBraid: High-Performance Layer-Parallel Training of Deep Neural Networks with MPI and GPU Acceleration

TorchBraid is a high-performance implementation of layer-parallel training for deep neural networks (DNNs) supporting MPI-based parallelism and GPU acceleration. Layer-parallel training has been developed to overcome the serialization inherent in forward and backward propagation of DNNs that limits utilization of computational resources in the strong scaling limit. To achieve this, TorchBraid integrates the PyTorch neural network framework with the state-of-the-art XBraid time-parallel library. Furthermore, this article presents the use and performance of TorchBraid, in addition to solutions for overcoming the algorithmic challenges inherent in combining automatic differentiation with layer-parallel. Results are presented with and without GPU acceleration for the Tiny ImageNet and MNIST image classification data sets, as well as recurrent neural networks. Overall, TorchBraid enables fast training of DNNs, both in a strong and weak scaling context. In addition to the TorchBraid software, several new advances in applying layer-parallel algorithms are detailed. Integration of layer-parallel with data-parallel algorithms is presented for the first time, showing the computational advantages of the combination. Standard deep learning techniques, like batch-normalization, are developed for layer-parallel training. Finally, a new approach combining layer-parallel with spatial coarsening in order to accelerate training for 3D image classification shows roughly a 10× speedup over serial execution.

Layer-parallel↗

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗

HQ-Tracer

HQ-Tracer is a software plugin related to the open-source project Ghidra. It is used for taking an emulation instruction trace from open source platforms and displaying inside of Ghidra. It allows for software debugging, reverse engineering, and enables interactive forward and backward stepping and searching through the execution traces. 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. SAND2021-4552 O

Wright, ChristopherM.↗