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

Mistake in PISO loop for interFoam solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 27, 2006, 21:04
Default Hello Friends I was wonderi
  #1
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hello Friends

I was wondering if there is a mistake in the PISO loop for the interFoam , rasInterFoam solvers.

In Dr. Rusche's PHD page 148,

(Del.V) = 0 , but at the interface density is rapidly changing so shouln't (Del.(rho*V)= 0);

Looking at the PISO loop in rasInterFoam we have laplacian(rUAf,pd)==div(phi) .

should this be changed to something like

laplacian(rUAF*rho,pd) == div(phi*rho)

Thanks a lot in advance

Kumar
kumar2 is offline   Reply With Quote

Old   June 28, 2006, 12:57
Default Dear Friends Could someone
  #2
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Dear Friends

Could someone comment on this post ?

thanks

kumar
kumar2 is offline   Reply With Quote

Old   June 28, 2006, 13:09
Default No, because the interFoam solv
  #3
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
No, because the interFoam solver conserves volume not mass (both fluid are assumed incompressible). This is done specifically to get around the problem of large gradients at the interface.
eugene is offline   Reply With Quote

Old   June 30, 2006, 18:26
Default Hi Eugene Thanks a lot for
  #4
Senior Member
 
kumar
Join Date: Mar 2009
Posts: 112
Rep Power: 8
kumar2 is on a distinguished road
Hi Eugene

Thanks a lot for your reply. I have just one more question.

In equation of 4.28 of henrik rusche's phd thesis we have phi=phi* - (1/A_D)_face * |S|*faceGradient(P)

But in 4.30 we have something like

div( (1/A_D)_face , gradient(P) ] = diver.(phi*)

How is |S|*faceGradient(P) converted to gradient(P) ??

In all openfoam solvers the equation 4.30 is implemented and they are dimensionally correct.

Another question if rUA = 1.0/UEqn.A()

is rUAf = rUA interpolated on the face and multiplied with surface area or it is just rUA interpolated to the faces

Thanks a lot once again

Kumar
kumar2 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
About interFoam solver zou_mo OpenFOAM Running, Solving & CFD 127 May 25, 2011 16:30
PEqnflux in compressible and incompressible PISO loop dbxmcf OpenFOAM Running, Solving & CFD 0 October 7, 2008 10:58
PISO loop 21kalee OpenFOAM Running, Solving & CFD 2 January 15, 2008 06:31
About interfoam solver qiu OpenFOAM Running, Solving & CFD 0 May 6, 2007 22:48
Dimensions of laplacian in PISO loop kumar2 OpenFOAM Running, Solving & CFD 2 July 3, 2006 14:34


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