Search NASA⌕ Search

SEARCH · Search NASA

Results for “self similarity”

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 235 records · Page 13

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗

Aircraft observations in a tropical supercluster over the equatorial Indian Ocean during MISO-BOB field campaign

The Monsoon Intra-Seasonal Oscillations in the Bay of Bengal (MISO-BOB) field campaign was conducted in the Indian Ocean during the 2018 and 2019 summer monsoon seasons. WC-130J aircraft of the 53rd Weather Reconnaissance Squadron of the US Air Force participated in the campaign in June 2018. The dropsonde observations across a tropical supercluster showed zonal wind variations in association with the structure of the convectively coupled Kelvin wave (CCKW). Within the supercluster, easterlies (westerlies) were observed in the upper (lower) troposphere; this transformation occurred just below the 0° C level. The cold pool had an easterly component throughout, and it was coldest (by 2.5° C) at the center of the supercluster, deepest (~ 1000 m) at its rear/western end, and shallowest (~ 300 m) at the front/eastern end. The level of free convection (LFC) at the front end was at 897 m altitude. At the eastern flank of the supercluster, zonal convergence in the lower troposphere occurred between 500-1500 m levels above the surface between the westerlies within the supercluster and opposing ambient easterlies. Thus, the uplifting of conditionally unstable air parcels above LFC to the east of the supercluster was likely to occur due to this convergence rather than the cold pool influence. Conversely, the western flank of the supercluster had low-level zonal divergence. These observations support the notion of ‘self-similarity’ among the mesoscale convective systems and large-scale waves.

54 ENVIRONMENTAL SCIENCES↗

Interplay of intercalation dynamics and lithium plating in monolithic and architectured graphite anodes during fast charging

Fast charging of high-capacity anodes is challenging due to lithium plating reactions, which lead to poor cycling performance and safety concerns. Thus, accurate predictions of plating onset and an understanding of this electrochemical process are crucial for robust battery design. However, the most commonly used models, based on porous electrode theory (e.g., the pseudo-2D model), are notoriously difficult to calibrate due to their complexity, limiting their predictive power. This work studies the process of lithium plating during fast charging of (small-particle) graphite half-cells by measuring local reaction progression and plating behavior using optical operando techniques. These experiments employ a realistic 1D graphite electrode geometry with commercially-relevant mass loading charged at fast charge rates. It is demonstrated that the local reaction progression and plating onset can not only be predicted accurately with a p2D numerical model, but that these processes follow a simple scaling law. Remarkably, the entire reaction histories of different electrodes charged at different rates (e.g., 160 μm thickness at 0.5C, 111 μm at 1C or 66 μm at 4C) were observed to have self-similar intercalation profiles. It is demonstrated that plating onset is in turn governed by the reaction profile which explains why both processes exhibit the same scaling behavior. Finally, operando measurements of local reaction dynamics are conducted for the first time in electrodes with channeled architectures, quantitatively determining how channels affect reaction uniformity and plating onset. Together, these results reveal underlying simplicity in the complex electrochemical environment of fast charging and lithium plating, improving understanding of this process. These fundamental insights are broadly applicable for design processes, modeling and experimental evaluation of lithium ion batteries.

25 ENERGY STORAGE↗

The SRG/eROSITA All-Sky Survey: Dark Energy Survey year 3 weak gravitational lensing by eRASS1 selected galaxy clusters

Context. Number counts of galaxy clusters across redshift are a powerful cosmological probe if a precise and accurate reconstruction of the underlying mass distribution is performed – a challenge called mass calibration. With the advent of wide and deep photometric surveys, weak gravitational lensing (WL) by clusters has become the method of choice for this measurement. Aims. We measured and validated the WL signature in the shape of galaxies observed in the first three years of the Dark Energy Survey (DES Y3) caused by galaxy clusters and groups selected in the first all-sky survey performed by SRG (Spectrum Roentgen Gamma)/eROSITA (eRASS1). These data were then used to determine the scaling between the X-ray photon count rate of the clusters and their halo mass and redshift. Methods. We empirically determined the degree of cluster member contamination in our background source sample. The individual cluster shear profiles were then analyzed with a Bayesian population model that self-consistently accounts for the lens sample selection and contamination and includes marginalization over a host of instrumental and astrophysical systematics. To quantify the accuracy of the mass extraction of that model, we performed mass measurements on mock cluster catalogs with realistic synthetic shear profiles. This allowed us to establish that hydrodynamical modeling uncertainties at low lens redshifts (z < 0.6) are the dominant systematic limitation. At high lens redshift, the uncertainties of the sources’ photometric redshift calibration dominate. Results. With regard to the X-ray count rate to halo mass relation, we determined its amplitude, its mass trend, the redshift evolution of the mass trend, the deviation from self-similar redshift evolution, and the intrinsic scatter around this relation. Conclusions. The mass calibration analysis performed here sets the stage for a joint analysis with the number counts of eRASS1 clusters to constrain a host of cosmological parameters. We demonstrate that WL mass calibration of galaxy clusters can be performed successfully with source galaxies whose calibration was performed primarily for cosmic shear experiments, opening the way for the cluster cosmological exploitation of future optical and NIR surveys like Euclid and LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

A conservative approach to scaling magneto-inertial fusion concepts to larger pulsed-power drivers

The Magnetized Liner Inertial Fusion (MagLIF) experimental platform [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)] represents the most successful demonstration of magneto-inertial fusion (MIF) techniques to date in pursuit of ignition and significant fusion yields. The pressing question remains regarding how to scale MIF concepts like MagLIF to more powerful pulsed-power drivers while avoiding significant changes in physical regimes that could adversely impact performance. In this work, we propose a conservative approach for scaling general MIF implosions, including MagLIF. Underpinning our scaling approach is a theoretical framework describing the evolution of the trajectory and thickness of a thin-walled, cylindrical, current-driven shell imploding on preheated, adiabatic fuel. By imposing that scaled implosions remain self-similar, we obtain a set of scaling rules expressing key target design parameters and performance metrics as functions of the maximum driver current I max . We identify several scaling paths offering unique, complementary benefits and trade-offs in terms of physics risks and driver requirements. Remarkably, when scaling present-day experiments to higher coupled energies, these paths are predicted to preserve or reduce the majority of known performance-degrading effects, including hydrodynamic instabilities, impurity mix, fuel energy losses, and laser-plasma interactions, with notable exceptions clearly delineated. In the absence of α heating, our scaling paths exhibit neutron yield per-unit-length scaling as $\tilde{Y}$∝[I$3\atop{max}$,I${4.14}\atop{max}$] and ignition parameter scaling as χ∝[I max ,I${2.14}\atop{max}$]. By considering the specific physics risks unique to each scaling path, we provide a roadmap for future investigations to evaluate different scaling options through detailed numerical studies and scaling-focused experiments on present-day facilities. Overall, these results highlight the potential of MIF as a key component of the national ignition effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗

A modal wave-packet model for the multi-mode Richtmyer–Meshkov instability

In this work, a model for multimode perturbations subject to the Richtmyer–Meshkov (RM) instability is presented and compared with simulations and experiments for conditions relevant to inertial confinement fusion. The model utilizes the single mode response to the RM impulse whereby its amplitude h(k, t) first grows with an initial velocity V 0 ∝ kh(k, 0) that eventually decays in time as 1/kV 0 t. Both the growth and saturation stages are subject to nonlinearities since they depend explicitly on the initial amplitude. However, rather than using the individual mode amplitude h(k, t), nonlinearity is taken to occur when the root-mean-square amplitude h rms (k, t) of a wave-packet within wavenumbers k ± δk becomes comparable to 1/k. This is done because nearby sidebands can act in unison for an auto-correlation distance 1/δk beyond nonlinearity as observed in the beam-plasma instability. Thus, the nonlinear saturation amplitude for each mode is reduced from the usual 1/k by a phase space factor that depends on the physical dimensionality, as in the Haan model for the Rayleigh–Taylor instability. In addition, for RM, the average value of kh rms for the initial spectrum is used to calculate a nonlinear factor F NL that reduces V 0 , as observed for single modes. For broadband perturbations, the model describes self-similar growth ∝t θ as successively longer wavelength modes reach saturation. The growing and saturated modes must be discerned because only the former promote θ and are enhanced by reshock and spherical convergence. All of these flows are described here by the model in good agreement with simulations and experiments.

42 ENGINEERING↗

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗

Simulating a pulsed-power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the magnetohydrodynamics (MHD) domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This much more affordable ideal MHD representation is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scaling patch analysis of planar turbulent mixing layers

Proper scales for the mean flow and Reynolds shear stress in planar turbulent mixing layers are determined from a scaling patch analysis of the mean continuity and momentum equations. By seeking an admissible scaling of the mean continuity equation, a proper scale for the mean transverse flow is determined as V ref =(dδ/dx)U ref , where dδ/dx is the growth rate of the mixing layer width and U ref =U h –U l is the difference between the velocity of the high speed stream U h and the velocity of the low speed stream U l . By seeking an admissible scaling for the mean momentum equation, a proper scale for the kinematic Reynolds shear stress is determined as R uv,ref =U avg V ref =[$\frac{1}{2A_u}$$\frac{dδ}{dx}$]$U^{2}_{ref}$ where A u $_{=}^{def}$(U h –U l )/(U h +U l ) is the normalized velocity difference that emerges naturally in the admissible scaling of the mean momentum equation. Self-similar equations for the scaled mean transverse flow V* and Reynolds shear stress $R^{*}_{uv}$=R uv /R uv,ref are derived from the mean continuity and mean momentum equations. Finally, approximate equations for V* and $R^{*}_{uv}$ are developed and found to agree well with experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental and theoretical comparison of ion properties from nanosecond laser-produced plasmas of metal targets

The ion emission properties of laser-produced plasmas as a function of laser intensities between 4–50 GW cm –2 and varying angles with respect to the target normal were investigated. The plasmas were produced by focusing 1064 nm, 6 ns pulses from an Nd:YAG laser on various metal targets. The targets used for this study include Ti, Mo, and Gd (Z = 22, 42, 64). It is noted that all ion profiles are composed of multiple peaks—a prompt emission peak trailed by three ion peaks (ultrafast, fast, and thermal). Experimentally, it is shown that each of these ion peaks follows a unique trend as a function of laser intensity, angle, and distance away from the target. Theoretically, it is shown that simple analytical models can be used to explain the properties of the ions. The variations in the ion velocity and density as a function of laser intensity are found to be in good agreement with theoretical models of sheath acceleration, isothermal self-similar expansion, and ablative plasma flow for various ion peaks.

42 ENGINEERING↗

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Focusing of laser-accelerated ions with spherically curved targets

A theoretical model for laser-accelerated ion focusing from spherically curved targets using a Gaussian self-similar solution is presented. This model describes the evolution of the focal location and focal radius with final ion energy and the initial radial boundary of the ions, the latter being a function of both target geometry and the ion acceleration radius for a given ion energy. The theory is supported by particle-in-cell simulations of a variety of target shapes with varying radii of curvature and target opening angles, as well as variations in the injected electron beam radius and energy spectrum. The theory and simulations suggest that the focal location varies linearly with the radius of curvature, with the square root of the ratio of ion energy to effective electron temperature, and monotonically increases with the energy-dependent ratio of the initial ion radial boundary to the radius of curvature. Considering ponderomotive acceleration, this suggests that the focal length should scale inversely with intensity, $(I_Lλ_L^2)$ −1/4 , suggesting a 10× increase in intensity will reduce the focal distance by ̃1.8× for a given ion energy.

Physics - Plasma physics↗

Out-of-time-order correlators bridge classical transport and quantum dynamics

The out-of-time-order correlator (OTOC) has emerged as a central tool for quantifying decoherence across wide-ranging physical platforms. Here, we demonstrate its direct measurement in a classical ensemble using nuclear magnetic resonance with a modulated gradient spin echo sequence and extend the method into a multidimensional correlation to track exchange phenomena. Position is encoded through magnetic field gradients and momentum through the velocity autocorrelation function, enabling experimental access to OTOCs for proton motion confined within the self-similar lattice of the metal–organic framework MOF-808. Here, water confined to specified geometries within the MOF pores gives rise to spatially distinct diffusive eigenmodes with characteristic relative entropies. We demonstrate that periodic radio frequency driving combined with gradient modulation yields entropy evolution through the selection of distinct diffusion modes. Frequency-resolved diffusion spectra connect these entropy dynamics to classical heat exchange laws, revealing how operational features of quantum systems are mirrored in confined, macroscopic spin ensembles.

Fricke, Sophia N. [University of California, Berke↗