[zhulianhua]

Hi! all,

I'm reading the source code of the van Leer limiter based surface interpolation scheme. I think I need your help to understand some of the code below:

template<class LimiterFunc>
class vanLeerLimiter
:
    public LimiterFunc
{

public:

    vanLeerLimiter(Istream&)
    {}
    scalar limiter
    (
        const scalar,
        const scalar faceFlux,
        const typename LimiterFunc::phiType& phiP,
        const typename LimiterFunc::phiType& phiN,
        const typename LimiterFunc::gradPhiType& gradcP,
        const typename LimiterFunc::gradPhiType& gradcN,
        const vector& d
    ) const
    {
        scalar r = LimiterFunc::r
        (
            faceFlux, phiP, phiN, gradcP, gradcN, d
        );

        return (r + mag(r))/(1 + mag(r));
    }
};

Where can I find the source code for the definition of function r() and the phiType/gradPhiType?

Thanks,

Lianhua

[zhulianhua]

The code in the post is from here:

https://github.com/OpenFOAM/OpenFOAM...Leer/vanLeer.H

Originally Posted by zhulianhua

Hi! all,

I'm reading the source code of the van Leer limiter based surface interpolation scheme. I think I need your help to understand some of the code below:

template<class LimiterFunc>
class vanLeerLimiter
:
    public LimiterFunc
{

public:

    vanLeerLimiter(Istream&)
    {}
    scalar limiter
    (
        const scalar,
        const scalar faceFlux,
        const typename LimiterFunc::phiType& phiP,
        const typename LimiterFunc::phiType& phiN,
        const typename LimiterFunc::gradPhiType& gradcP,
        const typename LimiterFunc::gradPhiType& gradcN,
        const vector& d
    ) const
    {
        scalar r = LimiterFunc::r
        (
            faceFlux, phiP, phiN, gradcP, gradcN, d
        );

        return (r + mag(r))/(1 + mag(r));
    }
};

Where can I find the source code for the definition of function r() and the phiType/gradPhiType?

Thanks,

Lianhua

[ngj]

Hallo Lianhua,

Note that the file you have shown is a templated class in the type name "LimiterFunc", why you cannot directly find that class. I believe that this is a template in e.g. the NVDTVD and NVDTVDV definitions, which you will find here:

https://github.com/OpenFOAM/OpenFOAM...eme/NVDVTVDV.H

and

https://github.com/OpenFOAM/OpenFOAM...cheme/NVDTVD.H

These are then invoked through the macros in e.g. vanLeer.C. I am not 100% sure that this is the correct understanding of the source code, as it borders my current understanding of the C++ language. I hope, nonetheless, that this will help you on your way.

Kind regards,

Niels

[zhulianhua]

Hello, Neils

Thanks for your reply. It really helps me. Yes, I have never seen such C++ usage on textbooks and web -- The base class is a template! I even don't know how to google this kind of usage.

Best,

Lianhua

[ngj]

Yes, I know. It is pretty advanced stuff

Kind regards,

Niels

[yourvahid]

I wanted to add a new limiter scheme such as "stacs" in openFOAM. I've already written the codes, however, I dont know how shall i compile and add it to openFOAM.
please note that I'm new in openFOAM environment and the interface is unfamiliar for me.
thanks a million.

[Jibran]

Originally Posted by yourvahid

I wanted to add a new limiter scheme such as "stacs" in openFOAM. I've already written the codes, however, I dont know how shall i compile and add it to openFOAM.
please note that I'm new in openFOAM environment and the interface is unfamiliar for me.
thanks a million.

Hi,

I suppose you have already added the .C and .H files in the relevant folder (for ex. \src\finiteVolume\interpolation\surfaceInterpolati on....). Then edit the 'files' and 'options' files in the 'Make' folder in \src\finiteVolume. Then compile it with 'wclean libso && wmake libso'. Finally make sure that your solver picks up the newly compiled library.

For more information please read the chapter 3.2 of the OF User Manual

Jibran