In this study we have developed a flexible and efficient numerical scheme for the simulation of three–dimensional incompressible flows in spherical coordinates. The main idea, inspired by a similar strategy as [1] for cylindrical coordinates, consists of a change of variables combined with a discretization on a staggered mesh and the special treatment of few discrete terms that remove the singularities of the Navier–Stokes equations at the sphere centre and along the polar axis. This new method alleviates also the time step restrictions introduced by the discretization around the polar axis while it still suffers from strong limitations if convection at the sphere centre dominates the flow. The scheme is second–order accurate in space and is verified and validated by computing numerical examples that are compared with similar results produced by other codes or available from the literature. The method can cope with flows evolving in the whole sphere, in a spherical shell and in a sector without any change and, thanks to the flexibility of finite–differences, it can employ generic mesh stretching (in two of the three directions) and complex boundary conditions.

A finite–difference scheme for three–dimensional incompressible flows in spherical coordinates

Santelli, L.;
2021-01-01

Abstract

In this study we have developed a flexible and efficient numerical scheme for the simulation of three–dimensional incompressible flows in spherical coordinates. The main idea, inspired by a similar strategy as [1] for cylindrical coordinates, consists of a change of variables combined with a discretization on a staggered mesh and the special treatment of few discrete terms that remove the singularities of the Navier–Stokes equations at the sphere centre and along the polar axis. This new method alleviates also the time step restrictions introduced by the discretization around the polar axis while it still suffers from strong limitations if convection at the sphere centre dominates the flow. The scheme is second–order accurate in space and is verified and validated by computing numerical examples that are compared with similar results produced by other codes or available from the literature. The method can cope with flows evolving in the whole sphere, in a spherical shell and in a sector without any change and, thanks to the flexibility of finite–differences, it can employ generic mesh stretching (in two of the three directions) and complex boundary conditions.
2021
Spherical coordinates, Singularity at the centre, Singularity at the polar axis
File in questo prodotto:
File Dimensione Formato  
2021_JComputPhys_424_ 109848_Santelli.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non pubblico
Dimensione 7.31 MB
Formato Adobe PDF
7.31 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.12571/34704
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 12
  • ???jsp.display-item.citation.isi??? 9
social impact