Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Baffle aperture design study of hollow cathode equipped ion thrusters

A simple theoretical model which can be used as an aid in the design of the baffle aperture region of a hollow cathode equipped ion thruster was developed. An analysis of the ion and electron currents in both the main and cathode discharge chambers is presented. From this analysis a model of current flow through the aperture, which is required as an input to the design model, was developed. This model was verified experimentally. The dominant force driving electrons through the aperture was the force due to the electrical potential gradient. The diffusion process was modeled according to the Bolm diffusion theory. A number of simplifications were made to limit the amount of detailed plasma information required as input to the model to facilitate the use of the model in thruster design. This simplified model gave remarkably consistant results with experimental results obtained with a given thruster geometry over substantial changes in operating conditions. The model was uncertain to about a factor of two for different thruster cathode region geometries. The design usefulness was limited by this factor of two uncertainty and by the accuracy to which the plasma parameters required as inputs to the model were specified.

Brophy, J. R., Jr.↗
NASA NTRS1984

Calculations of turbulent mass transport in a bluff-body diffusion-flame combustor

Experimental measurements of turbulent mass transit in a bluff-body diffusion-flame combustor (using CO2 instead of fuel) are analyzed to evaluate the accuracy of physical models used in flow computations. The data of Lightman and Magill (1981), Lightman et al. (1983), and Roquemore et al. (1983) are used to calculate apparent turbulent Schmidt numbers (TSN) for a series of flow conditions by a modified TEACH viscous-flow code. The modeling principles, calculation grid, and boundary conditions are discussed, and the results are presented in graphs comparing calculated and measured values. It is shown that models such as the two-equation (K - epsilon) model which use a single value of the TSN are inappropriate for conditions typical of gas-turbine combustors: local TSN variations due to turbulence behavior must be taken into account in improved models.

Sturgess, G. J.↗
NASA NTRS1985

The universal spectrum of AGNs and QSOs

The effects of the feedback of e(+)-e(-) pair reinjection in a plasma due to photon photon absorption of its own radiation are examined. A mechanism is presented which can produce an electron distribution function that can account for the overall spectral distribution of radiation of AGNs and QSOs and the specific slopes observed in the IR to UV and 2-50 keV bands. It is interesting to note that the necessary condition for this mechanism to work (i.e., most of energy injected at e(M sub e)(C sup 2) is realized in the accretion shock model of Kazanas and Ellison. This mechanism involves only one free parameter the compactness of the sources, L/R, whose mean value can also account for the diffuse gamma ray background in terms of AGNs.

Kazanas, D.↗
NASA NTRS1988

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

This paper will describe the experimental and analytical work that has been done to establish justification and feasibility for a Shuttle mid-deck experiment involving mass transfer between a gas bubble and a liquid. The experiment involves the observation and measurement of the dissolution of an isolated, immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble have been successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration has been accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model has been developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity. The results will yield information on transport under conditions of pure diffusion.

Dewitt, Kenneth J.↗
NASA NTRS2001

High Pressure, Transport Properties of Fluids: Theory and Data from Levitated Fluid-Drops at Combustion-Relevant Temperatures

The objective of this investigation is to derive a set of consistent mixing rules for calculating diffusivities and thermal diffusion factors over a thermodynamic regime encompassing the subcritical and supercritical ranges. These should serve for modeling purposes, and therefore for accurate simulations of high pressure phenomena such as fluid disintegration, turbulent flows and sprays. A particular consequence of this work will be the determination of effective Lewis numbers for supercritical conditions, thus enabling the examination of the relative importance of heat and mass transfer at supercritical pressures.

Bellan, J.↗
DOE OSTI2023

Efficient 3D Molecular Design with an E(3) Invariant Transformer VAE

This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. Here, we show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.

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

Multiscale framework for predicting the coupling between deformation and fluid diffusion in porous rocks

In this project, a predictive multiscale framework will be developed to simulate the strong coupling between solid deformations and fluid diffusion in porous rocks. We intend to improve macroscale modeling by incorporating fundamental physical modeling at the microscale in a computationally efficient way. This is an essential step toward further developments in multiphysics modeling, linking hydraulic, thermal, chemical, and geomechanical processes. This research will focus on areas where severe deformations are observed, such as deformation bands, where classical phenomenology breaks down. Multiscale geometric complexities and key geomechanical and hydraulic attributes of deformation bands (e.g., grain sliding and crushing, and pore collapse, causing interstitial fluid expulsion under saturated conditions), can significantly affect the constitutive response of the skeleton and the intrinsic permeability. Discrete mechanics (DEM) and the lattice Boltzmann method (LBM) will be used to probe the microstructure---under the current state---to extract the evolution of macroscopic constitutive parameters and the permeability tensor. These evolving macroscopic constitutive parameters are then directly used in continuum scale predictions using the finite element method (FEM) accounting for the coupled solid deformation and fluid diffusion. A particularly valuable aspect of this research is the thorough quantitative verification and validation program at different scales. The multiscale homogenization framework will be validated using X-ray computed tomography and 3D digital image correlation in situ at the Advanced Photon Source in Argonne National Laboratories. Also, the hierarchical computations at the specimen level will be validated using the aforementioned techniques in samples of sandstone undergoing deformation bands.

58 GEOSCIENCES↗
DOE OSTI2021

Scattering by whistler-mode waves during a quiet period perturbed by substorm activity

We study the dynamics of radiation belt electrons during a 10-day quiet period perturbed by substorm activity and preceding a high-speed stream (HSS), aiming at a global description of the radiation belts in L-shell, L in [2, 6], and energy [0.1, 10] MeV. We combine Van Allen Probes observations and Fokker-Planck numerical simulations of pitch-angle diffusion. The Fokker-Planck model uses event-driven pitch angle diffusion coefficients from whistler-mode waves, built from the wave properties and the ambient plasma density measurements from the Van Allen Probes. We first find this event has some similar characteristics to regular quiet times previously studied; a widely extended plasmasphere within which we observe strong and varying whistler-mode waves. These ambient conditions lead to strong pitch-angle scattering, which contributes to the creation of a wide slot region as well as a significant decay of the outer radiation belts, which are observed and qualitatively well simulated. In addition, we find the substorm activity causes short duration (within ± 4h) decay of the plasma density and a lowering amplitude of the whistler-mode waves within the plasmasphere, both causing opposite effects in terms of pitch angle diffusion. This leads to a diminution of pitch-angle diffusion at the time of the main substorm activity. Conversely, whistler-mode waves become enhanced in the time periods between the substorm injections. All effects cumulated, we find an enhancement of pitch angle diffusion by whistler-mode waves above L~4.7 during the 10-day period. This directly relates to the combination of quietness and substorm activity which allows pitch angle diffusing of up to 1 MeV electrons in the outer belt. Relativistic electrons of 1–2 MeV remain trapped in the outer belt, from L~4.7 to L~5.2, forming, in both the observations and the simulations, a distinct pocket of remnant electrons.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2010

Simplified Vicarious Radiometric Calibration

A measurement-based radiance estimation approach for vicarious radiometric calibration of spaceborne multispectral remote sensing systems has been developed. This simplified process eliminates the use of radiative transfer codes and reduces the number of atmospheric assumptions required to perform sensor calibrations. Like prior approaches, the simplified method involves the collection of ground truth data coincident with the overpass of the remote sensing system being calibrated, but this approach differs from the prior techniques in both the nature of the data collected and the manner in which the data are processed. In traditional vicarious radiometric calibration, ground truth data are gathered using ground-viewing spectroradiometers and one or more sun photometer( s), among other instruments, located at a ground target area. The measured data from the ground-based instruments are used in radiative transfer models to estimate the top-of-atmosphere (TOA) target radiances at the time of satellite overpass. These TOA radiances are compared with the satellite sensor readings to radiometrically calibrate the sensor. Traditional vicarious radiometric calibration methods require that an atmospheric model be defined such that the ground-based observations of solar transmission and diffuse-to-global ratios are in close agreement with the radiative transfer code estimation of these parameters. This process is labor-intensive and complex, and can be prone to errors. The errors can be compounded because of approximations in the model and inaccurate assumptions about the radiative coupling between the atmosphere and the terrain. The errors can increase the uncertainty of the TOA radiance estimates used to perform the radiometric calibration. In comparison, the simplified approach does not use atmospheric radiative transfer models and involves fewer assumptions concerning the radiative transfer properties of the atmosphere. This new technique uses two neighboring uniform ground target areas having different reflectance values. The target areas can be natural or artificial and must be large enough to minimize adjacent-pixel contamination effects. The radiative coupling between the atmosphere and the terrain needs to be approximately the same for the two targets. This condition can be met for relatively uniform backgrounds when the distance between the targets is within a few hundred meters. For each target area, the radiance leaving the ground in the direction of the satellite is measured with a radiometrically calibrated spectroradiometer. Using the radiance measurements from the two targets, atmospheric adjacency and atmospheric scattering effects can be subtracted, thereby eliminating many assumptions about the atmosphere and the radiative interaction between the atmosphere and the terrain. In addition, the radiometrically calibrated spectroradiometer can be used with a known reflectance target to estimate atmospheric transmission and diffuse- to-global ratios without the need for ancillary sun photometers. Several comparisons between the simplified method and traditional techniques were found to agree within a few percent. Hence, the simplified method reduces the overall complexity of performing vicarious calibrations and can serve as a method for validating traditional radiative transfer models

Stanley, Thomas↗
NASA NTRS1996

Annual Research Briefs - 1996, Center for Turbulence Research

Partial contents include: (1) Simulation of jet impingement heat transfer with the k-e-V2 model; (2) A new methodology for turbulence modelers using Direct Numerical Simulation (DNS) database analysis; (3) A particle representation model for the deformation of homogeneous turbulence; (4) Modeling of inhomogeneous compressible turbulence using a two-scale statistical theory; (5) Lewis number effects on partially premixed flames; (6) Conditional moment closure of mixing and reaction in turbulent non-premixed combustion; (7) Dynamic models for LES of turbulent front propagation with a spectral method; (8) Jet noise using large eddy simulation; (9) Towards numerical simulations of trailing-edge aeroacoustics; (10) A model for the onset of vortex breakdown; (11) Compressing an elliptic vortex: transition to turbulence by tumble breakdown; (12) A new method for the adaptive control of vortex-wall interactions; (13) Symmetries in turbulent boundary layer flows; (14) Near-wall models in large eddy simulations of flow behind a backward-facing step; (15) Assessment of the subgrid-scale models at low and high Reynolds numbers; (16) Large-eddy simulations of flow around a NACA 4412 airfoil using unstructured grids; (17) Progress on LES of flow past a circular cylinder; (18) Experimental investigation of flow through an asymmetric plane diffuser; (19) Progress in the large-eddy simulation of an asymmetric plane diffuser; (20) A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation; (21) On the generation of turbulent inflow conditions for boundary layer simulations, and (22) A new class of finite difference schemes.

COMPUTATIONAL FLUID DYNAMICS↗
DOE OSTI2015

Sorption Modeling and Verification for Off-Gas Treatment

The project has made progress toward developing a comprehensive modeling capability for the capture of target species in off gas evolved during the reprocessing of nuclear fuel. The effort has integrated experimentation, model development, and computer code development for adsorption and absorption processes. For adsorption, a modeling library has been initiated to include (a) equilibrium models for uptake of off-gas components by adsorbents, (b) mass transfer models to describe mass transfer to a particle, diffusion through the pores of the particle and adsorption on the active sites of the particle, and (c) interconnection of these models to fixed bed adsorption modeling which includes advection through the bed. For single-component equilibria, a Generalized Statistical Thermodynamic Adsorption (GSTA) code was developed to represent experimental data from a broad range of isotherm types; this is equivalent to a Langmuir isotherm in the two-parameter case, and was demonstrated for Kr on INL-engineered sorbent HZ PAN, water sorption on molecular sieve A sorbent material (MS3A), and Kr and Xe capture on metal-organic framework (MOF) materials. The GSTA isotherm was extended to multicomponent systems through application of a modified spreading pressure surface activity model and generalized predictive adsorbed solution theory; the result is the capability to estimate multicomponent adsorption equilibria from single-component isotherms. This advance, which enhances the capability to simulate systems related to off-gas treatment, has been demonstrated for a range of real-gas systems in the literature and is ready for testing with data currently being collected for multicomponent systems of interest, including iodine and water on MS3A. A diffusion kinetic model for sorbent pellets involving pore and surface diffusion as well as external mass transfer has been established, and a methodology was developed for determining unknown diffusivity parameters from transient uptake data. Two parallel approaches have been explored for integrating the kernels described above into a mass-transport model for adsorption in fixed beds. In one, the GSTA isotherm kernel has been incorporated into the MOOSE framework; in the other approach, a focused finite-difference framework and PDE kernels have been developed. Issues, including oscillatory behavior in MOOSE solutions to advection-diffusion problems, and opportunities have been identified for each approach, and a path forward has been identified toward developing a stronger modeling platform. Experimental systems were established for collection of microscopic kinetics and equilibria data for single and multicomponent uptake of gaseous species on solid sorbents. The systems, which can operate at ambient temperature to 250°C and dew points from -69 to 17°C, are useful for collecting data needed for modeling performance of sorbents of interest. Experiments were conducted to determine applicable models and parameters for isotherms and mass transfer for water and/or iodine adsorption on MS3A. Validation experiments were also conducted for water adsorption on fixed beds of MS3A. For absorption, work involved modeling with supportive experimentation. A dynamic model was developed to simulate CO 2 absorption with chemical reaction using high alkaline content water solutions. A computer code was developed to implement the model based upon transient mass and energy balances. Experiments were conducted in a laboratory-scale column to determine model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. This project has resulted in 7 publications, with 3 manuscripts in preparation. Also, 15 presentations were given at national meetings of ANS and AIChE and at Material Recovery and Waste Forms Campaign Working Group meetings.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Global numerical modeling of energetic proton acceleration in a coronal mass ejection traveling through the solar corona

The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE AGORA HIGH-RESOLUTION GALAXY SIMULATIONS COMPARISON PROJECT. II. ISOLATED DISK TEST

Using an isolated Milky Way-mass galaxy simulation, we compare results from nine state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt–Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly formed stellar clump mass functions show more significant variation (difference by up to a factor of ∼3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low-density region, and between more diffusive and less diffusive schemes in the high-density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Viscous organic aerosol particles in the upper troposphere: diffusivity-controlled water uptake and ice nucleation?

Here, new measurements of water diffusion in secondary organic aerosol (SOA) material produced by oxidation of α-pinene and in a number of organic/inorganic model mixtures (3-methylbutane-1,2,3-tricarboxylic acid (3-MBTCA), levoglucosan, levoglucosan/NH 4 HSO 4 , raffinose) are presented. These indicate that water diffusion coefficients are determined by several properties of the aerosol substance and cannot be inferred from the glass transition temperature or bouncing properties. Our results suggest that water diffusion in SOA particles is faster than often assumed and imposes no significant kinetic limitation on water uptake and release at temperatures above 220 K. The fast diffusion of water suggests that heterogeneous ice nucleation on a glassy core is very unlikely in these systems. At temperatures below 220 K, model simulations of SOA particles suggest that heterogeneous ice nucleation may occur in the immersion mode on glassy cores which remain embedded in a liquid shell when experiencing fast updraft velocities. The particles absorb significant quantities of water during these updrafts which plasticize their outer layers such that these layers equilibrate readily with the gas phase humidity before the homogeneous ice nucleation threshold is reached. Glass formation is thus unlikely to restrict homogeneous ice nucleation. Only under most extreme conditions near the very high tropical tropopause may the homogeneous ice nucleation rate coefficient be reduced as a consequence of slow condensed-phase water diffusion. Since the differences between the behavior limited or non limited by diffusion are small even at the very high tropical tropopause, condensed-phase water diffusivity is unlikely to have significant consequences on the direct climatic effects of SOA particles under tropospheric conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Denoising diffusion probabilistic models for generative alloy design

Inverse material design is an extremely challenging optimization task made difficult by, in part, the highly nonlinear relationship linking performance with composition. Quantitative approaches have improved significantly owing to advances in high throughput experimentation and computational thermodynamics. However, existing physics-based tools are mostly forward models; input a chemistry and obtain a prediction. More recently the materials community has leveraged advances in the machine learning community to establish novel inverse design frameworks. Very recently denoising diffusion probabilistic models have been shown to be extremely powerful generators producing synthetic data of various modalities e.g. images, text, audio, tables, etc.. In this work a novel framework for alloy design and optimization is proposed leveraging these class of models. Five key generative tasks are demonstrated (1) unconditional generation (2) composition conditioned generation (3) property conditioned generation (4) multi-feedstock conditioned generation and (5) generative optimization. These methods were tested on three case studies: high entropy alloy design, superalloy binder jet additive manufacturing, and in-situ dual-feedstock wire-arc additive manufacturing. Results indicate that the established models are extremely flexible, expressive, and robust. The architecture’s flexibility and training procedure empower the model to learn complex intra-compositional and composition-property relationships. Furthermore, the probabilistic nature of these models makes them well suited for addressing solution non-uniqueness and tackling uncertainty quantification tasks. While the fidelity and quantity of the underlying training data is paramount, we envision that future alloy design frameworks will make extensive use of these kinds of machine learning models as “search” tools bolstering the utility of experimental and computational approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Decoupling thermal and irradiation effects on grain boundary segregation

Radiation-induced segregation (RIS) is most often measured by peak solute concentration at a boundary. However, this may give an incomplete picture of segregation quantity and phenomena. Radiation-induced and thermal segregation at grain boundaries was investigated in Fe-9.6 at.% Cr after 9 MeV Fe 3+ ion irradiation at 400 °C. The experimental results were compared to kinetic Monte Carlo (KMC) simulations. The study revealed that Cr enrichment (peak segregation) at the grain boundaries was comparable in both the irradiated and non-irradiated conditions, although irradiation resulted in broader segregation profiles in both experiments and simulations, indicating an overall increase in grain boundary segregation due to irradiation. This broadening is attributed to back diffusion into the grain interior. While it is an established phenomenon, this study offers a quantitative evaluation using experimental data and KMC modeling. Further, these results emphasize the importance of analyzing the entire segregation profile and decoupling the thermal and irradiation contributions to solute segregation.

Grain boundary segregation↗
DOE OSTI2026

Digital quantum magnetism on a trapped-ion quantum computer

Digital quantum matter—realized when discrete quantum gates approximate continuous time evolution—is susceptible to heating into chaotic, structureless states. If digitization errors are adequately suppressed, a long-lived transient regime of approximately energy-conserving dynamics can be observed on gate-based quantum computers. Conservation of energy, in turn, enables the exploration of a wide variety of complex behaviours observed in equilibrium systems, ranging from the non-trivial microscopic origins of thermalization itself to the stabilization of effective models hosting exotic emergent properties. Here we use Quantinuum’s H2 quantum computer to simulate digitized dynamics of the quantum Ising model, suppressing digitization errors well enough to observe thermalization on timescales that severely challenge classical simulation methods. Relaxation of an inhomogeneous state reveals an emergent hydrodynamics owing to approximate energy conservation and we compute the associated diffusion constant. By reprogramming our simulations to take place on a triangular lattice with periodic boundary conditions, we observe thermalization consistent with emergent gauge and topological constraints resulting from lattice frustration. Furthermore, our results were enabled by continued advances in two-qubit gate quality (native partial entangler fidelities of 99.94(1)%) and establish digital quantum computers as powerful tools for studying (effectively) continuous-time dynamics.

Information theory and computation↗