[BenChaECN]

Hello everyone,


I am trying to code a new surface interpolation scheme for alpha in OpenFOAM.


For this new scheme, I need to calculate the gradient of alpha.


But I have a problem. Indeed, I use the existing blended scheme as starting point. But, the argumentat at the entrance is a GeometricField<Type, fvPatchField, volMesh>. And when I am trying the command fvc::grad() (I included fvc.H) it doesn’t work.


Does someone have a solution to this problem ?

[fumiya]

You might want to check the source code of the existing interpolation schemes that use the gradient operation such as linearUpwind:
https://github.com/OpenFOAM/OpenFOAM...s/linearUpwind

    tmp<fv::gradScheme<scalar>> gradScheme_
    (
        fv::gradScheme<scalar>::New
        (
            mesh,
            mesh.gradScheme(gradSchemeName_)
        )
    );

    for (direction cmpt = 0; cmpt < pTraits<Type>::nComponents; cmpt++)
    {
        tmp<volVectorField> tgradVf =
            gradScheme_().grad(vf.component(cmpt), gradSchemeName_);

        const volVectorField& gradVf = tgradVf();