CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam (http://www.cfd-online.com/Forums/openfoam-programming-development/115569-how-could-i-modify-peqn-h-adding-ajustphi-twophaseeulerfoam.html)

cheng1988sjtu April 2, 2013 11:32

How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam
 
Hi Foamers,

When I use zeroGradient pressure condition at the wall, I could see pressure sigularity at the pReff(), and I searched in the forum, and found that that's due to that fact that there is no adjustPhi after the buoyancy term.

So, I added the following lines before phi = alphaf*phia + betaf*phib;

adjustPhi(phia, Ua, p);
adjustPhi(phib, Ub, p);

however, after I added this lines, I could see that my Dispersed volume fraction is decreasing though slowly.

Could anyone give me a hint on how adjustPhi() should work?

Thanks

Zhen Cheng

cheng1988sjtu April 2, 2013 15:36

Never mind, I figured out that the mass reduction is not due to the adjustPhi, and FYI, this approach do eliminate the pressure singularity at the zeroGradient boundary (p).

Charlie

chegdan April 3, 2013 11:05

Greetings Zhen Cheng,

Just curious about your findings,

Quote:

Never mind, I figured out that the mass reduction is not due to the adjustPh
What was the source of mass reduction in your simulations? Just in case others run into the same issue.

Dan

Popmech March 3, 2014 10:13

Hi, I am using a modified version of the twoPhaseEulerFoam solver to model a mixing vessel with settling particles. I have a zero-gradient pressure condition on all surfaces including the free-surface which is represented by a symmertryPlane type boundary condition. I suffer from a similar problem with a singularity occurring at the pReffCell in the simulation which causes the volume fractions to be solved incorrectly. I have tried implementing the adjustPhi call as described in this thread, though this causes an error stating there is no outflow boundary to adjust. Does someone know how to work around this problem?

Kind Regards, JVZ


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