Foam::tmp<Foam::volScalarField> Foam::laminarModels::generalizedNewtonianViscosityModels::strainRateFunction:: nu ( const volScalarField& nu0, const volScalarField& strainRate ) const { tmp<volScalarField> tnu ( volScalarField::New ( IOobject::groupName(type() + ":nu", nu0.group()), nu0.mesh(), dimensionedScalar(dimViscosity, 0) ) ); tnu.ref().primitiveFieldRef() = strainRateFunction_->value(strainRate); volScalarField::Boundary& nuBf = tnu.ref().boundaryFieldRef(); const volScalarField::Boundary& sigmaBf = strainRate.boundaryField(); forAll(nuBf, patchi) { nuBf[patchi] = strainRateFunction_->value(sigmaBf[patchi]); }