SEARCH · Search NASA

Results for “operator exponential”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Observations of the magneto-Rayleigh-Taylor instability and shock dynamics in gas-puff Z-pinch experiments

Here, we describe a series of experiments performed to study the shock structure generated during the implosion of a gas-puff Z-pinch. The Z-pinch is produced by a double-annular gas-puff with a center jet driven by Cornell University's COBRA generator operating with a 1 MA, 200 ns current pulse. Using 532 nm laser interferometry and 100 MHz multi-frame cameras, a shock structure is observed to form early in the implosion. The shock appears to be created by a current layer at the outer radius of the imploding plasma which acts as a piston moving inward at several hundred km s –1 . The dynamics of the shock and its radial position ahead of the piston agree well with a simple uniform density model outlined in the study by Potter [Nucl. Fusion 18(6), 813 (1978)]. The outer surface of the current layer is observed to be Magneto-Rayleigh-Taylor unstable. The growth rate of this instability is found to depend on the radial density profile of the material within the layer of high-density fluid between the shock and the piston, as predicted by recent theoretical work [see, e.g., D. Livescu, Phys. Fluids 16(1), 118 (2004)]. Growth rates measured in krypton implosions, where the post-shock material is found to decay quasi-exponentially away from the piston, were more than ten times smaller than those recorded in otherwise identical implosions in argon plasmas, where the material between the shock and the piston was observed to maintain a uniform density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2000

Intelligence Applied to Air Vehicles

The exponential growth in information technology has provided the potential for air vehicle capabilities that were previously unavailable to mission and vehicle designers. The increasing capabilities of computer hardware and software, including new developments such as neural networks, provide a new balance of work between humans and machines. This paper will describe several NASA projects, and review results and conclusions from ground and flight investigations where vehicle intelligence was developed and applied to aeronautical and space systems. In the first example, flight results from a neural network flight control demonstration will be reviewed. Using, a highly-modified F-15 aircraft, a NASA/Dryden experimental flight test program has demonstrated how the neural network software can correctly identify and respond to changes in aircraft stability and control characteristics. Using its on-line learning capability, the neural net software would identify that something in the vehicle has changed, then reconfigure the flight control computer system to adapt to those changes. The results of the Remote Agent software project will be presented. This capability will reduce the cost of future spacecraft operations as computers become "thinking" partners along with humans. In addition, the paper will describe the objectives and plans for the autonomous airplane program and the autonomous rotorcraft project. Technologies will also be developed.

Rosen, Robert↗

DOE OSTI2016

Exascale Virtualized and Programmable Distributed Cyber Resource Control: Final Scientific Technical Report

Extreme-science drives the need for distributed exascale processing and communications that are carefully, yet flexibly, managed. Exponential growth of data for scientific simulations, experimental data, collaborative data analyses, remote visualization and GRID computing requirements of scientists in fields as diverse as high energy physics, climate change, genomics, fusion, synchrotron radiation, material science, medicine, and other scientific disciplines cannot be accommodated by simply applying existing transport protocols to faster pipes. Further, scientific challenges today demand diverse research teams, heightening the need for and increasing the complexity of collaboration. To address these issues within the network layer and physical layer, we have performed a number of research activities surrounding effective allocation and management of elastic optical network (EON) resources, particularly focusing on FlexGrid transponders. FlexGrid transponders support the opportunity to build Layer-1 connections at a wide range of bandwidths and to reconfigure them rapidly. The new flexibility supports complex new ways of using the physical layer that must be carefully managed and hidden from the scientist end-users. FlexGrid networks utilize flexible (or elastic) spectral bandwidths for each data link without using fixed wavelength grids. The flexibility in spectrum allocation brings many appealing features to network operations. Current networks are designed for the worst case impairments in transmission performance and the assigned spectrum is over-provisioned. In contrast, the FlexGrid networks can operate with the highest spectral efficiency and minimum bandwidth for the given traffic demand while meeting the minimum quality of transmission (QoT) requirement. Two primary focuses of our research are: (1) resource and spectrum allocation (RSA) for IP traffic over EONs, and (2) RSA for cross-domain optical networks. Previous work concentrates primarily on large file transfers within a single domain. Adding support for IP traffic changes the nature of the RSA problem: instead of choosing to accept or deny each request for network support, IP traffic is inherently elastic and thus lends itself to a bandwidth maximization formulation. We developed a number of algorithms that could be easily deployed within existing and new FlexGrid networks, leading to networks that better support scientific collaboration. Cross-domain RSA research is essential to support large-scale FlexGrid networks, since configuration information is generally not shared or coordinated across domains. The results presented here are in their early stages. They are technically feasible and practical, but still require coordination among organizations and equipment owners and a higher-layer framework for managing network requests.

42 ENGINEERING↗

NASA NTRS2004

Thermal and Nonthermal Contributions to the Solar Flare X-Ray Flux

The relative thermal and nonthermal contributions to the total energy budget of a solar flare are being determined through analysis of RHESSI X-ray imaging and spectral observations in the energy range from approx. 5 to approx. 50 keV. The classic ways of differentiating between the thermal and nonthermal components - exponential vs. sources - can now be combined for individual flares. In addition, RHESSI's sensitivity down to approx. 4 keV and energy resolution of approx. 1 keV FWHM allow the intensities and equivalent widths of the complex of highly ionized iron lines at approx. 6.7 keV and the complex of highly ionized iron and nickel lines at approx. 8 keV to be measured as a function of time. Using the spectral line and continuum intensities from the Chianti (version 4.2) atomic code, the thermal component of the total flare emission can be more reliably separated from the nonthermal component in the measured X-ray spectrum. The abundance of iron can also be determined from RHESSI line-to-continuum measurements as a function of time during larger flares. Results will be shown of the intensity and equivalent widths of these line complexes for several flares and the temperatures, emission measures, and iron abundances derived from them. Comparisons will be made with 6.7-keV Fe-line fluxes measured with the RESIK bent crystal spectrometer on the Coronas-F spacecraft operating in third order during the peak times of three flares (2002 May 31 at 00:12 UT, 2002 December 2 at 19:26 UT, and 2003 April 26 at 03:OO UT). During the rise and decay of these flares, RESIK was operating in first order allowing the continuum flux to be measured between 2.9 and 3.7 keV for comparison with RHESSI fluxes at its low-energy end.

Dennis, Brian R.↗

DOE OSTI2016

The offset and host light distributions of long gamma-ray bursts: A new view from HST observations of Swift bursts

We present the results of an extensive Hubble Space Telescope imaging study of 105, mostly Swift, long-duration gamma-ray bursts (LGRBs) spanning 0.03≲z≲9.4, which were localized using relative astrometry from ground- and space-based afterglow observations. We measure the distribution of LGRB offsets from their host centers and their relation to the underlying host light distribution. We find that the host-normalized offsets of LGRBs are more centrally concentrated than expected for an exponential disk profile, 〈R/R{sub h}〉 = 0.63, and in particular they are more concentrated than the underlying surface brightness profiles of their host galaxies and more concentrated than supernovae. The fractional flux distribution, with a median of 0.78, indicates that LGRBs prefer some of the brightest locations in their host galaxies but are not as strongly correlated as previous studies indicated. Importantly, we find a clear correlation between offset and fractional flux, where bursts at offsets R/R{sub h}≲0.5 exclusively occur at fractional fluxes ≳0.6, while bursts at R/R{sub h}≳0.5 have a uniform fractional flux distribution. This indicates that the spatial correlation of LGRBs with bright star-forming regions seen in the full sample is dominated by the contribution from bursts at small offset and that LGRBs in the outer parts of galaxies show no preference for unusually bright regions. We conclude that LGRBs strongly prefer the bright, inner regions of their hosts, indicating that the star formation taking place there is more favorable for LGRB progenitor production. This indicates that environmental factors beyond metallicity, such as binary interactions or IMF differences, may operate in the central regions of LGRB hosts.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Algorithms for Performance, Dependability, and Performability Evaluation using Stochastic Activity Networks

Modeling tools and technologies are important for aerospace development. At the University of Illinois, we have worked on advancing the state of the art in modeling by Markov reward models in two important areas: reducing the memory necessary to numerically solve systems represented as stochastic activity networks and other stochastic Petri net extensions while still obtaining solutions in a reasonable amount of time, and finding numerically stable and memory-efficient methods to solve for the reward accumulated during a finite mission time. A long standing problem when modeling with high level formalisms such as stochastic activity networks is the so-called state space explosion, where the number of states increases exponentially with size of the high level model. Thus, the corresponding Markov model becomes prohibitively large and solution is constrained by the the size of primary memory. To reduce the memory necessary to numerically solve complex systems, we propose new methods that can tolerate such large state spaces that do not require any special structure in the model (as many other techniques do). First, we develop methods that generate row and columns of the state transition-rate-matrix on-the-fly, eliminating the need to explicitly store the matrix at all. Next, we introduce a new iterative solution method, called modified adaptive Gauss-Seidel, that exhibits locality in its use of data from the state transition-rate-matrix, permitting us to cache portions of the matrix and hence reduce the solution time. Finally, we develop a new memory and computationally efficient technique for Gauss-Seidel based solvers that avoids the need for generating rows of A in order to solve Ax = b. This is a significant performance improvement for on-the-fly methods as well as other recent solution techniques based on Kronecker operators. Taken together, these new results show that one can solve very large models without any special structure.

Deavours, Daniel D.↗

DOE OSTI2022

Wettability and CHF limits of Accident-Tolerant nuclear fuel cladding materials in light water reactor conditions

Here, we present the results of experimental investigations aimed at evaluating the thermal–hydraulic performance of chromium-coated zircaloy, i.e., one of the most promising accident tolerant fuel (ATF) cladding material for light water nuclear reactors. Precisely, we investigate the wettability and critical heat flux (CHF) limits of chromium-coated and conventional zircaloy surfaces in prototypical reactor conditions. For both surface types, we measure the contact angle in a vapor-saturated atmosphere from ambient pressure to the operating pressure of pressurized water reactors (PWRs), i.e., ~15 MPa. We measure the ambient-pressure steady-state flow boiling CHF with a spatially uniform heat flux. We measure the high-pressure steady-state flow boiling CHF with a cosine shape heat flux (up to 20 MPa) and with a uniform heat flux (up to 15 MPa), also exploring the effect of CRUD deposits on the chromium-coated surface. Our results reveal that the chromium surface and the bare zircaloy surface have similar wettability and both become super-hydrophilic in PWR conditions, and that there is practically no difference in the steady-state CHF limits, both at low-pressure and high-pressure conditions, also when the chromium-coated surface is covered by a CRUD deposit. However, while the chromium-coating does not improve the CHF compared to the bare zircaloy surface, it improves the post-CHF behavior. The chromium coating prevents the reaction between zircaloy and steam, which results in the formation of a brittle zirconium oxide through the surface of the cladding. We also measure the transient CHF under exponentially escalating heat flux inputs of a nano-smooth and a rough surface mimicking a chromium-coated zircaloy cladding. Interestingly, the results of the transient heat flux tests suggest that the CHF limit for very short periods (i.e., fast transients) is independent of the surface finish, being the same for a rough chromium surface or a nano-smooth surface.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Theoretical Models of Optical Transients. I. A Broad Exploration of the Duration–Luminosity Phase Space

The duration–luminosity phase space (DLPS) of optical transients is used, mostly heuristically, to compare various classes of transient events, to explore the origin of new transients, and to influence optical survey observing strategies. For example, several observational searches have been guided by intriguing voids and gaps in this phase space. However, we should ask, do we expect to find transients in these voids given our understanding of the various heating sources operating in astrophysical transients? In this work, we explore a broad range of theoretical models and empirical relations to generate optical light curves and to populate the DLPS. We explore transients powered by adiabatic expansion, radioactive decay, magnetar spin-down, and circumstellar interaction. For each heating source, we provide a concise summary of the basic physical processes, a physically motivated choice of model parameter ranges, an overall summary of the resulting light curves and their occupied range in the DLPS, and how the various model input parameters affect the light curves. We specifically explore the key voids discussed in the literature: the intermediate-luminosity gap between classical novae and supernovae, and short-duration transients (≲10 days). We find that few physical models lead to transients that occupy these voids. Moreover, we find that only relativistic expansion can produce fast and luminous transients, while for all other heating sources events with durations ≲10 days are dim (M{sub R}≳−15 mag). Finally, we explore the detection potential of optical surveys (e.g., Large Synoptic Survey Telescope) in the DLPS and quantify the notion that short-duration and dim transients are exponentially more difficult to discover in untargeted surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2015

Locating the Most Energetic Electrons in Cassiopeia A

We present deep (>2.4 Ms) observations of the Cassiopeia A supernova remnant with NuSTAR, which operates in the 3-79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the >15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power laws in the 15-50 keV bandpass, though the central knots have a steeper (Gamma approx. −3.35) spectrum than the outer filaments (Gamma approx. −3.06). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of >15 keV emission does not follow that of the radio emission nor that of the low energy (<12 keV) X-rays, leaving the origin of the >15 keV emission an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the NuSTAR band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work.

acceleration of particles aEuro" ISM: individual o↗

DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗

NASA NTRS2009

Exploring Alternative Parameterizations for Snowfall with Validation from Satellite and Terrestrial Radars

Increases in computational resources have allowed operational forecast centers to pursue experimental, high resolution simulations that resolve the microphysical characteristics of clouds and precipitation. These experiments are motivated by a desire to improve the representation of weather and climate, but will also benefit current and future satellite campaigns, which often use forecast model output to guide the retrieval process. The combination of reliable cloud microphysics and radar reflectivity may constrain radiative transfer models used in satellite simulators during future missions, including EarthCARE and the NASA Global Precipitation Measurement. Aircraft, surface and radar data from the Canadian CloudSat/CALIPSO Validation Project are used to check the validity of size distribution and density characteristics for snowfall simulated by the NASA Goddard six-class, single moment bulk water microphysics scheme, currently available within the Weather Research and Forecast (WRF) Model. Widespread snowfall developed across the region on January 22, 2007, forced by the passing of a mid latitude cyclone, and was observed by the dual-polarimetric, C-band radar King City, Ontario, as well as the NASA 94 GHz CloudSat Cloud Profiling Radar. Combined, these data sets provide key metrics for validating model output: estimates of size distribution parameters fit to the inverse-exponential equations prescribed within the model, bulk density and crystal habit characteristics sampled by the aircraft, and representation of size characteristics as inferred by the radar reflectivity at C- and W-band. Specified constants for distribution intercept and density differ significantly from observations throughout much of the cloud depth. Alternate parameterizations are explored, using column-integrated values of vapor excess to avoid problems encountered with temperature-based parameterizations in an environment where inversions and isothermal layers are present. Simulation of CloudSat reflectivity is performed by adopting the discrete-dipole parameterizations and databases provided in literature, and demonstrate an improved capability in simulating radar reflectivity at W-band versus Mie scattering assumptions.

Molthan, Andrew L.↗

NASA NTRS2007

Loss-of-Control-Inhibitor Systems for Aircraft

Systems to provide improved tactile feedback to aircraft pilots are being developed to help the pilots maintain harmony between their control actions and the positions of aircraft control surfaces, thereby helping to prevent loss of control. A system of this type, denoted a loss-of-control-inhibitor system (LOCIS) can be implemented as a relatively simple addition to almost any pre-existing flight-control system. The LOCIS concept offers at least a partial solution to the problem of (1) keeping a pilot aware of the state of the control system and the aircraft and (2) maintaining sufficient control under conditions that, as described below, have been known to lead to loss of control. Current commercial aircraft exhibit uneven responses of primary flight-control surfaces to aggressive pilot control commands, leading to deterioration of pilots ability to control their aircraft. In severe cases, this phenomenon can result in loss of control and consequent loss of aircraft. For an older aircraft equipped with a purely mechanical control system, the loss of harmony between a pilot s command action and the control- surface response can be attributed to compliance in the control system (caused, for example, by stretching of control cables, flexing of push rods, or servo-valve distortion). In a newer aircraft equipped with a fly-by-wire control system, the major contributions to loss of harmony between the pilot and the control surfaces are delays attributable to computer cycle time, control shaping, filtering, aliasing, servo-valve distortion, and actuator rate limiting. In addition, a fly-by-wire control system provides no tactile feedback that would enable the pilot to sense such features of the control state as surface flutter, surface jam, position limiting, actuator rate limiting, and control limiting imposed by the aircraft operational envelope. Hence, for example, when a pilot is involved in aggressive closed-loop maneuvering, as when encountering a wake-vortex upset on final landing approach, the control-surface delay can lead to loss of control. Aggressive piloting can be triggered and exacerbated by control-system anomalies, which the pilot cannot diagnose because of the lack of symptoms caused by the absence of feedback through the controls. The purpose served by a LOCIS is to counteract these adverse effects by providing real-time feedback that notifies the pilot that the aircraft is tending to lag the pilot s commands. A LOCIS (see figure) includes cockpit control input-position sensors, control-surface output-position sensors, variable dampers (for example, shock absorbers containing magneto-rheological fluids such that the damping forces can be varied within times of the order of milliseconds by varying applied magnetic fields) attached to the cockpit control levers, electromagnet coils to apply the magnetic fields, and feedback control circuits to drive the electromagnet coils. The feedback control gains are chosen so that the current applied to each electromagnet coil results in a damping force that increases in a suitable nonlinear manner (e.g., exponentially) with the difference between the actual and commanded positions of the affected control surface. The increasing damping force both alerts the pilot to the onset of a potentially dangerous situation and resists the pilot s effort to command a control surface to change position at an excessive rate

AHarrah, Ralph C.↗

DOE OSTI2022

Classically verifiable quantum advantage from a computational Bell test

Existing experimental demonstrations of quantum computational advantage have had the limitation that verifying the correctness of the quantum device requires exponentially costly classical computations. Here we propose and analyse an interactive protocol for demonstrating quantum computational advantage, which is efficiently classically verifiable. Our protocol relies on a class of cryptographic tools called trapdoor claw-free functions. Although this type of function has been applied to quantum advantage protocols before, our protocol employs a surprising connection to Bell’s inequality to avoid the need for a demanding cryptographic property called the adaptive hardcore bit, while maintaining essentially no increase in the quantum circuit complexity and no extra assumptions. Leveraging the relaxed cryptographic requirements of the protocol, we present two trapdoor claw-free function constructions, based on Rabin’s function and the Diffie–Hellman problem, which have not been used in this context before. We also present two independent innovations that improve the efficiency of our implementation and can be applied to other quantum cryptographic protocols. First, we give a scheme to discard so-called garbage bits, removing the need for reversibility in the quantum circuits. Second, we show a natural way of performing postselection that reduces the fidelity needed to demonstrate quantum advantage. Combining these results, we describe a blueprint for implementing our protocol on Rydberg atom-based quantum devices, using hardware-native operations that have already been demonstrated experimentally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Simulating scalar field theories on quantum computers with limited resources

Here we present a quantum algorithm for implementing φ 4 lattice scalar field theory on qubit computers. The field is represented in the discretized field amplitude basis. The number of qubits and elementary gates required by the implementation of the evolution operator is proportional to the lattice size. The algorithm allows efficient φ 4 state preparation for a large range of input parameters in both the normal and broken-symmetry phases. The states are prepared using a combination of variational and adiabatic evolution methods. First, the ground state of a local Hamiltonian, which includes the φ 4 self-interaction, is prepared using short variational circuits. Next, this state is evolved by switching on the coupling between the lattice sites adiabatically. The parameters defining the local Hamiltonian are adjustable and constitute the input of our algorithm. We present a method to optimize these parameters in order to reduce the adiabatic time required for state preparation. For preparing broken-symmetry states, the adiabatic evolution problems caused by crossing the phase transition critical line and by the degeneracy of the broken-symmetry ground state can be addressed using an auxiliary external field which gradually turns off during the adiabatic process. We show that the time dependence of the external field during the adiabatic evolution is important for addressing the broken-symmetry ground state degeneracy. The adiabatic time dependence on the inverse error tolerance can be reduced from quadratic to linear by using a field strength that decreases exponentially in time relative to one that decreases linearly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Inferring the scrape-off layer heat flux width in a divertor with a low degree of axisymmetry

Plasma facing components (PFCs) in the next generation of tokamak devices will operate in challenging environments, with heat loads predicted to exceed 10 MWm -2 . The magnitude of these heat loads is set by the width of the channel, the ‘scrape-off layer’ (SOL), into which heat is exhausted, and can be characterised by an e-folding length scale for the decay of heat flux across the channel. It is expected this channel will narrow as tokamaks move towards reactor relevant conditions. Understanding the processes involved in setting the SOL heat flux width is imperative to be able to predict the heat loads PFCs must handle in future devices. Measurements of the SOL width are performed on the high-field spherical tokamak, ST40, using a newly commissioned infrared thermography system. With its high on-axis toroidal magnetic field (≥1.5 T) ST40 is uniquely positioned to investigate the influence of toroidal field on the heat flux width in spherical tokamaks, whilst also extending measurements of the SOL width in spherical tokamaks to increased poloidal field (≥0.3 T). Due to the divertor on ST40 having a low degree of axisymmetry, it is necessary for a set of radial measurements of the heat flux to be taken across the divertor, made possible using an automated toolchain that fully incorporates its 3D geometry. These radial profiles are combined with the magnetic geometry of the plasma to infer the width of the SOL, with both Eich and double exponential profiles of heat flux observed. A reduction in the heat flux is observed toroidally across part of the divertor, along with increased heat loads observed locally around the edges of the tiles. Future work in characterising the impact of tile misalignment and uncertainties in the reconstructed divertor magnetic geometry is required in order to further understand the observed heat flux patterns, as are additional investigations into the role potentially being played by an inhomogeneous sheath electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Simulating Effective QED on Quantum Computers

In recent years simulations of chemistry and condensed materials has emerged as one of the preeminent applications of quantum computing, offering an exponential speedup for the solution of the electronic structure for certain strongly correlated electronic systems. To date, most treatments have ignored the question of whether relativistic effects, which are described most generally by quantum electrodynamics (QED), can also be simulated on a quantum computer in polynomial time. Here we show that effective QED, which is equivalent to QED to second order in perturbation theory, can be simulated in polynomial time under reasonable assumptions while properly treating all four components of the wavefunction of the fermionic field. In particular, we provide a detailed analysis of such simulations in position and momentum basis using Trotter-Suzuki formulas. We find that the number of T -gates needed to perform such simulations on a 3 D lattice of n s sites scales at worst as O ( n s 3 / ϵ ) 1 + o ( 1 ) in the thermodynamic limit for position basis simulations and O ( n s 4 + 2 / 3 / ϵ ) 1 + o ( 1 ) in momentum basis. We also find that qubitization scales slightly better with a worst case scaling of O ~ ( n s 2 + 2 / 3 / ϵ ) for lattice eQED and complications in the prepare circuit leads to a slightly worse scaling in momentum basis of O ~ ( n s 5 + 2 / 3 / ϵ ) . We further provide concrete gate counts for simulating a relativistic version of the uniform electron gas that show challenging problems can be simulated using fewer than 10 13 non-Clifford operations and also provide a detailed discussion of how to prepare multi-reference configuration interaction states in effective QED which can provide a reasonable initial guess for the ground state. Finally, we estimate the planewave cutoffs needed to accurately simulate heavy elements such as gold.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2012

Greenhouse Gas Concentration Data Recovery Algorithm for a Low Cost, Laser Heterodyne Radiometer

The goal of a coordinated effort between groups at GWU and NASA GSFC is the development of a low-cost, global, surface instrument network that continuously monitors three key carbon cycle gases in the atmospheric column: carbon dioxide (CO2), methane (CH4), carbon monoxide (CO), as well as oxygen (O2) for atmospheric pressure profiles. The network will implement a low-cost, miniaturized, laser heterodyne radiometer (mini-LHR) that has recently been developed at NASA Goddard Space Flight Center. This mini-LHR is designed to operate in tandem with the passive aerosol sensor currently used in AERONET (a well established network of more than 450 ground aerosol monitoring instruments worldwide), and could be rapidly deployed into this established global network. Laser heterodyne radiometry is a well-established technique for detecting weak signals that was adapted from radio receiver technology. Here, a weak light signal, that has undergone absorption by atmospheric components, is mixed with light from a distributed feedback (DFB) telecommunications laser on a single-mode optical fiber. The RF component of the signal is detected on a fast photoreceiver. Scanning the laser through an absorption feature in the infrared, results in a scanned heterodyne signal io the RF. Deconvolution of this signal through the retrieval algorithm allows for the extraction of altitude contributions to the column signal. The retrieval algorithm is based on a spectral simulation program, SpecSyn, developed at GWU for high-resolution infrared spectroscopies. Variations io pressure, temperature, composition, and refractive index through the atmosphere; that are all functions of latitude, longitude, time of day, altitude, etc.; are modeled using algorithms developed in the MODTRAN program developed in part by the US Air Force Research Laboratory. In these calculations the atmosphere is modeled as a series of spherically symmetric shells with boundaries specified at defined altitudes. Temperature, pressure, and species mixing ratios are defined at these boundaries. Between the boundaries, temperature is assumed to vary linearly with altitude while pressure (and thus gas density) vary exponentially. The observed spectrum at the LHR instrument will be the integration of the contributions along this light path. For any absorption measurement the signal at a particular spectral frequency is a linear combination of spectral line contributions from several species. For each species that might absorb in a spectral region, we have pre-calculated its contribution as a function of temperature and pressure. The integrated path absorption spectrum can then by calculated using the initial sun angle (from location, date, and time) and assumptions about pressure and temperature profiles from an atmospheric model. The modeled spectrum is iterated to match the experimental observation using standard multilinear regression techniques. In addition to the layer concentrations, the numerical technique also provides uncertainty estimates for these quantities as well as dependencies on assumptions inherent in the atmospheric models.

Miller, J. Houston↗