[sur4j]

Hi,

I am using the OpenFOAM version 2.2.1 and have been following the tutorial for adding thermal diffusion to the existing icoFoam solver (http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam) and have encountered some problems.

The solver had compiled fine (some warning messages appeared but I read that this was because of some unused parts of the header files and should not cause problems). I continued and made the changes to the existing cavity tutorial as shown in the wiki page. I ran icoFoam and everything seemed fine apart from a time step continuity error appearing at every time:

time step continuity errors : sum local = 5.25344e-09, global = 1.40373e-18, cumulative = -2.29381e-18

I presumed that this was something to do with the corrector methods used within the solvers, is this the case?

I then viewed the results in paraview and found that after time 0.1 the temperature information just disappear's, the pressure (p) and velocity (u) information is visible and works fine, it is just the temperature that can not be viewed.

I am really confused why this is happening, I have tried redoing this many times and searched the forum but cant find any answers. I would really appreciate help with this.

Thank you

[sur4j]

This is what appeared when compiling the solver:

sk@sk-VirtualBox:~/OpenFOAM/sk-2.2.1/applications/solvers/my_icoFoam$ wmake
Making dependency list for source file my_icoFoam.C
SOURCE=my_icoFoam.C ;  g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam221/src/finiteVolume/lnInclude -I/opt/openfoam221/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam221/src/OpenFOAM/lnInclude -I/opt/openfoam221/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linuxGccDPOpt/my_icoFoam.o
/opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H: In function ‘int main(int, char**)’:
/opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
/opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H:12:16: warning: unused variable ‘momentumPredictor’ [-Wunused-variable]
/opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3  -DNoRepository -ftemplate-depth-100 -I/opt/openfoam221/src/finiteVolume/lnInclude -I/opt/openfoam221/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam221/src/OpenFOAM/lnInclude -I/opt/openfoam221/src/OSspecific/POSIX/lnInclude   -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linuxGccDPOpt/my_icoFoam.o -L/opt/openfoam221/platforms/linuxGccDPOpt/lib \
	     -lfiniteVolume -lsampling -lOpenFOAM -ldl   -lm -o /home/sk/OpenFOAM/sk-2.2.1/platforms/linuxGccDPOpt/bin/my_icoFoam

This is the final part of what I got when running the tutorial with the icoFoam command:

Time = 0.48

Courant Number mean: 0.222158 max: 0.852134
DILUPBiCG:  Solving for Ux, Initial residual = 2.3555e-07, Final residual = 2.3555e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 5.14355e-07, Final residual = 5.14355e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1.11022e-06, Final residual = 3.67162e-07, No Iterations 1
time step continuity errors : sum local = 5.49626e-09, global = -9.21479e-19, cumulative = -2.38786e-18
DICPCG:  Solving for p, Initial residual = 7.06599e-07, Final residual = 7.06599e-07, No Iterations 0
time step continuity errors : sum local = 8.27603e-09, global = -4.60574e-19, cumulative = -2.84844e-18
ExecutionTime = 0.22 s  ClockTime = 0 s

Time = 0.485

Courant Number mean: 0.222158 max: 0.852134
DILUPBiCG:  Solving for Ux, Initial residual = 2.22829e-07, Final residual = 2.22829e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 4.87534e-07, Final residual = 4.87534e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1.10464e-06, Final residual = 3.50083e-07, No Iterations 1
time step continuity errors : sum local = 5.3505e-09, global = -1.90665e-19, cumulative = -3.0391e-18
DICPCG:  Solving for p, Initial residual = 5.52457e-07, Final residual = 5.52457e-07, No Iterations 0
time step continuity errors : sum local = 7.00941e-09, global = 9.60605e-19, cumulative = -2.0785e-18
ExecutionTime = 0.22 s  ClockTime = 0 s

Time = 0.49

Courant Number mean: 0.222158 max: 0.852134
DILUPBiCG:  Solving for Ux, Initial residual = 2.09588e-07, Final residual = 2.09588e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 4.59868e-07, Final residual = 4.59868e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 8.08884e-07, Final residual = 8.08884e-07, No Iterations 0
time step continuity errors : sum local = 9.1113e-09, global = -1.24822e-19, cumulative = -2.20332e-18
DICPCG:  Solving for p, Initial residual = 9.46436e-07, Final residual = 9.46436e-07, No Iterations 0
time step continuity errors : sum local = 1.02383e-08, global = -3.21112e-19, cumulative = -2.52443e-18
ExecutionTime = 0.22 s  ClockTime = 0 s

Time = 0.495

Courant Number mean: 0.222158 max: 0.852134
DILUPBiCG:  Solving for Ux, Initial residual = 1.99665e-07, Final residual = 1.99665e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 4.36311e-07, Final residual = 4.36311e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1.0746e-06, Final residual = 3.53797e-07, No Iterations 1
time step continuity errors : sum local = 5.37651e-09, global = -1.16781e-18, cumulative = -3.69224e-18
DICPCG:  Solving for p, Initial residual = 6.81574e-07, Final residual = 6.81574e-07, No Iterations 0
time step continuity errors : sum local = 8.06059e-09, global = -5.29395e-21, cumulative = -3.69754e-18
ExecutionTime = 0.22 s  ClockTime = 0 s

Time = 0.5

Courant Number mean: 0.222158 max: 0.852134
DILUPBiCG:  Solving for Ux, Initial residual = 1.89493e-07, Final residual = 1.89493e-07, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 4.14522e-07, Final residual = 4.14522e-07, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1.06665e-06, Final residual = 3.39604e-07, No Iterations 1
time step continuity errors : sum local = 5.25344e-09, global = 1.40373e-18, cumulative = -2.29381e-18
DICPCG:  Solving for p, Initial residual = 5.36118e-07, Final residual = 5.36118e-07, No Iterations 0
time step continuity errors : sum local = 6.86432e-09, global = 5.09378e-19, cumulative = -1.78443e-18
ExecutionTime = 0.23 s  ClockTime = 1 s

End

Also, in my time folders I only have a temperature (T) file in the 0 folder, it does not appear in any of the others (0.1, 0.2, 0.3, 0.4, 0.5) whilst U and P appear in them all.

If any further information is required please let me know.

Thank you.

[wyldckat]

Greetings sur4j,

Although the tutorial does say "add temperature to icoFoam", the idea is that you're meant to run the customized version of icoFoam, not icoFoam itself. According to the tutorial, you're meant to run my_icoFoam:

my_icoFoam

Best regards,
Bruno

[sur4j]

Thank you so much for your help Bruno, the simulation is working fine now.

My simulation is a static simulation in which the fluid is not moving, I am trying to have a variable viscosity in this and as far as I am aware the viscosity of a stationary fluid does not have any effect on the heat transfer, please correct me on this if I am wrong.

I would like this viscosity to be visible in para view and it is variable such that when it gets to a certain temperature it starts increasing at a set rate and then at some other higher temperature it stops increases and stays at the final value. How would you suggest I start this simulation, should I create some sort of dummy variable in the same sort of way temperature was added to the simulation?

Again, thank you so much for your help I really appreciated it.

[wyldckat]

Hi sur4j,

I'll try to answer you on the thread you have dedicated for that question: http://www.cfd-online.com/Forums/ope...viscosity.html

[Hamzeh_Mirab]

Hello every one,
I have a same problem in adding temperatture to icoFoam solver.
I have done all steps in manual in openWiki.net and wmake and modify initial and boundary condition and scheme and solution but after running of case by typing my_interFoam command only u and p in paraFoam and in time step output folders appear and T output doesnt appear.please guide me

BestWishes

[Bahram]

Originally Posted by Hamzeh_Mirab

Hello every one,
I have a same problem in adding temperatture to icoFoam solver.
I have done all steps in manual in openWiki.net and wmake and modify initial and boundary condition and scheme and solution but after running of case by typing my_interFoam command only u and p in paraFoam and in time step output folders appear and T output doesnt appear.please guide me

BestWishes

Hi

Can you post your changes and also your example? It is really hard to judge without seeing the files!

Bests
Bahram

[demichie]

Originally Posted by Hamzeh_Mirab

Hello every one,
I have a same problem in adding temperatture to icoFoam solver.
I have done all steps in manual in openWiki.net and wmake and modify initial and boundary condition and scheme and solution but after running of case by typing my_interFoam command only u and p in paraFoam and in time step output folders appear and T output doesnt appear.please guide me

BestWishes

Are you modifying icoFoam or interFoam?

Mattia

[Hamzeh_Mirab]

It fixed after restarting of my system???????