initial commit for version 1.6.x patch release

[OpenFOAM-1.6.x.git] / applications / solvers / multiphase / settlingFoam / plasticViscosity.H

volScalarField plasticViscosity
(
    const dimensionedScalar& plasticViscosityCoeff,
    const dimensionedScalar& plasticViscosityExponent,
    const volScalarField& alpha
)
{
    return
    (
         plasticViscosityCoeff*
         (
             pow(10.0, plasticViscosityExponent*alpha + SMALL) - scalar(1)
         )
    );
}