关于可压缩高速流求解sonicFoam求解器

[vivian]

刚刚开始接触openfoam，要做高超声速的模拟，看了求解器，认为sonicFoam比较合适，看这个求解器有很多疑问，向大神们请教。

1. sonicFoam求解器

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

        #include "compressibleCourantNo.H"

        #include "rhoEqn.H"

        // --- Pressure-velocity PIMPLE corrector loop
        while (pimple.loop())
        {
            #include "UEqn.H"
            #include "EEqn.H"

            // --- Pressure corrector loop
            while (pimple.correct())
            {
                #include "pEqn.H"
            }

            if (pimple.turbCorr())
            {
                turbulence->correct();
            }
        }

        rho = thermo.rho();

        runTime.write();

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

Q1. 基于压力基求解，我的理解是"rhoEqn.H"是基于质量守恒方程求得密度，这个方程加fvOptions是基于源项中的密度修正吗？主要来源是哪里呢？

{
    fvScalarMatrix rhoEqn
    (
        fvm::ddt(rho)
      + fvc::div(phi)
      ==
        fvOptions(rho)
    );

Q2. 求解能量方程，为什么有这一项fvc::ddt(rho, K) + fvc::div(phi, K)，求解的能量方程是什么样子的那？感觉跟平时推导的能量方程有很多差距？
能指点一下这个求解器的推导吗？不胜感激。

[李东岳]

粗略看了一下，针对Q2，希望这个有帮助：CFD中的能量方程
fvOptions主要来源于多重参考系，多孔介质，附加的能量产生等，OpenFOAM把所有的对不连续的变化都放在了fvOptions里面。他是个非常普适性的概念。

[vivian]

谢谢东岳大神的回复，看了您推荐的，感觉有点明白了，还在挣扎中。