[randon]

Hello all,

I am working on cell center to face center implementation in an inhouse code. I am trying out different schemes. Currently, a limiter is defined in the similar approach as limitedLinear.

where is the face center value. Then the limiter involves a cell gradient calculation in the process. So, for this gradient, I am using Gauss linear algorithm.

1. My primary question is, how to compute the in the Gauss linear algorithm at the boundary faces ? Should I simply use the owner cell center value or use weights from makeWeights() function. If the face is internal computing weights based on owner and cell centre values is clear.

2. If the weights are to be used and facei is at boundary, owner[facei] exists but not neighbour[facei]. Then what is the value of C[neighbour[facei]] ?

Reference code : https://github.com/OpenFOAM/OpenFOAM...nterpolation.C

 
    scalarField& w = weights.primitiveFieldRef();

    forAll(owner, facei)
    {
        scalar SfdOwn = mag(Sf[facei] & (Cf[facei] - C[owner[facei]]));
        scalar SfdNei = mag(Sf[facei] & (C[neighbour[face]] - Cf[facei]));
        w[facei] = SfdNei/(SfdOwn + SfdNei);
    }

    surfaceScalarField::Boundary& wBf = weights.boundaryFieldRef();

    forAll(mesh_.boundary(), patchi)
    {
        mesh_.boundary()[patchi].makeWeights(wBf[patchi]);
    }

3. Also, could you please help me understand the makeWeights function implementaion at boundary ? I am improving my C++ and OpenFOAM programming skills. Please don't mind if the question is preliminary.

I can provide any further information that is needed. Based on your suggestions, I can also implement the value for which uses linear scheme. Thanks for your help.

Best regards.