How could I modify the pEqn.H by adding ajustPhi() in twoPhaseEulerFoam
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?
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).
Greetings Zhen Cheng,
Just curious about your findings,
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 11:51.|