MULES, CMULES and IMULES 测试

[李东岳]

对比一下MULES（FCT）和vanLeer格式。下面是某点的M0，这个点均有界。

另外一点的M0，vanLeer越界了

整个流场的M0最小值，vanLeer越界了

[影川风]

东岳老师，能否提个不情之请，可以分享下算例吗？我目前也在看对流项的差值格式，想学习下这方面的知识。

[李东岳]

算例非常简单，就是1维对流问题。不过可以给你个思路制作求解器。改编scalarTransportFoam：

#include "fvCFD.H"
#include "simpleControl.H"
#include "MULEStest.H"
#include "FCT.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    #include "createTime.H"
    #include "createMesh.H"

    simpleControl simple(mesh);

    #include "createFields.H"

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    Info<< "\nCalculating scalar transport\n" << endl;

    #include "CourantNo.H"

    while (simple.loop())
    {
        Info<< "\n\nTime = " << runTime.timeName() << nl << endl;

		scalarField phiTout;

        while (simple.correctNonOrthogonal())
        {
			phiT = phi * fvc::interpolate(T);

			phiTout = phiT.ref();
			Info << "phiT before MULES " << phiTout << nl << nl;

			//MULEStest::explicitSolve(T, phi, phiT, 9, 0);
			FCT(T, phi, phiT, 9, 0);

            //fvScalarMatrix TEqn
            //(
            //    fvm::ddt(T)
            //  + fvc::div(phiT)
            //);
            //TEqn.relax();
            //TEqn.solve();
        }

		scalarField TOUTPUT = T.ref();

		phiTout = phiT.ref();

		Info << "\n\nT is " << TOUTPUT << nl;
		Info << "phiT after MULES " << phiTout << nl;

        runTime.write();
    }

    //Info<< "End\n" << endl;

    return 0;
}

通过MULES学习插值格式太复杂了

[影川风]

@李东岳 好的，非常感谢东岳老师:happy:

[asininno]

@李东岳 MULES(FCT)格式是什么意思，用mules 的vanleer 和不用的对比吗

[李东岳]

MULES(FCT)

MULES就是FCT的变种，区别在于限制器的求解和原始文献有点不同。使用FCT的时候，需要指定一种高阶格式，比如linear或者vanLeer

[yhdthu]

前辈，弱弱的问一句，这个分布图是用paraview做的么？像这种不同时刻的怎么做呢😂

[李东岳]

不，是gnuplot做的，目前在小岛等汽车，只能帮你到这里了:happy:

[yhdthu]

@李东岳 多谢前辈 我自己研究一下～

[金石为开]

@李东岳 老师 我也想知道测试格式算例是啥样的？您写的上面代码解标量方程那个Matrix代码是不是不该注释掉啊

[李东岳]

@金石为开
等回德国我回复你。MULES这个2018年应该可以写个解析出来。2016年就在计划中，但是有别的事情要搞。