[Voulet]

Hi there. I use openfoam V8.


I've my own solver with some custom volVectorField.


According to the banana method and https://cfd.direct/openfoam/user-guide/v8-fvSchemes/ for divergences schemes I can use some schemes with some limitors in order to make the divergence not to numerically explode. I'm interested with these schemes.


In my simulation the classical velocity advection terme use this scheme, which is fine:

div(phi,U)      bounded Gauss upwind;

When i copy/paste this scheme with my field :

 div((gradPt*1))  bounded Gauss upwind;

--> FOAM FATAL IO ERROR: 
unknown div scheme bounded

Valid div schemes are :
1(Gauss)

I also try

div((gradPt*1))  Gauss limitedLinear 1;

and I got

wrong token type - expected word, found on line 84 the doubleScalar 1

But I'm sure this is the good grammar since when I grep -ri Gauss I can found for exemple:

div(phi,k)      Gauss limitedLinear 1;

in combustion/fireFoam/LES/smallPoolFire2D/system/fvSchemes


What is the issue ? Are some schemes which are not available when explicitly calculating fvc::div(volVectorfield) instead of the more classical constructor fvc::div(surfacePhield,volVectorfield) ?


For information my volVectorField is defined like this :

 volVectorField gradPt(fvc::grad(p)+0.5*fvc::grad(U&U));


 
volVectorField directionGradPt(
    gradPt*dimensionedScalar(pow(dimLength,-1)*pow(dimTime,2),1)
);
forAll(directionGradPt,celli){
    directionGradPt[celli] = (mag(gradPt[celli]) == 0)? vector::zero : directionGradPt[celli]/mag(gradPt[celli]);
};

Thanks a lot.

[Voulet]

In fact, looking in the whole tutorial folder with :

tut
grep -ri div > tmp
sed -i '/fvSchemes/!d' tmp
sed -i '/divSchemes/d' tmp
sed -i '/default/d' tmp
sed -i '/,/d' tmp

I'm looking for all the schemes with div() withouth any comma inside in order to remove all the div(whatever,whatever) and have only the constructor i'm interested in.



I got a 264 lines file with containes 264 times the Gauss linear scheme according to grep -ri "Gauss linear;" tmp | wc -l


Thus why is there only this scheme available for a fvc::div() wich is not convective ? i.e not fvc::div(phi,thing) ?


How can i limit such a divergence using another scheme ?