# (chtMultiRegionFoam) reducing gravity increase the simulation time !

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

 February 8, 2010, 02:00 (chtMultiRegionFoam) reducing gravity increase the simulation time ! #1 Member   Join Date: Dec 2009 Posts: 46 Rep Power: 8 Hi foamers; when i simulate a simple conjugate heat transfer problem using chtMultiRegionFoam solver and i don't want any buoyancy effects appear (plug flow) ,, i changed the g file from (0 0 -9.81) to (0 0 0) , when i did this the simulation time increase horribly i don't know what is the reason why g effects the simulation greatly ,, and i tried it with other simulation and also the tutorial case and did the same thing (great increase in the simulation time for every step ! ) any one can help thanks very much ..

 February 9, 2010, 17:10 #2 Member   Join Date: Dec 2009 Posts: 46 Rep Power: 8 can any one help me to get red of that bad behavior of chtMultiRegionFoam solver (it may be a bug) thanks in advance..

 February 12, 2010, 03:42 effects of g=0 #3 New Member   Join Date: Nov 2009 Posts: 3 Rep Power: 8 Hi openofam1, hi all trying to run a heat transfer simulation with chtMultiRegionFoam takes very long till I reach steady conditions (if I reach it!?). Then I read your post and tried to run the simulation faster without any gravity. Analysing the written log.file I had to realize that there is nearly no temperature change but a significant change in execution time as you already mentioned. Here you can see the difference: Temperature at the beginning: 298,15K As a heat source I use q=606[W/(m*m)] in changeDictionaryDict. According to the used area it matches P=1,5W without g: Solving for solid region cube DICPCG: Solving for T, Initial residual = 0.009375165, Final residual = 8.318404e-07, No Iterations 24 DICPCG: Solving for T, Initial residual = 0.001351302, Final residual = 7.010276e-07, No Iterations 17 Min/max T:min(T) [0 0 0 1 0 0 0] 298.2413 max(T) [0 0 0 1 0 0 0] 298.3844 ExecutionTime = 1849.95 s ClockTime = 2964 s Region: domain0 Courant Number mean: 3.232212e-05 max: 0.2909719 deltaT = 1.428571 Time = 150.1 Moreover after 2522sec it only reaches max(T) [0 0 0 1 0 0 0] 299.2234 with g: Solving for solid region cube DICPCG: Solving for T, Initial residual = 0.00038313, Final residual = 3.640629e-08, No Iterations 2 DICPCG: Solving for T, Initial residual = 2.75657e-06, Final residual = 4.424209e-08, No Iterations 1 Min/max T:min(T) [0 0 0 1 0 0 0] 299.3279 max(T) [0 0 0 1 0 0 0] 299.4493 ExecutionTime = 27342.29 s ClockTime = 40150 s Region: domain0 Courant Number mean: 0.001197062 max: 0.2999681 deltaT = 0.002828811 Time = 150.445 -> more than 1K difference and nearly 15 times faster! In my Opinion without g the convection is not considered and thats the reason for the lower executionTime Does anybody know whether it is possible to run the case faster and how?? I can't imagine that it takes that long but I do not know what's to change!? overall the domain0 has 92520 cells cube has 15424 cells plate 1540 cells but I think this number is necessary because a coarsen mesh has bad effects on the stl-files. If anybody knows a solution I would be very grateful! best regards Bene

 February 12, 2010, 08:55 #4 New Member   J H Join Date: Mar 2009 Location: Germany Posts: 20 Rep Power: 9 Hi Bene85, how did you use the heat source flux? What have i to type in the changeDictionaryDict file to use it? It would be very helpful for me.

 March 6, 2010, 04:44 #5 New Member   Join Date: Nov 2009 Posts: 3 Rep Power: 8 Hi kawuppdich, sorry for the late answer, but I have been quite busy. Using the heat-flux I used the following in the changeDictionaryDict: cube_to_plate { type solidWallHeatFluxTemperature; K K; q uniform 6302; // the value depends on Power and face: here e.g. P=15,6W and face=0,05m x 0,04951m -> q=6302W/(m*m) value uniform 298.15; } }This entry has to be made in both changeDictionaryDicts of cube (cube_to_plate) and plate (plate_to_cube) to specify both sides, but one has to be 0 that the heat flux only takes place in one direction. I hope you understood what I meant.

 March 10, 2010, 05:06 #6 New Member   J H Join Date: Mar 2009 Location: Germany Posts: 20 Rep Power: 9 thanks for reply, it works fine

 March 10, 2010, 07:40 #7 New Member   Join Date: Nov 2009 Posts: 3 Rep Power: 8 Hi kawuppdich and OpenFOAMers, I'm glad that it works! What can you say about your ExecutionTime? Does it take very long? Do you run your simulation in parallel? I have another question. Maybe you can help me!? OpenFOAM has libraries for so called RAS turbulence models. Do you know whether RAS is conform to the abbreviation RANS? Would be very grateful for any help, but it's quite urgent. Thanks Bene85

 March 10, 2010, 10:41 #8 New Member   J H Join Date: Mar 2009 Location: Germany Posts: 20 Rep Power: 9 in my case there is no change of calculation time. I run my testcase in parallel on 2 cpus RAS stands for Reynold Averaged Stress and RANS for Reynold Averaged Navier Stokes and thats all what i know about it. I´m just a beginner in cfd so sorry. greetings

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Whyman OpenFOAM Programming & Development 36 March 30, 2015 13:35 teguhtf FLUENT 12 October 7, 2009 22:54 Janendra Telisinghe FLUENT 0 December 11, 2007 07:14 msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58 Colin Main CFD Forum 6 November 2, 2003 15:13

All times are GMT -4. The time now is 04:29.