## Reconstructing Richtmyer-Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks

DANIEL A. SERINO, 1,* MARC L. KLASKY, 1 BALASUBRAMANYA T. NADIGA, 1 XIAOJIAN XU, 2 AND TREVOR WILCOX 1

1 Los Alamos National Laboratory, Los Alamos, NM 87545, USA

- 2 Dept. of Electrical and Computer Engineering, University of Michigan, Ann Arbor, MI 48109
- * dserino@lanl.gov

Abstract: A trained attention-based transformer network can robustly recover the complex topologies given by the Richtmyer-Meshkoff instability from a sequence of hydrodynamic features derived from radiographic images corrupted with blur, scatter, and noise. This approach is demonstrated on ICF-like double shell hydrodynamic simulations. The key component of this network is a transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model. This approach is demonstrated to exhibit an excellent ability to accurately recover the RichtmyerMeshkov instability growth rates, even despite the gas-metal interface being greatly obscured by radiographic noise.

## © 2024 Optica Publishing Group

## 1. Introduction

The ability to isolate precise material interfaces from radiographic images is of paramount importance in describing a broad array of physical phenomena including a number of hydrodynamic instabilities including the Richtmyer Meshkoff instability (RMI) [1,2], the Rayleigh Taylor instability [3-5], and the Kelvin-Helmholtz instability [6]. These phenomena arise in a broad range of situations from astronomical size events like supernova collapse [7], to microscopic events such as gas bubble sonoluminescence [8] and the supersonic combustion of ramjets [9]. Indeed, extraction of the peaks and troughs associated with these phenomena at the material interfaces, RMI in particular, is essential in capturing the growth rates and optimizing designs to mitigate these instabilities. Here we focus on the RMI that originates from the interaction of a shock wave with an interface separating two materials with different densities. In most cases an initial perturbation will be amplified following the passage of the shock. A deposition of baroclinic vorticity that is practically-instantaneous and results from a misalignment between the pressure gradient across the shock and the local density gradient across the interface leads to a subsequent (slower) growth of the perturbation.

While extensive examinations and reviews of the RMI have been performed over the past sixty years, most of of the experimental work to validate the RMI has been performed within shock tubes using rectangular geometry with gaseous or fluid materials, although some recent work has been performed in cylindrical geometry [10-16]. Additionally, some recent works have also examined the RMI in solid-solid settings to estimate the strength in ductile materials [17, 18].

The emergence of Inertial Confinement Fusion (ICF) as a potential power source has been a major impetus for the continued examination of the RMI [11, 12]. However, the spherical convergent geometry as well as large attenuation present in these experiments precludes previous analysis methods and diagnostics. For example, planar laser induced fluorescence, planar Rayleigh scattering, and particle image velocimetry. Consequently, examination of the RMI

Fig. 1. Example double shell capsule specification based on the 1.06 MJ yield design from Ref. [19].

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-2bcc9cb5-eee0-4de2-8578-45544dd6e2b3.jpg)

Fig. 2. Sample ( 𝑟, 𝑧 ) projection of the density (1 st column), zoomed-in view of the Richtmyer-Meshkov interface (2 nd column), synthetic radiograph (3 rd columns), and a zoomed-in view of the radiograph (4 th columns) labeled with the RMI interface (left half) and Canny edge labels (right half).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-a386c793-41f3-4b41-9812-e8caf315eadb.jpg)

instability in promising ICF concepts, such as the double shell capsules, depicted in Figure 1 demand the development of new experimental and analysis techniques to capture details of the RMI whose behavior and control are thought to be crucial for continued progress in ICF [11,20]. That is, the RMI may induce mixing due on the interface between the outer shell and the inner fuel and in so doing inhibit the fusion reaction and can consequently be a limiting factor in the energy produced [21-24].

As an illustration in the difficulty in locating the interfaces necessary to quantitatively determine the RMI growth rates in spherical geometries, we present the result of a ICF hydrodynamic simulation in which a shock wave generated during the implosion of a nearly-spherical double shell capsule into a gas medium interacts with the irregular surface at the gas-metal interface during rebound to generate an RMI. It should be noted that this configuration is of particular interest in the proposed double shell ICF simulations as illustrated in Figure 1.

While the material interface in the hydrodynamic simulation is clearly observable and easily captured by applying traditional image processing methods, the observation from the noisy radiographic image makes the identification of the peaks and troughs (which are necessary to quantitatively characterize the growth rates of the RMI) problematic. Indeed, as may be observed from Figure 2, traditional edge detection algorithms are generally insufficient to accurately find these features. Furthermore, these methods (Sobel, Canny, and Laplacian of Gaussian edge filters) generally require fine-tuning of hyper-parameters associated with the method, which may lead to significant uncertainty/error in the results. Finally, the sensitivity of the RMI to the initial characteristics of the perturbed surface further complicate the quantitative determination of growth rates.

This paper reports on the development of a new physics-based density reconstruction method capable of capturing, through a subsequent interface finding algorithm, RMI growth rates of ICF double shell implosions in a spherically convergent geometry using low dimensional features that are robustly identifiable from dynamic radiographic images. This approach may enable experimental validation in a regime that, to our knowledge, been largely unexplored.

## 2. Density Reconstruction via Dynamic Radiography

To examine the dynamic evolution of materials undergoing strong deformation in material science, shock physics, and ICF, dynamic radiography is routinely used as an experimental diagnostic [25-28]. Indeed, capturing the dynamic evolution of the material interface is essential in quantifying the growth rates of the RMI. As previously discussed, the limitations of traditional image processing algorithms require consideration of alternative techniques for capturing this interface with sufficient accuracy to infer the growth rates from noisy radiographs. One such technique for capturing the material interface is to use traditional density reconstruction techniques followed by applying gradient-based approaches to locate the material interfaces associated with the RMI growth rates. That is, an inversion of the projection data may first be performed to obtain the density field. Traditional methods of performing density inversions date back to Radon [29]. Over the past several decades, image reconstruction methods have evolved from simple analytical methods such filtered-back projection (FBP) methods, X-ray CT (e.g. Feldkamp-Davis-Kress or FDK methods), and the Inverse Abel Transform for axisymmetric systems [30-32] to models that enable inclusion of poly-energetic sources and also attempt to address the noise field (e.g. statistical noise as well as scatter). The models are complex nonlinear, non-convex forward models and employ iterative reconstruction techniques [33]. Iterative reconstruction (IR) algorithms are based on more sophisticated models for the imaging system's physics and models for sensor and noise statistics. These methods are often called model based image reconstruction (MBIR) methods or statistical image reconstruction (SIR) methods [34]. Generally, density reconstruction algorithms using dynamic experimental radiographic data with complex noise fields (e.g., non-Gaussian noise, scattered radiation, complex beam dynamics, etc.) still encounter difficulty when sufficient accuracy is required to extract the intricate details necessary to characterize the RMI phenomena [28]. Furthermore, the nature of these methods introduces a certain degree of regularization that may impact the ability to accurately preserve edges. Research continues in this area to improve the ability to preserve edges in settings where uncertainty in the characteristics of scatter and noise are present [35-37].

More recently, machine learning (ML) approaches have been applied to radiographic reconstruction. Many of these ML architectures have outperformed IR methods by a large margin at a specified degradation level [38, 39]. Similar to traditional non-ML-based algorithms, ML-based approaches experience difficulty with complex noise fields. ML architectures that are trained using a specific noise distribution typically fail to obtain favorable results in the presence of out-of-sample noise [40,41]. This problem can be partially addressed by training a number of models targeting each level of degradation. However, in many cases, the noise depends on the object and other imaging artifacts and therefore it is impractical to produce a general noise model to cover all unknown eventualities. In recent work [41], a machine learning-based Wasserstein generative adversarial network (WGAN) was developed to reconstruct density fields from noisy radiographic projections. This approach demonstrated excellent performance in reconstructing density for testing cases that had the same noise properties as the training set. However, when testing noise levels outside of the training population, the reconstructions exhibited a rapid degradation in performance. Interestingly, degradation in performance was observed even when the noise level was reduced from the level used in the training.

In light of the inability to train a general denoising algorithm to address the issues associated with out-of-population noise and other imaging artifacts, we instead use a class of new ML-based methods that rely on features that are robustly identifiable from a series of dynamic radiographs. These methods are composed of two major components; one component is responsible for feature extraction, in this case the robust out-going shock, and another component performs density reconstructions from the extracted shock. This approach has been previously demonstrated on a relatively simple 1D system, where robust features were combined with the underlying hydrodynamic equations of motion to produce density reconstructions [28]. This methodology outperformed a traditional, direct radiograph-to-density reconstruction method in the presence of scatter. Furthermore, the method was also capable of generating families of solutions consistent with the observed features and a methodology for examining the uncertainty in the predictions was proposed.

To investigate the RMI we consider a test problem of shock propagation due to the implosion of a gas-filled ICF double shell capsule with a sinusoidal perturbation on the interior surface. Numerous simulations were performed with variable initial conditions and material parameters to produce a data set of density field sequences in time. Synthetic radiographs were produced for each time snapshot of density and a noise field was generated to mimic a realistic experimental setup. Ground truth features were generated using the density fields corresponding to coefficients of a curve fit for the shocks. We focused on a features-to-density network, for which we investigated two approaches including a generative variational autoencoder (VAE) approach called the ShockDecoderViT based on the vision transformer [42], and a deterministic structurepreserving model called the Mass-Conserving Transformer, based on the original transformer [43] to reconstruct the density fields to enable the subsequent extraction of the gas/metal interface and the RMI growth rates. Each architecture transforms a temporal sequence of features into their corresponding density fields using networks that incorporate attention layers. In previous work, attention blocks increased the expressive power of the networks by correlating temporal dependencies in the sequences [43]. In this paper, we demonstrate that incorporating attention in the features-to-density networks improves density reconstruction errors relative to a similar network without attention. The density fields are then processed to obtain the gas/metal interface to obtain the RMI growth rates.

To test the effectiveness of the features-to-density network in the full radiograph-to-density pipeline, we developed a noise model to characterize the errors of a feature extractor network. For the noise model, we used the results from a feature extractor network developed in [44] that consists of a convolutional neural network (CNN) with an image Fourier feature encoding (IFFE) layer. This network was trained on the noisy radiograph-to-features problem and is capable of recovering the shock features with sub-pixel accuracy. We trained the features-to-density network using noiseless features and applied the noise model during testing. We show that the features-todensity network remains robust in the presence of significant amounts of out-of-sample noise. Using the results of these ML methods, we investigated the limits of the radiograph-to-density pipeline in reconstructing the much higher dimensional phenomena present. Using the density fields we then assessed the ability of the ML method to reconstruct the gas metal interface attributed to the RMI with lower dimensional features, such as the outgoing shock captured from the radiographs. This pipeline provides a possible means to address the difficulty in validating simulations of 3D ICF double shell RMI experiments.

In the remainder of this work we present details of our new approach for performing density reconstructions and subsequent extraction of the gas-metal interface necessary for characterizing

RMI phenomena using the robust features extracted from the dynamic radiographic images. Section 3 introduces the dynamic radiography problem and model formulation. Section 4 introduces details of the hydrodynamic ICF test problem and data generation. Section 5 presents two machine learning approaches for determining the density fields from robust features identified from a dynamic radiographic sequence. Section 6 presents the results of the density reconstructions and ability to extract RMI growth rates using the ML architectures. Finally section 7 presents discussion of the simulation results and conclusions drawn from the numerical experiments.

## 3. Dynamic Radiography Problem Setup

Our objective is to recover the intricate details of the gas-metal interface by first estimating a time-series of densities GLYPH&lt;8&gt; 𝜌 𝑡,𝑖 ( 𝑥 ) GLYPH&lt;9&gt; 𝑡 𝑁 𝑡 = 𝑡 1 , with 𝑥 ∈ R 3 and 𝑖 representing a material index, from their corrupted radiographs { 𝑚 𝑡 } 𝑡 𝑁 𝑡 = 𝑡 1 , where 𝑚 𝑡 ∈ R 𝑁 1 × 𝑁 2 is a monochrome image with 𝑁 1 × 𝑁 2 pixels. The times 𝑡 1 , 𝑡 2 , ..., 𝑡 𝑁 denote 𝑁 time points at which the radiographic measurements are collected. In this investigation we consider the axisymmetric problem, where density objects are fully characterized by the central slice through the near-spherical object. Following the determination of the density fields we then extract the gas-metal interface as described in Section 6.

Our imaging model may be described as follows: The areal density of the object along a ray 𝑟 connecting the source and detector is denoted

$$\rho _ { A _ { i } } ( r ) = \mathcal { A } [ \rho _ { i } ] = \int _ { - \infty } ^ { \infty } \rho _ { i } ( r _ { x } ( u ), r _ { y } ( u ), r _ { z } ( u ) ) \, \mathrm d u, & & ( 1 )$$

where A is the forward operator corresponding to the Abel transform, ( 𝑟 𝑥 ( 𝑢 , 𝑟 ) 𝑦 ( 𝑢 , 𝑟 ) 𝑧 ( 𝑢 )) is a parameterization of the ray 𝑟 , 𝜌 𝑖 (·) denotes a time snapshot of a spatially varying density, and 𝑖 is an index corresponding to the material being imaged. Using a simple measurement model of a mono-energetic X-ray source, the number density of photons reaching the detector along ray 𝑟 is approximately given by

$$I ( r ) = I _ { 0 } \exp \left ( - \sum _ { i } \xi _ { i } \rho _ { A _ { i } } ( r ) \right ),$$

where 𝐼 0 is the number density of the incident beam and 𝜉 𝑖 is the mass attenuation coefficient of material 𝑖 [45]. A discrete, finite radiograph is measured in practice representing a finite grid of detectors. This direct radiograph (without scatter) 𝑑 at each detector pixel is approximated as

$$d [ m, n ] = \int _ { R _ { m, n } } I ( r ) \, \mathrm d r \approx C I ( r _ { m, n } )$$

where 𝑅 𝑚,𝑛 denotes the rays impinging pixel ( 𝑚, 𝑛 ) , 𝑟 𝑚,𝑛 is the ray through the pixel center, and 𝐶 is a constant that depends on factors such as the detector pixel area.

The transmission (noisy radiograph) includes contamination from several noise terms, and is given by

$$T _ { t } & = d _ { t } + n _ { t },$$

where 𝑑 𝑡 and 𝑛 𝑡 are the direct radiograph and noise at time 𝑡 . At each time, we model the noise as

$$n = D _ { \text{dsh} } + D _ { s } + B _ { s } + \eta$$

where 𝐷 dsb represents the blur from the source and detector, given by

$$D _ { d s b } = D _ { s b } * \phi _ { \text{db} }, \quad D _ { s b } = d * G _ { \text{blur} } ( \sigma _ { \text{blur} } ).$$

Here, 𝜙 db is a custom kernel associated with the detector and 𝐺 blur ( 𝜎 blur ) is a 2D Gaussian kernel with standard deviation 𝜎 blur . In addition, the noise model includes correlated scattered radiation, 𝐷 s , and an uncorrelated linear tilted background scatter field, 𝐵 s , given by

$$D _ { s } = \kappa d * G _ { \text{scatter} } ( \sigma _ { \text{scatter} } ), \quad B _ { s } = a x + b y,$$

where 𝐺 scatter ( 𝜎 blur ) is a 2D Gaussian kernel with standard deviation 𝜎 scatter , 𝑎 and 𝑏 are scalar constants, and 𝑥 and 𝑦 are projections of the 2D functions 𝑥 and 𝑦 to the image plane. Finally, both gamma and photon noise are added as

$$\eta = \eta ^ { \text{Po} } ( \gamma _ { \text{g} } ) * \phi _ { \text{g} } + \eta ^ { \text{Po} } ( \gamma _ { \text{p} } ) * \phi _ { \text{p} }.$$

where 𝜂 Po is Poisson Noise and 𝛾 g and 𝛾 p are the rates and 𝜙 g and 𝜙 p are the gamma and photon kernels, respectively.

The choice to model scatter as a kernel convolved with the direct signal is common in the scatter correction literature [46, 47]. This approach provides a fast scatter model that is representative of models used in practice. The above model can be readily extended to polyenergetic X-ray sources [47,48]. This work focuses on material objects that are of interest in double shell ICF capsules, and we work with monoenergetic X-ray sources for simplicity.

Great difficulty in performing the density time-series reconstruction arises due to the presence of scatter, noise, and additional factors including the variability of the beam spot, energy spectra, and model mismatch in the forward operator. Incorporating a time-series for the evolution of the underlying density fields over time could provide improved dynamic reconstructions [28,49-52]. To this end we adopt a machine learning approach to incorporate hydrodynamic priors to encourage the estimated density time-series to be consistent with the equations of motion that govern their evolution.

The evolution of the density over time in a dynamic experiment can be modeled by a system of partial differential equations (PDEs) describing radiation hydrodynamics. To facilitate the analysis of the radiation-hydrodynamic system we utilize the Euler equations in a manner analogous to Bello-Maldonado. [53,54] The system of PDEs comprise a continuum model governing density evolution and is applicable in many realistic scenarios where dynamic radiography is applied. A closure of the hydrodynamic variables in the form of a material-dependent equation of state (EOS) is required, along with appropriate initial and boundary conditions, to uniquely prescribe the time evolution of any hydrodynamic variable, including the density. In practice, the EOS and other conditions or parameters for a specific experiment or test case are usually unknown. In this case, directly using the PDEs to enable dynamic reconstruction is not feasible. In our study, we numerically integrate the Euler equations with varying choices of EOS model parameters to generate a large dataset of density time-series and learn an EOS-agnostic model for reconstructing densities.

## 4. Hydrodynamic Test Problem

## 4.1. Generation of Density Time Series

As a test problem, we study shock propagation in a time-dependent (3D) density profile, created by an implosion of a nearly-spherical ICF double shell configuration. The study is limited to azimuthal symmetry so that the density at any time can be described in cylindrical coordinates ( 𝑟, 𝑧 ) . Additionally we restrict our attention to the Mie-Grüneisen (MG) EOS model. Simulations are performed on a 440 × 440 uniform Cartesian grid on a computational domain given by the quarter-plane [ 0 , 𝐿 ] × [ 0 , 𝐿 ] , where 𝐿 = 341 𝜇 m. The uniform grid cell size is Δ 𝑟 = Δ 𝑧 = 440 𝐿 . The metallic shell is made of Tantalum and its density is initially uniform at a value of 16.65 g/cc. The inner surface of the Tantalum can be described as the set of coordinates

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-1888be96-013f-4365-98c2-4575c0b8ca87.jpg)

Fig. 3. Left: 3D mock-up of a Tantalum shell (green) with a perturbation on the interior surface and an outer ablator layer (gray). Middle: projection of the Tantalum shell onto ( 𝑟, 𝑧 ) coordinates. The inner radius is parameterized by the angle, 𝑢 , between the white dotted line and the 𝑟 axis. The drive from the ablator is modelled as an initial velocity on the Tantalum shell ( 𝑣 impl ). Right: Plot of the 20 separate profiles for radius of the perturbed inner surface verses angle 𝑢 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-3276a8a4-96e0-4808-9408-284964ba7029.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-97f383e4-0497-432c-83b0-d8e7242c1e8a.jpg)

( 𝑟 in ( 𝑢 , 𝑧 ) in ( 𝑢 )) satisfying

$$( r _ { \text{in} } ( u ) ^ { 2 } + z _ { \text{in} } ( u ) ^ { 2 } ) ^ { 1 / 2 } = R _ { \text{in} } + \sum _ { k = 1 } ^ { 8 } F _ { k } \cos ( 2 k u ), \quad u \in [ 0, \pi / 2 ],$$

where 𝑅 in = 248 𝜇 m, 𝐹 𝑘 , 𝑘 = 1 , . . . , 8, are coefficients of the perturbation corresponding to the 𝑘 th cosine harmonic. The outer surface of the shell is a sphere with radius 𝑅 out = 310 𝜇 m. There are 20 different inner surface perturbation profiles considered in our dataset. The corresponding coefficients are recorded in Table 2. Figure 3 presents an initial perturbation given to the interior shell. As an initial condition, the shell is given a uniform implosion velocity, 𝑣 impl, in the direction of the origin to initiate an implosion.

In addition to varying the profile of the inner radius and the initial implosion velocity, our dataset consists of simulations covering parameters characterizing the MG equation of state [55],

$$p \left ( \chi, T \right ) = \frac { \rho _ { 0 } c _ { s } ^ { 2 } \chi \left ( 1 - \frac { 1 } { 2 } \Gamma _ { 0 } \chi \right ) } { ( 1 - s _ { 1 } \chi ) ^ { 2 } } + \Gamma _ { 0 } \rho _ { 0 } c _ { V } ( T - T _ { 0 } ),$$

where 𝜒 = 1 -𝜌 0 𝜌 , 𝜌 0 and 𝑇 0 are the reference density and temperature, respectively, 𝑐 𝑠 is the speed of sound, Γ 0 is the Grüneisen parameter at the reference state, 𝑠 1 is the slope of the linear shock Hugoniot curve, and 𝑐 𝑉 is the specific heat capacity at constant volume. Out of these parameters, we keep the reference density 𝜌 0 fixed at 16.65 g/cc and the reference temperature 𝑇 0 fixed at 0.0253 eV. The parameter set { 𝑐 𝑠 , 𝑠 1 , Γ 0 , 𝑐 𝑉 } is varied as shown in Table 1.

Table 1. Matrix of parameter values used to develop the simulated dataset. All combinations of above parameters are used to simulate our data.

| Options                  | 1             | 2              | 3              | 4      | 5     |
|--------------------------|---------------|----------------|----------------|--------|-------|
| Γ 0                      | 1.6           | 1.7            | 1.76           | 1.568  | 1.472 |
| 𝑠 1                      | 1.22          | 1.464          | 1.342          |        |       |
| 𝑐 𝑠 [m/s]                | 339000        | 372900         | 305100         | 355000 |       |
| 𝑐 𝑉 [ erg g - 1 eV - 1 ] | 1 . 6 × 10 10 | 1 . 76 × 10 10 | 1 . 44 × 10 10 |        |       |

Altogether, the dataset realizes every unique parameter combination in a 7-dimensional parameter cube with 28 800 total simulations. Each hydrodynamic simulation is comprised of ,

(a) Example of the density evolution for an initial condition with inner surface perturbation profile 1 for time indices 0 , 5 , . . . , 35 . The images are 440x440 pixels representing the domain [ 0 , 𝐿 ] × [ 0 , 𝐿 ] .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-37e66fe2-ff65-4144-80f0-0d24d88fbce6.jpg)

Example Densities for Each Profile

(b) Examples of the density at time index 40 for each inner surface perturbation profile. The images are 150x150 pixels representing the domain h 0 , 15 44 𝐿 i × h 0 , 15 44 𝐿 i .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-d000be84-ce55-4ae1-b2f1-dacc52f86577.jpg)

Fig. 4. Example plots of the density evolution (a) and the various RMI profiles representing each inner surface perturbation profile (b).

density field snapshots at later times when the instability is present. We label these times as 𝑛 = 0 1 , , . . . , 40. An example of a density time series is shown in Figure 4. Once the Tantalum shell has collapsed, a shock is formed and reflected from the axis. The shock then interacts with the perturbed inner Tantalum edge. This creates an RMI. The topology of this interior evolves as depicted in Figure 4. The expanding shock proceeds to propagate into the non-constant dynamic density background. We chose frames corresponding to the time instants at 𝑛 = 25 30 35 40 to , , , train the network in our studies.

## 4.2. Generation of Synthetic Radiographs

Synthetic radiographs are produced at each time step. The create\_sino\_3d function from the ASTRAToolbox [56] is used to evaluate the forward operator in (1). Equation (2) is used to obtain the direct radiograph signal from the areal mass using 𝐼 0 = 3 201 . · 10 -4 , 𝜉 ( gas ) = 9 40 cm 2 . / g, 𝜉 ( Tantalum ) = 𝜉 ( coll ) = 13 03 cm 2 . / g. The source blur kernel is a 2D Gaussian kernel with 𝜎 blur chosen randomly between 1 and 3 pixels with a random orientation between 5 and 26 degrees. The correlated scatter kernel is also a 2D Gaussian kernel with 𝜎 scatter chosen randomly between 10 and 30 pixels and scatter level, 𝜅 , chosen randomly between 10 and 30. The coefficients of the

Fig. 5. Upper left: radiograph before applying noise. Upper right: radiograph after applying noise. Bottom left: horizontal line out through the center of the radiographs. Bottom right: vertical line out through the center of the radiographs. In both line plots, the orange line corresponds the noisy radiograph and the blue line corresponds to the radiograph without noise.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-e20f7b67-39b5-4058-ad85-c19fe79fa211.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-4b39b9e9-96df-43db-ad90-ab58a74666d3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-36966b23-3086-4a47-b96c-b283ada4d02e.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-f2cd8b3f-8d6a-4d8a-b996-5ee0ddaeef0b.jpg)

background scatter field are chosen such that the level is randomly between 0.5 and 1.5 times the mean signal level in the center of the image and the tilt is between -10% and 10%. The level of the gamma noise is randomly set in the range (39,000, 50,000) and the level of the photon noise is randomly set in the range (350, 450). Each random number is generated independently for each simulation and time step.

Figure 2 shows an example of a density field at time index 40 and a synthetic radiograph, 𝑚 , generated using the above method. For the same example, Figure 5 shows the direct signal, 𝑑 , radiograph 𝑚 , and profiles across the center of the two images.

## 4.3. Generation of Shock Features

One of the primary aspects of the ICF shell dynamics is the evolution of the inner gas-metal interface, i.e., the growth of the instability. This is because the passage of the incoming and outgoing shocks through this interface renders it unstable to the RMI. Considering temporally evolving simulations, we are interested in times when the instability on this interface has permitted the growth of perturbations to the extent that the inner gas-metal interface displays significant asymmetry. As such, we assume that the interface as identified by the feature extraction procedure is not robust. That is, we expect that the interface as identified by the imaging and feature extraction procedures is sensitive to the dependent the measurement model. This is in contrast with the shock and outer edge features that we assume are robust. The robustness of these latter features is due to the symmetry of the setup and dynamics and the stable nature of their evolution. Nevertheless, because of its passage across the unstable inner gas-metal interface, we expect the stably evolving shock to be imprinted with a decaying set of perturbations that can be reliably identified.

Accordingly, we have extracted shock and edge features at each time for each sequence of density fields. An edge-detection algorithm was utilized so as to enable the determination of a parametric representation of the shock and edge as a function of polar angle [28]. These features are subsequently compressed into a low-dimensional representation in terms of cosine harmonic coefficients,

$$r ^ { ( i ) } ( \theta ) = \sum _ { j = 0 } ^ { N ^ { ( i ) } } F _ { j } ^ { ( i ) } \cos ( 2 j \theta ),$$

for 𝑖 = shock edge. We found that , 𝑁 ( shock ) = 8 and 𝑁 ( edge ) = 5 can represent the shock and edge features with sufficient accuracy across the dataset.

## 4.4. Noise Model for Feature Extraction

The focus of this paper is to develop a features-to-density network to be used as the second component in a radiograph-to-features-to-density pipeline for the purpose of quantifying the growth rates of RMI in double shell ICF environments. For simplicity and modularity, we built each component independently, where the feature extractor was trained on radiograph and feature data pairs and the features-to-density network was trained on feature and density data pairs. To test the performance of the trained features-to-density network, we developed a noise model to simulate the errors that would arise from the feature extraction model.

The feature extraction model is a CNN consisting of an image Fourier feature encoding (IFFE) layer, convolution layers, and fully connected layers. The feature extractor is applied to a subset of the synthetic radiographs. Errors between the reconstructed cosine harmonic coefficients and that of the ground truth were calculated and used to compute a sample mean and covariance matrix for the shock. Using these computed statistics, we simulated the error produced from reconstructing the cosine harmonic coefficients from synthetic radiographs using a multivariate Gaussian noise model. During model testing, we added noise consistent with this distribution to the inputs to better characterize the performance of the model.

## 5. Description of Features-to-Density Architectures

We compared two different architectures for solving the features-to-density problem. These approaches involve a generative variational autoencoder (VAE) network based on the vision transformer [42], which we refer to as the ShockDecoderViT, and a deterministic structurepreserving network based on the original transformer [43], which we call the Mass-Conserving Transformer. Due in part to the presence of radiographic noise, the features identified by the radiograph-to-features network can only be determined up to a certain level of precision. The two networks have different approaches for handling this uncertainty. The generative approach builds uncertainty into the model and is trained to minimize degeneracy of the density reconstructions. The deterministic network uses mass-conservation to improve robustness of predictions.

Both architectures are trained on the features-to-density problem using data at time steps 𝑛 = 25 30 35 40. The training set consists of 80% ( , , , 𝑁 = 23,040) randomly selected density time series for each inner surface perturbation profile and the testing set consists of the remaining 20% ( 𝑁 = 5,760) of data. Since most of the action is downstream of the outgoing shock, we primarily restrict our attention to reconstructing a smaller region encapsulating the gas metal interface. For purposes of comparison, the architectures were trained to output a 150 × 150 square pixel image representing the density field in the cylindrical domain, GLYPH&lt;2&gt; 0 , 15 44 𝐿 GLYPH&lt;3&gt; × GLYPH&lt;2&gt; 0 , 15 44 𝐿 GLYPH&lt;3&gt; . This smaller domain window captures the most interesting physics. Additionally, the Mass-Conserving Transformer architecture was trained to output a 440 × 440 square pixel image representing the full domain of the simulation data, [ 0 , 𝐿 ] × [ 0 , 𝐿 ] .

## 5.1. ShockDecoderViT: A Vision Transformer-based Shock Decoder

We propose a new generative architecture to decode the density field from the shock features based on the vision transformer [42] and a conditional variational autoencoder (cVAE) learning framework [57-59]. In particular, the data consists of the density field 𝜌 𝑡 ( ) and a corresponding small number of numeric features characterizing the shock F( ) 𝑡 at four times in a dynamically evolving flow. The aim is to obtain a trained model (in this case just the decoder part of the cV AE) that, given a set of time evolving shock features alone-the conditions, generates a sequence of 2D density fields at the corresponding times. We call the decoder a ShockDecoderViT because the details of the flow can be reconstructed given a minimal set of shock characteristics. In common with most current ML approaches, the utility of this model is largely restricted to the training data distribution.

The cVAE consists of an encoder and a decoder. First, the encoder embeds the input images and conditions and processes them using the transformer. Here ideas of the vision transformer are used to process the density field. That is, the density field is considered as an image that is first split into patches which serve as tokens. In the present setting, only the region interior of the outermost extent of the outgoing shock is considered for the reason that most of the rapid and large dynamical changes are confined to this region. In particular, the interior most 150 × 150 region is partitioned into 10 patches of size 15 × 15. The patches and the conditions, which are taken to be the zeroth through seventh cosine harmonics of the outgoing shock at the four times, are then projected onto the tokens' embedding space. After adding positional encoding to the projections, dot-product attention-based transformer blocks are applied to obtain a compact variational latent space representation of the temporal sequence of densities and shock features at the bottleneck. That is, the encoder produces the mean 𝜇 𝜌 𝑡 , ( ( ) F( )) 𝑡 and variance Σ ( 𝜌 𝑡 ( ) , F( )) 𝑡 of the latent space features as functions of density and shock features. Thereafter, the decoder takes inputs of the latent variables realizations and the shock features and embeds them into a series of tokens that are processed using attention-based transformer blocks before they are projected back into density space.

The transformer-based encoder and decoder (ShockDecoderViT) of the cVAE are trained simultaneously with the Evidence Lower Bound (ELBo) of variational inference as the loss function. This loss function may be thought of as having two components: the reconstruction loss of the autoencoder and the KL divergence between the latent distribution produced by the encoder and a standard multivariate normal. In the testing phase, only the decoder is used; it takes as input the shock features. The shock features are combined with random realizations of the latent space variables to then generate the density fields.

In summary, our new architecture, the ShockDecoderViT, combines the strengths of the dot product attention-based Vision Transformer and the conditional variational autoencoder to generate temporally coherent sequences of density fields from dynamically evolving sequences of shock features. While the patch-based approach allows it to handle large images efficiently, the transformer provides the ability to capture long range dependencies in both the spatial and temporal domains. Results are presented for an architecture that uses two transformer blocks in the encoder and eight transformer blocks for the decoder. Increasing the number of encoder transformer blocks and small variations of the number of transformer blocks in the decoder led to minor changes in results (not shown).

## 5.2. Mass-Conserving Transformer Architecture

In addition to considering the ShockDecoderViT, which is a generative network, we also investigated a purely deterministic architecture. The Mass-Conserving Transformer network is a structure-preserving architecture that uses the transformer blocks introduced in [43]. While the ShockDecoderViT uses a vision transformer in both the encoder and decoder to go between latent representations and density field image patches, the Mass-Conserving Transformer treats the input shock features using a series of transformer blocks before projecting to density. The network uses a mass-conservation layer to enforce that the dynamics conserve mass between time steps.

The input to the Mass-Conserving Transformer is a matrix 𝑥 ∈ R 𝑁 𝑡 × 𝑁 𝑓 , and the output is 𝑦 ∈ R 𝑁 𝑡 × 𝑁 𝑟 × 𝑁 𝑧 . 𝑁 𝑡 is the number of times in the sequence, 𝑁 𝑓 is the total number of features, which consists of shock features a temporal encoding vector which represents polynomials in 𝑛 , 𝑛 𝑖 , for 𝑖 = 0 1 , , . . . 𝑁 𝑝 , and 𝑁 𝑟 × 𝑁 𝑧 is the resolution of the output image. The shock features are each centered and normalized using the mean and standard deviation of the training set.

The features are first inputted into a network consisting of multiple transformer blocks,

$$\chi = T _ { B } \circ \cdots \circ T _ { 1 } ( x ),$$

where 𝜒 ∈ R 𝑁 𝑡 × 𝑁 𝑓 is a matrix representing a latent representation of the density fields at each time in the sequence. The transformer blocks have the same structure as the blocks proposed in [43]. Each block applies dot-product attention to incorporate temporal dependencies.

After the series of transformer blocks, the latent variables are projected into density fields. A critical component of projection step is the learning of characteristic eigenfunctions, Z ( ) 𝑖 ∈ R 𝑁 𝑟 × 𝑁 𝑧 , for 𝑖 = 1 , . . . , 𝑁 Z on the image plane. The first set of transformations is given by

$$\mathcal { Y } _ { i j } ^ { n } = h \left ( \mathcal { Z } _ { i j } ^ { ( 1 ) }, \dots, \mathcal { Z } _ { i j } ^ { ( N z ) } ; \, \theta _ { n } \right ), \quad i = 1, \dots, N _ { r }, \ \ j = 1, \dots, N _ { z }, \ \ n = 1, \dots, N _ { t }, \ ( 1 3 )$$

where ℎ is a feedforward neural network parameterized by weights, 𝜃 𝑛 , which are in turn determined by

$$\theta _ { n } = g ( \chi _ { n } ), \quad n = 1, \dots, N _ { t },$$

where 𝑔 is also a feedforward neural network. The intermediate output, Y ∈ R 𝑁 𝑡 × 𝑁 𝑟 × 𝑁 𝑧 , represents a sequence of images with the same resolution at the final output 𝑦 . The reconstructed density is obtained by passing the sequence, Y 𝑛 , 𝑛 = 1 , . . . , 𝑁 𝑡 , through a mass conservation layer. Let M be the total mass of the object and let

$$M ( \gamma ) = \sum _ { i, j } 4 \pi r _ { i } \gamma _ { i j } \Delta r \Delta z,$$

be the calculation of mass of in the image 𝛾 , where 𝑟 𝑖 is the radius at index 𝑖 and 𝛾 𝑖 𝑗 is the density value in pixel ( 𝑖, 𝑗 ) . The output density reconstruction is the sequence of images

$$y ^ { n } = \frac { \mathcal { M } } { M ( | \mathcal { Y } ^ { n } | ) } | \mathcal { Y } ^ { n } |, \quad n = 1, \dots, N _ { t }.$$

Our architecture used 6 consecutive transformer blocks, 𝐵 = 6, each with 8 heads, 𝐻 = 8, a latent dimension of 64, 𝑘 = 64, and a feedforward neural network with inner dimension 2048 and tanh activation function. For the projection layer, we chose 𝑁 Z = 40 eigenfunctions and represented both 𝑔 and ℎ using feedforward neural networks with 2 hidden layers and latent dimension of 100. The architecture was trained to minimize the root-mean-squared error (RMSE) of the density reconstructions on the training set.

## 6. Results

This section presents results of the trained ShockDecoderViT and Mass-Conserving Transformer networks. First we compare the two approaches using ensemble metrics such as root-meansquared error (RMSE) and structural similarity. Section 6.1 examines the generative nature of ShockDecoderViT. Section 6.2 analyzes the density reconstructions of each inner surface perturbation profile and its corruptions due to errors arising from the radiograph-to-features network. Additionally, we use an example reconstruction to demonstrate how peak-to-trough evolution can be calculated accurately, as indicated by capturing the RMI growth rates from the reconstruction and comparing to those obtained from the ground truth simulations.

|                             | RMSE (g/cc)   | RMSE (g/cc)   | RMSE (g/cc)   | Structural Similarity   | Structural Similarity   | Structural Similarity   |
|-----------------------------|---------------|---------------|---------------|-------------------------|-------------------------|-------------------------|
| Architecture                | min           | average       | max           | min                     | average                 | max                     |
| ShockDecoderViT             | 0.023         | 0.044         | 0.148         | 0.58                    | 0.95                    | 0.98                    |
| Mass Conserving Transformer | 0.017         | 0.043         | 0.147         | 0.76                    | 0.97                    | 1.00                    |

Fig. 6. Histogram of root-mean-squared errors (RMSE) (top left) and structural similarity (top right) between the density reconstruction and ground truth for the testing set for the ShockDecoderViT and Mass-conserving transformer. Bottom: table of summary statistics for the above histograms.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-b335b1e3-b8e7-4f07-9c45-04a39db5aa12.jpg)

Figure 6 shows histograms of RMSE and structural similarity evaluated on the testing set for both architectures along with a table of summary statistics. Both networks exhibit similar performance in RMSE, with the Mass-Conserving Transformer having better minimum, average, and maximum errors compared to the ShockDecoderViT. The Mass-Conserving Transformer significantly outperforms the ShockDecoderViT in structural similarity, with a mode shifted to the right and better minimum, average, and maximum structural similarity values.

## 6.1. ShockDecoderViT Density Reconstruction

This section examines the generative nature of the ShockDecoderViT network. For each set of shock features in the testing set, we generated 𝑁 = 7 realizations of density reconstructions using the ShockDecoderViT and computed the standard deviation of the density, RSME, and structural similarity. Figure 7 summarizes the results of this study, showing histograms of the standard deviations of both RMSE and structural similarity, along with summary statistics in a table. The variation in the RMSE and Similarity metric are very small relative to the ensemble averages, indicating that the variations due to the degeneracy are small. We reiterate that the standard deviation of deterministic methods, including the Mass-Conserving Transformer, on a fixed input is zero by definition.

Figure 8 shows examples of generations from six sets of shock features. These six features were chosen to correspond to the extremities of the standard deviation histogram, minimums, averages, and maximums for the standard deviation of RMSE and structure similarity. For each set of features, we chose the best, average, and worst performing reconstruction to display. Additionally, the standard deviation of density is shown in the last column. For the best and

Fig. 7. Histogram of the standard deviation of RMSE (left) and SS (right) for 𝑁 = 7 generations of each example in the testing set. Bottom: table of summary statistics for the above histograms.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-2fb1c339-28bc-4d65-b8d2-e2cc8eafd55f.jpg)

| Architecture                | Std of RMSE ( 𝑁 = 7 ) (g/cc)   | Std of RMSE ( 𝑁 = 7 ) (g/cc)   | Std of RMSE ( 𝑁 = 7 ) (g/cc)   | Std of Structural Similarity ( 𝑁 = 7 )   | Std of Structural Similarity ( 𝑁 = 7 )   | Std of Structural Similarity ( 𝑁 = 7 )   |
|-----------------------------|--------------------------------|--------------------------------|--------------------------------|------------------------------------------|------------------------------------------|------------------------------------------|
|                             | min                            | average                        | max                            | min                                      | average                                  | max                                      |
| ShockDecoderViT             | 4.82e-06                       | 2.56e-04                       | 2.27e-02                       | 2.09e-05                                 | 5.83e-04                                 | 4.14e-02                                 |
| Mass-Conserving Transformer | 0                              | 0                              | 0                              | 0                                        | 0                                        | 0                                        |

average examples, the reconstructions are nearly identical under visual inspection, while the worst performing reconstruction exhibits significant errors from the ground truth and are visually different from the other reconstructions in the example. This is also confirmed by examining the large values present in the density standard deviation heat map. This extreme example is however one example of a number of outliers represented by the narrow wide tail of the histogram in Figure 7.

## 6.2. Mass-Conserving Transformer Density Reconstruction

This section analyzes the density reconstructions produced by the Mass-Conserving Transformer network. Figure 9 shows the training loss, representing RMSE over the training set, verses the training epoch. The trained model was used to produce the density reconstructions and their corresponding errors shown in Figure 10. Each example reconstruction in this figure, labeled 1 through 20, corresponds to a different set of shock features, chosen randomly, corresponding to each inner-surface perturbation profile. The labeled structural similarity and RMSE were computed for the ensemble of testing data for each profile. The top plot (a) displays the density reconstructions sorted from lowest to highest structural similarity. The bottom plot shows the corresponding errors of the above plot, labeled by profile ensemble average RMSE in g/cc. Despite the low dimensional feature space, the density, including the complex details of the RMI, can be reconstructed to a high level of accuracy.

Next, we investigated the effect of errors from the radiograph-to-features network on the Mass-Conserving Transformer results. For each set of shock features in the testing set, a random error is generated using the noise model described in Section 4.4. Reconstructions were produced for this generation of random error multiplied by factors of 0 , 1 , . . . , 19, respectively, and the corresponding RMSE was computed. Figure 11 shows box plots summarizing the error statistics of this study for each noise multiplier. A clear degradation in accuracy is observed as the noise multiplier is increased. For the case of 1 × multiplier, which corresponds to the expected error produced from the radiograph-to-features network, the boxplot represents a small shift from the clean features. The trend towards greater noise remains gradual and bounded, demonstrating the

(a) Multiple generative reconstructions for examples with the lowest (top row), average (middle row), and highest (bottom row) standard deviation in RMSE.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-cb74845f-3de6-48a8-b8d5-130cf18a4c50.jpg)

(b) Multiple generative reconstructions for examples with the lowest (top row), average (middle row), and highest (bottom row) standard deviation in SS.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-edbd6a75-47ff-4740-8d32-55d934f28f39.jpg)

Fig. 8. Examples of generative reconstructions produced by the ShockDecoderViT corresponding to the lowest, average, and highest standard deviation in RMSE (a) and Structural Similarity (b). In both sub-figures, the left-most column is the ground truth, the middle three columns correspond to the best, average, and worst errors, and the right-most column shows the standard deviation in density for the 𝑁 = 7 generative reconstructions for the corresponding case. Each plot shows quarter planes corresponding data at each of the four times in the sequence counter-clockwise from the first quadrant.

robustness of the model, which may be attributed to it's structure-preserving properties.

Next, we studied the expressive power of the attention mechanism used in the Mass-Conserving transformer. By default, the dot product attention represents a fully connected network between each of the 4 terms of the sequence. For comparison, we also considered a sequence length 1 network, formed by breaking all the cross-temporal connections of the network, and a sequence length 2 network, formed by breaking the cross-temporal connections between the first two times and the latter two times. Figure 12 shows the training loss of the three approaches. The results indicate that incorporating more temporal correlations led to smaller training losses.

While our previous results focused on reconstructing density in a smaller region around the RMI, we also demonstrate that our approach can be applied to reconstruct the density in the entirety of the domain using the Mass-Conserving Transformer. A new model was trained using the full 440 × 440 domain available in the data set to produce a new set of architecture weights. The two plots in Figure 13 show histograms of RMSE and structural similarity evaluated on the

Fig. 9. Root mean squared error verses epoch during training.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-f3de67fd-2a02-4891-bad5-7c6d9a74fef0.jpg)

testing set. We observe that the complexities of the gas metal interface are still reconstructed by the full domain model. Figure 15 shows an example of density reconstructions and corresponding horizontal and vertical line-outs through the center and Figure 14 shows a zoomed in view of the reconstruction and its error. In these examples, the RMI and shock location are still captured accurately.

Next, we demonstrate that the Mass-Conserving Transformer can accurately measure properties of the RMI such as the growth of the peak-to-trough radial distance. We chose a set of shock features, applied a random perturbation generated from the noise model of the radiographto-features network, and performed a density reconstruction. First we identified pixels that are contained in the gas domain, which is contained in a region near ( 𝑟, 𝑧 ) = ( 0 0 , which , ) corresponds to a region of significantly lower density compared to that of the surrounding metal. A discontinuity separates the interface between the gas and metal. On both sides of the discontinuity front, the densities are approximately constant. In the reconstruction, there is a small region consisting of less than a few pixels where the density transitions smoothly between the nominal fluid density and the nominal solid density. Therefore, we attempt to split the difference. We find the maximum density of the metal, 𝜌 max and take the fluid domain to be the region where the density is less than 1 2 𝜌 max. From this, we identify the pixel locations on the edge of the domain with the furthest and closest distances to the origin. Figure 16 shows the result of this identification algorithm applied to both the reconstruction and the ground truth for two examples of RMI profiles, representing high and low frequency modes. Additionally, we apply a Canny edge filter for various choices of the smoothing and threshold parameters to the noisy radiograph to identify the peak and trough of the RMI. Figure 16 shows a line plot comparing the peak-to-trough evolution between the ground truth, the reconstruction, and Canny edge filter. The density reconstruction approach is successful in accurately identifying the peak and trough locations and consequently their radial distance. The Canny edge filter however is unable to obtain an acceptable level of accuracy and suffers from large variations due to choices in meta-parameters. Figure 2 shows the resulting edges that are detected for the high frequency profile (profile 1). This demonstrates that our density reconstruction algorithm is a viable approach for making growth-rate estimates in the evolution of a spherically-symmetric RMI. We cannot overemphasize the importance of this discovery. That is, this method of reconstruction may allow the RMI growth rates to be experimentally verified in a spherically convergent geometry for the first time.

## 7. Conclusions

This work presents a new density reconstruction approach that uses a trained attention-based transformer network to enable accurate density reconstructions from a series of noisy radiographic images. The key component of this network is the transformer encoder that acts on a sequence of features extracted from noisy radiographs. This encoder includes numerous self-attention layers that act to learn temporal dependencies in the input sequences and increase the expressiveness of the model.

The two architectures had comparable performance, with the Mass-conserving Transformer slightly outperforming the ShockDecoderViT. Furthermore, the shock harmonic features provided sufficient constraints to limit the variability in the density predictions, with variations that are insignificant compared to the mean reconstruction errors.

Our examinations also indicated the benefit of the dynamics relative to the traditional approach of reconstruction of single images. Finally, we demonstrated the ability to use our sparse features to accurately reconstruct the fine features of the image, enabling, for the first time, the determination of the RMI growth rates. This work opens the potential to experimentally capture RMI growth rates in spherical geometries.

## A. Cosine Coefficients for Inner Surface Perturbation Profile

The coefficients of the cosine harmonic series of the initial inner surface perturbation profile is scaled according to 𝐹 𝑖 = 𝑅 in ¯ 𝐹 𝑖 / 8 , for 𝑖 = 0 , . . . , 8, where ¯ 𝐹 0 = 8, ¯ 𝐹 5 = ¯ 𝐹 7 = 0, and the rest of the coefficients are provided by Table 2.

Funding. This work was supported by the U.S. Department of Energy through the Los Alamos National Laboratory (LANL) and the Laboratory Directed Research and Development program of LANL.

Acknowledgments. The authors thank Jeff Fessler for helpful comments on the draft and Oleg Korobkin for extracting features for the data set.

Disclosures. The authors declare no conflicts of interest.

Data Availability Statement. Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

## References

- 1. R. D. Richtmyer, 'Taylor instability in shock acceleration of compressible fluids,' Tech. rep., Los Alamos Scientific Lab., N. Mex. (1954).
- 2. E. E. Meshkov, 'Instability of the interface of two gases accelerated by a shock wave,' Fluid Dyn. 4 , 101-104 (1969).
- 3. Rayleigh, 'Investigation of the character of the equilibrium of an incompressible heavy fluid of variable density,' Proc. Lond. mathematical society 1 , 170-177 (1882).
- 4. G. I. Taylor, 'The instability of liquid surfaces when accelerated in a direction perpendicular to their planes. I,' Proc. Royal Soc. London. Ser. A. Math. Phys. Sci. 201 , 192-196 (1950).
- 5. D. H. Sharp, 'An overview of Rayleigh-Taylor instability,' Phys. D: Nonlinear Phenom. 12 , 3-18 (1984).
- 6. P. Drazin, 'Kelvin-Helmholtz instability of finite amplitude,' J. Fluid Mech. 42 , 321-335 (1970).
- 7. J. Kane, R. Drake, and B. Remington, 'An evaluation of the Richtmyer-Meshkov instability in supernova remnant formation,' The Astrophys. J. 511 , 335 (1999).
- 8. J. Holzfuss, 'Surface-wave instabilities, period doubling, and an approximate universal boundary of bubble stability at the upper threshold of sonoluminescence,' Phys. Rev. E 77 , 066309 (2008).
- 9. G. B. Goodwin and E. S. Oran, 'Premixed flame stability and transition to detonation in a supersonic combustor,' Combust. Flame 197 , 145-160 (2018).
- 10. V. Rupert, 'Shock-interface interaction: current research on the Richtmyer-Meshkov problem,' in Shock Waves: Proceedings of the 18th International Symposium on Shock Waves, Held at Sendai, Japan 21-26 July 1991, (Springer, 1992), pp. 83-94.
- 11. M. Brouillette, 'The Richtmyer-Meshkov instability,' Annu. Rev. Fluid Mech. 34 , 445-468 (2002).
- 12. Y. Zhou, R. J. Williams, P. Ramaprabhu, et al. , 'Rayleigh-Taylor and Richtmyer-Meshkov instabilities: A journey through scales,' Phys. D: Nonlinear Phenom. 423 , 132838 (2021).
- 13. E. Leinov, G. Malamud, Y. Elbaz, et al. , 'Experimental and numerical investigation of the Richtmyer-Meshkov instability under re-shock conditions,' J. Fluid Mech. 626 , 449-475 (2009).
- 14. R. L. Holmes, G. Dimonte, B. Fryxell, et al. , 'Richtmyer-Meshkov instability growth: experiment, simulation and theory,' J. Fluid Mech. 389 , 55-79 (1999).
- 15. Y. Zhou, 'Rayleigh-Taylor and Richtmyer-Meshkov instability induced flow, turbulence, and mixing. II,' Phys. Reports 723 , 1-160 (2017).
- 16. Q. Zhang and M. J. Graham, 'A numerical study of Richtmyer-Meshkov instability driven by cylindrical shocks,' Phys. Fluids 10 , 974-992 (1998).
- 17. M. B. Prime, W. T. Buttler, S. J. Fensin, et al. , 'Tantalum strength at extreme strain rates from impact-driven Richtmyer-Meshkov instabilities,' Phys. Rev. E 100 , 053002 (2019).
- 18. M. B. Prime, W. T. Buttler, M. A. Buechler, et al. , 'Estimation of metal strength at very high rates using free-surface Richtmyer-Meshkov instabilities,' J. Dyn. Behav. Mater. 3 , 189-202 (2017).
- 19. E. C. Merritt, J. P. Sauppe, E. N. Loomis, et al. , 'Experimental study of energy transfer in double shell implosions,' Phys. Plasmas 26 , 052702 (2019).
- 20. K. Prestridge, P. Vorobieff, P. Rightley, and R. Benjamin, 'Validation of an instability growth model using particle image velocimetry measurements,' Phys. review letters 84 , 4353 (2000).
- 21. M. Emery, J. Gardner, R. Lehmberg, and S. Obenschain, 'Hydrodynamic target response to an induced spatial incoherence-smoothed laser beam,' Phys. Fluids B: Plasma Phys. 3 , 2640-2651 (1991).
- 22. R. Ishizaki and K. Nishihara, 'Propagation of a rippled shock wave driven by nonuniform laser ablation,' Phys. review letters 78 , 1920 (1997).
- 23. J. D. Lindl, R. L. McCrory, and E. M. Campbell, 'Progress toward ignition and burn propagation in inertial confinement fusion,' Phys. Today 45 , 32-40 (1992).
- 24. R. Taylor, A. Velikovich, J. Dahlburg, and J. Gardner, 'Saturation of laser imprint on ablatively driven plastic targets,' Phys. review letters 79 , 1861 (1997).
- 25. E. Lifshin, X-ray Characterization of Materials (John Wiley &amp; Sons, 2008).
- 26. D. C. Swift, A. L. Kritcher, J. A. Hawreliak, et al. , 'Absolute Hugoniot measurements from a spherically convergent shock using X-ray radiography,' Rev. Sci. Instruments 89 (2018).
- 27. J. Rygg, O. Jones, J. Field, et al. , '2D x-ray radiography of imploding capsules at the National Ignition Facility,' Phys. Rev. Lett. 112 , 195001 (2014).
- 28. M. Hossain, B. T. Nadiga, O. Korobkin, et al. , 'High-precision inversion of dynamic radiography using hydrodynamic features,' Opt. Express 30 , 14432-14452 (2022).
- 29. J. Radon, 'Uber die bestimmung von funktionen durch ihre integralwerte längs gewisser mannigfaltigkeiten (on the determination of function from their integral value along certain manifolds),' Berichte Sächsische Akadamie der Wissenschaften, Leipzig, Math.-Phys. Kl 69 , 262-267 (1917).
- 30. A. M. Cormack, 'Representation of a function by its line integrals, with some radiological applications,' J. Appl. Phys. 34 , 2722-2727 (1963).
- 31. L. A. Feldkamp, L. C. Davis, and J. W. Kress, 'Practical cone-beam algorithm,' JOSA A 1 , 612-619 (1984).
- 32. R. N. Bracewell and R. N. Bracewell, The Fourier transform and its applications , vol. 31999 (McGraw-Hill New York, 1986).
- 33. I. A. Elbakri and J. A. Fessler, 'Statistical image reconstruction for polyenergetic x-ray computed tomography,' IEEE transactions on medical imaging 21 , 89-99 (2002).
- 34. S. Ravishankar, J. C. Ye, and J. A. Fessler, 'Image reconstruction: From sparsity to data-adaptive methods and machine learning,' Proc. IEEE 108 , 86-109 (2019).
- 35. M. Pasha, A. K. Saibaba, S. Gazzola, et al. , 'A computational framework for edge-preserving regularization in dynamic inverse problems,' Electron. Trans. on Numer. Anal. 58 , 486-516 (2023).
- 36. D. F. Yu and J. A. Fessler, 'Edge-preserving tomographic reconstruction with nonlocal regularization,' IEEE transactions on medical imaging 21 , 159-173 (2002).

Table 2. Scaled cosine coefficients for the initial shell profile used for each profile.

|   profile |    𝐹 1 |    𝐹 2 |    𝐹 3 |   𝐹 4 |    𝐹 6 |   𝐹 8 |
|-----------|--------|--------|--------|-------|--------|-------|
|         1 | 0      | 0      | 0      | 0     | 0      | 0.08  |
|         2 | 0      | 0      | 0      | 0.08  | 0      | 0     |
|         3 | 0      | 0.08   | 0      | 0     | 0      | 0     |
|         4 | 0      | 0      | 0      | 0     | 0      | 0.075 |
|         5 | 0      | 0      | 0      | 0.075 | 0      | 0     |
|         6 | 0      | 0.075  | 0      | 0     | 0      | 0     |
|         7 | 0      | 0.0075 | 0      | 0     | 0.0025 | 0.065 |
|         8 | 0.0075 | 0      | 0.0025 | 0.065 | 0      | 0     |
|         9 | 0.005  | 0.0657 | 0      | 0     | 0      | 0     |
|        10 | 0      | 0      | 0      | 0     | 0      | 0.06  |
|        11 | 0      | 0      | 0      | 0.06  | 0      | 0     |
|        12 | 0      | 0.06   | 0      | 0     | 0      | 0     |
|        13 | 0      | 0      | 0      | 0     | 0      | 0.055 |
|        14 | 0      | 0      | 0      | 0.055 | 0      | 0     |
|        15 | 0      | 0.055  | 0      | 0     | 0      | 0     |
|        16 | 0      | 0.0075 | 0      | 0     | 0.0025 | 0.045 |
|        17 | 0.0075 | 0      | 0.0025 | 0.045 | 0      | 0     |
|        18 | 0.0051 | 0.0457 | 0      | 0     | 0      | 0     |
|        19 | 0      | 0      | 0      | 0.04  | 0      | 0     |
|        20 | 0      | 0.04   | 0      | 0     | 0      | 0     |

37. A. H. Delaney and Y. Bresler, 'Globally convergent edge-preserving regularized reconstruction: an application to limited-angle tomography,' IEEE Trans. on Image Process.

7

, 204-221 (1998).

- 38. M. Unberath, J.-N. Zaech, C. Gao, et al. , 'Enabling machine learning in x-ray-based procedures via realistic simulation of image formation,' Int. J. Comput. Assist. Radiol. Ssurgery 14 , 1517-1528 (2019).

39. S. Sajed, A. Sanati, J. E. Garcia, et al.

,

'The effectiveness of deep learning vs. traditional methods for lung disease diagnosis using chest x-ray images: A systematic review,' Appl. Soft Comput. p. 110817 (2023).

- 40. C. Sun, Y. Liu, and H. Yang, 'Degradation-aware deep learning framework for sparse-view CT reconstruction,'
- Tomography 7 , 932-949 (2021).
- 41. Z. Huang, M. Klasky, T. Wilcox, and S. Ravishankar, 'Physics-driven learning of Wasserstein GAN for density reconstruction in dynamic tomography,' Appl. Opt. 61 , 2805-2817 (2022).
- 42. A. Dosovitskiy, L. Beyer, A. Kolesnikov, et al. , 'An image is worth 16x16 words: Transformers for image recognition at scale,' in International Conference on Learning Representations, (2021).
- 43. A. Vaswani, N. Shazeer, N. Parmar, et al. , 'Attention is all you need,' in Advances in Neural Information Processing Systems, vol. 30 I. Guyon, U. V. Luxburg, S. Bengio, et al. , eds. (Curran Associates, Inc., 2017).
- 44. X. Xu, M. Klasky, J. A. Fessler, et al. , 'An end-to-end learning approach for subpixel edge feature extraction,' Radiogr. Imaging Tomogr. Summer Symp. (2023).
- 45. M. Berger, J. Hubbell, S. Seltzer, et al. , 'XCOM: Photon cross sections database (version 1.5),' NIST Stand. Ref. Database 8 , 87-3597 (2010).
- 46. M. Sun and J. M. Star-Lack, 'Improved scatter correction using adaptive scatter kernel superposition,' Phys. Med. Biol. 55 , 6695-6720 (2010).
- 47. M. T. McCann, M. L. Klasky, J. L. Schei, and S. Ravishankar, 'Local models for scatter estimation and descattering in polyenergetic X-ray tomography,' Opt. Express 29 , 29423-29438 (2021).
- 48. I. Elbakri and J. Fessler, 'Statistical image reconstruction for polyenergetic X-ray computed tomography,' IEEE Trans. on Med. Imaging 21 , 89-99 (2002).
- 49. G. R. Myers, A. M. Kingston, T. K. Varslot, et al. , 'Dynamic tomography with a priori information,' Appl. Opt. 50 , 3685-3690 (2011).
- 50. S. Bonnet, A. Koenig, S. Roux, et al. , 'Dynamic X-ray computed tomography,' Proc. IEEE 91 , 1574-1587 (2003).
- 51. L. Desbat, S. Rit, R. Clackdoyle, et al. , 'Algebraic and analytic reconstruction methods for dynamic tomography,' in 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, (IEEE, 2007), pp. 726-730.
- 52. C. Jailin and S. Roux, 'Dynamic tomographic reconstruction of deforming volumes,' Materials 11 , 1395 (2018).
- 53. E. F. Toro, Riemann solvers and numerical methods for fluid dynamics: a practical introduction (Springer Science &amp; Business Media, 2013).
- 54. P. D. Bello-Maldonado, T. V. Kolev, R. N. Rieben, and V. Z. Tomov, 'A matrix-free hyperviscosity formulation for high-order ale hydrodynamics,' Comput. &amp; Fluids 205 , 104577 (2020).
- 55. Hertel, E. S. Jr. and Kerley, G. I., 'CTH Reference Manual: The Equation of State Package,' report SAND98-0947, Sandia National Laboratories (1998).
- 56. W. van Aarle, W. J. Palenstĳn, J. De Beenhouwer, et al. , 'The ASTRA toolbox: A platform for advanced algorithm development in electron tomography,' Ultramicroscopy 157 , 35-47 (2015).
- 57. D. P. Kingma and M. Welling, 'Auto-encoding variational Bayes,' CoRR abs/1312.6114 (2013).
- 58. L. P. Cinelli, M. A. Marins, E. A. B. Da Silva, and S. L. Netto, Variational methods for machine learning with applications to deep networks (Springer, 2021).
- 59. D. P. Kingma, S. Mohamed, D. Jimenez Rezende, and M. Welling, 'Semi-supervised learning with deep generative models,' Adv. Neural Inf. Process. Syst. 27 (2014).

Reconstructions and Testing Set Structural Similarity for each Profile

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-a32b65ea-5d3a-4d1f-8bdc-6b4eb7973df5.jpg)

(a) Mass-Conserving Transformer: Examples of density reconstructions for each inner surface perturbation profile sorted by the profile ensemble average Structural Similarity in the testing set. The plots are labeled by their profile number and the corresponding profile ensemble average Structural Similarity.

Reconstruction Error and Testing Set RSME for each Profile

(b) Mass-Conserving Transformer: Corresponding errors between the model and the ground truth on a log scale for the above reconstructions, labeled by their profile number and the corresponding profile ensemble average RSME in g/cc.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-ca58812e-5208-4612-9333-74f86e0611b1.jpg)

Fig. 10. Mass-Conserving Transformer: (a) Density reconstructions and (b) their errors for a random choice of simulation parameters corresponding to each inner surface perturbation profile. Each plot shows quarter planes corresponding data at each of the four times in the sequence counter-clockwise from the first quadrant.

Fig. 11. Box-plots of RSME (g/cc) evaluated on the testing set corresponding to noise multipliers 0 , 1 , . . . , 19. Each box extends from the first quartile to the third quartile of the data with an orange line representing the median. The whiskers extend from the box to the farthest data point lying within 1.5x the inter-quartile range from the box. The red dots represent points outside of the whiskers.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-5d1b242b-41f6-4fb0-95ca-79dbc5b3e2e5.jpg)

Fig. 12. Comparison of training loss when the attention blocks of the mass-conserving transformer are fully connected (sequence length 4), connected only between the first two and latter two times (sequence length 2), and fully disconnected (sequence length 1). For consistency between methods, RMSE is evaluated on the entire sequence of 4 density reconstructions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-472ad8b2-6904-49a6-b83e-5dc474b46c7b.jpg)

120

100

80

60

40

20

RMSE

Fig. 13. Histogram of root-mean-squared errors (RMSE) (top left) and structural similarity (top right) between the density reconstruction and ground truth for the testing set for the Mass-conserving transformer trained on the entire domain. Bottom: table of summary statistics for the above histograms.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-b88618dd-10a8-441e-9b9b-bb49faa63230.jpg)

|                             | RMSE (g/cc)   | RMSE (g/cc)   | RMSE (g/cc)   | Structural Similarity   | Structural Similarity   | Structural Similarity   |
|-----------------------------|---------------|---------------|---------------|-------------------------|-------------------------|-------------------------|
| Architecture                | min           | average       | max           | min                     | average                 | max                     |
| Mass Conserving Transformer | 0.010         | 0.019         | 0.049         | 0.92                    | 0.98                    | 0.99                    |

Fig. 14. Example of full density reconstructions and their corresponding errors produced by the Mass-Conserving Transformer. The input shock features are polluted with simulated noise from the radiograph-to-feature network.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-4402a20a-bb61-4f8a-ad19-bb15e94042e9.jpg)

RMSE on

Testing Set (440x440 Domain)

Mass-conserving Transformer (FD)

Fig. 15. Plot of reconstruction and horizontal and vertical density line plots produced by the Mass-Conserving Transformer trained on the entire domain. The input shock features are polluted with simulated noise from the radiograph-to-feature network.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-9f56f968-8d4c-41e1-9050-16f8da9d8b44.jpg)

(a) High frequency RMI mode example (profile 1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-3c9348a2-7729-4684-b773-69a65bb633c1.jpg)

RMI Peak and Trough Identification

(b) Low frequency RMI mode example (profile 9).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00985v1-figure-92610e76-ba51-4277-b01b-143de54d9779.jpg)

Fig. 16. Top plots: density fields and identified peak and trough points of the RMI (green markers) for the ground truth and reconstructions using inputs polluted with simulated noise from the radiograph-to-feature network. Bottom plots: evolution of the maximum RMI peak-to-trough radial distance for the ground truth and reconstructions corresponding to the above plots.
Daniel A. Serino Marc L. Klasky Balasubramanya T. Nadiga Xiaojian Xu Trevor Wilcox
Reconstructing Richtmyer-Meshkov instabilities from noisy radiographs using low dimensional features and attention-based neural networks
A trained attention-based transformer network can robustly recover the
complex topologies given by the Richtmyer-Meshkoff instability from a sequence
of hydrodynamic features derived from radiographic images corrupted with blur,
scatter, and noise. This approach is demonstrated on ICF-like double shell
hydrodynamic simulations. The key component of this network is a transformer
encoder that acts on a sequence of features extracted from noisy radiographs.
This encoder includes numerous self-attention layers that act to learn temporal
dependencies in the input sequences and increase the expressiveness of the
model. This approach is demonstrated to exhibit an excellent ability to
accurately recover the Richtmyer-Meshkov instability growth rates, even despite
the gas-metal interface being greatly obscured by radiographic noise.