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

@李东岳

老师，我的求解器在计算时设置为层流，然后湍流的方程和粘性在求解器中解算。不知道是不是因为设置为层流的原因，wallshearstress在后处理的时候提示unable to find turbulence model in the database，所以没法用postProcess算

各位老师打扰了，本人才疏学浅学艺不精，最近在后处理，花了好久终于看懂了网上解释壁面剪切应力的帖子。
由于使用的湍流模型得到的湍流粘性系数为张量，在求解器中编写后处理程序时 ，遇到了不知道该如何输出应力在边界上的值的问题：

先前看到，在OpenFOAM源码的wallShearStress.C中，壁面剪切应力由以下方法求得：

wallShearStressBf[patchi] = (-Sfp/magSfp) & tau.boundaryField()[patchi];

其中的tau，其实就是有效雷诺应力：

tau = model.devSigma();

但由于我的后处理程序植入在solver中（直接写bool不太会），因此自己在createFields里定义了一个应力张量场stress，然后把层流的粘性应力：

dimensionedScalar("nu", dimViscosity, 1.388e-5) * dev(twosymm(fvc::grad(U)))

和湍流粘性应力相加得到，经检验可以输出内部场（通过stress[cellID]）的应力值。

正准备和法向单位向量相乘，发现输不出来边界上的stress.boundaryField()[faceID]

我知道我这样的处理肯定是过于简单和理想化了（简直就是瞎搞，瞎搞都算不上）只能盼望各位老师赐教了，该如何得到边界上的值呢

@stardust
请问楼主写完了嘛，如果可以的话可不可以贴上来让大家学习一下 （最近正好在做一个类似的东西）

@cccrrryyy
当然老师您可能会问为什么非定常问题还要用稳态求解器，这个的话可能主要是这几点原因：

一是原有的经典模型如k-epsilon等对这些问题可以用稳态算法进行计算，通过模型理论中的平均方法得到稳态解（相关问题没有做足够的调研，应该也有使用瞬态求解器的）；

二是我现在使用的模型的提出者基本上都使用稳态算法进行求解的；

三是我初步使用了瞬态求解器进行计算，但不稳定，会发散（当然也有可能是算法和设置有问题导致的）

@cccrrryyy
老师说的是，我没有说太清楚。

毕竟所有的实际问题在严格的意义上都不是定常的。如果广义上这个流场的脉动比较小，那么理论上我们用一个稳态的求解器的话，是可以算出来一个残差很小的收敛解的，老师您看我这样理解是对的吗。

然而同时对于一些问题，虽然实际上并不能认为是一个严格的定常问题，但是我们现在使用的很多湍流模型用稳态求解器计算却能得到一个稳态的解，比如平面射流。

那么这个时候，对这样的非定常问题，如果一个湍流模型在稳态求解器的基础上，算得的结果呈现出了涡的发展变化过程，老师您是怎样看这个结果的呢

想请教各位大佬一个有关稳态求解器与湍流的问题

就是理论上，稳态的求解器应该算出来一个最终残差很小的稳定解对吗

但是如果使用的湍流模型说，自己会解算出湍流的拟序结构（according to 它的提出者），所以使用稳态求解器算出来的解，会呈现出一个脉动的特征（比如涡的变化发展等过程），以至于没有一个收敛的解...

请问各位大佬是怎么看这样的结果的呢，这样是可以理解或接受的吗，又或请问有没有相关的资料可以学习一下这方面内容的呢（本人菜鸡，跪求指教）

@李东岳 老师，公式中的p是坐标，然后是共顶点的两条边构成的向量叉乘吗？

如果是这样，叉乘不是两边的模长相乘，乘以两边夹角的正弦么，就相当于底乘高了，再乘上0.5就是面积了？（基于首句的假设）

$0.5\ast \vec{l_{AB}} \times \vec{l_{AC}} = 0.5 \ast l_{AB} \ast l_{AC} \ast \mathbf{sinA} = 0.5 \ast l_{AB} \ast H = S$

不知道是不是这样，求各位大佬轻喷

@李东岳
老师，请问对于定义在体心上的物理量，这里correct的边界场，是指边界上的一层网格吗

@wangfei9088
请问老师，有没有可能，是因为我们将p在pressureReference类下的对象命名为pressureReference后，再写pf时：

pressureReference pressureReference(pf, simple.dict());

编译器会认为这行第一个pressureReference就是一个对象了，不再识别为一个类的名称进而创建新的对象呢