[uli]

hi

I am trying to get a turbulence model (for heat flux) running that was implemented in simplefoam the following way:

• introducing a volScalarField T (temperature) in createfields
• introducing a TEqn.H file that includes the energy equation AND the turbulence model for the heat flux, which itself includes wallfunctions for some coefficients
• adding "include TEqn.H" to mysimpleFoam.C
• making other *.H files (like wallFvPatch.H and nearWallDist.H) available by adding "include <filename.H>"

After some minor changes due to the newer OF version I use (2.1.1) I am able to compile and run the modified simpleFoam solver.

But for pimpleFoam it does not work (I applied the above steps on pimpleFoam). Simplefoam uses RASModel.H whereas pimpleFoam uses turbulenceModel.H. What is the difference?

Without replacing

include "turbulenceModel.H"

by

include "RASModel.H"

in the mypimpleFoam.C file I get

/home/uli/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude/nearWallDist.H: In function 'int main(int, char**)':
/home/uli/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude/nearWallDist.H:44:1: error: 'namespace' definition is not allowed here
In file included from pimpleFoamT_YSC.C:75:0:
TEqn.H:50:77: error: 'nearWallDist' was not declared in this scope

when I try to compile the solver.

So I replaced turbulenceModel.H by RASModel.H and it compiled without error message, but running a case I get:

Starting time loop
Courant Number mean: 0.0219895 max: 0.686747
Time = 0.002059
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.15401e-07, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.46016e-07, No Iterations 4
　
--> FOAM FATAL ERROR: 
hanging pointer, cannot dereference
From function PtrList::operator[] const
in file /home/uli/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/PtrListI.H at line 134.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/home/uli/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/uli/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 
at pimpleFoamT_YSC.C:0
#3 
in "/home/uli/OpenFOAM/uli-2.1.1/platforms/linux64GccDPOpt/bin/pimpleFoamT_YSC"
#4 __libc_start_main in "/lib64/libc.so.6"
#5 
at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116

Any help is appreciated.

I use RAS and normally pisoFoam but since the model was implemented in simpleFoam I thought it might be easier to use pimpleFoam for this case.

Uli