36行Burgers方程的OpenFOAM源代码

[李东岳]

#include "fvCFD.H"

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

    volVectorField U
    (
        IOobject
        (
            "U",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );
    
    while (runTime.loop())
    {
        Info<< "Time = " << runTime.userTimeName() << nl << endl;
        surfaceScalarField phi(fvc::flux(U));

        solve
        (
            fvm::ddt(U) 
          + 0.5*fvm::div(phi, U, "div(phi,U)")
        );

        runTime.write();
    }
    return 0;
}

测试算例： burgers.tar.xz

相应的结果：一个好一个坏

[JamieZ]

李老师，请问对流项的系数0.5有何说法？

[李东岳]

写成的守恒形式。$u\frac{\p u}{\p t}\neq\frac{\p uu}{\p t}$

[JamieZ]

明白，不过这样似乎只对一维情况成立吧？如果对二维或三维问题，应该通过

 - fvm::Sp(fvc::div(phi), U)

转换成非守恒形式？不知道是否正确~

[李东岳]

你说的对。我那个只是在1D成立。

[JamieZ]

感谢李老师

[李东岳]

http://dyfluid.com/class.html#burgers

针对你说的1D情况。我修改了一下我写的，防止歧义。同时课堂算例增加一个，这个更具有挑战性