Mathematical modeling of the incompressible fluid flows

Abstract

The aim of the present paper is the demonstration of the opportunities of the mathematical modeling of the separated flows of the river and the sea water around blunt bodies on the basis of the Navier-Stokes equations (NSE) in the Boussinesq approximation. In the river water only the wakes of the obstacles are observed. In the sea water along with wakes the internal waves (IWs) are generated. The 3D homogeneous (river) and density stratified (sea) incompressible viscous fluid flows around a sphere and a square cylinder have been investigated by means of the direct numerical simulation on supercomputers and the visualization of the 3D vortex structures in the wake. For solving of NSE the Splitting on physical factors Method for Incompressible Fluid flows (SMIF) with hybrid explicit finite difference scheme (second-order accuracy in space, minimum scheme viscosity and dispersion, capable for work in wide range of the Reynolds (Re) and the internal Froude (Fr) numbers and monotonous) has been developed and successfully applied. The different transitions in sphere wakes (2D-3D, laminar-turbulent and others) with increasing of Re (1 < Re < 500000) and decreasing of Fr (0.005 < Fr < 100) have been investigated in details. Thus the classifications of the viscous fluid flow regimes around a sphere have been refined. The original classification of the 2D stratified viscous fluid flow regimes around a square cylinder has been obtained at Re < 200, 0.1 < Fr < 100. At Fr = 0.1, Re = 50 the formation process of the two symmetric wavy hanging (soaring) sheets of density (starting from two hanging vortices) is demonstrated. It was found also that the values of the maximum phase difference along the circular cylinder axis are approximately equal to 0.1-0.2 Tf (for mode A, 191 < Re ≤ 300) and 0.015-0.030 Tf (for mode B, 300 ≤ Re ≤ 400), where the time Tf is the period of the flow.