void omegaWallFunctionFvPatchScalarField::setMaster() { if (master_ != -1) { return; } const volScalarField& omega = static_cast<const volScalarField&>(this->dimensionedInternalField()); const volScalarField::GeometricBoundaryField& bf = omega.boundaryField(); label master = -1; forAll(bf, patchI) { if (isA<omegaWallFunctionFvPatchScalarField>(bf[patchI])) { omegaWallFunctionFvPatchScalarField& opf = omegaPatch(patchI); if (master == -1) { master = patchI; } opf.master() = master; } } }
// apply zero-gradient condition for omega forAll(cornerWeights_, patchI) { if (!cornerWeights_[patchI].empty()) { omegaWallFunctionFvPatchScalarField& opf = omegaPatch(patchI); opf == scalargpuField(omega0, opf.patch().faceCells()); } }
// accumulate all of the G and omega contributions forAll(cornerWeights_, patchI) { if (!cornerWeights_[patchI].empty()) { omegaWallFunctionFvPatchScalarField& opf = omegaPatch(patchI); const gpuList<scalar>& w = cornerWeights_[patchI]; opf.calculate(turbulence, w, opf.patch(), G0, omega0); } }