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

rhoSimplecFoam - which algorithm is implemented - references

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 22, 2013, 13:22
Default rhoSimplecFoam - which algorithm is implemented - references
  #1
New Member
 
Manuel
Join Date: Jan 2013
Posts: 21
Rep Power: 4
Manuel CFD is on a distinguished road
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
Manuel CFD is offline   Reply With Quote

Old   April 23, 2013, 06:47
Default
  #2
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 8
fredo490 is on a distinguished road
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.
fredo490 is offline   Reply With Quote

Old   April 23, 2013, 08:55
Default
  #3
New Member
 
Manuel
Join Date: Jan 2013
Posts: 21
Rep Power: 4
Manuel CFD is on a distinguished road
Dear Frederic,

I thank you a lot for your attention and for your quick answer.

I had already read Peric-Ferziger book concerning the SIMPLE, SIMPLEC and PISO methods, and things seem to work well: I have implemented some different versions of these algorithms inside OpenFOAM, and they seem to perform well for incompressible flows.

But problems arise when I try to simulate all-Mach flows. In particular I would like to implement a pressure-based all-Mach flow solver. I read the Peric-Ferziger book concerning this subject (see chapter 10), and I tried to implement the corresponding solver based on SIMPLE algorithm. Unfortunately, when I tested it on the shock-reflecting-wall test case, the code diverges. So I tried to use the codes already implemented in OpenFOAM, that is, "rhoSimpleFoam" and "rhoSimplecFoam": the first one diverges too, the second one seems to perform well.
This is why I had a look at the algorithms implemented therein. But I found an algorithm which does not seem to be the one described in Ferziger-Peric book. In the following I post a short part of the "pEqn.H" file, which I am not able to understand:

surfaceScalarField phic
(
"phic",
fvc::interpolate(rho/AtU - rho/AU)*fvc::snGrad(p)*mesh.magSf()
+ phid*(fvc::interpolate(p) - fvc::interpolate(p, "UD"))
);

fvScalarMatrix pEqn
(
fvm::div(phid, p)
+ fvc::div(phic)
- fvm::Sp(fvc::div(phid), p)
+ fvc::div(phid)*p
- fvm::laplacian(rho/AtU, p)
);

What does the variable "phic" represent? Where does the difference between interpolated pressures come from? And, moreover, where does this kind of pressure correction equation arise from?

These are the main problems I am not able to solve. Do you know something about those?

I thank you a lot in advance.

Bye
Manuel CFD is offline   Reply With Quote

Old   April 23, 2013, 11:07
Default
  #4
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 8
fredo490 is on a distinguished road
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.
fredo490 is offline   Reply With Quote

Old   April 23, 2013, 11:25
Default
  #5
New Member
 
Manuel
Join Date: Jan 2013
Posts: 21
Rep Power: 4
Manuel CFD is on a distinguished road
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
Manuel CFD is offline   Reply With Quote

Old   April 23, 2013, 12:33
Default
  #6
Senior Member
 
HECKMANN Frédéric
Join Date: Jul 2010
Posts: 236
Rep Power: 8
fredo490 is on a distinguished road
I've never seen this exact code and I don't know where you can find the reference... Sorry
fredo490 is offline   Reply With Quote

Old   April 24, 2013, 03:39
Default
  #7
New Member
 
Manuel
Join Date: Jan 2013
Posts: 21
Rep Power: 4
Manuel CFD is on a distinguished road
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
Manuel CFD is offline   Reply With Quote

Old   September 8, 2014, 07:28
Default
  #8
Senior Member
 
Tushar Chourushi
Join Date: Jul 2009
Location: IIT-Indore, India
Posts: 319
Blog Entries: 1
Rep Power: 9
Tushar@cfd is on a distinguished road
Quote:
Originally Posted by Manuel CFD View Post
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?
Tushar@cfd is offline   Reply With Quote

Reply

Tags
all-mach, compressible, rhosimplecfoam, rhosimplefoam

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
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai Main CFD Forum 0 April 18, 2007 03:48
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai FLUENT 0 April 13, 2007 23:17
CFD algorithm design with reconfigurable computing Beatríz Navarro NUMECA 0 June 24, 2006 18:23
references about the fan/radiator model Mihai ARGHIR Main CFD Forum 1 December 17, 2000 08:01
SIMPLE algorithm Jonathan Castro Main CFD Forum 3 December 10, 1999 05:59


All times are GMT -4. The time now is 17:32.