initial commit for version 1.6.x patch release

[OpenFOAM-1.6.x.git] / applications / solvers / multiphase / interFoam / UBlendingFactor.H

    surfaceScalarField alpha1f = fvc::interpolate(alpha1);
    surfaceScalarField UBlendingFactor
    (
        "UBlendingFactor",
        sqrt(max(min(4*alpha1f*(1.0 - alpha1f), 1.0), 0.0))
    );