Question

fluent 模拟泵

在管道中设置一段类似于泵的功能的边界条件，有文献说用动量源项，不用进出口边界条件。自己编的动量源项只在源项附近有流体，整个计算域都没有流体。而且残差曲线的连续性方程精度到10-2次方就不下降了，请求高手帮助！！
如果能解决问题，可以追加悬赏分。
这是我自己编的动量源程序，流场中只有该动量源附近有流体流动，请问我编写的程序存在哪些问题，请高手解答！
#include "udf.h"
#define CON 20.0
DEFINE_SOURCE(mou_source,cell,thread,ds,equ)
{
real vel=0.45;
real x[ND_ND];
real y;
real z;
real source=0;
C_CENTROID(x,cell,thread);
y=x[1];
z=x[2];
if((((x[0]+0.0305621)*(x[0]+0.0305621)+(y-0.1733262)*(y-0.1733262))<=0.000049) && (z<=0.24 && z>=0))
source=CON*vel;
ds[equ]=0;
return source;
}

Answer 1

我做过循环回路里加风扇的模拟
可以用fan 边界去模拟， 中间含泵区域用Fan, 将泵的性能曲线拟合多项式系数定义进去。
具体设置可参照 HELP

另外，你的UDF源项正确加在泵区域的流体吗？其余流体部分设置有没有做好？这种方法完全可以实现的，我没细看你的函数，但我用FLUENT现成的例子加动量项没有你的这个问题。

追问

风扇模型是不是必须添加入口和出口呢？我的管路中间没有截面，所以用动量源代替泵了。还请您不吝赐教！

追答

可以不要入口出口。 中间必须有截面啊。
你要用动量源也要有截面把那部分流体区分出来

追问

我的动量源判断语句如下：if((((x[0]+0.0305621)*(x[0]+0.0305621)+(y-0.1733262)*(y-0.1733262))=0))
也就是在一片x，y方向上半径为7mm的圆，z方向为0—0.24m区域内添加动量源，方向向上。
请问这样对吗？多谢了！

追答

这个语句看起来没问题。 关键是动量方向怎么加呢？

为什么不网格中把这个区域划分出来？， 可以省去很多事，

追问

不会您说的网格划分方法，我加动量源方向是z正方向，因此source=CON*vel;
不知您说的区域划网格的方法在什么地方可以查询到，还请赐教。谢谢！

追答

有点无语。
这是网格的基本操作啊， GAMBIT中，把那块体分割出来，和其它区域交界设置interior，这个区域定义Fluid 单独命个名字，比如"Pump"。

这样导入到FLUENT里面就会找到那个PUMP区域了，然后UDF动量加载到这个区域。

我只能解释到这个程度了，见谅！

Answer 2

请问你的问题解决了吗？我也解决了这样的问题，想请教一下，qq932888020

Answer 3

模拟泵的流场，现在比较常用的都是设压力入口，压力出口，给定压力边界条件。连续性方程精度只到10-2次，可能是你的试验次数到设定值就停止了，这样的话增加迭代次数就行。如果不是上述原因，那你就得重新划分网格，得把网格划分的再密些。

追问

因为我必须在流场内添加示踪剂，不知压力入口与压力出口是否能模拟混合过程？
具体设置可否有参考资料？多谢！

追答

你可以参考两相流模型这一章，看下FLUENT流体工程仿真计算实例与应用（北京理工大学出版社）或者FLUENT帮助即可

追问

我翻阅了您说的这本书。但是还是没有解决泵的模拟，我是想将一个循环回路中间添加一段类似泵的动量源，我没有办法在管路中添加截面，也没办法添加出口和入口，请问您能否提供建议？