Non-oscillatory forward-in-time integrators for viscous incompressible flows past a sphere

Joanna Szmelter, Piotr K. Smolarkiewicz, Zhao Zhang, Zhixin Cao

Research output: Contribution to journalArticlepeer-review

7 Citations (Scopus)
45 Downloads (Pure)


A non-oscillatory forward-in-time (NFT) integrator is developed to provide solutions of the Navier–Stokes equations for incompressible flows. Simulations of flows past a sphere are chosen as a benchmark representative of a class of engineering flows past obstacles. The methodology is further extended to moderate Reynolds number, stably stratified flows under gravity, for Froude numbers that typify the characteristic regimes of natural flows past distinct isolated features of topography in weather and climate models. The key elements of the proposed method consist of the Multidimensional Positive Definite Advection Transport Algorithm (MPDATA) and a robust non-symmetric Krylov-subspace elliptic solver. The solutions employ a finite volume spatial discretisation on unstructured and hybrid meshes and benefit from a collocated arrangement of all flow variables while being inherently stable. The development includes the implementation of viscous terms with the detached-eddy simulation (DES) approach employed for turbulent flows. Results demonstrate that the proposed methodology enables direct comparisons of the numerical solutions with corresponding laboratory studies of viscous and stratified flows while illustrating accuracy, robustness and flexibility of the NFT schemes. The presented simulations also offer a better insight into stably stratified flows past a sphere.
Original languageEnglish
Pages (from-to)365-383
Number of pages19
JournalJournal of Computational Physics
Early online date27 Feb 2019
Publication statusPublished - 1 Jun 2019


Dive into the research topics of 'Non-oscillatory forward-in-time integrators for viscous incompressible flows past a sphere'. Together they form a unique fingerprint.

Cite this