CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Slow run with rhoCentralFoam

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

Like Tree4Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 18, 2013, 10:23
Default Slow run with rhoCentralFoam
  #1
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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 ?
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 10:38
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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:
  • The total number of cells.
  • The minimum cell volume.
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   August 18, 2013, 10:43
Default
  #3
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Old   August 18, 2013, 11:00
Default
  #4
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 11:08
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Old   August 18, 2013, 11:10
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Well, reading doesn't mean that you studied it

OK then, a few questions:
  1. In the output of the rhoCentralFoam solver, what Courant number values is it giving you for each time iteration?
  2. Are you running the case in serial or in parallel?
  3. How long does it take to run each time iteration?
  4. How exactly did you run the case? More specifically, what steps did you take to run it?
  5. What settings do you have inside the files "system/fvSchemes" and "system/fvSolution"?
    • (Attach the files or follow the instructions from the second link on my signature for pasting code with "[CODE]".)
  6. On which tutorial case did you base your case on?
wyldckat is offline   Reply With Quote

Old   August 18, 2013, 11:12
Default
  #7
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 11:34
Default
  #8
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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 ....
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 11:46
Default
  #9
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Old   August 18, 2013, 12:36
Default
  #10
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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 .
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 12:41
Default
  #11
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Old   August 18, 2013, 12:48
Default
  #12
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by immortality View Post
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.
Sonu already posted that info in post #4.

@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
It should give you a new field named "Co" that you can see in ParaView. Use that field and the Slice to search for where you have the highest Courant number values.


Beyond this, I suggest that you also read this post: Strange Results at Tank Outlet with InterFoam post #17 - it should give you some ideas on how to diagnose cases.

Best regards,
Bruno
Endel and immortality like this.
wyldckat is offline   Reply With Quote

Old   August 18, 2013, 13:26
Default
  #13
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
ok i will do the above said and check my mesh again .
thanks alot Bruno and Ehsan for helping ...
immortality likes this.
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 13:53
Default
  #14
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 14:11
Default
  #15
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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;
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   August 18, 2013, 14:58
Default
  #16
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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 ....
yash.aesi is offline   Reply With Quote

Old   August 18, 2013, 15:49
Default
  #17
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
Quote:
Originally Posted by yash.aesi View Post
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
Hi Suno
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.
yash.aesi likes this.
__________________
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.
immortality is offline   Reply With Quote

Old   August 21, 2013, 08:33
Default
  #18
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,186
Rep Power: 16
immortality is on a distinguished road
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.
immortality is offline   Reply With Quote

Old   August 31, 2013, 08:55
Default
  #19
Member
 
sonu
Join Date: Jul 2013
Location: delhi
Posts: 79
Rep Power: 3
yash.aesi is on a distinguished road
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 .
yash.aesi is offline   Reply With Quote

Old   April 21, 2015, 05:42
Default Maximum velocity in rhoCentralFoam
  #20
New Member
 
Amolkumar
Join Date: Mar 2015
Posts: 2
Rep Power: 0
amolkumar is on a distinguished road
Please tell me what is the maximum velocity for operation in rhoCentralFoam, Since I can't exceed above 50 m/s.
amolkumar 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
SnappyHexMesh OF-1.6-ext crashes on a parallel run norman1981 OpenFOAM Bugs 5 December 7, 2011 13:48
Slow very localized divergence kyle OpenFOAM Running, Solving & CFD 1 November 30, 2010 04: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 12:54


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