Appearance of deterministic mixing behavior from ensembles of fluctuating hydrodynamics simulations of the Richtmyer-Meshkov instability

Kiran Narayanan, Ravi Samtaney
2018 Physical review. E  
We obtain numerical solutions of the two-fluid fluctuating compressible Navier-Stokes (FCNS) equations, which consistently account for thermal fluctuations from meso-to macroscales, in order to study the effect of such fluctuations on the mixing behavior in the Richtmyer-Meshkov instability (RMI). The numerical method used was successfully verified in two stages: for the deterministic fluxes by comparison against air-SF 6 RMI experiment, and for the stochastic terms by comparison against the
more » ... ect simulation Monte Carlo results for He-Ar RMI. We present results from fluctuating hydrodynamic RMI simulations for three He-Ar systems having length scales with decreasing order of magnitude that span from macroscopic to mesoscopic, with different levels of thermal fluctuations characterized by a nondimensional Boltzmann number (Bo). For a multidimensional FCNS system on a regular Cartesian grid, when using a discretization of a space-time stochastic flux Z(x,t) of the form Z(x,t) → 1/ √ h tN(ih,n t) for spatial interval h, time interval t, h, and Gaussian noise N should be greater than h 0 , with h 0 corresponding to a cell volume that contains a sufficient number of molecules of the fluid such that the fluctuations are physically meaningful and produce the right equilibrium spectrum. For the mesoscale RMI systems simulated, it was desirable to use a cell size smaller than this limit in order to resolve the viscous shock. This was achieved by using a modified regularization of the noise term via Z(x,t) → 1/ t max (h 3 ,h 3 0 )N (ih,n t), with h 0 = ξh ∀h < h 0 ,ξ ∈ Z + . Our simulations show that for systems with Bo 1 deterministic mixing behavior emerges as the ensemble-averaged behavior of several fluctuating instances, whereas when Bo ≈ 1, a deviation from deterministic behavior is observed. For all cases, the FCNS solution provides bounds on the growth rate of the amplitude of the mixing layer.
doi:10.1103/physreve.97.043111 pmid:29758761 fatcat:euwiwyg3dzee7bwjxnmuaky75u