[TCS]

Hi All,

OpenFOAM 4.1

I am adding a scalar transport equation to the InterFOAM solver
$\frac{\partial c}{\partial t} + \nabla\cdot(c{\bf u})$
to track the progress of a scalar concentration field (c) which I define in setFields to be 0,1,2 initially in certain locations.

I have successfully added a standard equation in the InterFoam.C file:

fvScalarMatrix CEqn
        (
              fvm::ddt(c)
            + fvm::div(phi, c)
        );

CEqn.solve();

However, this is affected by a large amount of numerical diffusion.

I understand I need to use the MULES FCT scheme to try and solve this, as noted in Using MULES to solve transport equation with sharp interface but am unsure of how to do this.

I've tried adding the following to the alphaEqn file:

        tmp<surfaceScalarField> phiCure
        (
            fvc::flux
            (
                phi,
                c,
                alphaScheme
            )
          + fvc::flux
            (
               -fvc::flux(-phir, alpha2, alpharScheme),
                c,
                alpharScheme
            )
        );

together with

MULES::explicitSolve(c, phiCN, phiCure, 2, 0);

but this does not compile with wmake:
error: no matching function for call to ‘explicitSolve(Foam::volScalarField&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >&, int, int)’
among many errors.

Any help appreciated! Cheers.