OpenFOAM-11，OpenFOAM-v2306

[李东岳]

最近一星期。ESI以及基金会纷纷推出OpenFOAM-v2306以及OpenFOAM-v11。后者改动尤其巨大。求解器模块化，比如之前的icoFoam是个求解器，但目前被处理成了一个库。OpenFOAM-v11比之前，在多物理场耦合方面，更加友好。但更加增加了使用难度。

[星星星星晴]

所以以后查资料得看年份了，年份越久，可能越没啥鸟用

[李东岳]

@星星星星晴 是的。尤其是对于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!