CFD中文网

@李东岳 这个案例我用的商软，发现空隙率0.01和0.5差别不大，影响最大的是阻力系数，我试了这两个空隙率，水都能很顺畅的流出来

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

@bestucan 好吧，谢谢。我再琢磨琢磨

@搬运工 好的，谢谢回复:yes:

最近计算一个方腔流动遇到一些问题，这个流动沿流向（x方向）是周期性边界，yz方向是壁面，流动固定x方向的平均速度，自动求解压力梯度，在采用rans计算时是正常的，可以计算出一个合理的结果

pressure gradient = 0.179097

但是当采用PANS（自己编译的湍流模型）计算的时候，压力梯度会变得非常小，导致计算结果很差，

pressure gradient = 0.0903267

自己编译的模型可以在channel395中正常计算，但是在该算例中出现压力的降低，不知道是什么原因引起的

@李东岳 感谢东岳老师的回复！

@coolhhh 感谢您的回复 ，确实是这个样，没发现那个字母大写了。

按网格类型分类的软件列表：
Mesh Generation & Grid Generation on the Web

你可以读取openFOAM处理好的网格文件啊，它的网格文件还是很有规律的，至少我现在已经实现读取和处理openFOAM的网格数据了，自己在编写求解器。

这个是不是用到了基于物理量+动网格技术。。。

@李东岳 老师，最近我也遇到了这个问题，单相流计算的时候会提醒，但是结果也可以收敛，但是计算完初场的时候，改为多相流空化计算，马上就发散了，不知道怎么可以解决这个问题，已经重新画了好几次网格了

@东岳 谢谢老师

@Gordonaero 在 两相流含气率计算 中说：

两相混合密度=气相密度乘含气率+液相密度乘（1-含气率）

并且，气体密度也是定值

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

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

感谢分享！！！:huahua:

@队长别开枪

谢谢您的建议，感觉这个方案可行。 我尽早试试，给您答复。。。

谢谢

谢谢版主。

IMG_20231127_193223.jpg
在计算单相流作为初场，可以收敛，但是湍流粘度比超限，继续设置为多相流空化计算的时候，开始发散，不知道有哪些解决办法，已经重新画过网格。
边界条件是速度入口，压力出口，湍流sst k-w，稳态计算