OpenFOAM-11，OpenFOAM-v2306

李东岳

@星星星星晴是的。尤其是对于OpenFOAM-11这面。OpenFOAM这面几个大改动，一个是OpenFOAM-3.0，一个是OpenFOAM-11。这个OpenFOAM-11改的，改的相当大。

星星星星晴

@李东岳当时看到 9 把lagrangian改成那个样子，就已经停留在8了，现在这么快都到11了。看起来也是有好多修修补补

星星星星晴

@李东岳

刚才看了一下gay hub，发现真的是完全不一样了啊。。
还臭不要脸的添加了legacy
就这更新速度，生产队的驴也没有这么能干啊。。

李东岳

@星星星星晴是，拉格朗日库是他们那面很重视的项目，更新很大。以及一些多相流模型。他们都更新都很大。

不过有一些纯粹是代码上的改动，跟公式算法没啥关系。比如openfoam11，代码重写占据80%，算法我并不觉得会有那么天翻地覆，因为很多都很成熟了也没有大bug。

Weller在open day的时候就说过这个想法。同时由于目前的商业竞争关系。openfoam11就这么出来了。其实openfoam-dev很久就这样了。他们一直也没当做官方版发出来。

李东岳

下面这个代码是foamRun的主程序。这个在之前都是处理速度、压力、各种离散的代码，目前里面的核心代码都处理成了类，全部动过库来处理。

本质上算法都是一样的。只不过代码全部推翻了重写。


int main(int argc, char *argv[])
{
    argList::addOption
    (
        "solver",
        "name",
        "Solver name"
    );

    #include "setRootCase.H"
    #include "createTime.H"

    // Read the solverName from the optional solver entry in controlDict
    word solverName
    (
        runTime.controlDict().lookupOrDefault("solver", word::null)
    );

    // Optionally reset the solver name from the -solver command-line argument
    args.optionReadIfPresent("solver", solverName);

    // Check the solverName has been set
    if (solverName == word::null)
    {
        args.printUsage();

        FatalErrorIn(args.executable())
            << "solver not specified in the controlDict or on the command-line"
            << exit(FatalError);
    }
    else
    {
        // Load the solver library
        solver::load(solverName);
    }

    // Create the default single region mesh
    #include "createMesh.H"

    // Instantiate the selected solver
    autoPtr<solver> solverPtr(solver::New(solverName, mesh));
    solver& solver = solverPtr();

    // Create the outer PIMPLE loop and control structure
    pimpleSingleRegionControl pimple(solver.pimple);

    // Set the initial time-step
    setDeltaT(runTime, solver);

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

    Info<< nl << "Starting time loop\n" << endl;

    while (pimple.run(runTime))
    {
        solver.preSolve();

        // Adjust the time-step according to the solver maxDeltaT
        adjustDeltaT(runTime, solver);

        runTime++;

        Info<< "Time = " << runTime.userTimeName() << nl << endl;

        // PIMPLE corrector loop
        while (pimple.loop())
        {
            solver.moveMesh();
            solver.fvModels().correct();
            solver.prePredictor();
            solver.momentumPredictor();
            solver.thermophysicalPredictor();
            solver.pressureCorrector();
            solver.postCorrector();
        }

        solver.postSolve();

        runTime.write();

        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
            << nl << endl;
    }

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

    return 0;
}

wangfei9088

都瞅了瞅，基金会的版本改动大，学习成本高了。习惯了老版本，就算要换方向，必须用新版本，也会选择ESI的版本。个人观点。

李东岳

@wangfei9088 你说的有道理。OpenFOAM基金会那面犟得很。哎。这肯定会有一拨人跑ESI版本那面去了

wangfei9088

@李东岳李老师，那您的CFD课会紧跟潮流普及新版本不？

李东岳

@wangfei9088 是的。要的。要持续学习哎。要不然就被后浪拍死了。

wangfei9088

@李东岳 oh, cool!

CFD中文网

OpenFOAM-11，OpenFOAM-v2306