|
[Sponsors] |
August 18, 2013, 10:23 |
Slow run with rhoCentralFoam
|
#1 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo ,
As a beginner in OF , i want to simulate the flow through my geometry after conversion . I have fuel and air inlet with different velocity . I am trying to use the rhoCentralFoam solver . but after giving a run my simulation is goin on from last 24 hours and i had given a end time of 10 and delta t as 0.0025 . any anybody please guide me why this is happening ? |
|
August 18, 2013, 10:38 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Greetings Sonu,
Since your a beginner with OpenFOAM, I have to ask: have you studied the first tutorial in the OpenFOAM User Guide? Namely this one: http://www.openfoam.org/docs/user/cavity.php#x5-40002.1 In addition, if you run checkMesh, you'll see some important information about your mesh, such as:
Bruno
__________________
|
|
August 18, 2013, 10:43 |
|
#3 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
I didn't understand your problem,you mean it has taken a long time?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
August 18, 2013, 11:00 |
|
#4 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo Bruno ,
i had read the tutorials and the problem is that its taking long time to finish simulation . and the mesh info is as follows : P { margin-bottom: 0.08in; } Mesh stats points: 298650 internal points: 0 faces: 594092 internal faces: 295444 cells: 148256 faces per cell: 6 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 0 Min volume = 0.030014. Max volume = 26.4353. Total volume = 329080. Cell volumes ok |
|
August 18, 2013, 11:08 |
|
#5 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
whats your case,maybe its the nature of your problem,my run takes 6 hours for .000933 s !
its hard to answer without further information,if its not the case,maybe your settings are not correct.what you have set maxCo if use?
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
August 18, 2013, 11:10 |
|
#6 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Well, reading doesn't mean that you studied it
OK then, a few questions:
__________________
|
|
August 18, 2013, 11:12 |
|
#7 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo Ehsan ,
actually my problenm is just simple i am just trying to simulate flow through my geometry which has 2 inlets one for fuel and other is for air . I had given end time as 10 s and its around 20 hrs and only 4 s simulation is over till now ....i had given only p,u and t variables in the 0 file . so can you please guide me is there anything wrong ?? thanks alot |
|
August 18, 2013, 11:34 |
|
#8 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo bruno ,
I am running it in series . and for running it i used following commands : fluentMeshToFoam //for mesh conversion) P { margin-bottom: 0.08in; } transformPoints -scale '(0.01 0.01 0.01)' ///scaling of geometry rhoCentralFoam //solver it is based on the tutorial of the oblique shock under the compressible flow in rhocentralfoam solver . this is what its showing in the output : Mean and max Courant Numbers = 0.0229302 0.200008 deltaT = 4.3472e-14 Time = 4.45301e-09 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 ExecutionTime = 70775.2 s ClockTime = 70801 s and i am using the files of system of oblique shock tutorials only . thanks for helping .... |
|
August 18, 2013, 11:46 |
|
#9 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
deltaT = 4.3472e-14 ???????
your BC's has some trouble,how much is the maximum velocity?when velocity be too high deltaT becomes too low so that courant number stay constant.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
August 18, 2013, 12:36 |
|
#10 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo Ehsan ,
The maximum velocity is 30 m/s i.e for fuel inlet and the maximum CO no in control dict file is .0.2 . |
|
August 18, 2013, 12:41 |
|
#11 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
then the problem may be with the mesh,maybe it has some too small cells.do a checkMesh as Bruno said and submit it here in the forum.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
August 18, 2013, 12:48 |
|
#12 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Quote:
@Sonu: My guess is that you have a cell somewhere that does not have a good-enough shape for the massive speed you have. I suggest that you stop the run and do some post-processing with ParaView, namely use the "Slice" filter and check where the speeds are the highest (by representing the "U" field). In addition, you can calculate the Courant number onto a field by running the following command, after the case has stopped: Code:
Co Beyond this, I suggest that you also read this post: http://www.cfd-online.com/Forums/ope...tml#post446350 post #17 - it should give you some ideas on how to diagnose cases. Best regards, Bruno
__________________
|
||
August 18, 2013, 13:26 |
|
#13 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
ok i will do the above said and check my mesh again .
thanks alot Bruno and Ehsan for helping ... |
|
August 18, 2013, 13:53 |
|
#14 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo ,
sori but i tried to see what happen when i wl give the run time as 1 s in control dict . but i think now also its gonna take long time for it as its showing in o/p : Mean and max Courant Numbers = 0.0212188 0.20001 deltaT = 4.3461e-16 Time = 1.16148e-12 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0 diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0 ExecutionTime = 1882.22 s ClockTime = 1882 s i.e large delta T so what it mean? please tell me |
|
August 18, 2013, 14:11 |
|
#15 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,975
Blog Entries: 45
Rep Power: 128 |
Hi Sonu,
The post on the other is just an example. On that other case, setting the total time to 1s was enough. But for you, given that it quickly reduces to such a small delta T, I advise you to change the keywords to the following ones in your "system/controlDict" file: Code:
stopAt nextWrite; writeControl timeStep; Bruno
__________________
|
|
August 18, 2013, 14:58 |
|
#16 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
helo Bruno ,
yes it worked and it stops at very next time . but i am wondering whats wrong with my mesh if everything else is rite , i had done simulation with this mesh only in fluent and there its working fine .... |
|
August 18, 2013, 15:49 |
|
#17 | |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
Quote:
I can have a look into the case if you want and send it to me to see if I can do something for the case.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
||
August 21, 2013, 08:33 |
|
#18 |
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26 |
Hi Sonu
1)mesh seems not be correct,it has lower resolution near walls,is it what you want? 2)p in outlet should not be zeroGradient,you have to set a back pressure in outlet. check these and let me know.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King. To Be or Not To Be,Thats the Question! The Only Stupid Question Is the One that Goes Unasked. |
|
August 31, 2013, 08:55 |
|
#19 |
Member
sonu
Join Date: Jul 2013
Location: delhi
Posts: 92
Rep Power: 12 |
greetings Ehsan ,
sori for late reply . mesh resolution is kept low near wall purposely and when i am giving the back pressure at outlet it is not accepting it , suggesting to choose other options . but same case when i tried in rhoPimpleFoam its working fine . |
|
April 21, 2015, 05:42 |
Maximum velocity in rhoCentralFoam
|
#20 |
New Member
Amolkumar
Join Date: Mar 2015
Posts: 2
Rep Power: 0 |
Please tell me what is the maximum velocity for operation in rhoCentralFoam, Since I can't exceed above 50 m/s.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
SnappyHexMesh OF-1.6-ext crashes on a parallel run | norman1981 | OpenFOAM Bugs | 5 | December 7, 2011 12:48 |
Slow very localized divergence | kyle | OpenFOAM Running, Solving & CFD | 1 | November 30, 2010 03:27 |
How to run and save tranisent case on Cray XD1 | Leon | FLUENT | 0 | October 3, 2006 21:59 |
Windows 64-bit, Distributed Parallel Run Issues... | Erich | CFX | 3 | March 28, 2006 16:36 |
Restatring run | Tim | CFX | 2 | March 3, 2006 11:54 |