CFD中文网

我觉得STARCCM+的网格功能挺6的，笛卡尔/多面体等，控制的参数有限，基本上零基础跟着教程做几遍就差不多了，网格生成速度也很快，比较复杂的几何还是很容易解决的；
其他的用point wise多一点，复杂的几何自己划分块有点麻烦，如果不能分层很容易混乱，简单几何还是非常便利的；ICEM复杂几何切起来还是比较舒心的；^_^

@东岳 我估计是抖音看上了罗老师的流量:yes:

@Albert-Lee 我跟你遇到了相同的问题，你的测试和我的结果也相差不大。目前我也在寻找一个解决方法，我看到的在UDF中获取Level set function 是使用 C_LSF(c,t)宏，其中c = cell index, t = mixture phase thread，原链接：https://www.cfd-online.com/Forums/fluent/144991-access-level-set-function-udf.html

注意：需要额外include "sg_ls.h" 因为他没有包含在udf.h里面

之前也困惑问过类似的问题，李老师说复杂的情况用ICEM 之后网格转换一下，这样确实是省时又省力哦。

https://www.cfd-online.com/Forums/openfoam-solving/230250-buoyantpimplefoam-boussinesq-incompressible.html

好像就是不可压的:141:

vim保存窗口：

:mksession ~/.session.vim vim -S ~/.session.vim

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

@Gordonaero 感谢，这两篇论文之前没有看到过。
这个动态填充的研究是想预测填充以及降温时间吗？关于这方面管道预冷的研究，之前看过国内的西交大有做这方面的研究，包括实验与数值计算。
按照我的理解，这个均相模型将混合物作为一种赝流体，封闭的关键其实是混合物的状态方程以及如何将混合物拆解成气液两相的比例。

如题，这段时间在做一个分子模拟的案例，需要应用到N-S方程反求粘度，想知道这里边的一些基本思路。

谢谢，已更新

大神如何让stl 导入 运行第一步blockMesh？
我一运行就报错

哈哈哈哈有点搞笑

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

OpenFOAM-3？？为啥装这么老的版本。另外，CGAL不安装也可以，不影响openfoam安装

@random_ran 你好 获取了速度场随时间变化的数据后，是用啥软件计算可以求出雷诺应力呀？

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

这本确实是神书，里面的示例代码给的真的相当赞。

两相流传质有专门的宏，不需要申请UDM，找到Euler两相流的传质宏，挂载到mass transfer那个选项里就行

@星星星星晴 真是厉害的，厉害的啊！又多一技能，哈哈哈。感谢分享！

@李东岳 妈耶，羞煞我也