[pablodecastillo]

Hello,

Why in alphaEqn.H is it not added alpha1 = max(min(alpha1, scalar(1)), scalar(0)); i mean this artificial bounded?

Thanks

[Phicau]

Hi

MULES solver ensures boundedness, that is why the last 2 inputs are 0 and 1, the bounds of alpha1. This is achieved by the use of limiters.

The final result may show values greater than 1 or lower than 0, but with differences of the order of 10^-5, which is negligible.

Regards

[santiagomarquezd]

Using an artificial bounding of this kinds leads to a non-phase-preserving algorithm due it 'cuts' with having in account the conservative properties of the alpha equation. In other words you will lose or gain mass.

Regards.

[pablodecastillo]

But looking for steady state loss or gain mass is not going to be relevant or yes?

Thanks

[santiagomarquezd]

Hhmm, I don't know, but since the solver is unsteady probably it will blow up before reach the steady state due mass preservation problems.

Regards.