Some aspects in Kelvin-Helmholtz instability with and without Boussinesq approximation

The topic of this paper is the Kelvin-Helmholtz instability, a phenomenon which occurs on the interface of a stratified fluid, in the presence of a parallel shear flow, when there is a velocity and density difference across the interface of two adjacent layers. This paper focuses on a numerical simulation modelled by the Taylor-Goldstein equation, which represents a more realistic case compared to the basic Kelvin-Helmholtz shear flow. The Euler system is solved with new modelled smooth velocity and density profiles at the interface. The flux at cell boundaries is reconstructed by implementing a third order WENO (Weighted Essentially Non-Oscillatory) method. Next, a Riemann solver builds the fluxes at cell interfaces. The use of both Rusanov and HLLC solvers is investigated. Temporal discretization is done by applying the second order TVD (total variation diminishing) Runge-Kutta method on a uniform grid. Numerical simulations are performed comparatively for both Kelvin-Helmholtz and Taylor-Goldstein instabilities, on the same simulation domains. We find that increasing the number of grid points leads to a better accuracy in shear layer vortices visualization. Thus, we can conclude that applying the Taylor-Goldstein equation improves the realism in the general fluid instability modelling.

New Closed-Form Thermal Boundary Layer Solutions in Shear Flow With Power-Law Velocity

The forced convection problem for a developing thermal boundary layer in a parallel shear flow is studied. If the shear flow has a power-law velocity profile, exact similarity thermal boundary layer solutions in terms of Gamma functions can be found. Specifically, three types of thermal boundary conditions are considered: a step temperature change, a step flux change, and a concentrated heat source. The latter is also analogous to mass diffusion form an isolated source. The mixing index for mass diffusion is found exactly.