OpenFOAM中的时间步

李东岳

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

        #include "CourantNo.H"
        // Pressure-velocity PISO corrector
        {
            fvModels.correct();
           
            #include "UEqn.H"
            如果不进行动量预测 U都是当前时间步
            // --- PISO loop
            while (piso.correct())
            {
                #include "pEqn.H"
            }
            在这面之后U是下一个时间步的
        }

        viscosity->correct();
        turbulence->correct();

        runTime.write();

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

你那个我看不太懂，用上面的代码举例是这样

kcol

我想实现上述算法，我如果用

    volScalarField alpha10("alpha10", 2*alpha1-alpha1.oldTime());
    volScalarField alpha11("alpha11", 2*alpha1-0.5*alpha1.oldTime());

来表示算出来的结果就有问题，正常算一个周期的迭代步数只能算三分之二周期

李东岳

你写在了alpha.H的前面还是后面

kcol

写在组建方程的前面，开始时间步循环计算alpha的开头，

李东岳

开始时间步循环计算alpha的开头，

alpha方程之前的alpha表示t时刻已知的alpha，alpha方程求解之后表示t+dt时刻的alpha

kcol

这样说我的代码写的就没问题，但是实际计算的时候这样写就是不对，同时我要是写成
2*alpha1.oldTime()-alpha1.oldTime().oldTime()
这样计算结果就是没问题的，

李东岳

那就奇怪了，不知道咋回事。你这个方程式没问题的？

kcol

按照上述算法来执行，也就是第一二步写成如下形成：

    volScalarField alpha10("alpha10", 2*alpha1-alpha1.oldTime());
    volScalarField alpha11("alpha11", 2*alpha1-0.5*alpha1.oldTime());

是不对的，如果写成

    volScalarField alpha10("alpha10", 2*alpha1.oldTime()-alpha1.oldTime().oldTime());
    volScalarField alpha11("alpha11", 2*alpha1.oldTime()-0.5*alpha1.oldTime().oldTime());

算法就变成了就可以求解，只是求解精度比matlab上写的精度差一点，应该是第一二步有点问题。图片中的算法在matlab中实施起来是没问题的

李东岳

alpha的时间离散格式是什么

kcol

是向后差分，但是我没有用到ddt，是自己写的函数推进时间步

CFD中文网

OpenFOAM中的时间步