@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就是一个对象了，不再识别为一个类的名称进而创建新的对象呢

@wangfei9088
谢谢老师，实在是醍醐灌顶，编译也通过了；

C++的类一直是一知半解，通过这次的改写理解进一步加深了；

另外老师我补充一下，如果p在pressureReference类下的对象命名为pressureReference，那么pf的对象怎么命名都通过不了，比如PRf，pressureRenferencef之类的，后来我尝试了这样的形式：

pressureReference PR(p, simple.dict());

mesh.setFluxRequired(p.name());

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

mesh.setFluxRequired(pf.name());

即两个对象分别命名为PR和PRf，就能通过了，这是为什么

另外老师，还有一个问题想请教一下：

如果我理解没错的话，我们这里命名的PR和PRf，是在pressureReference这个类下的两个对象对吗；
那么，之前p在pressureReference类下的对象命名为pressureReference后，pf对象创建不通过的问题，是由于构造函数名称和类名称相同导致的错误吗

我暂时的解决方法是在createFields.H写

pressureReference pressureReference(p,pf, simple.dict());
mesh.setFluxRequired(p.name());
mesh.setFluxRequired(pf.name());

因为有在$FOAM_SRC中看到类似的操作，把pcorr还是rho与p放在一起。
但还是恳求各位老师指点，若有相似经历或操作

向各位老师请教一个问题，我的求解器里有两组速度压力场，压力分别是p和pf，但是在求解器代码里我不知道怎么设置参考压力 ；

在最新的OF9中，设置参考压力的代码如下（以simpleFoam为例）：

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

mesh.setFluxRequired(p.name());

—createFields.H

pEqn.setReference
(
    pressureReference.refCell(),
    pressureReference.refValue()
);

—pEqn.H
如果我在createFields.H中再加一段

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

mesh.setFluxRequired(pf.name());

编译就会不通过，提示“ expected ';' before‘pressureReference’ ”
pressureReference的C文件我找到了，但读不太懂，想请教一下各位老师求解器代码里到底应该怎么设置

class pressureReference
{
    // Private Data

        //- Optional cell in which the reference pressure is set
        label refCell_;

        //- Optional pressure reference level
        scalar refValue_;


public:

    // Constructors

        //- Construct from the simple/pimple sub-dictionary
        //  supporting the pressure and derived pressure fields (e.g. p_rgh)
        pressureReference
        (
            const volScalarField& p,
            const volScalarField& pRef,
            const dictionary& dict,
            const bool pRefRequired = true
        );

        //- Construct from the simple/pimple sub-dictionary
        pressureReference
        (
            const volScalarField& p,
            const dictionary& dict,
            const bool pRefRequired = true
        );


    // Member Functions

        //- Return the cell in which the reference pressure is set
        inline label refCell() const;

        //- Return the pressure reference level
        inline scalar refValue() const;
};