CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   How is the turbulence model called in openfoam? (

hz283 November 17, 2012 19:10

How is the turbulence model called in openfoam?

I am looking into how the RAS models are called (for example how k and epsilon equations are called to be solved in pisofoam). I had some questions about this.

Take PISOFOAM as an example, I know that turbulence->correct(); is used to solve the k and epsilon equations and correct the turbulent viscosity. And then turbulence is defined in creatFields. H as follows:

00039 autoPtr<incompressible::turbulenceModel> turbulence
00040 (
00041 incompressible::turbulenceModel::New(U, phi, laminarTransport)
00042 );

in the class turbulenceModel, the function correct appears as follows:

00127 void turbulenceModel::correct()
00128 {
00129 transportModel_.correct();
00130 }

In fact, in transportModel, there is no function defination for correction(). Thus I do not know how the main solver, PISOFOAM, continue to call the k and epsilon equations. I really appreciate it if anyone can give me some help with this problem.


jona November 18, 2012 01:11

Hi hz283,

Foam::incompressible::turbulenceModel is an abstract class where the function correct() is virtual, so the child class, e.g. RASmodel, has to redefine it.

RASmodel declares it also as a virtual method, so the next child, e.g. kEpsilon, has to redefine it as well.

In the definition of correct() in kEpsilon finally the additional two equations for k and epsilon are solved.


explanation of virtual functions:
(see Virtual Members)

Have fun,


hz283 November 18, 2012 06:08

Hi Jorn,

Thank you so much for your help.


All times are GMT -4. The time now is 08:15.