Effect of gas viscosity on the interfacial instability development in a two-phase mixing layer

TL;DR

This study investigates how gas viscosity influences interfacial instability in two-phase mixing layers, revealing a viscosity-dependent transition from absolute to convective instability through simulations and stability analysis.

Contribution

It provides new insights into the Re-dependent transition between absolute and convective regimes, supported by both linear stability analysis and interface-resolved simulations.

Findings

Higher Re reduces temporal growth rate of instability.

Reynolds number threshold triggers absolute to convective transition.

Wave speed varies spatially and aligns with theoretical modes.

Abstract

The interfacial instability in a two-phase mixing layers between parallel gas and liquid streams is important to two-phase atomization. Depending on the inflow conditions and fluid properties, interfacial instability can be convective or absolute. The goal of the present study is to investigate the impact of gas viscosity on the interfacial instability. Both interface-resolved simulations and linear stability analysis (LSA) have been conducted. In LSA, the Orr-Sommerfeld equation is solved to analyze the spatio-temporal viscous modes. When the gas viscosity decreases, the Reynold number ($\text{Re}$) increases accordingly. The LSA demonstrates that when $\text{Re}$ is higher than a critical threshold, the instability transitions from the absolute to the convective (A/C) regimes. Such a $\text{Re}$-induced A/C transition is also observed in the numerical simulations, though the critical Re observed in simulations is significantly lower than that predicted by LSA. The LSA results indicate that the temporal growth rate decreases with Re. When the growth rate reaches zero, the A/C transition will occur. The $\text{Re}$-induced A/C transition is observed in both confined and unconfined mixing layers and also in cases with low and high gas-to-liquid density ratios. In the transition from typical absolute and convective regimes, a weak absolute regime is identified in the simulations, for which \tcr{the spectrograms} show both the absolute and convective modes. The dominant frequency in the weak absolute regime can be influenced by the perturbation introduced at the inlet. The simulation results also show that the wave propagation speed can vary in space. In the absolute instability regime, the wave propagation speed agrees well with the absolute mode celerity near the inlet and increases to the Dimotakis speed further downstream.