# rhoSimplecFoam - which algorithm is implemented - references

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

 April 22, 2013, 13:22 rhoSimplecFoam - which algorithm is implemented - references #1 New Member   Manuel Join Date: Jan 2013 Posts: 21 Rep Power: 5 Dear Foamers, I have a problem in implementing classical all mach flow solver. Mainly I would like to implement the all-Mach flow solver described in Ferziger and Peric book, but I have some problem in running it. As test case I chose the classical shock reflecting wall 2D domain, with both left and down inlet being supersonic. My algorithm diverges. I have tested on the same test case also the rhoSimpleFoam and rhoSimplecFoam solvers: the first one diverges too, the second one seems to work well. But, having a look at the algorithm implemented therein, I do not understand what kind of algorithm is implemented in the "pEqn.H" file. Could anybody please help me in understanding this algorithm? Or could anybody please indicate me some reference which I have to look for in order to understand the algorithm implemented therein? I thank you in advance for your help and attention. Manuel

 April 23, 2013, 06:47 #2 Senior Member   HECKMANN Frédéric Join Date: Jul 2010 Posts: 237 Rep Power: 9 rhoSimplecFoam is either SIMPLE, SIMPLEC or PISO depending of your settings in fvsolution. You can have a look into to Ferziger and Perics book for the algorithms. Ps. The algorithm depends of the number of correction iterations you set in fvsolution.

 April 23, 2013, 11:07 #4 Senior Member   HECKMANN Frédéric Join Date: Jul 2010 Posts: 237 Rep Power: 9 what you quote is only used in the case of the transonic = true. I guess it is a specific correction to avoid problems in the transonic region. From my old knowledge, the equations in the transonic region switch from parabolic to hyperbolic. This might be the reason of this expression.

 April 23, 2013, 11:25 #5 New Member   Manuel Join Date: Jan 2013 Posts: 21 Rep Power: 5 Yes, exactly, the part of code I posted refers to the transonic case. And, as you wrote, when the problem becomes supersonic, in the Navier-Stokes equations the hyperbolic part is not negligible. So, concerning the part of code I posted, do you know some reference in literature where I can find the algorithms and the adjustments implemented in "pEqn.H"? Because I am not able to understand the reason why the modifying terms are implemented in this way. I thank you a lot again. Bye

 April 23, 2013, 12:33 #6 Senior Member   HECKMANN Frédéric Join Date: Jul 2010 Posts: 237 Rep Power: 9 I've never seen this exact code and I don't know where you can find the reference... Sorry

 April 24, 2013, 03:39 #7 New Member   Manuel Join Date: Jan 2013 Posts: 21 Rep Power: 5 Ok. No problem. I'll continue to look for any helpful references in literature. If anybody will be able to help me, please do not hesitate to contact me..!! Thanks Bye

September 8, 2014, 07:28
#8
Senior Member

Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 10
Quote:
 Originally Posted by Manuel CFD Ok. No problem. I'll continue to look for any helpful references in literature. If anybody will be able to help me, please do not hesitate to contact me..!! Thanks Bye
Hello Manuel,

I would like to know, Were you able to solve the stated problem in your post?

 Tags all-mach, compressible, rhosimplecfoam, rhosimplefoam

 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 Yan Kai Main CFD Forum 0 April 18, 2007 03:48 Yan Kai FLUENT 0 April 13, 2007 23:17 BeatrÃ­z Navarro NUMECA 0 June 24, 2006 18:23 Mihai ARGHIR Main CFD Forum 1 December 17, 2000 08:01 Jonathan Castro Main CFD Forum 3 December 10, 1999 05:59

All times are GMT -4. The time now is 03:21.