CFD中文网

@李东岳 东岳老师，个人原因，麻烦删掉这个帖子。

@lixin19981013 可以联系我，我教你。qq，1259451345

@dr-yiqi-wang 求解器计算过程中肯定求解了dc/dx，但是没有储存下来，自己写个udf把dc/dx储存下来就可以了

@dr-yiqi-wang 可以将dc/dx储存到user define memory

@木子可 这篇文献我看过。作者的意思是，如果你要模拟颗粒之间孔隙内的流动，就把颗粒的几何画出来，抽取流道。这样必然会遇到一个问题，就是如何处理颗粒之间的接触点，以避免产生畸形网格（因为接触点附近的缝隙是非常狭窄的）。

如果局部网格尺寸为1e-5m时，由于vof中库朗数的限制，时间步长也需要设置的很小，为1e-6或者1e-7，此时计算时间过长，该如何解决呢？

@sungda 可能是。加密试试。

一楼图放错了。网格和一楼一致，结果如下图所示。

@sungda 我也遇到一个VOF的问题，不知道您有没有经验。
https://cfd-china.com/topic/5207/fluent-vof模拟液滴悬挂-液滴中心出现气泡?_=1634704563012

使用Fluent中的VOF方法模拟液滴悬挂。模型是轴对称的，红色为液体，蓝色为气体。在运行一段时间后会发现，对称轴中心会产生气泡，并将一部分液体挤出液滴。

目前的猜想是数值误差，请问各位老师遇到过吗，该如何避免。

网格如下图，

@cccrrryyy 和温度是不一样的。两个不同温度的物体碰到一起，温度在边界上不会间断，应该是连续的。但溶解，气液界面的浓度是间断的。比如下图，

@李东岳 感觉这个模型和我的想法不一样。感觉它是一个模型同时适用于气体域、液体域和气液边界，并没有分开处理，气液界面也不需要边界条件。

@李东岳 是静止的计算域，不用VOF。用多计算域求解。气体溶解到水中，如果压力梯度大，可以用简易的微分方程计算速度，但暂时不考虑压力梯度产生的速度。

至于PDE，就是存粹的扩散方程，$\frac{\partial c}{\partial t}=D\frac{\partial^2 c}{\partial x^2}$。

主要想问，如果按我的想法，那么气相的气液界面上边界条件应该怎么设置？

@wwzhao @bestucan @cccrrryyy
抱歉打扰几位老师，您看是否有经验，求指导一下。

@tens 感谢。这些文献应该很有用！

@tens 如果去结合 VOF，我感觉太复杂了。我可以完全不考虑液体流动、相变等等。**只考虑组分传递的事情，用对流扩散方程就好。**有没有这样简单的模型。

@李东岳

模拟气体溶解，我认为可能主要看气液界面。如果气液界面非常复杂，肯定就没法直接模拟了，这时候就用到模型。比如说，鼓泡，这个气液界面非常多，可以把气体溶解的过程看成一个源项，利用传质系数$h_m$、气液比表面积$a_gl$、Henry系数$H^{cp}$、气体压力$p_g$和水中的气体浓度$c_l$计算气体溶解的量。

但是，如果气液界面单一，是不是可以直接模拟呢？（我不确定）
比如下图，气液界面就是一条直线，如果要模拟溶解，把模型分成两个域，气体域（物理量$c_g$）、液体域（$c_l$）。液体域在气液界面上的液相浓度$c_l$可以利用Henry定律计算，并且和气体域的局部气体浓度相关。那么气体域在气液界面上的通量怎么计算呢？

如图，是 Fluent 模型的收敛残差，连续性和三个动量方程的残差都是 $1\times10^{-3}$，只有能量方程是 $1\times10^{-6}$，这有什么特别的考量吗？

按我的经验来看，自定义源项后发散，一般都是源项太大了。你可以这么测试，

1. 设定源项大小为较小的常数，比如 1e-6，如果不发散，基本就确定是源项太大了；
2. 用 UDM 把源项储存下来，看看大小。