A study on temporal accuracy of OpenFOAM

  • Lee, Sang Bong (Department of Naval Architecture and Offshore Engineering, Dong-A University)
  • Received : 2016.08.01
  • Accepted : 2016.11.26
  • Published : 2017.07.31


Cranke-Nicolson scheme in native OpenFOAM source libraries was not able to provide 2nd order temporal accuracy of velocity and pressure since the volume flux of convective nonlinear terms was 1st accurate in time. In the present study the simplest way of getting the volume flux with 2nd order accuracy was proposed by using old fluxes. A possible numerical instability originated from an explicit estimation of volume fluxes could be handled by introducing a weighting factor which was determined by observing the ratio of the finally corrected volume flux to the intermediate volume flux at the previous step. The new calculation of volume fluxes was able to provide temporally accurate velocity and pressure with 2nd order. The improvement of temporal accuracy was validated by performing numerical simulations of 2D Taylor-Green vortex of which an exact solution was known and 2D vortex shedding from a circular cylinder.


Supported by : National Research Foundation of Korea (NRF)


  1. Cao, S., Ge, Y., Tamura, Y., 2010. Mechanism of the lift force on the circular and square cylinders in shear flows. In: 5th International Symposium on Computational Wind Engineering, May 23-27, North Carolina, USA.
  2. Kim, K., Baek, S.-J., Sung, H.J., 2002. An implicit velocity decoupling procedure for the incompressible Navier-Stokes equations. Int. J. Numer. Methods Fluids 38, 125-138.
  3. Lee, H.B., Rhee, S.H., 2015. A dynamic interface compression method for VOF simulations of high-speed planing watercraft. J. Mech. Sci. Technol. 29, 1849-1857.
  4. Lee, S.B., 2014. Application of OpenFOAM to prediction of hull resistance. In: 9th International OpenFOAM Workshop, June 23-26, Zagreb, Croatia.
  5. Leonard, B.P., 1979. A stable and accurate convective modelling procedure based on quadratic upstream interpolation. Comput. Methods Appl. Mech. Eng. 19 (1), 59-98.
  6. OpenFOAM User Guide by OpenFOAM Foundation Ltd. (founded in
  7. Park, S., Rhee, S.H., 2012. Computational analysis of turbulent supercavitating flow around a two-dimensional wedge-shaped cavitator geometry. Comput. Fluids 70, 73-85.
  8. Park, S., Rhee, S.H., 2013. Numerical analysis of the three-dimensional cloud cavitating flow around a twisted hydrofoil. Fluid Dyn. Res. 45, 1-20.
  9. Rajani, B.N., Kandasamy, A., Majumdar, S., 2009. Numerical simulation of laminar flow past a circular cylinder. Appl. Math. Model. 33, 1228-1247.
  10. van Leer, B., 1979. Towards the ultimate conservative difference scheme, V. A second order sequel to Godunov's method. J. Comput. Phys. 32, 101-136.
  11. Williamson, C.H.K., 1989. Oblique and parallel modes of vortex shedding in the wake of a circular cylinder at low Reynolds numbers. J. Fluid Mech. 206, 579-627.