CFD中文网

缩小模型只能保证一个无量纲数相同，一般的流动考虑雷诺数就行，缩放会带来一些非真实的物理效应，但是计算量不够就不要讲究了

$\epsilon \rightarrow 0$ 并不是$\epsilon = 0$，否则不会有...if $u$ is smooth at $x_1$ and $x_2$...。在粘度趋向于很小的时候，不连续变成具备一定厚度的光滑解，同样承认有厚度的激波。所以
$$
\epsilon\int_\Omega\frac{\p}{\p u}\left(\frac{\p \eta}{\p u}\right)\left(\frac{\p u}{\p x}\right)^2\rd x\rd t \geq 0
$$
另外，
\begin{equation}
\int_{x_2}^{x_1}\left(\epsilon(\eta_q q_x)_x\rd x -\epsilon\eta(\eta_q)_q q_x^2\right)\rd x=\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)-\epsilon\eta(\eta_q)_q q_x^2\Delta x
\end{equation}
考虑一个非常小的$\epsilon=1e-10$，在控制体内$\epsilon\left(\eta_q q_x|_{x=x_1}-\eta_q q_x|_{x=x_2}\right)\rightarrow 0$，$\epsilon\eta(\eta_q)_q q_x^2\Delta x$还是大于0.

缺这个包：binutils

或者路径没配置对，可以用echo $PATH看看有没有/usr/bin这一项

@tens 好的，谢谢你，我试一下

@东岳 总压减小是因为速度变小，还是静压变小，或二者都有？无粘不可压的时候就是动能与压力势能的转化，是守恒的。

可以在 OpenFOAM cavity 算例里执行 blockMesh 后，查看 constant/plymMesh 文件夹。就是这些东西。
不过记录形式向会有区别。但都是典型的“非结构网格”。

网格系统，自己开发的通用性不强，处理起来全都得自己手写一套。当然好处就是不用去学成熟的网格系统。便于上手，但越往后，想要的功能越多，越麻烦。

如果是 OpenFOAM 的用的网格系统，paraview 就可以看，也可以用 blender 结合 blockMesh 生成。snapHexMesh 也可以。对网格、流场的各种处理也可以调用 VTK 库直接在程序里进行。

看了一下你上传的数据，没有像样的文件头，估计是刚接手一套祖传代码:136:

已经解决了！谢谢各位

@东岳 感谢回复。建议原文适当增加配图，这将提高可阅读性。

牛顿内摩擦定律与stokes定律间有什么关系？二者计算壁面应力有何区别？

:136: :136: :136: :136: :136: :136:

只能发图了..

网格数量太多了。15510 * 16000 * 22个网格？

@bestucan 在别人电脑上是正常的 ，很奇怪，我最后重装电脑解决了

@ycx1042341217 @huningdong 请问两位，壁面第一层网格高度可以做到什么程度？我发现最小只能设置到大约5e-5m:136:

$sdsfa$

$sdsfa$

sdsfa

Info << patchCentre.size() << endl; // 若注释掉上一条Info语句，在这里输出patchCentre的size得到的值变为了0

换成这个呢？

Info << patch.faceCentres();

更新一下
发现了问题的源头，是网格的问题。
我使用ICEM画网格，或者OpenFOAM的blockMesh工具画网格，都是先画背景网格，然后对近壁面进行多次加密。加密是八叉树类型的加密，所以在网格尺寸变化的地方，大网格的尺寸直接是小网格的二倍。这样的网格确实不好，可能一些拓扑结构上的不好导致了singleGraph程序进行插值时找不到插值的基点。
解决办法：
举例说距离壁面0~1的网格尺寸是1,1~2的网格尺寸是2，那么就使用singleGraph在0~1这个范围内提取一条线，在1~2再提取一条线，这时的type可以是uniform，也可以是midPoint或者midPointAndFace。一句话，就是分段提取信息。
我是重新画了网格， 老师说之前画的这网格根本不能用:mihu:

@东岳 谢谢李老师的回复！
我看了下跟公式（7）到（8）的过程很相似，但是为什么公式7到公式8的过程中将公式7左边第四项给略去呢？

改重心 或者自定义坐标系

$\alpha$

第一步，先算稳态的， 用simple算算看看能不能产生非对称尾迹。然后用这个作为初值计算。
第二步，换成其他的湍流模型试试， SA模型我用的不多，fluent的SA模型的话，里面的参数设置可能跟标准模型有差别， 可以参考NACA turbulence resources 里面关于SA模型的说明。
还有一点，默认的收敛依据10-3可能不是特别好，可以试试把收敛依据降到1-10以下， 固定迭代步数，例如可以设置成10， 即每个时间步迭代10次收敛。