On the development of OpenFOAM solvers based on explicit and implicit high-order Runge-Kutta schemes for incompressible flows with heat transfer

Abstract Open-source CFD codes provide suitable environments for implementing and testing low-dissipative algorithms typically used to simulate turbulence. In this research work we developed CFD solvers for incompressible flows based on high-order explicit and diagonally implicit Runge–Kutta (RK) schemes for time integration. In particular, an iterated PISO-like procedure based on Rhie–Chow correction was used to handle pressure–velocity coupling within each implicit RK stage. For the explicit approach, a projected scheme was used to avoid the “checker-board” effect. The above-mentioned approaches were also extended to flow problems involving heat transfer. It is worth noting that the numerical technology available in the OpenFOAM library was used for space discretization. In this work, we additionally explore the reliability and effectiveness of the proposed implementations by computing several unsteady flow benchmarks; we also show that the numerical diffusion due to the time integration approach is completely canceled using the solution techniques proposed here.

[1]  L. Schouveiler,et al.  Self-sustained oscillations in the wake of a sphere , 2002 .

[2]  Vinayak Eswaran,et al.  HEAT AND FLUID FLOW ACROSS A SQUARE CYLINDER IN THE TWO-DIMENSIONAL LAMINAR FLOW REGIME , 2004 .

[3]  Hassan M. Badr,et al.  Forced convection from a rotationally oscillating cylinder placed in a uniform stream , 2000 .

[4]  Jungwoo Kim,et al.  An immersed-boundary finite-volume method for simulations of flow in complex geometries , 2001 .

[5]  D. Lilly,et al.  A proposed modification of the Germano subgrid‐scale closure method , 1992 .

[6]  William W. Willmarth,et al.  Some experimental results on sphere and disk drag , 1971 .

[7]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[8]  Chun Shen,et al.  Implementation of density-based solver for all speeds in the framework of OpenFOAM , 2014, Comput. Phys. Commun..

[9]  R. Chhabra,et al.  Two-dimensional unsteady forced convection heat transfer in power-law fluids from a cylinder , 2010 .

[10]  E. Achenbach,et al.  Vortex shedding from spheres , 1974, Journal of Fluid Mechanics.

[11]  Iman Harimi,et al.  Numerical simulation of fluid flow and forced convection heat transfer from tandem circular cylinders using overset grid method , 2012 .

[12]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[13]  P. Moin,et al.  A dynamic subgrid‐scale eddy viscosity model , 1990 .

[14]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[15]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[16]  S. Churchill,et al.  A Correlating Equation for Forced Convection From Gases and Liquids to a Circular Cylinder in Crossflow , 1977 .

[17]  Osamu Inoue,et al.  Sound generation by a two-dimensional circular cylinder in a uniform flow , 2002, Journal of Fluid Mechanics.

[18]  Daniele Vitali,et al.  Experimental study on a Savonius wind rotor for street lighting systems , 2016 .

[19]  S. A. Isaev,et al.  Comparative analysis of the calculation data on an unsteady flow around a circular cylinder obtained using the VP2/3 and fluent packages and the Spalart-Allmaras and Menter turbulence models , 2005 .

[20]  C. Williamson Vortex Dynamics in the Cylinder Wake , 1996 .

[21]  R. Ricci,et al.  Unsteady Aerodynamics of a Savonius wind rotor: a new computational approach for the simulation of energy performance , 2010 .

[22]  Sergio Pirozzoli,et al.  A low-dissipative solver for turbulent compressible flows on unstructured meshes, with OpenFOAM implementation , 2017 .

[23]  U. Ascher,et al.  Projected implicit Runge-Kutta methods for differential-algebraic equations , 1990 .

[24]  Barry Koren,et al.  Accuracy analysis of explicit Runge-Kutta methods applied to the incompressible Navier-Stokes equations , 2012, J. Comput. Phys..

[25]  Ville Vuorinen,et al.  DNSLab: A gateway to turbulent flow simulation in Matlab , 2016, Comput. Phys. Commun..

[26]  B. Schoenung,et al.  NUMERICAL CALCULATION OF LAMINAR VORTEX-SHEDDING FLOW PAST CYLINDERS , 1990 .

[27]  M. Liou,et al.  A New Flux Splitting Scheme , 1993 .

[28]  S. Orszag,et al.  Numerical investigation of transitional and weak turbulent flow past a sphere , 2000, Journal of Fluid Mechanics.

[29]  I. S. Ertesvåg,et al.  Towards Simulation of Far-Field Aerodynamic Sound from a Circular Cylinder Using OpenFoam , 2014 .

[30]  Necati Mahir,et al.  Numerical investigation of convective heat transfer in unsteady flow past two cylinders in tandem arrangements , 2008 .

[31]  Wolfgang Schröder,et al.  A lattice-Boltzmann method with hierarchically refined meshes , 2013 .

[32]  P. Moin,et al.  Suitability of upwind-biased finite difference schemes for large-eddy simulation of turbulent flows , 1997 .

[33]  Yoshimichi Tanida,et al.  On the Fluid Forces Acting on Rectangular Sectional Cylinders , 1978 .

[34]  H. Sakamoto,et al.  A STUDY ON VORTEX SHEDDING FROM SPHERES IN A UNIFORM FLOW , 1990 .

[35]  K. S. Yeo,et al.  A SVD-GFD scheme for computing 3D incompressible viscous fluid flows , 2008 .

[36]  Christophe Duwig,et al.  On the implementation of low-dissipative Runge–Kutta projection methods for time dependent flows using OpenFOAM® , 2014 .

[37]  Bernhard Müller,et al.  High order numerical simulation of Aeolian tones , 2008 .

[38]  V. C. Patel,et al.  Flow past a sphere up to a Reynolds number of 300 , 1999, Journal of Fluid Mechanics.

[39]  R. Alexander Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[40]  Michael S. Warren,et al.  Vortex Methods for Direct Numerical Simulation of Three-Dimensional Bluff Body Flows , 2002 .