Direct numerical simulation (DNS) of Rayleigh-Taylor instability (RTI) has been reported in Sengupta, et al. (2022) [35] using more than 4 billion points with initial time steps of 7.69 ×10−8s. The set-up consists of quiescent heavy (cold) fluid atop lighter (hot) fluid initially separated by an insulated partition. Removal of this partition creates acoustic pulses in the directions normal to the interface, which provides the receptivity route of RTI to transient acoustic pulses spread over several mega-Hz frequencies, far in excess of ultrasonic frequencies. The DNS and its analysis reveal the pulses to be severely attenuated, which cannot be explained by the classical wave equation. In the present research, after demonstrating the features of the DNS results, we report in detail, the theoretical development of a wave equation incorporating losses based on the Navier-Stokes equation without Stokes' hypothesis for a quiescent ambience. Apart from the physical properties of this altered wave equation, the numerical solution of the same is obtained. The governing partial differential equation (PDE) for the propagation of disturbances in the spectral plane, provides the dispersion relation between wavenumber and circular frequency in the dissipative medium accounting for viscous losses. The presented analysis provides physical properties using the global spectral analysis (GSA). This shows the far-field perturbation to propagate either as attenuated waves or strictly in a diffusive manner depending upon the wavenumber. The PDE is solved numerically by a high-accuracy compact scheme and the four-stage, Runge-Kutta scheme for time advancement. The computed solution is shown to match not only with the developed theoretical analysis but also explains the DNS results for RTI at early times when the computed flow field truly represents the quiescent ambience.
Read full abstract