PEqnflux in compressible and incompressible PISO loop

dbxmcf · October 7, 2008, 10:58

Hi, all:

I have a question about the pEqn.flux() in the PISO loop, for the incompressible(icoFoam) case, after solving the pressure correction equation:

phi -= pEqn.flux();

while in the compressible case (sonicFoam),

phi = pEqn.flux();

My understanding is the incompressible flow is just a special case for the compressible flow (rho=const), no matter the pressure correction equation looks like. Can anyone explain the difference for the two cases?

Thanks a lot!

David1 · July 11, 2015, 15:21

Hi Roy,

Did you eventually found your answer? Because I have the same question you have/had.

Thank you!

AshwaniAssam · September 27, 2019, 01:14

Did it get resolved? Because I do not see anywhere else in OpenFOAM except in the sonicFOAM this term

phi = pEqn.flux();

without any addition or subtraction.

Is it right?

AshwaniAssam · September 27, 2019, 06:13

I found the answer here: PISO Details

In sonicFoam
the ".flux()" function returns the influence of all implicit terms on the flux. In this case everything is implicit, so you can assign it like that
phi = pEqn.flux();

So, what I understand is that .flux() operator gives the flux terms of the equation which are handled implicitly.

October 7, 2008, 10:58	Hi, all: I have a question	#1
dbxmcf Member roy fokker Join Date: Mar 2009 Posts: 44 Rep Power: 17	Hi, all: I have a question about the pEqn.flux() in the PISO loop, for the incompressible(icoFoam) case, after solving the pressure correction equation: phi -= pEqn.flux(); while in the compressible case (sonicFoam), phi = pEqn.flux(); My understanding is the incompressible flow is just a special case for the compressible flow (rho=const), no matter the pressure correction equation looks like. Can anyone explain the difference for the two cases? Thanks a lot!

September 27, 2019, 01:14		#3
AshwaniAssam Senior Member Ashwani Join Date: Sep 2013 Location: Hyderabad Posts: 154 Rep Power: 12	Did it get resolved? Because I do not see anywhere else in OpenFOAM except in the sonicFOAM this term phi = pEqn.flux(); without any addition or subtraction. Is it right? Last edited by AshwaniAssam; September 27, 2019 at 01:15. Reason: Typos correction

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Incompressible boussiesq flow pressure BC piso and continuity problems	tehache	OpenFOAM Running, Solving & CFD	4	November 6, 2008 03:08
PISO loop	21kalee	OpenFOAM Running, Solving & CFD	2	January 15, 2008 05:31
Compressible -> incompressible.	Jinwon	Main CFD Forum	6	November 23, 2007 21:07
Dimensions of laplacian in PISO loop	kumar2	OpenFOAM Running, Solving & CFD	2	July 3, 2006 14:34
Mistake in PISO loop for interFoam solver	kumar2	OpenFOAM Running, Solving & CFD	3	June 30, 2006 18:26

July 11, 2015, 15:21		#2
David1 New Member België/Belgique Join Date: Jul 2015 Posts: 3 Rep Power: 10	Hi Roy, Did you eventually found your answer? Because I have the same question you have/had. Thank you!

September 27, 2019, 06:13		#4
AshwaniAssam Senior Member Ashwani Join Date: Sep 2013 Location: Hyderabad Posts: 154 Rep Power: 12	I found the answer here: PISO Details In sonicFoam the ".flux()" function returns the influence of all implicit terms on the flux. In this case everything is implicit, so you can assign it like that phi = pEqn.flux(); So, what I understand is that .flux() operator gives the flux terms of the equation which are handled implicitly.