各位大佬大家好,我最近准备在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"))
{}