在rhoSimpleFoam基础上加入一组U、p、T场

evensun

各位大佬大家好，我最近准备在rhoSimpleFoam的基础上加入自己的算法，算法中比较关键的点是有另一组速度、压力和温度场（暂且叫它们Uf，pf，Tf）。

之前在simpleFoam中顺利地植入了，但是在可压缩这里我发现：thermo.rho()这些函数在调用时，读取的物理量是求解器原来的U、p，如果我再另外定义Uf，pf的话，好像不能直接使用thermo.rho()函数进行另一组速度压力场的密度修正？

附上一段src/thermophysicalModels/basic/fluidThermo/fluidThermo.C中的代码，麻烦大佬们看下是不是这样的

Foam::fluidThermo::implementation::implementation
(
    const fvMesh& mesh,
    const word& phaseName
)
:
    p_(lookupOrConstruct(mesh, "p"))
{}


Foam::fluidThermo::implementation::implementation
(
    const fvMesh& mesh,
    const dictionary& dict,
    const word& phaseName
)
:
    p_(lookupOrConstruct(mesh, "p"))
{}