Mass source term model is working incorrectly for condensation

[Johannes-Kind]

TwoPhaseFlow/src/phaseChange/phaseChangeModels/massSourceTermModel/hardtWondra/hardtWondra.C

Lines 148 to 155 in fefa227

	if (intPsiLiquid.value() > 1e-99)
	{
	Nl = intPsi0/intPsiLiquid;
	}
	if (intPsiVapor.value() > 1e-99)
	{
	Nv = intPsi0/intPsiVapor;
	}

TwoPhaseFlow/src/phaseChange/phaseChangeModels/massSourceTermModel/hardtWondraGasPhase/hardtWondraGasPhase.C

Lines 149 to 156 in fefa227

	if (intPsiLiquid.value() > 1e-99)
	{
	Nl = intPsi0/intPsiLiquid;
	}
	if (intPsiVapor.value() > 1e-99)
	{
	Nv = intPsi0/intPsiVapor;
	}

In both Hardt & Wondra models for phase change (see above), the scaling of the smeared mass source terms is not activated, if the total mass source is smaller than 1e-99 (i.e., if condensation occurs). The code should be changed to

if (mag(intPsiLiquid.value()) > 1e-99)
{
    Nl = intPsi0/intPsiLiquid;
}
if (mag(intPsiVapor.value()) > 1e-99)
{
    Nv = intPsi0/intPsiVapor;
} 

Now both evaporation and condensation can be modeled, even simultaneously in one system.