CFD中文网

是没添加fvc的头文件吗？

你有在etc\caseDicts\postProcessing\fields下添加文件吗

怎么用openfoam提供的工具实现我不知道，我是利用paraview提取出乱序的数据后，利用python脚本重新排序，再作图

@wenke
不知道你的问题解决了没，你可以看下，可视化的计算数据应该在那个列表的下面

@宝丁 我的意思是比如case 1背压是p1，case 2背压是p2，case 1计算收敛后，将case 1的结果mapFields到case 2，再进行计算，case 2的边界条件是可以修改的

mapFields如果不添加-consistent，应该就不会映射边界条件，可以先修改边界条件，再mapFields；
也可以mapFields后，直接修改边界条件

@wwzhao 在 关于tmp<volScalarField>&用法的疑问 中说：

tmp是智能指针，不应该用引用的方式返回。

请问是指这句

const tmp<volScalarField>& tmuEff = turbModel.muEff();

还是指这句

const volScalarField& muEff = tmuEff();

我是照着下面写的，把三句并成了两句

const volVectorField& U = lookupObject<volVectorField>(fieldName_);
const tmp<volTensorField> tgradU(fvc::grad(U));
const volTensorField& gradU = tgradU();

@东岳
那没用tmp时，有时出错可能是内存占完了？
我是在functionobject里用，只是后处理时算一下，倒是不太在意速度，只求不要时不时出错。

在编写涉及到muEff的后处理函数时发现，使用下面方法得到muEff时，可以编译，但是计算结果是错的，或者浮点数溢出

        const volScalarField& muEff = turbModel.muEff();

而改成以下形式时，计算结果是正确的

        const tmp<volScalarField>& tmuEff = turbModel.muEff();
        const volScalarField& muEff = tmuEff();

而对于其他变量，如T，U等，两种形式的写法没发现对计算结果有影响，请问tmp<volScalarField>&有什么意义吗？为什么有的变量要用这种写法？谢谢！

@东岳
另一个alphatJayatillekeWallFunction和Fluent中的壁面函数一样，根据$y^+$使用不同的公式计算，打算试一下这个

标准壁面函数是下面这样的

这样$\alpha_t$在粘性层应该是0，在惯性层应该用公式（17.109）计算，但是alphatWallFunction在边界处给定$\alpha_t = \frac{\mu_t}{Pr_t}$，这应该是流场内部的计算方法吧，这样的话壁面$y^+$应该放到多大？

@队长别开枪
十分感谢您的回复，我之前主要是对

SIMPLEC改进了通量修正方法，没有为压力设置松弛因子。

有疑问，因为在rhoSimpleFoam的pEqn.H与pcEqn.H的代码中都有p.relax();，结合您最新的回复，OpenFOAM中SIMPLEC应该也可以在fields中设置压力松弛。

另外，请问OpenFOAM中哪个求解器支持SIMPLER，没有找到。

@d-benjamin

我们在修改求解器的时候是否还需要删除该依赖文件，如果需要，它在哪个文件夹中。

我使用of4、5，修改求解器时，直接复制原求解器文件夹、重命名、修改求解器及make文件，然后直接编译就行了

@izumi

看了一下rhoSimpleFoam的代码，发现transonic循环中才有以下压力方程松弛的设置，因此调节equations中的压力松弛因子应该只对transonic的情况有用

// Relax the pressure equation to maintain diagonal dominance
pEqn.relax();

@队长别开枪

你好，我看rhoSimpleFoam的pEqn.H与pcEqn.H中均有

// Explicitly relax pressure for momentum corrector
p.relax();

请问这个是不是压力松弛？

@东岳

这几天在算复杂模型，发现对于我的算例初场很差时，将equations中压力松弛因子先调低再调到1是有必要的，可以保证前期不发散，后期time step error的sum local持续降低。
值得注意的是我的算例中包含10个圆通道，采用相同的fvSolution时（equations中压力松弛因子为1），8个通道的sum可以降到0.001，2个通道的sum维持在1左右，但是检查进出口质量流量相差不大，不知道为什么。
这个可以当作一个trick吧，如果time step error降不下去，可以试着调一下equations中的压力松弛因子。

@东岳
time step error指的是sum local；
phi=0.00751（0.00779）中，0.00751是出口质量流量，括号中的0.00779是进口质量流量，边界条件是给定的进口质量流量为0.00779

使用rhoSimpleFoam计算圆管内流动时，发现收敛时time step error很大，进出口质量流量也相差较大，修改网格、数值格式、边界条件均无效，后对照tutorial，发现问题出在压力松弛因子。原本为避免发散，将松弛因子均设为了0.3，然后逐渐增加至0.7~0.9，这样设置就会出现上述问题，参照tutorial将压力松弛给到1，此时进出口流量守恒。
针对压力松弛因子，从第5200步开始，做了如下测试：

• 将field中的压力松弛因子设置为1，其余0.9，5951收敛：time step error=118.87，phi=0.00751（0.00779）

• 将equations中的压力松弛因子设置为1，其余0.9，5275收敛：time step error=0.00196268，phi=0.00779（0.00779）

• 将equations与field中的压力松弛因子均设置为1，其余0.9，5268收敛：time step error=0.00283641，phi=0.00779（0.00779）

• 将equations与field中的压力松弛因子均设置为0.9，其余0.9，5786收敛：time step error=130.74，phi=0.00748（0.00779）

可以看出：equations中的压力松弛因子的设定是进出口量守恒的关键，请问这是为什么？是否所有计算中都要将equations中的压力松弛因子最终给到1？ 谢谢！

@forward
应该可以独立划分网格，我做过简单测试，两个固体区域夹着一个流体区域，界面处固体网格与流体网格疏密不同，可以计算，没有比较过精度，不知道有无影响。

@东岳
谢谢！用Surface Vectors filter就直接画出来了