Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “balance”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A rigorous physics-based enhanced parameter estimation (EPE) methodology for calibration of building energy simulations

Buildings rarely perform as designed/simulated and there are numerous tangible benefits if this gap is reconciled. A new scientifically rigorous yet pragmatic methodology for calibrating building energy simulations - called Enhanced Parameter Estimation (EPE) - is proposed that allows physically relevant parameter estimation rather than a blind force-fit to energy use data. Starting with a rapidly created simulation model, calibration is performed in two stages: (a) building shell calibration with the HVAC system replaced by an ideal system that meets the loads (b) HVAC system calibration with the building shell and all internal loads replaced by a box with only process loads. In the first stage, EPE identifies a small number of high-level heat flows in the energy balance, calculates them with specifically tailored individual driving functions, introduces physically significant parameters to best accomplish energy balance, and, estimates the parameters and their uncertainty bounds. Calibration is thus done with corrective heat flows without any arbitrary tuning of input parameters. Calibration accuracy is enhanced by machine learning of the residual errors. The EPE methodology is demonstrated by means of: a synthetic building and an actual 75,000 sq. ft. building in Pennsylvania. Future work needed for widespread application is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

On the non-existence of stepped-pressure equilibria far from symmetry

The Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502) has been successful in the construction of equilibria in 3D configurations that contain a mixture of flux surfaces, islands and chaotic magnetic field lines. In this model, the plasma is sliced into sub-volumes separated by ideal interfaces, and in each volume the magnetic field is a Beltrami field. In the cases where the system is far from possessing a continuous symmetry, such as in stellarators, the existence of solutions to a stepped-pressure equilibrium with given constraints, such as a multi-region relaxed MHD minimum energy state, is not guaranteed but is often taken for granted. Using SPEC, we have studied two different scenarios in which a solution fails to exist in a slab with analytic boundary perturbations. Here we found that with a large boundary perturbation, a certain interface becomes fractal, corresponding to the break up of a Kolmogorov–Arnold–Moser (KAM) surface. Moreover, an interface can only support a maximum pressure jump while a solution of the magnetic field consistent with the force balance condition can be found. An interface closer to break-up can support a smaller pressure jump. We discovered that the pressure jump can push the interface closer to being non-smooth through force balance, thus significantly decreasing the maximum pressure it can support. Our work shows that a convergence study must be performed on a SPEC equilibrium with interfaces close to break-up. These results may also provide insights into the choice of interfaces and have applications in finding out the maximum pressure a machine can support.

3D equilibrium↗
DOE OSTI2024

Core inductive electric field during sawtooth crashes on DIII-D

Sawtooth crashes on tokamak plasmas exhibit relaxation much faster than resistive time scales via a mechanism not fully understood. Using core magnetic measurements from the Radial Interferometer-Polarimeter (RIP) diagnostic on the DIII-D tokamak, Grad–Shafranov equilibria constrained by internal magnetic measurements that have high time resolution (<1μs) can be computed, allowing analysis of how equilibrium parameters such as safety factor q, current density J, and parallel electric field E ∥ , particularly on-axis, evolve. At the sawtooth crash, on-axis safety factor q 0 is observed to rise by 5% but remain below 1 throughout the cycle, and on-axis current density J 0 is observed to drop by 5%. On-axis parallel electric field E ∥ (0) is found to be balanced by ηJ 0 (resistivity times on-axis current density) except during the 200 µs crash period, where E ∥ (0) reaches 22 V m –1 , exceeding ηJ 0 by a factor of more than 2000. These first measurements in tokamak plasmas verify that generalized Ohm's law is not balanced during the crash by resistive effects alone; this is a finding expected due to the relaxation being much faster than resistive timescales. As a result, measurement of the electric field during the tokamak sawtooth serves to illuminate the physical mechanisms at work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of collisional ion orbit loss on neoclassical tokamak radial electric fields

We report that ion orbit loss is considered important for generating the radially inward electric field E r in a tokamak edge plasma. In particular, this effect is emphasized in diverted tokamaks with a magnetic X point. In neoclassical equilibria, Coulomb collisions can scatter ions onto loss orbits and generate a radially outward current, which in steady state is balanced by the radially inward current from viscosity. To quantitatively measure this loss-orbit current in an edge pedestal, an ion-orbit-flux diagnostic has been implemented in the axisymmetric version of the gyrokinetic particle-in-cell code XGC. As the first application of this diagnostic, a neoclassical DIII-D H-mode plasma is studied using gyrokinetic ions and adiabatic electrons. The validity of the diagnostic is demonstrated by studying the collisional relaxation of E r in the core. After this demonstration, the loss-orbit current is numerically measured in the edge pedestal in quasisteady state. In this plasma, it is found that the radial electric force on ions from E r approximately balances the ion radial pressure gradient in the edge pedestal, with the radial force from the plasma flow term being a minor component. The effect of orbit loss on E r is found to be only mild.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The role of ETG modes in JET–ILW pedestals with varying levels of power and fuelling

Abstract We present the results of GENE gyrokinetic calculations based on a series of JET–ITER-like-wall (ILW) type I ELMy H-mode discharges operating with similar experimental inputs but at different levels of power and gas fuelling. We show that turbulence due to electron-temperature-gradient (ETGs) modes produces a significant amount of heat flux in four JET–ILW discharges, and, when combined with neoclassical simulations, is able to reproduce the experimental heat flux for the two low gas pulses. The simulations plausibly reproduce the high-gas heat fluxes as well, although power balance analysis is complicated by short ELM cycles. By independently varying the normalised temperature gradients ( ω T e ) and normalised density gradients ( ω n e ) around their experimental values, we demonstrate that it is the ratio of these two quantities η e = ω T e / ω n e that determines the location of the peak in the ETG growth rate and heat flux spectra. The heat flux increases rapidly as η e increases above the experimental point, suggesting that ETGs limit the temperature gradient in these pulses. When quantities are normalised using the minor radius, only increases in ω T e produce appreciable increases in the ETG growth rates, as well as the largest increases in turbulent heat flux which follow scalings similar to that of critical balance theory. However, when the heat flux is normalised to the electron gyro-Bohm heat flux using the temperature gradient scale length L T e , it follows a linear trend in correspondence with previous work by different authors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of collisional ion orbit loss on tokamak radial electric field and toroidal rotation in an L-mode plasma

Abstract Ion orbit loss has been used to model the formation of a strong negative radial electric field E r in the tokamak edge, as well as edge momentum transport and toroidal rotation. To quantitatively measure ion orbit loss, an orbit-flux formulation has been developed and numerically applied to the gyrokinetic particle-in-cell code XGC. We study collisional ion orbit loss in an axisymmetric DIII-D L-mode plasma using gyrokinetic ions and drift-kinetic electrons. Numerical simulations, where the plasma density and temperature profiles are maintained through neutral ionization and heating, show the formation of a quasisteady negative E r in the edge. We have measured a radially outgoing ion gyrocenter flux due to collisional scattering of ions into the loss orbits, which is balanced by the radially incoming ion gyrocenter flux from confined orbits on the collisional time scale. This suggests that collisional ion orbit loss can shift E r in the negative direction compared to that in plasmas without orbit loss. It is also found that collisional ion orbit loss can contribute to a radially outgoing (counter-current) toroidal-angular-momentum flux, which is not balanced by the toroidal-angular-momentum flux carried by ions on the confined orbits. Therefore, the edge toroidal rotation shifts in the co-current direction on the collisional time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal quench induced by a composite pellet-produced plasmoid

Injecting shattered pellets is the critical concept of the envisaged ITER disruption mitigation system (DMS). Rapid deposition of large amounts of material should presumably result in controlled cooling of the entire plasma. A considerable transfer of thermal energy from the electrons of the background plasma to the ions accompanies a localized material injection due to the ambipolar expansion along the magnetic field line of the cold and dense plasmoid produced by the ablated pellet. Radiation initially plays the dominant role in the energy balance of a composite plasmoid containing high-Z impurities. A competition between the ambipolar expansion and the radiative losses defines the Thermal Quench scenario, including the amount of pre-quench thermal energy radiated on a short collisional timescale—possibly detrimental for the plasma-facing components. The present work quantifies plasmoid energy balance for disruption mitigation parameters. For pure hydrogen injection, up to 90% of the pre-pellet electron thermal energy may go to the newly injected ions. We also demonstrate that a moderate high-Z impurity content within the plasmoid can reduce highly localized radiation at the beginning of the expansion. The thermal energy will then dissipate on the much longer ion collisional timescale, which would be attractive for ITER DMS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigations of atomic and molecular processes of NBI-heated discharges in the MAST Upgrade Super-X divertor with implications for reactors

This experimental study presents an in-depth investigation of the performance of the MAST-U Super-X divertor during NBI-heated operation (up to 2.5 MW) focussing on volumetric ion sources and sinks as well as power losses during detachment. The particle balance and power loss analysis revealed the crucial role of Molecular Activated Recombination and Dissociation (MAR and MAD) ion sinks in divertor particle and power balance, which remain pronounced in the change from ohmic to higher power (NBI heated) L-mode conditions. The importance of MAR and MAD remains with double the absorbed NBI heating. MAD results in significant power dissipation (up to ${\sim} 20\%$ of $P_\textrm{SOL}$), mostly in the cold ($T_e \lt 5$ eV) detached region. Theoretical and experimental evidence is found for the potential contribution of $D^-$ to MAR and MAD, which warrants further study. These results suggest that MAR and MAD can be relevant in higher power conditions than the ohmic conditions studied previously. Post-processing reactor-scale simulations suggests that MAR and MAD can play a significant role in divertor physics and synthetic diagnostic signals of reactor-scale devices, which are currently underestimated in exhaust simulations. This raises implications for the accuracy of reactor-scale divertor simulations of particularly tightly baffled (alternative) divertor configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗
DOE OSTI2021

Albedo-induced global warming impact of Conservation Reserve Program grasslands converted to annual and perennial bioenergy crops

Abstract Climate benefit assessments of bioenergy crops often focus on biogeochemical impacts, paying little if any attention to biogeophysical impacts. However, land conversions required for large-scale bioenergy crop production are substantial and may directly affect the climate by altering surface energy balance. In the US, such land conversions are likely to be met in part by converting Conservation Reserve Program (CRP) grassland to bioenergy crops. Here, we converted three 22 year old CRP smooth brome grass fields into no-till corn, switchgrass, or restored prairie bioenergy crops. We assessed the biogeophysical climate impact of the conversions using albedo changes relative to unconverted reference CRP grassland. The corn and perennial fields had higher annual albedo than the grassland they replaced—causing cooling of the local climate. The cooling of the corn field occurred solely during the non-growing season—especially when surfaces were snow-covered, whereas the cooling of the perennial fields was more prominent during the growing season. Compared to biogeochemical impacts with fossil fuel offsets for the same land conversions over eight years, the annual albedo-induced climate benefits add ∼35% and ∼78% to the annual biogeochemical benefits provided from the switchgrass and restored prairie fields, respectively, and offset ∼3.3% of the annual greenhouse gas (GHG) emissions from the corn field. We conclude that albedo-induced climate mitigation from conversion of CRP lands to perennial but not annual bioenergy crops can be substantial, and future climate impact assessments of bioenergy crops should include albedo changes in addition to GHG balances in order to better inform climate policies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Strong influence of convective heat transfer efficiency on the cooling benefits of green roof irrigation

Abstract By enhancing evapotranspiration (ET), green roofs provide cooling benefits for the urban environment and are recognized as a promising heat mitigation strategy. The evaporative cooling effects of green roofs strongly depend on the soil moisture conditions and thus irrigation may be needed to sustain the cooling benefits. It has been shown that the magnitude of cooling benefits offered by green roof irrigation varies spatially, but its controlling factors remain elusive. In this study, we combine a surface energy balance (SEB) model with global simulations generated by an improved Earth System Model to illustrate the key factors controlling the cooling benefits of green roof irrigation. We employ a simple irrigation scheme, which is only active when there is no ice in soil layers and when the soil moisture is below field capacity. As a result, most of the irrigation water leaves the green roof system via ET. We find that the magnitude and also the spatial variability of the cooling benefits of green roof irrigation are controlled by the irrigation amount, and a surface energy redistribution factor that encodes the efficiencies of different SEB components in transferring heat. Further analysis indicates that the enhancement of latent heat flux due to irrigation is largely balanced by the reduction of sensible heat flux on green roofs. Therefore, the amount of irrigation needed per unit decrease of green roof surface temperature is mainly controlled by the convective heat transfer efficiency. A lower convective heat transfer efficiency (e.g. under a lower wind speed) helps reduce the amount of irrigation needed per unit decrease of green roof surface temperature. This study highlights the importance of SEB in constraining the cooling benefits of green roof irrigation and provides valuable guidance for urban planners and policy makers in terms of heat mitigation and sustainable water management.

Wang, Linying (ORCID:0000000308769220)↗
DOE OSTI2022

Rapidly changing high-latitude seasonality: implications for the 21st century carbon cycle in Alaska

Seasonal variations in high-latitude terrestrial carbon (C) fluxes are predominantly driven by air temperature and radiation. At present, high-latitude net C uptake is largest during the summer. Recent observations and modeling studies have demonstrated that ongoing and projected climate change will increase plant productivity, microbial respiration, and growing season lengths at high-latitudes, but impacts on high-latitude C cycle seasonality (and potential feedbacks to the climate system) remain uncertain. Here we use ecosys, a well-tested and process-rich mechanistic ecosystem model that we evaluate further in this study, to explore how climate warming under an RCP8.5 scenario will shift C cycle seasonality in Alaska throughout the 21st century. The model successfully reproduced recently reported large high-latitude C losses during the fall and winter and yet still predicts a high-latitude C sink, pointing to a resolution of the current conflict between process-model and observation-based estimates of high-latitude C balance. We find that warming will result in surprisingly large changes in net ecosystem exchange (NEE; defined as negative for uptake) seasonality, with spring net C uptake overtaking summer net C uptake by year 2100. This shift is driven by a factor of 3 relaxation of spring temperature limitation to plant productivity that results in earlier C uptake and a corresponding increase in magnitude of spring NEE from –19 to –144 gC m –2 season –1 by the end of the century. Although a similar relaxation of temperature limitation will occur in the fall, radiation limitation during those months will limit increases in C fixation. Additionally, warmer soil temperatures and increased carbon inputs from plants lead to combined fall and winter C losses (163 gC m –2 ) that are larger than summer net uptake (123 gC m –2 season –1 ) by year 2100. However, this increase in microbial activity leads to more rapid N cycling and increased plant N uptake during the fall and winter months that supports large increases in spring NPP. Due to the large increases in spring net C uptake, the high-latitude atmospheric C sink is projected to sustain throughout this century. Our analysis disentangles the effects of key environmental drivers of high-latitude seasonal C balances as climate changes over the 21st century.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluating power grid model hydropower feasibility with a river operations model

Production cost models (PCMs) simulate dispatch of generators across a large power grid and are used widely by planners to study the reliability of electricity supply. As energy systems transition away from the thermoelectric technologies that have traditionally balanced electricity supply and demand, hydropower and its representation in PCMs is of increasing importance. A limitation of PCMs applied to continental power grids with diverse generation portfolios is that hydropower generation is simulated without full consideration of complex river dynamics, leading to possible misrepresentation of grid flexibility and performance. In addition, data used in PCMs may reflect outdated operating policies. In this paper we propose a hydropower generation feasibility test for PCMs. The approach uses a detailed hydropower model to determine whether the hourly hydropower schedule from a PCM with simplified monthly parameterization can be attained after accounting for realistic river dynamics and operating policies, such as spill requirements and general water movement and balance through a cascade reservoir system. We perform this hydropower generation test for the ‘Big 10’ hydropower system on the Columbia River (part of the Western Interconnect of the United States), revealing 9% overestimation of available hydropower generation in a PCM solution in an average hydrologic year. Our evaluation provides insight on the cost and opportunities for better representing hydropower in PCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cryogenic aspects of a 20 MW class low-temperature superconducting generator for the renewables industry

In this paper we give a progress update on the cryogenic design for cooling a 20 MW class, partially superconducting generator with stationary field coils for offshore wind renewables industry [1-2]. This is a continuation of an earlier program on a 10 MW system dating back ten years. Whereas this new power rating increase leads to a radial diameter expansion of the superconducting field coils from 4 to > 9.5 m, it maintains the axial length. We show the design based on this scaled up with enlarged diameter. The field coil size increase asks for higher cooling power that leads to a bigger cold box size to accommodate the cryocoolers. In the design, as many as 8 cryocoolers can be accessed and serviced from the nacelle that houses the main cryogenic components. A typical thermal load balance sheet with all components is given and compared against the available cryocooler cooling power at different operating conditions. Due to the cryocoolers’ local point of contact cooling within the nacelle, the temperature gradient of the thermal shield that fully encloses the cold mass with its embedded field coils needs to be balanced out to minimize the heat burden on the field coils. This requires additional analytical efforts and implementation of further design features. The extended nacelle houses the cold box with its cryogenic infrastructure. The interface for the envisaged cryogenic pushbutton closed-loop circulating system remains invisible, requires no handling of cryogenic liquids and is hermetically closed. The field coil diameter increase also leads to a greater initial helium gas storage volume. In addition, it requires higher initial room temperature fill pressure within the toroidal helium vapor storage tanks and in cooling tubes connecting to the cryocoolers. The toroidal helium gas tanks are thermally coupled with the thermal shield so that helium convection inside the storage tank can improve heat transfer and reduce the temperature gradient of the thermal shield. Besides those heat transfer challenges, additional mechanical strain within this large structure is exerted on the torque tubes during initial cooldown and when energizing field coils. Some of those design challenges are quite unexpected, leading to novel workarounds in order to maintain the chosen cooling strategy. Finally, we assess those design limitations in view of further cryogenic scalability with emphasis on manufacturability and assembly.

cryogenics↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Physics makes the difference: Bayesian optimization and active learning via augmented Gaussian process

Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.

42 ENGINEERING↗