Influence of Density Ratios on Richtmyer–Meshkov Instability with Non-Equilibrium Effects in the Reshock Process

: The Richtmyer–Meshkov instability in a two-component system during the reshock process for various density ratios is studied through the discrete Boltzmann method. Detailed investigations are conducted on both hydrodynamic and thermodynamic non-equilibrium behaviors. Specifically, the analysis focuses on the density gradient, viscous stress tensor, heat flux strength, thermodynamic non-equilibrium intensity, and thermodynamic non-equilibrium area. It is interesting to observe the complex variations to non-equilibrium quantities with the changing shock front, rarefaction wave, transverse wave, and material interface. Physically, the non-equilibrium area is extended as the perturbed material interface grows after the passing of the shock wave or secondary impact. Moreover, the global non-equilibrium manifestation decreases when the transmitted shock front and transverse waves leave or when the reflected rarefaction wave weakens. Additionally, the global thermodynamic non-equilibrium effect is enhanced as the physical gradients or non-equilibrium area increase. Finally, the local non-equilibrium effect decreases when the fluid structure gradually disappears under the action of dissipation/diffusion.


Introduction
The Richtmyer-Meshkov (RM) [1,2] instability is induced by a shock wave impacting a wrinkled material interface between two fluids.It is mainly due to the baroclinic effect caused by the non-coincidence of the density gradient and pressure gradient across the material interface [3].The RM instability occurs frequently in nature and engineering applications, including supernova explosion [4], weapon implosion [5], aerospace [6], and inertial confinement fusion (ICF) [7,8].Particularly, in the ICF process, the development of the RM instability is enhanced by multiple interactions between the incident shock wave and the material interface, which thereby limit the energy output [9].In supersonic combustion engines, the RM instability accelerates contact between the fuel and oxidizer, thus enhancing the combustion rate [10].Indeed, understanding the physical phenomenon of the RM instability is of great significance in the engineering field.
Recently, the RM instability induced by the secondary impact has become an active research topic and has produced fruitful results in theoretical [11][12][13] and experimental [14][15][16][17] studies.In particular, numerical simulation has promoted relative research due to its convenience and efficiency, and it can enrich our understanding of the RM instability from multiple perspectives.These perspectives are mainly attributed to the following three levels: macroscopic [18][19][20][21][22][23][24], microscopic [25,26], and mesoscopic aspects [27][28][29].In the macroscopic view, various codes based on hydrodynamic governing equations, such as Euler or Navier-Stokes equations, have been widely used to simulate fluid instability.For example, in 2011, Ukai et al. [18] employed the compressible Navier-Stokes equations to study the growth of the RM instability during the reshock process for various interface shapes.In 2014, Tritschler et al. [19] explored the RM instability at the deterministic multimode interface and found that the vorticity deposition produced by the secondary impact material interface is greater than that of the primary impact.In 2019, Li et al. [21] studied the RM instability with reshock in the interaction between the mixing zone and multiple waves by using a direct numerical simulation and found that the mixing zone is accelerated by rarefaction and compression waves alternatively with decaying strength after reshock.In 2021, Bin et al. [23] successfully captured the growth of the RM instability mixed width under reshock and spatial structure in three dimensions by constrained largeeddy simulations.In 2022, Mohaghar et al. [24] performed three-dimensional simulation of the reshocked inclined RM instability turbulent mixing layer under single-mode and multimode configurations through Flash code and discovered that the growth of the mixing layer depends on the initial perturbation.Moreover, in the microscopic view, Wu et al. [25] utilized classical molecular dynamics to investigate the microscopic-scale RM instability of a single-mode Cu-He interface subjected to a cylindrically converging shock in 2018.Later, in 2020, Liu et al. [26] investigated the circulation deposition processes under different initial conditions by a direct simulation Monte Carlo method; they found the viscosity gradient in the region of interaction increases with an increase in the shock strength, and the incident shock Mach number determines the relative contribution of the viscous source.Although these microscopic approaches can provide insight into the RM instability, they are limited to small-scale problems due to the high cost of computational resources.
Fortunately, the mesoscopic method can provide a bridge between the macroscopic and microscopic approaches and can resolve the contradiction between physical accuracy and calculation speed.Recently, the discrete Boltzmann method (DBM) has become an important tool for studying complex fluid systems at the mesoscopic level.Based upon non-equilibrium statistical physics, the DBM not only recovers fluid equations at the hydrodynamic limit, it also provides thermodynamic non-equilibrium (TNE) information beyond that of traditional macroscopic models [30,31].In fact, the DBM is a variant version of the lattice Boltzmann method, which has been widely used as a solver of the Navier-Stokes equations and other partial differential equations [32][33][34][35][36][37][38].As a coarse-grained physical method, the DBM can capture abundant hydrodynamic and TNE behaviors and describe the evolution of sharp interfaces, small-scale structures, and rapidly changing modes in fluid systems.For instance, in 2014, Lin et al. [27] presented a polar-coordinate DBM for compressible flows, simulated the RM instability induced by a shock wave traveling outwards through a heavier (lighter) medium to a lighter (heavier) medium, and analyzed the non-equilibrium manifestations around the material and mechanical interfaces.In 2018, Chen et al. [28] used the multiple-relaxation-time DBM to investigate the non-equilibrium characteristics of a coexistence system involving both RM and Rayleigh-Taylor instability and found that the non-equilibrium effect of the system increases as the Mach number increases, while the effect of gravity acceleration on the non-equilibrium characteristics is different at various phases.In 2023, Shan et al. [29] employed the two-fluid DBM to investigate the non-equilibrium kinetic effects on the RM instability and reshock process and found that lighter fluid has a higher entropy production rate than heavier fluid.Moreover, the DBM also provides new insights for complex fluid systems, such as hydrodynamic instability [39][40][41][42][43], multiphase flow [44,45], reactive flow [46], etc.
It is noteworthy that the Atwood number has attracted great attention because it plays a crucial role in the RM instability.Early in 1960, under the assumptions of a small disturbance and incompressible fluid, Richtmyer proposed a linear growth model that revealed a positive correlation between the growth rate of a perturbation and the Atwood number [1].In 2011, Lombardini et al. [47] found that the evolution of mixed regions showed irregularities with an increasing Atwood number in the process of the RM instability.In 2016, Zhou et al. [48] investigated the RM instability under single and secondary shocks for various Atwood numbers and found that for the secondary shock, the larger the Atwood number, the smaller the normalized mixed mass.In 2019, Chen et al. [49] systematically analyzed the influence of the Atwood number on the RM instability for slightly perturbed interfaces in elastic-plastic media.In the same year, Liao et al. [50] studied the RM instability driven by a perturbed shock wave and found that the amplitude of the material interface decreases with the increasing Atwood number.In 2021, Tang et al. [51] adopted the Euler equations to study the single-mode RM instability induced by the convergent shock wave and observed that the amplitude perturbation grows more rapidly and the nonlinearity becomes more pronounced as the Atwood number increases.In 2022, Ren et al. [52] used the hydrocode Pagosa to explore the multimode RM instability with a variety of initial conditions, including differing the number of modes, Atwood number, and perturbation type.Up to now, the study of the Atwood number or density ratio has mainly involved the macroscopic aspect.In fact, in the RM instability process, the interaction between the shock wave and material interface also results in abundant and pronounced TNE behavior.Unfortunately, these behaviors are inappropriately ignored due to the limitations of computational models.
To further study the hydrodynamic non-equilibrium and TNE effects of the RM instability in the reshock process for various density ratios, the two-component DBM [53] is employed in this work.The structure of this paper is as follows: In Section 2, the DBM is briefly introduced.In Section 3, we simulate the RM instability and discuss its evolution mechanism in detail.Finally, a conclusion is given in Section 4.

Two-Component Discrete Boltzmann Model
The discrete Boltzmann equation for a two-component fluid system takes the following form [53]: where σ(= A, B) represents chemical species, i(= 0, 1, . . ., N) is the index of discrete velocities, f i σ ( f i σeq ) is the discrete (equilibrium) distribution function of each component, t is the time, v i is the discrete velocity, and r is the space coordinate.The collision term Ω σ i is calculated as where the relaxation time τ σ = (n A /θ A + n B /θ B ) −1 , and θ A and θ B are the relaxation parameters.The number density n σ , mass density ρ σ , and flow velocity u σ of each component are expressed by where m σ represents the molecular mass of σ.The mixing number density n, mass density ρ, and flow velocity u are The internal energy density of each component and the internal energy density of the system are where the symbol η i is used to describe the internal energy density for extra degrees of freedom.Let us introduce D and I, which represent the spatial dimension and extra degrees of freedom, respectively, and choose D = 2 and I = 3 in this paper.Then the individual temperature T σ and mixing temperature T can be calculated as Via the Chapman-Enskog expansion, it is proved that the Navier-Stokes equations can be recovered from Equation (1) in the continuum limit: see details in Appendix A. To this end, f i σeq should satisfy seven kinetic moments: see Appendix B. Those relationships can be written in the following matrix form: in terms of with an N × N matrix C = [C 1 , C 2 , . . ., C N ] T that links moment space to velocity space, and the elements within C are determined by the discrete velocity model.From Equation ( 13), the discrete equilibrium distribution function of each component can be obtained where C −1 is the inverse matrix of C.
It is noteworthy that the first three kinetic moments in Equations (A6)-(A8) are related to the conservation of mass, momentum, and energy, respectively.Consequently, f σeq i can be replaced by f σ i in those three formulas.However, in Equations (A9)-(A12), if f σeq i is replaced by f σ i , the values of the left and right sides may deviate.In fact, the deviation of the non-conservation moments of f σ i and f σeq i can be used to measure the degree of departure from the equilibrium state of the system, which is the core and purpose of the DBM.To be specific, the symbols M σ * m,n ( f σ i ) and M σ * m,n ( f σeq i ) are introduced to stand for the central kinetic moments of f σ i and f σeq i , respectively.Here "m, n" denotes that the m-order tensor is contracted to the n-order tensor.Both Then, the non-equilibrium quantities are defined as To better describe the global TNE of the fluid system, several TNE quantities are defined below: where ∆ σ * 2 denotes the non-organized momentum flux and is related to the viscosity, ∆ σ * 3,1 and ∆ σ * 3 represent the non-organized energy flux and are associated with the heat flux, and ∆ σ * 4,2 indicates the flux of the non-organized energy flux.By integrating several TNE quantities and then calculating the average over the whole computational region L x × L y , where L x and L y indicate the length and width of the fluid system, the average TNE intensities can be obtained as below: These average TNE intensities describe the overall non-equilibrium effects of the system from different perspectives.
In addition, for the discretization of the velocity space, we adopt the discrete velocity model D2V16 that takes the following form [46]: where v a , v b , v c , and v d are flexible parameters; and η a , η b , η c , η d are tunable parameters.Furthermore, the first-order forward Euler scheme is adopted for the time derivative, and the second-order nonoscillatory and nonfree-parameter dissipation difference scheme is adopted for the space derivatives [54].Additionally, all physical quantities are dimensionless in this paper [40].In fact, dimensionless quantities can be converted from/to dimensional physical quantities, which is a technique that is widely used in both the discrete Boltzmann method [40,55] and the lattice Boltzmann method [56][57][58][59].

Results and Discussions
In this section, let us perform two-dimensional simulations of the RM instability in the reshock process within a region L x × L y = 0.5 × 0.1.Figure 1 delineates the initial configuration of the RM instability as follows: where the subscripts L and M represent the left and middle regions, respectively, which are filled with species A, and R stands for the right area, which is full of B. The shock wave is initially located at x = 0.025 between the left and middle regions, and it propagates rightwards with a Mach number Ma = 1.2.(The current research is based on a relatively small Mach number, and other cases for larger Mach numbers may be considered in the future.)The physical quantities across the shock wave satisfy the Rankine-Hugoniot relations [60].Further, a perturbation in the material interface at around x = 0.125 separates the M part from the R part.To stimulate the growth of the RM instability, a sinusoidal perturbation function is applied to the material interface, y = 0.125 + A 0 cos(ky), in terms of the perturbation amplitude A 0 = 0.01 and the wave number k = 2π/L y .Moreover, to ensure the fluid system is in mechanical and thermal equilibrium initially, the pressures and temperatures of components A and B are the same alongside the material interface, and the concentrations of the two species are chosen as n A = n B = 1 here.Consequently, the density ratio between them is In addition, the Atwood number is defined as It can be found that there is a relationship between the density ratio and the Atwood number.In this paper, without loss of generality, we give a fixed m A = 1 and a tunable m B .Accordingly, the density ρ R in Equation ( 29) and the density ratio r in Equation ( 30) are flexible.In addition, we choose the relaxation parameters θ A = θ B = 2.0 × 10 −5 , and other parameters (v a , v b , v c , v d ) = (0.6, 1.6, 2.9, 5.9) and (η a , η b , η c , η d ) = (0.0, 2.9, 0.0, 0.0).The periodic boundary condition is used on the top and bottom, the inlet boundary condition is used on the left, and the symmetrical boundary condition is used on the right.

Verification and Validation
First of all, the grid-independent test is performed: see details in Appendix C. For the sake of effective and accurate simulation, we choose mesh grids of Nx × Ny = 2000 × 400, a space step ∆x = ∆y = 2.5 × 10 −4 , and a time step ∆t = 5.0 × 10 −6 .Then, to further validate the DBM for the simulation of the RM instability, a comparison is made between the DBM and the Zhang-Sohn model [61].Mathematically, the theoretical solution of the overall amplitude growth from the early to the late stage is predicted by the Zhang-Sohn model as follows [61]: where a 1 is the post-shock amplitude, A 1 is the post-shock Atwood number, and v 0 = k∆uA 1 a 1 , with ∆u denoting the change in the interface velocity caused by the shock wave.Figure 2 delineates the perturbation growth rates under the condition of ρ R = 2.5.The solid and dashed lines stand for the numerical and theoretical results, respectively.The shock front encounters the material interface at the time instant t = 0.0 and compresses it subsequently, so the simulated growth rate is less than zero in the early stage.After the time t = 0.08, the shock wave departs from the material interface, the perturbed interface increases, and its growth rate is greater than zero.As time goes on, the growth rate decreases gradually.From Figure 2, it can be seen that the DBM results agree with the theoretical solutions of the Zhang-Sohn model [61] for the evolution of the RM instability, indicating the reliability of the present numerical method.For an intuitive understanding of the evolution of the RM instability, Figure 3 illustrates the snapshots for the case r = 2.5 at different moments: the contours of the density, non-equilibrium strength, and vorticity are plotted in the three columns from left to right, respectively.Clearly, in the first column, the shock wave propagates from light to heavy fluids in the early state.At t = 0.1, the shock wave first passes through the material interface, producing a leftward reflected shock wave and a rightward transmitted shock wave.The transmitted shock wave reaches the right boundary then is reflected and propagates leftwards from t = 0.3∼0.7.At about t = 0.7, the secondary shock wave interacts with the material interface; then, the reflected rarefaction wave enters the heavier fluid and the transmitted shock wave enters the lighter fluid.Meanwhile, transverse waves are generated as a result of the impact of the shock wave on the wrinkled material interface and are continuously elongated, forming remarkable crisscrossing patterns.Afterward, the material interface starts to deform, and both the peak and valley shift into each other.This is also one of the classical phenomena of the RM instability.When t > 1.3, the two components penetrate each other, creating "bubble" and "spike" structures, and the "spike" structure gradually evolves into the "mushroom" structure.Finally, as the fluid is fully mixed, the fine fluid structure slowly disappears.These phenomena are qualitatively consistent with previous experimental and numerical studies [14,22].
Further, it is apparent in the second column of Figure 3 that the fine fluid structures can be vividly depicted by the non-equilibrium contours.Here, the non-equilibrium strength is calculated by summing the TNE intensity of each component, i.When the (reflected transmitted) shock wave passes through the material interface, it leads to a sudden surge in the TNE strength.Subsequently, the RM instability continues to evolve, the material interface is elongated, and the nonequilibrium region increases.As the two components are fully mixed, the material interface is gradually blurred, and the fluid system tends to equilibrium eventually.In addition, the last column in Figure 3 displays the vorticity in the evolution of the RM instability.Here, the vorticity is mathematically expressed by and physically it is closely related to the flow-field structure and can reflect the deformation of the material interface.In the initial stage, symmetry and a small vortex structure begin to appear around the material interface due to the velocity difference between two sides of the material interface.With the evolution of the RM instability, the vortex structure develops, which further promotes the deformation and elongation of the material interface.These simulation results are qualitatively consistent with previous studies [22,62].
To gain insight into the evolution of the RM instability with different density ratios, we show the case of r = 5.0 in Figure 4. Similar to the process in Figure 3, the linear, non-linear, deformation, and mixing stages of the RM stability can also be observed in Figure 4. Different from the case in Figure 3, here, the transmitted shock wave travels slower in component B in the period of t = 0.1∼0.8.As shown in Figure 4, after t > 1.0, the evolution speed of the material interface is accelerated significantly, leading to a more remarkable structural change and a more complex flow field.All the above numerical results agree with previous research [51,63].

Hydrodynamic Non-Equilibrium Effects
Next, let us study the hydrodynamic non-equilibrium behavior in the evolution of the RM instability.For this purpose, we introduce the average density gradients as below: where the integral is over the whole computational region.Physically, |∇ x ρ| and |∇ y ρ| stand for the average density gradients in the x and y directions, respectively.The average density gradient |∇ρ| increases with increasing |∇ x ρ| or |∇ y ρ|. Figure 5a-c show the evolution of the average density gradients |∇ x ρ|, |∇ y ρ|, and |∇ρ|, respectively.The lines with different symbols denote the results for various density ratios, i.e., r = 2.5, 3.0, . . ., 6.0, respectively.It can be found in Figure 5a that the tendency of |∇ x ρ| is similar in each case.Take the case with r = 5.0: the shock wave impacts the material interface from t = 0.04 to 0.1, leading to a rapid upward trend of |∇ x ρ|.Then, |∇ x ρ| continues to rise and reaches the peak for the first time at t = 1.2, when the density gradient across the material interface, shock wave, rarefaction wave, and transverse waves all together contribute to the maximum of |∇ x ρ|.From t = 1.2 to 1.5, |∇ x ρ| shows a downward trend because the transmitted shock wave and transverse waves begin to leave the computational domain, although the material interface undergoes a reversal.Subsequently, from t = 1.5 to 2.4, the RM instability continues to evolve, the material interface is elongated, the physical field becomes complex, and |∇ x ρ| rises again.Finally (t > 2.4), the material interface gradually becomes blurred due to the effects of dissipation/diffusion; the physical field becomes smoother, leading to a decrease in |∇ x ρ|.From Figure 5b, it can be observed that the evolution of the average density gradient in the y direction |∇ y ρ| is different from that of |∇ x ρ|.Before the time t = 1.0, under the influence of the RM instability, the sum of the material surface and mechanical interface increases in the y direction; hence, the global density gradient in the y direction rises.Subsequently, |∇ y ρ| decreases due to the interface reversal when the length of the contact interface decreases.During the stage of t = 1.1∼2.7, the length of the material interface increases because of the RM instability, so |∇ y ρ| rises again.Finally (t > 2.7), as the fluids are fully mixed, the fine structure gradually disappears, the system tends to reach equilibrium, and |∇ y ρ| decreases.
Moreover, as shown in Figure 5c, it is easy to understand the tendency of |∇ρ| through analysis of the variation to |∇ x ρ| and |∇ y ρ|.Moreover, to have a deeper insight into the law of density gradients, the second column of Figure 5 shows the results of t s 0 |∇ x ρ|dt/t s , t s 0 |∇ y ρ|dt/t s , and t s 0 |∇ρ|dt/t s , where the integral is conducted over the period from t = 0 to t s , and the final instant is chosen at t s = 4 in this paper.Physically, the three quantities reflect the time-integration results of |∇ x ρ|, |∇ y ρ|, and |∇ρ|, respectively.The fitting functions are 0 |∇ y ρ|dt/t s ) show a linear relationship with the density ratio.Physically, as the density ratio increases, the material interface evolves faster and the fluid structure changes are more complex, resulting in a more significant change to the system density gradients.

Thermodynamic Non-Equilibrium Effects
In this subsection, let us study the essential TNE behaviors in the process of the RM instability.First of all, the viscous stress tensor D σ * 2 is explored in Figure 6.Subsequently, the average heat flux strength D σ * 3,1 is investigated in Figures 7 and 8.Then, the average TNE strength D σ * is probed in Figure 9. Lastly, the proportion of the non-equilibrium region Sr σ is discussed in Figure 10.3,1x shows a trend of increasing, afterward decreasing, then increasing, and finally decreasing overall in each case.Let us take the case of r = 5.0 for example.To the material interface is impacted by the shock wave before t = 0.1 and by the secondary shock wave during the phase of t = 0.9∼1.18.During the two stages, the mixing of two components accelerates, which promotes the heat exchange simultaneously, resulting in a rapid increase in D A * 3,1x .During the period t = 0.1∼0.5, the transmitted shock wave propagates into the heavier fluid and the change in D A * 3,1x is not obvious.For the phase of t = 0.5∼0.9, the transmitted shock wave is reflected back and moves leftwards, and D A * 3,1x rises gradually.Afterwards, as the transmitted shock wave departs from the physical system, the value of

n s i t y r a t i o
increases because the vortex structure on the interface becomes more complex, the contact area of the fluid increases, and the heat exchange of the fluid system increases remarkably.After t = 2.3, the two components are fully mixed, the fluid structure gradually disappears under the action of dissipation/diffusion, and D A * 3,1x shows a gradual declining trend.Figure 7b shows the evolution of the heat flux intensity D A * 3,1y over time.It can be seen that for all cases, D A * 3,1y firstly rises, afterward falls, then increases, and eventually decreases.For the specific case of r = 5.0, before t = 0.9, the material interface is widened in the x direction continuously and the heat exchange in the y direction is enhanced due to the growth of the RM instability; hence, D A * 3,1y increases.During t = 0.9∼1.1, the reshock wave impacts the material interface; the material interface is compressed, resulting in a slight decrease in D A * 3,1y .Afterwards, as the RM instability evolves, the material interface is elongated constantly, the contact area between the two components increases, and the heat exchange is enhanced.Finally (t > 3.0), under the action of dissipation and heat conduction, the fluid structure slowly disappears, the mixing degree of the two components is deeper, the physical field becomes smooth, and thereby, D A * 3,1y declines.The evolution of the heat flux strength 3,1 first shows a rising and then a falling tendency in each case on the whole.That is to say, the heat exchange of the system is elevated in the early stage and then diminished in the evolution of the RM instability.In fact, the variation of D Next, to gain a deeper understanding of the TNE behavior in the RM instability system, the average TNE strength D σ * is probed below.As seen in Figure 9a, for all examples, the profile of D A * over time has a similar pattern and is positively related to the density ratio.Specifically, take r = 5.0 as an example.The TNE intensity of component A shows a fluctuating increase before t = 1.2 due to the impact of the shock wave on the material interface and the increase of the non-equilibrium area around the transmitted shock wave, reflected shock or rarefaction wave, and transverse waves.During t = 1.2∼1.5, the transverse waves and transmitted shock wave leave the computational domain, which causes a decrease in the physical gradient of the system.And when t > 1.5, D A * increases again and decreases eventually.This phenomenon results from competition between the following two mechanisms.On the one hand, the RM instability induced by the secondary impact continues to evolve, the width of the material interface increases, and the non-equilibrium region also increases: thereby, the TNE intensity is enhanced.On the other hand, under the dissipation/diffusion effect, the vortex structure gradually disappears and the physical field becomes smooth.The former (latter) plays a prominent role in the increase (decrease) stage.
In addition, Figure 9b shows the evolution of D B * for various density ratios.For more detailed analysis of the evolution trend of D B * , take r = 5 as an example.In the very early phase t < 0.1, the shock wave rapidly passes through the material interface and enters the heavier medium.Sharp physical gradients start to emerge in the B medium, causing D B * to surge rapidly.During the stage of t = 0.1∼0.9,D B * rises with a small disturbance as the transmitted shock wave continues to propagate rightwards until it reaches the right boundary, at which point the transmitted shock wave is reflected and then propagates leftwards.Subsequently, the physical gradients decrease at t = 0.9∼1.1 as the transmitted shock wave passes through the material interface and leaves component B, and a reflected rarefaction wave is generated in the right region and then dissipates gradually.When t > 1.1, it can be observed that D B * shows a trend of increasing and then decreasing.The reason for the rise is that the material interface is reversed and elongated, the physical gradients change to become more complex, and the TNE is enhanced.Meanwhile, the decrease in the TNE strength D B * results from the fact that the physical gradients become smooth under the effect of dissipation/diffusion.Moreover, Figure 9c,d portray the relationship between the time integration of D σ * and the density ratios.The fitting functions are , where the integral is performed from t = 0 to t s = 4. Physically, for a larger density ratio, the material interface becomes more complex and the increase to the physical gradient is larger, leading to the growth of the TNE intensity of each component.
Finally, for a more detailed description of the evolution law of the global TNE intensity, let us introduce the proportion of the non-equilibrium region Sr σ , which is defined as the ratio of the non-equilibrium area to the total area of the fluid system [41].Here, the non-equilibrium region is where the TNE intensity of species σ is greater than a threshold θ.In this work, the threshold value is chosen as θ = 0.02.
Figure 10a shows that for any case, Sr A firstly rises, afterwards falls, then rises again, and finally falls.Take r = 5.0 as an example: it can be seen that there is an increasing trend before t = 1.1.The rise is attributed to the impact of the shock wave on the material interface as well as the increase in physical gradients near the transmitted shock wave, reflected shock or rarefaction wave, and transverse waves, enhancing the TNE intensity.From t = 1.1 to 1.3, the TNE intensity is weakened as the transverse waves and transmitted shock wave depart the calculation region.During the period of t = 1.3∼2.5, the RM instability continues to evolve, the interface is continuously elongated, and the non-equilibrium region increases, and so Sr A rises slowly.When t > 2.5, the TNE intensity is decreased as the fine structure gradually disappears under the effect of dissipation/diffusion.
In the same way, Figure 10b describes the evolution of Sr B for various density ratios r.The evolutionary trend of Sr B is similar to that of Sr A , and the physical mechanisms can be obtained by the analysis of Sr A .Further, the time integration results of Sr σ for various density ratios are depicted in Figure 10c,d.Their relationships are t s 0 Sr A dt/t s = −0.10exp(−0.27r)+ 0.09 and t s 0 Sr B dt/t s = −0.12exp(−0.67r)+ 0.06, respectively.As a result, for a larger density ratio, the interface evolves faster and the macroscopic physical gradient changes are more significant, and thus, the TNE becomes stronger.

Conclusions
In this paper, the RM instability in the reshock process for various density ratios is investigated through the two-component DBM.Both hydrodynamic and TNE effects are analyzed in detail.On the whole, the average density gradients demonstrate a trend of rising, falling, rising again, and then falling eventually.A linearly increasing relationship is obtained between the time integration of the average density gradient and the density ratio.The main physical mechanisms are as follows.Firstly, the areas of density gradients are extended as the mixing of the two components accelerates under the action of the shock wave.Secondly, the global physical gradients decrease while the transmitted shock wave and transverse waves leave the calculational domain.Thirdly, the local macroscopic physical gradients decrease when the fluid structure gradually disappears due to dissipation/diffusion.Finally, the TNE effects are studied from various perspectives, including the average viscous stress tensor, heat flux intensity, TNE intensity, and proportion of the non-equilibrium region.It is found that the time integration of the TNE quantities and the proportion of the non-equilibrium region of each species are higher for a larger density ratio.However, the heat flux intensity of the heavier medium shows the opposite trend.Physically, there are several competitive mechanisms associated with these TNE quantities.(i) The TNE strength is significantly enhanced under the impact of the shock wave.(ii) The non-equilibrium region increases with the elongation of the material interface.(iii) The global non-equilibrium manifestation decreases as the transmitted shock wave and transverse waves leave or when the reflected rarefaction wave weakens.(iv) The global TNE intensity is enhanced (weakened) with an increase (decrease) in the macroscopic physical gradients.These findings can deepen our understanding of the RM instability with non-equilibrium effects during the reshock process from mesoscopic to macroscopic levels and are useful for understanding and mitigating the influence of hydrodynamic instabilities on indirect-drive ICF implosions [7,[64][65][66].Meanwhile, the research can also contribute to the knowledge of the interaction mechanism between shock waves and material interfaces, explosions, supersonic combustors [5,10] where p σ = n σ T σ , e σ = (D + I)T σ /(2m σ ), µ σ = p σ τ σ , κ = γµ σ , γ = (D + I + 2)/(D + I), and e = ∑ σ ρ σ (e σ + u σ2 /2)/ρ − u 2 /2 denote the pressure, internal energy per unit mass, dynamic viscosity coefficient, heat conductivity, specific heat ratio, and internal energy of the mixing system per unit mass, respectively.The term δ αβ represents the Kronecker function, and the subscripts α, β, χ stand for x or y.

Appendix B. Expressions of the Kinetic Moments
The equilibrium distribution function f σeq i satisfies the following seven kinetic moments: M 2,0 f m σ e α e β δ αβ + u • u , (A9)

Appendix C. The Verification of Grid Convergence
To ensure an accurate and effective simulation of the RM instability, the grid independence is verified.Figure A1 depicts the profile of the density along the center line L y /2 in the case of r = 2.5 at the time instant t = 2.0.We choose four different grid sizes: N x × N y = 250 × 50, 500 × 100, 1000 × 200, and 2000 × 400.It can be found that the simulation results gradually converge as the number of mesh grids increases.To be specific, the simulation results for 1000 × 200 and 2000 × 400 are close to each other.This means that the simulation accuracy is high enough for 2000 × 400.Consequently, a grid size of N x × N y = 2000 × 400 with a corresponding spatial step of ∆x = ∆y = 2.5 × 10 −4 is adopted for all simulations in this paper.

Figure 1 .
Figure 1.The initial configuration of the RM instability.

Figure 2 .
Figure 2.Comparison of the growth rates of the RM instabilities given by the DBM and the Zhang-Sohn model[61].

Figure 3 .
Figure 3.The first to third columns show contours of the density, non-equilibrium strength and vorticity in the RM process in the case of r = 2.5, respectively.

Figure 4 .
Figure 4.The first to third columns show contours of the density, non-equilibrium strength, and vorticity, respectively, for the RM process in the case of r = 5.0.

Figure 5 .
Figure 5. (a) Average density gradient in the x direction versus time.(b) Average density gradient in the y direction versus time.(c) Average density gradient versus time.(d) t s 0 |∇ x ρ|dt/t s versus density ratios.(e) t s 0 |∇ y ρ|dt/t s versus density ratios.(f) t s 0 |∇ρ|dt/t s versus density ratios.

Figure 6 .2
Figure 6.(a) Viscous stress tensor of component A versus time.(b) Viscous stress tensor of component B versus time.(c) t s 0 D A * 2 dt/t s versus density ratios.(d) t s 0 D B * 2 dt/t s versus density ratios.

Figure 6a delineates theFigure 7 .22
Figure 6a delineates the evolution of the viscous stress tensor D A * 2 with different density ratios.It can be found that the whole level of D A * 2 is higher for a larger density ratio.To better analyze the evolutionary mechanism of the RM instability, let us take r = 5.0 for example.Clearly, D A * 2 decreases before t = 0.1 as the shock wave passes through the material interface.Then, D A * 2 goes up during t = 0.1∼0.2owing to the RM instability when the amplitude of the material interface rises.Subsequently, it goes down until about t = 0.4 because the reflected shock wave leaves the left boundary.Afterwards, there are few changes to D A * 2 during 0.4 < t ≤ 0.9 while the RM instability has a weak effect on D A * 2 .Later, at 0.9 < t ≤ 1.2, D A * 2 soars sharply because of the secondary impact of the reflected transmitted shock wave.During t = 1.2∼1.5,D A * 2 decreases with the departure of the transmitted shock wave and transverse waves.After that, the material interface
s i t y r a t i o

d e n s i t y r a t iFigure 9 .
Figure 9. (a) Average TNE strength of component A versus time.(b) Average TNE strength of component B versus time.(c) t s 0 D A * dt/t s versus density ratios.(d) t s 0 D B * dt/t s versus density ratios.

Figure 10 .00
Figure 10.(a) The proportion of the non-equilibrium region of component A versus time.(b) The proportion of the non-equilibrium region of component B versus time.(c) t s 0 Sr A dt/t s versus density ratios.(d) t s 0 Sr B dt/t s versus density ratios.Next, we investigate the heat flux intensity in the RM process.Figure 7 exhibits the evolution of the heat flux intensities D A * 3,1x , D A * 3,1y , and D A * 3,1 over time.It is observed in Figure 7a that D A * in Figure7b.To have a quantitative study, Figure7ddepicts the time integration of D A * 3,1 versus the density ratio.The squares denote the numerical results and the line represents the fitting function dt/t s = [−2.85exp(−0.18r)+ 3.00] × 10 −3 .It is evident that the larger the density ratio is, the faster the material interface evolves and the more significant the temperature gradient changes are.Subsequently, the evolutions of the heat flux strengths D various density ratios are discussed: see Figure8a-c.Their overall evolutionary trends are similar to that of component A. It is noteworthy that the heat flux strength of component B is negatively correlated with the density ratio.Specifically, for the heavier media, as r increases, the corresponding molecular mass becomes larger and the heat exchange of component B becomes weaker.It is further found that the relationship between dt/t s = [4.21exp(−0.26r)+ 2.49] × 10 −4 .Obviously, with the increasing density ratio, D B * 3,1x decreases monotonically.

Figure A1 .
Figure A1.The verification of grid convergence: the density profiles are plotted along the center line L y /2 at the time instant t = 2.0 for the RM process.