Foam::tmp<Foam::volScalarField> Foam::radiation::cloudScatter::sigmaEff() const { tmp<volScalarField> tsigma ( new volScalarField ( IOobject ( "sigma", mesh_.time().timeName(), mesh_, IOobject::NO_READ, IOobject::NO_WRITE, false ), mesh_, dimensionedScalar("zero", dimless/dimLength, 0.0) ) ); forAll(cloudNames_, i) { const thermoCloud& tc ( mesh_.objectRegistry::lookupObject<thermoCloud>(cloudNames_[i]) ); tsigma() += tc.sigmap(); } return 3.0*tsigma; }
void Foam::viscosityModels::strainRateFunction::correct() { tmp<volScalarField> tsigma = strainRate(); const volScalarField& sigma = tsigma(); nu_.primitiveFieldRef() = strainRateFunction_->value(sigma()); volScalarField::Boundary& nuBf = nu_.boundaryFieldRef(); const volScalarField::Boundary& sigmaBf = sigma.boundaryField(); forAll(nuBf, patchi) { nuBf[patchi] = strainRateFunction_->value(sigmaBf[patchi]); }
Foam::tmp<Foam::volScalarField> Foam::surfaceTensionModels::temperatureDependent::sigma() const { tmp<volScalarField> tsigma ( volScalarField::New("sigma", mesh_, dimSigma) ); volScalarField& sigma = tsigma.ref(); const volScalarField& T = mesh_.lookupObject<volScalarField>(TName_); sigma.field() = sigma_->value(T.field()); volScalarField::Boundary& sigmaBf = sigma.boundaryFieldRef(); const volScalarField::Boundary& TBf = T.boundaryField(); forAll(sigmaBf, patchi) { sigmaBf[patchi] = sigma_->value(TBf[patchi]); }