[rajinthss]

Hello all,

I am having some trouble understanding the implementation of the convection scheme of the alpha transport equation (in the alphaEqn.H file). In the interFoam/ VOF solver. The original code looks like this:

fvScalarMatrix alpha1Eqn
(
    (
        LTS
      ? fv::localEulerDdtScheme<scalar>(mesh).fvmDdt(alpha1)
      : fv::EulerDdtScheme<scalar>(mesh).fvmDdt(alpha1)
    )
  + fv::gaussConvectionScheme<scalar>(mesh, phiCN, upwind<scalar>(mesh, phiCN)).fvmDiv(phiCN, alpha1)

 ==
    Su + fvm::Sp(Sp + divU, alpha1)

So it seems that an upwind convection scheme has been used. If wanted to change this to a higher-order scheme, say SuperBee, how would I go about doing that? The documentation implies that it should be implemented as:

LimitedScheme< Type, Limiter, LimitFunc >

But I cannot seem to figure out the arguments for the limiter and limitfunc.

Any guidance would be much appreciated.

Thanks