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;
        }
    }
}
Example #2
0
    // apply zero-gradient condition for omega
    forAll(cornerWeights_, patchI)
    {
        if (!cornerWeights_[patchI].empty())
        {
            omegaWallFunctionFvPatchScalarField& opf = omegaPatch(patchI);

            opf == scalargpuField(omega0, opf.patch().faceCells());
        }
    }
Example #3
0
    // 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);
        }
    }