[Fauster]

Dear Foamers,

I am working on the implementation of a new limited scheme in OpenFOAM. I am not very good in C++ and I almost try to work by adpating existing files of OpenFOAM. In my case I started with the QUICK.H file located in

/opt/OpenFOAM-v1806/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/QUICK

The limiter is calculated using the centre and gradient values of neighboring node P and N of the face f.

scalar phiCD = cdWeight*phiP + (1 - cdWeight)*phiN;

scalar phiU, phif;

if (faceFlux > 0)
{
phiU = phiP;
phif = 0.5*(phiCD + phiP + (1 - cdWeight)*(d & gradcP));
}
else
{
phiU = phiN;
phif = 0.5*(phiCD + phiN - cdWeight*(d & gradcN));
}
// Calculate the effective limiter for the QUICK interpolation
scalar QLimiter = (phif - phiU)/stabilise(phiCD - phiU, SMALL);
// Limit the limiter between upwind and downwind
return max(min(QLimiter, 2), 0);

I would like to add a Courant number dependency in the calculation of the returned limiter. If I am correct I need to call :

const fvMesh& mesh = this->mesh();

because I will need to access to the time step mesh.time().deltaT(), the cell volume mesh.V() and the flux magnitude fvc::surfaceSum(mag(this->faceFlux_))()(). I also guess that I need to add #include "fvMesh.H" and #include "surfaceFields.H".

The problem is that I obtain various errors when I try to use
const fvMesh& mesh = this->mesh() in the scalar limiter function :

‘const class Foam::mHRIC2Limiter<Foam::NVDTVD>’ has no member named ‘mesh’
const fvMesh& mesh = this->mesh();

Is there any smart way to do such a thing without changing all the limiter function and/or the NVDTVD.H file ?
Thanks in advance for helping me !

[mAlletto]

Maybe the courant number blended scheme can give you some hints of how to access the variables you want:

https://www.openfoam.com/documentati...8H_source.html