CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Problems adding temperature to icoFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By wyldckat

Reply
 
LinkBack Thread Tools Display Modes
Old   November 23, 2013, 13:39
Default Problems adding temperature to icoFoam
  #1
Member
 
Join Date: Aug 2013
Posts: 60
Rep Power: 3
sur4j is on a distinguished road
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 is offline   Reply With Quote

Old   November 23, 2013, 15:48
Default
  #2
Member
 
Join Date: Aug 2013
Posts: 60
Rep Power: 3
sur4j is on a distinguished road
This is what appeared when compiling the solver:

Code:
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:

Code:
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.
sur4j is offline   Reply With Quote

Old   November 23, 2013, 16:23
Default
  #3
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
Best regards,
Bruno
Annier likes this.
wyldckat is offline   Reply With Quote

Old   November 23, 2013, 17:56
Default
  #4
Member
 
Join Date: Aug 2013
Posts: 60
Rep Power: 3
sur4j is on a distinguished road
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.
sur4j is offline   Reply With Quote

Old   November 24, 2013, 07:37
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi sur4j,

I'll try to answer you on the thread you have dedicated for that question: Time and temperature dependent viscosity
wyldckat is offline   Reply With Quote

Old   January 1, 2015, 10:49
Default dont apperar T value in adding temperature to icoFoam solver
  #6
New Member
 
Hamzeh
Join Date: Oct 2014
Posts: 11
Rep Power: 2
Hamzeh_Mirab is on a distinguished road
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
Hamzeh_Mirab is offline   Reply With Quote

Old   January 2, 2015, 06:49
Default
  #7
New Member
 
Bahram's Avatar
 
Bahram Haddadi
Join Date: Feb 2014
Location: Vienna, Austria
Posts: 17
Rep Power: 3
Bahram is on a distinguished road
Quote:
Originally Posted by Hamzeh_Mirab View Post
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
Bahram is offline   Reply With Quote

Old   January 2, 2015, 16:12
Default
  #8
New Member
 
Mattia de\' Michieli Vitturi
Join Date: Mar 2009
Posts: 23
Rep Power: 8
demichie is on a distinguished road
Quote:
Originally Posted by Hamzeh_Mirab View Post
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
demichie is offline   Reply With Quote

Old   January 11, 2015, 02:14
Default
  #9
New Member
 
Hamzeh
Join Date: Oct 2014
Posts: 11
Rep Power: 2
Hamzeh_Mirab is on a distinguished road
It fixed after restarting of my system???????
Hamzeh_Mirab is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
libOpenSMOKE Tobi OpenFOAM Programming & Development 496 Yesterday 02:41
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


All times are GMT -4. The time now is 00:20.