Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical schemes”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Performance-portable Binary Neutron Star Mergers with AthenaK

We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume methods to evolve the fluid and a constrained-transport scheme to preserve the divergence-free constraint for the magnetic field. We also utilize a first-order flux correction (FOFC) scheme to reduce the need for an artificial atmosphere and optionally enforce a maximum principle to improve robustness. We demonstrate the accuracy of AthenaK using a set of standard tests in flat and curved spacetimes. Using a SANE accretion disk around a Kerr black hole, we compare the new solver to the existing solver for stationary spacetimes using the so-called "HARM-like" formulation. We find that both formulations converge to similar results. We also include the first published binary neutron star (BNS) mergers performed on graphical processing units (GPUs). Thanks to the FOFC scheme, our BNS mergers maintain a relative error of $\mathcal{O}$(10 –11 ) or better in baryon mass conservation up to collapse. Finally, we perform scaling tests of AthenaK on OLCF Frontier, where we show excellent weak scaling of ≥80% efficiency up to 32,768 GPUs and 74% up to 65,536 GPUs for a GRMHD problem in dynamical spacetimes with six levels of mesh refinement. AthenaK achieves an order-of-magnitude speedup using GPUs compared to CPUs, demonstrating that it is suitable for performing numerical relativity problems on modern exascale resources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Effect of dust on rainfall over the Red Sea coast based on WRF-Chem model simulations

Water is the single most important element of life. Rainfall plays an important role in the spatial and temporal distribution of this precious natural resource, and it has a direct impact on agricultural production, daily life activities, and human health. One of the important elements that govern rainfall formation and distribution is atmospheric aerosol, which also affects the Earth's radiation balance and climate. Therefore, understanding how dust compositions and distributions affect the regional rainfall pattern is crucial, particularly in regions with high atmospheric dust loads such as the Middle East. Although aerosol and rainfall research has garnered increasing attention as both an independent and interdisciplinary topic in the last few decades, the details of various direct and indirect pathways by which dust affects rainfall are not yet fully understood. Here, we explored the effects of dust on rainfall formation and distribution as well as the physical mechanisms that govern these phenomena, using high-resolution WRF-Chem simulations (~1.5 km × 1.5 km) configured with an advanced double-moment cloud microphysics scheme coupled with a sectional eight-bin aerosol scheme. Our model-simulated results were realistic, as evaluated from multiple perspectives including vertical profiles of aerosol concentrations, aerosol size distributions, vertical profiles of air temperature, diurnal wind cycles, and spatio-temporal rainfall patterns. Rainfall over the Red Sea coast is mainly caused by warm rain processes, which are typically confined within a height of ~6 km over the Sarawat mountains and exhibit a strong diurnal cycle that peaks in the evening at approximately 18:00 local time under the influence of sea breezes. Numerical experiments indicated that dust could both suppress or enhance rainfall. The effect of dust on rainfall was calculated as total, indirect, and direct effects, based on 10-year August-average daily-accumulated rainfall over the study domain covering the eastern Red Sea coast. For extreme rainfall events (domain-average daily-accumulated rainfall of ≥ 1.33 mm), the net effect of dust on rainfall was positive or enhancement (6.05 %), with the indirect effect (4.54 %) and direct effect (1.51 %) both causing rainfall increase. At a 5 % significance level, the total and indirect effects were statistically significant whereas the direct effect was not. For normal rainfall events (domain-average daily-accumulated rainfall < 1.33 mm), the indirect effect enhanced rainfall (4.76 %) whereas the direct effect suppressed rainfall (-5.78 %), resulting in a negative net suppressing effect (-1.02 %), all of which were statistically significant. We investigated the possible physical mechanisms of the effects and found that the rainfall suppression by dust direct effects was mainly caused by the scattering of solar radiation by dust. The surface cooling induced by dust weakens the sea breeze circulation, which decreases the associated landward moisture transport, ultimately suppressing rainfall. For extreme rainfall events, dust causes net rainfall enhancement through indirect effects as the high dust concentration facilitates raindrops to grow when the water vapor is sufficiently available. Our results have broader scientific and environmental implications. Specifically, although dust is considered a problem from an air quality perspective, our results highlight the important role of dust on sea breeze circulation and associated rainfall over the Red Sea coastal regions. Our results also have implications for cloud seeding and water resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Near-zero-index ultra-fast pulse characterization

Transparent conducting oxides exhibit giant optical nonlinearities in the near-infrared window where their linear index approaches zero. Despite the magnitude and speed of these nonlinearities, a “killer” optical application for these compounds has yet to be found. Because of the absorptive nature of the typically used intraband transitions, out-of-plane configurations with short optical paths should be considered. In this direction, we propose an alternative frequency-resolved optical gating scheme for the characterization of ultra-fast optical pulses that exploits near-zero-index aluminium zinc oxide thin films. Besides the technological advantages in terms of manufacturability and cost, our system outperforms commercial modules in key metrics, such as operational bandwidth, sensitivity, and robustness. The performance enhancement comes with the additional benefit of simultaneous self-phase-matched second and third harmonic generation. Because of the fundamental importance of novel methodologies to characterise ultra-fast events, our solution could be of fundamental use for numerous research labs and industries.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Qubit-Oscillator Concatenated Codes: Decoding Formalism and Code Comparison

Concatenating bosonic error-correcting codes with qubit codes can substantially boost the errorcorrecting power of the original qubit codes. It is not clear how to concatenate optimally, given that there are several bosonic codes and concatenation schemes to choose from, including the recently discovered Gottesman-Kitaev-Preskill (GKP) – stabilizer codes [Phys. Rev. Lett. 125, 080503 (2020)] that allow protection of a logical bosonic mode from fluctuations of the conjugate variables of the mode. We develop efficient maximum-likelihood decoders for and analyze the performance of three different concatenations of codes taken from the following set: qubit stabilizer codes, analog or Gaussian stabilizer codes, GKP codes, and GKP-stabilizer codes. We benchmark decoder performance against additive Gaussian white noise, corroborating our numerics with analytical calculations. We observe that the concatenation involving GKP-stabilizer codes outperforms the more conventional concatenation of a qubit stabilizer code with a GKP code in some cases. We also propose a GKP-stabilizer code that suppresses fluctuations in both conjugate variables without extra quadrature squeezing and formulate qudit versions of GKP-stabilizer codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Effect of heating/cooling dynamics in the hysteresis loop and tunable IR emissivity of VO 2 thin films

We experimentally investigate the semiconductor-to-metal transition (SMT) in vanadium dioxide thin films using an infrared thermographic technique. During the semiconductor to metal phase change process, VO 2 optical properties dynamically change and infrared emission undergoes a hysteresis loop due to differences between heating and cooling stages. The shape of the hysteresis loop was accurately monitored under different dynamic heating/cooling rates. In order to quantify and understand the effects of different rates, we used a numerical modelling approach in which a VO 2 thin layer was modeled as metamaterial. The main experimental findings are interpreted assuming that both the rate of formation and shape of metallic inclusions are tuned with the heating/cooling rate. The structural transition from monoclinic to tetragonal phases is the main mechanism for controlling the global properties of the phase transition. However, our experimental results reveal that the dynamics of the heating/cooling process can become a useful parameter for further tuning options and lays out a macroscopic optical sensing scheme for the microscopic phase change dynamics of VO 2 . Our study sheds light on phase-transition dynamics and their effect on the infrared emission spectra of VO 2 thin films, therefore enabling the heating/cooling rate to be an additional parameter to control infrared emission characteristics of thermal emitters. The hysteresis loop represents the phase coexistence region, thus being of fundamental importance for several applications, such as the operation of radiative thermal logic elements based on phase transition materials. For such applications, the phase transition region is shifted for heating and cooling processes. We also show that, depending on the way the phase change elements are heated, the temperature operation range will be slightly modified.

42 ENGINEERING↗
DOE OSTI2020

Coupling reduced‐order modeling and coarse‐grained CFD‐DEM to accelerate coal gasifier simulation and optimization

Abstract High fidelity three‐dimensional (3‐D) numerical simulations of commercial‐scale coal gasifiers are very time consuming and expensive. This article proposes a reduced‐order modeling approach that uses quasi one‐dimensional (1‐D) CFD–DEM simulation results to serve as an accurate initial condition for the 3‐D simulation, which significantly shortens the physical time from several hours to minutes to achieve steady states. The 3‐D simulation employs a previously validated coarse‐grained CFD‐DEM approach, which further accelerates the simulation. Both 1‐D and 3‐D simulations adopt a particle shrinkage scheme to correctly account for the particle volume change due to chemical reaction and thus provides an accurate bed height. Comparison with available experimental measurements of the bed temperature and synthesis gas composition is used to validate the simulation. The final syngas composition and flowrates are strongly affected by the fuel or gasification agent rate and process operating conditions as expected. This approach is then used to optimize two important operating conditions (air and steam flows) to achieve specific gasifier performance goals for the production of Fischer–Tropsch synthesis gas and production of fuel gas for gas turbines.

Yu, Jia↗
DOE OSTI2020

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗
DOE OSTI2020

Application of Numerical Atmospheric Models

In this study, a numerical modeling framework for simulating extreme storm events was established using the Weather Research and Forecasting (WRF) model. Such a framework is necessary for the derivation of engineering parameters such as probable maximum precipitation that are the cornerstone of large water-management infrastructure design. Here, this framework was built based on a heavy storm that occurred in Nashville, Tennessee (U.S.), in 2010, and verified using two other extreme storms. To achieve the optimal setup, several combinations of model resolutions, initial/boundary conditions (IC/BC), cloud microphysics, and cumulus parameterization schemes were evaluated using multiple metrics of precipitation characteristics. The evaluation suggests that WRF is most sensitive to the IC/BC option. Simulation generally benefits from finer resolutions up to 5 km. At the 15 km level, NCEP2 IC/BC produces better results, whereas NAM IC/BC performs best at the 5 km level. The recommended model configuration from this study is: NAM or NCEP2 IC/BC (depending on data availability), 15 km or 15–5 km nested grids, Morrison microphysics, and Kain-Fritsch cumulus schemes. Validation of the optimal framework suggests that these options are good starting choices for modeling extreme events similar to the test cases. This optimal framework is proposed in response to emerging engineering demands of extreme storm event forecasting and analyses for design, operations, and risk assessment of large water infrastructures.

Chen, Xiaodong↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2021

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗
DOE OSTI2026

Dynamical Downscaling of Earth System Model Data for Energy System Analysis

Assessing energy resources (e.g., solar, wind, and hydro) under future scenarios requires datasets with sufficient spatial and temporal detail to capture variability and extreme events. While global-scale Earth System Model (ESM) projections are widely used, their coarse resolution limits direct application to regional energy system analyses. Dynamical downscaling offers a robust approach to generate physically consistent, fine-scale datasets that better represent local atmospheric processes impacting energy resources. In this work, we present a two-stage approach for producing high-resolution historical and future projections over the contiguous United States (CONUS). First, we optimize the Weather Research and Forecasting (WRF) model configuration for energy-relevant variables - solar irradiance, wind speed, and precipitation - by conducting ERA5-driven simulations at 8-km and 28-km resolution. Multiple physics schemes and model configurations within the WRF are evaluated against observational datasets including the National Solar Radiation Database (NSRDB), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the Stage IV multi-radar/multi-sensor precipitation product for the CONUS domain. Using the best-performing configuration, we dynamically downscale MPI-ESM1-2-HR simulations for 2000-2060 under SSP2-4.5 and SSP5-8.5 scenarios at 4-km spatial and hourly temporal resolution. This presentation will provide a comprehensive analysis of the results from multiple numerical experiments and high-resolution ESM projections. In addition, we will discuss potential applications of our high-resolution datasets within the energy sector and outline future research avenues dedicated to evaluating how extreme weather events influence system performance and resilience.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Reflectivity and spectrum of relativistic flying plasma mirrors

Flying plasma mirrors induced by intense lasers via laser wake field acceleration scheme have been proposed as a promising way to generate few-cycle EUV or x-ray lasers. In addition, if such a relativistic plasma mirror can accelerate, then it would serve as an analog black hole to investigate the information loss paradox associated with the black hole Hawking evaporation. Among these applications, the reflectivity, which is usually frequency-dependent, would affect the outgoing photon spectrum and, therefore, impact on the analysis of the physics under investigation. In this paper, these two issues are investigated analytically and numerically with one-dimensional particle-in-cell simulations. Based on our simulation results, we propose a new model that provides a better estimate of the reflectivity than those studied previously. Moreover, we found that the peak frequency of the reflected spectrum of a Gaussian incident wave deviates from the expected value, 4γ 2 ω, due to the dependence of reflectivity on the frequency of the incident wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Strictly Enforcing Invertibility and Conservation in CNN-Based Super Resolution for Scientific Datasets

Abstract Recently, deep convolutional neural networks (CNNs) have revolutionized image “super resolution” (SR), dramatically outperforming past methods for enhancing image resolution. They could be a boon for the many scientific fields that involve imaging or any regularly gridded datasets: satellite remote sensing, radar meteorology, medical imaging, numerical modeling, and so on. Unfortunately, while SR-CNNs produce visually compelling results, they do not necessarily conserve physical quantities between their low-resolution inputs and high-resolution outputs when applied to scientific datasets. Here, a method for “downsampling enforcement” in SR-CNNs is proposed. A differentiable operator is derived that, when applied as the final transfer function of a CNN, ensures the high-resolution outputs exactly reproduce the low-resolution inputs under 2D-average downsampling while improving performance of the SR schemes. The method is demonstrated across seven modern CNN-based SR schemes on several benchmark image datasets, and applications to weather radar, satellite imager, and climate model data are shown. The approach improves training time and performance while ensuring physical consistency between the super-resolved and low-resolution data. Significance Statement Recent advancements in using deep learning to increase the resolution of images have substantial potential across the many scientific fields that use images and image-like data. Most image super-resolution research has focused on the visual quality of outputs, however, and is not necessarily well suited for use with scientific data where known physics constraints may need to be enforced. Here, we introduce a method to modify existing deep neural network architectures so that they strictly conserve physical quantities in the input field when “super resolving” scientific data and find that the method can improve performance across a wide range of datasets and neural networks. Integration of known physics and adherence to established physical constraints into deep neural networks will be a critical step before their potential can be fully realized in the physical sciences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗
DOE OSTI2026

Ion Trapping Studies and Mitigation Strategies for the EIC ERL-Based Strong Hadron Cooler

An Energy Recovery Linac based strong hadron cooler was previously considered for the Electron-Ion Collider. The required electron beam parameters for variable-energy strong hadron cooling place significant constraints on ion trapping and collective effects. This paper presents initial studies of these constraints through a combination of analytical modelling and numerical simulations of ion production, trapping behaviour, and mitigation strategies. A multi-bunch tracking framework based on ELEGANT with the ionEffects module is used to simulate machine operation over millisecond time scales, corresponding to more than 3 × 10^5 electron bunches. The simulations include modelling of ionisation processes together with transverse electron–ion dynamics, allowing the evolution and accumulation of ions to be investigated. Analytical expressions based on Gaussian beam distributions are used to estimate ion trapping conditions and benchmark the simulation results. A bi-periodic bunch spacing scheme is also investigated as a possible mitigation method by detuning the ion oscillation frequency. These studies provide an initial assessment of ion trapping in the strong hadron cooler and demonstrate possible approaches for reducing beam–ion effects.

Bi, R. [Lancaster University, Cockcroft Institute]↗
DOE OSTI2020

Engineered thermalization and cooling of quantum many-body systems

We develop a scheme for engineering genuine thermal states in analog quantum simulation platforms by coupling local degrees of freedom to driven, dissipative ancilla pseudospins. We demonstrate the scheme in a many-body quantum spin lattice simulation setting. A Born-Markov master equation describing the dynamics of the many-body system is developed, and we show that if the ancilla energies are periodically modulated, with a carefully chosen hierarchy of timescales, one can effectively thermalize the many-body system. Through analysis of the time-dependent dynamical generator, we determine the conditions under which the true thermal state is an approximate dynamical fixed point for general system Hamiltonians. Finally, we evaluate the thermalization protocol through numerical simulation and discuss prospects for implementation on current quantum simulation hardware.

36 MATERIALS SCIENCE↗