CFD中文网

@浪迹天大 谢啦 什么是CFDNotepad啊，给我解答问题还要送给我吗，哈哈，太好了

老帖重新挤上来了，有办法解决么？

@程迪 回答地太清晰了！非常感谢！

@李东岳 多谢啊！我去看看

@cfd-china
在非正交修正的循环内即，

fvScalarMatrix pEqn
                (
                    fvm::laplacian(rAU, p) == fvc::div(phiHbyA)//公式(26)，有关div(phiHbyA)请参考：更新3
                );
   pEqn.setReference(pRefCell, pRefValue);
                pEqn.solve();//公式(26)

反复计算pEqn，rAU与phiHbyA 都不变，每次更新的p也就不变啊，看不出来非正交修正体现在哪个地方，请问这里怎么理解？

请问fluent里面，如何将我已经time-average 的速度场和压力场作为初始值初始化（initialize）流场？

是啊，苦恼中。。。

曾经试过cutplane， 但是，该方法截取的某一平面是无限大的，而我只想要某一个有限矩形范围内的数据，多谢啦:happy:

@cfd-china 多谢啦，我试试！

多谢啊，我是小白，不知道在哪个位置加入这个代码，以icoFoam为例，能给我指点一二吗？

@cfd-china 用途这么广，就是有点疑惑。是不是所有的湍流粘度的表达都是用流场的一些不变量表示（LES的Smargorinsky-Lily亚格子粘性需要应变张量的模表示），与坐标系无关？这样就可以直接用旋转坐标系下的相对量直接得到湍流粘度。

边界条件、实验与湍流模型没问题的话，再增加网格。你的网格貌似还没有收敛，网格数量少与实验一致很可能就是碰巧的结果。

openfoam，fluent中都有将参考系固定在旋转机械上的求解器。湍流模型包括大涡模拟里面的很多经验参数都是针对静止坐标系的，把他们用于旋转坐标系下合适吗？

在N-S微分方程中加入狄拉克函数（脉冲函数）源项，在OF中如何实现？请教一下各位有什么思路吗？这个函数怎么用编程语言实现？

@李东岳 想明白了，多谢东岳大神！

@李东岳 多谢啦，从公式14，总感觉系数矩阵是对角阵，如下式
2u1 =-u2 （此时的u2及其系数都是当前已知时间步的值）
2u2 =-u1-u3（此时的u1、u3及其系数都是当前已知时间步的值）
2u3=-u2 (同上)
提取矩阵
A=
2 0 0
0 2 0
0 0 2
b={-u2 -u1-u3 -u2}'
因为“仔细观察公式(14)，除了$u^r_p$均为已知。”可能我看的不够仔细，我再仔细揣摩一下，多谢大神解答。

@李东岳 也就是说只需要知道k和epsilon或omega，都可以求解得到湍流粘度，边界上也是如此吗？那还需要设置边界上的初始值吗？也就是0时刻

@李东岳 临点不都是已知的量吗（用了已知时间步的量）？正如我第二个问题所问。在“icoFoam解析”中，公式14下面一行字也写了：“仔细观察公式(14)，除了$u^r_p$均为已知。”

@李东岳 多谢啊，这个我有所了解。可能我没有说清楚。仔细观察公式（14），系数矩阵A是个对角阵啊，直接除过去就好，x=b/A，感觉没有必要用各种复杂的矩阵求解器啊。

多谢东岳大神啊，回到公式（14），如你所说，我若是10万网格，那么就有十万个公式，组成10万阶矩阵，Ax=b。但是仔细观察公式（14），这10万阶系数矩阵A是个对角阵啊，直接除过去就好，x=b/A，感觉没有必要用各种复杂的矩阵求解器啊。请您指点一下。