A numerical study of stably stratified flows past spheres at moderate Reynolds numbers is presented. The resolved flows can adequately describe a wide class of geophysical, environmental, and engineering flows characterised by the density stratification of the terrestrial atmosphere and oceanic thermocline. The range of physical phenomena developing when stratified flows impact single and multiple spheres constitute a convenient benchmark for complex geometry applications, e.g. mountains, islands, wind turbines, and buildings. Solutions of Navier-Stokes equations, in the incompressible Boussinesq limit, are obtained by applying a semi-implicit finite volume (FV) non-oscillatory forward-in-time (NFT) integration scheme enhanced by MPI parallelization. The developed model is applied for a systematic investigation of stratified flow patterns arising for a range of Froude numbers Fr ∈ [0.1,∞] at Reynolds numbers Re = 200 and Re = 300, for which the neutrally stratified flows induces distinctly different near-wake features
