[bala_gk1988]

Hi guys,
I am trying to create a convective boundary condition for my LES jet flow problem. Although other user written boundary condition are available in this forum, I am trying to implement a BC which uses the average velocity of the patch as the convective velocity.

I decided to modify the advective boundary condition, which uses the local flux to calculate the convective velocity. I just modified the portion which calculates the advection speed(convective velocity). The below code is compiling without errors. Is it correct.

Foam::convectiveFvPatchField<Type>::advectionSpeed () const
{
const surfaceScalarField& phi =
this->db().objectRegistry::template lookupObject<surfaceScalarField>
(phiName_);

fvsPatchField<scalar> phip =
this->patch().template lookupPatchField<surfaceScalarField, scalar>
(
phiName_
);
fvsPatchField<scalar> phipt(phip);
phipt = Foam::gSum(phip); //****** Line one ******//

const scalarField& areap = this->patch().magSf();
scalarField area(areap.size(),Foam::gSum(areap)); //*** Line two ***/

return phipt/area;
}


In line one can I equate a surface scalar patch field to a surface scalar ?
Is line two correct ?

[calf.Z]

Do you test your BC successfully?