CFD中文网

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

最近在看k-e湍流模型壁面函数的知识，公式推导过程中假定了壁面附近脉动动能与产生相平衡，
1.JPG
有没有人知道这个式子怎么来的？湍动能耗散与产生相平衡的条件是什么？

$\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.

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

顶一下，如果有同学想组队团报，可以联系QQ:876440251
我问了客服小姐姐，说优惠可以享受hhh

各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？

同问,我也正在做类似的东西,但是算到最后连续性方程都是收敛的,但湍流方程总不收敛

@李东岳 感谢李老师，我试试看能不能画出来

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

提示信息如下：
...
rt -lm -o /myelasticThermalSolidFoam
/usr/bin/ld: cannot open output file /myelasticThermalSolidFoam: 权限不够
collect2: error: ld returned 1 exit status
make: *** [/home/zhangec/foam/foam-extend-5.0/wmake/Makefile:160：/myelasticThermalSolidFoam] 错误 1
请教各位老师同学怎么解决？

他们测试了10步、60步啥的，还有30步好像，上周看的，现在忘了

以及对于入口气态氧化剂超声速射入燃烧室，入口液态燃料低速进入燃烧室进行非预混燃烧，最后速度为低马赫数的用sprayFoam进行燃烧仿真合适吗
请各位给小白点建议

@装深沉同学 请问，最后您解决了吗？

多谢老铁，已更正，将在下一版本更新

这是用离散法求解PBE的增长过程。如果没有增长过程，如果不是离散法，就不需要这么做。a是波速，Minmod是limiter。如果Minmod目前你还不熟悉，我大体知道你的CFD基础在哪，你需要补CFD基本理论，再来看PBE:138:

@guanghuilineu 那看起来发散是个很正常的结果，如果相变区速度不为0时，也变速也变密度。

如果只变速，和周围空间交换能量，即动能转换为其他形式的能量，没问题。
如果速度和密度一起变，动能是依赖于质量的。这时候，质量也变了，不守恒了，就崩了。

消失的那部分质量，带有的动能如何传递出去。

压力就是传递能量的桥梁。在传递的时候一头突变了，就断了。可不可以先使相变区域停止流动再变密度。如果这样，不会出现动能突变的状况。或者在变密度操作后立刻更新一下速度；变前计算动能；变后用动能，根据变后的质量，计算新的速度。

或者最粗鲁的方法，用极小的时间步，这样突变也可以稳住。

多谢老师解答，问题已经解决了

多面体明显质量好过关 过基础计算关 相对正六面体 同样参数控制 数量变小了

@李东岳 对哦，我审题不仔细，周期边界的话，哪来的入流。

@东岳 Thin-Layer Navier—Stokes Approximation