
[Sponsors] 
November 23, 2013, 13:39 
Problems adding temperature to icoFoam

#1 
Member
Join Date: Aug 2013
Posts: 60
Rep Power: 4 
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.25344e09, global = 1.40373e18, cumulative = 2.29381e18 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 

November 23, 2013, 15:48 

#2 
Member
Join Date: Aug 2013
Posts: 60
Rep Power: 4 
This is what appeared when compiling the solver:
Code:
sk@skVirtualBox:~/OpenFOAM/sk2.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 Wnounusedparameter Woldstylecast O3 DNoRepository ftemplatedepth100 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’ [Wunusedvariable] /opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H:12:16: warning: unused variable ‘momentumPredictor’ [Wunusedvariable] /opt/openfoam221/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [Wunusedvariable] g++ m32 Dlinux DWM_DP Wall Wextra Wnounusedparameter Woldstylecast O3 DNoRepository ftemplatedepth100 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 addneeded Xlinker noasneeded Make/linuxGccDPOpt/my_icoFoam.o L/opt/openfoam221/platforms/linuxGccDPOpt/lib \ lfiniteVolume lsampling lOpenFOAM ldl lm o /home/sk/OpenFOAM/sk2.2.1/platforms/linuxGccDPOpt/bin/my_icoFoam Code:
Time = 0.48 Courant Number mean: 0.222158 max: 0.852134 DILUPBiCG: Solving for Ux, Initial residual = 2.3555e07, Final residual = 2.3555e07, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 5.14355e07, Final residual = 5.14355e07, No Iterations 0 DICPCG: Solving for p, Initial residual = 1.11022e06, Final residual = 3.67162e07, No Iterations 1 time step continuity errors : sum local = 5.49626e09, global = 9.21479e19, cumulative = 2.38786e18 DICPCG: Solving for p, Initial residual = 7.06599e07, Final residual = 7.06599e07, No Iterations 0 time step continuity errors : sum local = 8.27603e09, global = 4.60574e19, cumulative = 2.84844e18 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.22829e07, Final residual = 2.22829e07, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 4.87534e07, Final residual = 4.87534e07, No Iterations 0 DICPCG: Solving for p, Initial residual = 1.10464e06, Final residual = 3.50083e07, No Iterations 1 time step continuity errors : sum local = 5.3505e09, global = 1.90665e19, cumulative = 3.0391e18 DICPCG: Solving for p, Initial residual = 5.52457e07, Final residual = 5.52457e07, No Iterations 0 time step continuity errors : sum local = 7.00941e09, global = 9.60605e19, cumulative = 2.0785e18 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.09588e07, Final residual = 2.09588e07, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 4.59868e07, Final residual = 4.59868e07, No Iterations 0 DICPCG: Solving for p, Initial residual = 8.08884e07, Final residual = 8.08884e07, No Iterations 0 time step continuity errors : sum local = 9.1113e09, global = 1.24822e19, cumulative = 2.20332e18 DICPCG: Solving for p, Initial residual = 9.46436e07, Final residual = 9.46436e07, No Iterations 0 time step continuity errors : sum local = 1.02383e08, global = 3.21112e19, cumulative = 2.52443e18 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.99665e07, Final residual = 1.99665e07, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 4.36311e07, Final residual = 4.36311e07, No Iterations 0 DICPCG: Solving for p, Initial residual = 1.0746e06, Final residual = 3.53797e07, No Iterations 1 time step continuity errors : sum local = 5.37651e09, global = 1.16781e18, cumulative = 3.69224e18 DICPCG: Solving for p, Initial residual = 6.81574e07, Final residual = 6.81574e07, No Iterations 0 time step continuity errors : sum local = 8.06059e09, global = 5.29395e21, cumulative = 3.69754e18 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.89493e07, Final residual = 1.89493e07, No Iterations 0 DILUPBiCG: Solving for Uy, Initial residual = 4.14522e07, Final residual = 4.14522e07, No Iterations 0 DICPCG: Solving for p, Initial residual = 1.06665e06, Final residual = 3.39604e07, No Iterations 1 time step continuity errors : sum local = 5.25344e09, global = 1.40373e18, cumulative = 2.29381e18 DICPCG: Solving for p, Initial residual = 5.36118e07, Final residual = 5.36118e07, No Iterations 0 time step continuity errors : sum local = 6.86432e09, global = 5.09378e19, cumulative = 1.78443e18 ExecutionTime = 0.23 s ClockTime = 1 s End If any further information is required please let me know. Thank you. 

November 23, 2013, 16:23 

#3 
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,215
Blog Entries: 35
Rep Power: 94 
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: Code:
my_icoFoam Bruno
__________________
I'll be at OFW11 in Portugal 

November 23, 2013, 17:56 

#4 
Member
Join Date: Aug 2013
Posts: 60
Rep Power: 4 
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. 

November 24, 2013, 07:37 

#5 
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,215
Blog Entries: 35
Rep Power: 94 
Hi sur4j,
I'll try to answer you on the thread you have dedicated for that question: Time and temperature dependent viscosity
__________________
I'll be at OFW11 in Portugal 

January 1, 2015, 10:49 
dont apperar T value in adding temperature to icoFoam solver

#6 
New Member
Hamzeh
Join Date: Oct 2014
Posts: 11
Rep Power: 3 
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 

January 2, 2015, 06:49 

#7  
New Member
Bahram Haddadi
Join Date: Feb 2014
Location: Vienna, Austria
Posts: 19
Rep Power: 4 
Quote:
Can you post your changes and also your example? It is really hard to judge without seeing the files! Bests Bahram 

January 2, 2015, 16:12 

#8  
Member
Mattia de\' Michieli Vitturi
Join Date: Mar 2009
Posts: 42
Rep Power: 8 
Quote:
Mattia 

January 11, 2015, 02:14 

#9 
New Member
Hamzeh
Join Date: Oct 2014
Posts: 11
Rep Power: 3 
It fixed after restarting of my system???????


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
libOpenSMOKE  Tobi  OpenFOAM Programming & Development  506  September 28, 2015 06:07 
abnormal temperature near interface when adding energy equation to interFoam  houkensjtu  OpenFOAM  4  June 26, 2013 13:40 
Problems With courant number Icofoam  basilwatson  OpenFOAM Running, Solving & CFD  0  May 21, 2011 09:39 
Problems with icoFoam  bego  OpenFOAM  1  November 16, 2009 06:01 
Problems with adding two meshes together in CFX Pre  peterputer1  CFX  2  September 23, 2009 08:08 