[jason]

Hi,
I'm trying to use the PhaseLimitStabilization source in multiphaseEulerFoam to stop my alpha.water, alpha.oil and alpha.gas dropping too low and ultimately crashing the k-e model.

I currently have this format in my fvOptions file:


stabilize_epsilon.water
{
type symmTensorPhaseLimitStabilization;

field epsilon.water;
rate rLambda.water;
residualAlpha 1e-3;
}
stabilize_epsilon.oil
{
type symmTensorPhaseLimitStabilization;

field epsilon.oil;
rate rLambda.oil;
residualAlpha 1e-3;
}
stabilize_epsilon.gas
{
type symmTensorPhaseLimitStabilization;

field epsilon.gas;
rate rLambda.gas;
residualAlpha 1e-3;
}

stabilize_k.water
{
type symmTensorPhaseLimitStabilization;

field k.water;
rate rLambda.water;
residualAlpha 1e-3;
}
stabilize_k.oil
{
type symmTensorPhaseLimitStabilization;

field k.oil;
rate rLambda.oil;
residualAlpha 1e-3;
}
stabilize_k.gas
{
type symmTensorPhaseLimitStabilization;

field k.gas;
rate rLambda.gas;
residualAlpha 1e-3;
}

stabilize_U.water
{
type symmTensorPhaseLimitStabilization;

field U.water;
rate rLambda.water;
residualAlpha 1e-3;
}
stabilize_U.oil
{
type symmTensorPhaseLimitStabilization;

field U.oil;
rate rLambda.oil;
residualAlpha 1e-3;
}
stabilize_U.gas
{
type symmTensorPhaseLimitStabilization;

field U.gas;
rate rLambda.gas;
residualAlpha 1e-3;
}


The solver seems to give the correct response


Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_epsilon.water
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_epsilon.oil
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_epsilon.gas
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_k.water
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_k.oil
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_k.gas
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_U.water
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_U.oil
Selecting finite volume options model type symmTensorPhaseLimitStabilization
Source: stabilize_U.gas


My problem is still occurring, the alpha drives to very low values and crashes the k-e model, would anyone know if this is the correct way to do this? I've tested with and without and this doesn't seem to make a difference.

Thanks

Jason