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



  • 刚刚开始接触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里面。他是个非常普适性的概念。



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


登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待