Checkerboard pressure is one of the most significant problems arising from the use of col- located grids for fluid dynamic simulations using the finite-volume method. The original Rhie–Chow momentum interpolation technique, termed Original Momentum Interpola- tion Method (OMIM) was proposed to eliminate the non-physical saw-tooth pressure os- cillations. However, it was soon proved that the steady-state solutions obtained with this technique were under-relaxation factor dependent. Nevertheless, standard OMIM is still R commonly used in several CFD codes like OpenFOAM, the software used for this work. In this paper the OMIM and a possible correction for under-relaxation dependency, which R has been implemented in OpenFOAM ⃝ , are discussed in detail. The proposed methodology is compared and contrasted with OMIM in terms of accuracy of the solution and speed of convergence for several classical pressure–velocity segregated algorithms for steady state solvers; namely SIMPLE, SIMPLE-C, SIMPLE-R and PISO. A classical laminar 2D cavity is used as the base test-case. The study is then extended to a more complex 2D airfoil profile (NACA0012). In the cases considered mesh uniformity and orthogonality are progressively reduced and turbulence starts playing an important role, limiting therefore the conver- gence of the cases and the performance of the correction.

Influence of momentum interpolation methods on the accuracy and convergence of pressure–velocity coupling algorithms in OpenFOAM

PISCAGLIA, FEDERICO;MONTORFANO, ANDREA;ONORATI, ANGELO;
2017

Abstract

Checkerboard pressure is one of the most significant problems arising from the use of col- located grids for fluid dynamic simulations using the finite-volume method. The original Rhie–Chow momentum interpolation technique, termed Original Momentum Interpola- tion Method (OMIM) was proposed to eliminate the non-physical saw-tooth pressure os- cillations. However, it was soon proved that the steady-state solutions obtained with this technique were under-relaxation factor dependent. Nevertheless, standard OMIM is still R commonly used in several CFD codes like OpenFOAM, the software used for this work. In this paper the OMIM and a possible correction for under-relaxation dependency, which R has been implemented in OpenFOAM ⃝ , are discussed in detail. The proposed methodology is compared and contrasted with OMIM in terms of accuracy of the solution and speed of convergence for several classical pressure–velocity segregated algorithms for steady state solvers; namely SIMPLE, SIMPLE-C, SIMPLE-R and PISO. A classical laminar 2D cavity is used as the base test-case. The study is then extended to a more complex 2D airfoil profile (NACA0012). In the cases considered mesh uniformity and orthogonality are progressively reduced and turbulence starts playing an important role, limiting therefore the conver- gence of the cases and the performance of the correction.
File in questo prodotto:
File Dimensione Formato  
cam2016.pdf

Accesso riservato

Descrizione: articolo
: Publisher’s version
Dimensione 1.34 MB
Formato Adobe PDF
1.34 MB Adobe PDF   Visualizza/Apri
MARTJ_OA_01-17.pdf

embargo fino al 20/04/2018

Descrizione: Paper open access
: Post-Print (DRAFT o Author’s Accepted Manuscript-AAM)
Dimensione 1.45 MB
Formato Adobe PDF
1.45 MB Adobe PDF Visualizza/Apri

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/11311/1010574
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 13
  • ???jsp.display-item.citation.isi??? 9
social impact