[kaifu]

Hi Foamers,

In alphaEqn.H (like in the solver bubbleFoam or twoPhaseEulerFoam) I have to make sure alpha is strictly in [0, 0.523] after each

alphaEqn.solve();

So I simply make

alpha = max(0,alpha);
alpha = min(0.523,alpha);

in the code.
So far I have realized it may break the mass continuity if I do so. So I considered if I should create a "packingLimiter" as it is in solver twoPhaseEulerFoam. Thanks

// Kai