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,
|All times are GMT -4. The time now is 15:51.|