applications/solvers/multiphase/settlingFoam/yieldStress.H

   1 volScalarField yieldStress
   2 (
   3     const dimensionedScalar& yieldStressCoeff,
   4     const dimensionedScalar& yieldStressExponent,
   5     const dimensionedScalar& yieldStressOffset,
   6     const volScalarField& alpha
   7 )
   8 {
   9     tmp<volScalarField> tfld
  10     (
  11         yieldStressCoeff*
  12         (
  13             pow(10.0, yieldStressExponent*(alpha + yieldStressOffset))
  14           - pow(10.0, yieldStressExponent*yieldStressOffset)
  15         )
  16     );
  17
  18     return tfld();
  19 }