[Shibi]

Hello to all,

I am writing a codedMixed boundary condition to emulate the behavior of a boundary implemented into OpenFOAM.

As an example let us consider the following:

rightWall
     {
        type            codedMixed;

        refValue       uniform (0 0 0);
        refGradient      uniform (0 0 0);
        valueFraction   uniform 1;

        name            velocityMixed;   

        code
        #{

               refValue() = some expression
               refGrad() = some expression
               valueFraction() = some expression
         #};
}

Afterwards the code in the boundary of OpenFOAM overloads the = operator as:

void Foam::boundaryFvPatchVectorField
::operator=
(
    const fvPatchField<vector>& pvf
)
{
    tmp<vectorField> n = patch().nf();
 
    fvPatchField<vector>::operator=
    (
        valueFraction()*(n()*(n() & pvf))
      + (1 - valueFraction())*pvf
    );
}

How can I input this into the codedMixed boundary condition? What is the value of pvf? Is it the value of the mixedBoundary? Is overleading the operator even possible?



Something like: [?]

(
                     valueFraction()*refValue()
                      +
                       (1.0 - valueFraction())*
                       (
                        this->patchInternalField()
                        + refGrad()/this->patch().deltaCoeffs()
                        )
                        );

Would really appreciate some help with this.

Best Regards