OpenFOAM——气泡上升

计算域的顶部为大气,其余部分为壁面

 

流体的物性参数为:

 

首先进行建模操作,任何建模软件均可,本算例采用ICEM直接建模,生成网格,然后利用OpenFOAM下转化网格,划分完成的网格如下: 

 

网格比较密集,高度方向设置了960个节点,宽度方向上设置了480个节点,在导出网格的同时,我们对网格进行了缩放。建议在导出网格的时候就缩放网格,避免在后面的操作中,由于忘记缩放网格,导致计算出错

接下来转入OpenFOAM的操作:

首先新建一个文件夹,名字任取,用来作为算例文件夹,本算例中我将该文件夹命名为:bubble 

 

然后进入OpenFOAM的安装目录,将安装目录下的damBreak算例(我的目录为/opt/openfoam5/tutorials/multiphase/interFoam/laminar/damBreak/damBreak)下的0文件夹、constant文件夹和system文件夹拷贝到bubble文件夹下,然后删除system目录下的blockMeshDict文件,因为我们利用OpenFOAM的命令转化.msh文件为OpenFOAM能接受的网格文件。

然后我们将刚才我们生成的.msh网格拷贝到bubble文件夹下。在算例文件夹下打开终端,输入fluentMeshToFoam命令:

 

我们打开constant文件夹下的transportproperties文件,在此处设置空气和水的物性参数,内容如下: 

OpenFOAM——气泡上升

turbulenceProperties文件和g文件无需修改

然后对初始边界条件进行设置,下面转入0文件夹下进行操作:

p_rgh文件当中的内容如下: 

OpenFOAM——气泡上升

U文件当中的内容如下: 

OpenFOAM——气泡上升

alpha.water文件当中的内容如下:

OpenFOAM——气泡上升

相关边界的具体说明,可以参考《OpenFOAM用户指南》的溃坝算例。以上设置完成以后,我们可以用下面的python代码设置气泡:

OpenFOAM——气泡上升
    以上python代码的详细解释,请见《OpenFOAM——设置自定义非均匀场区域》,链接:http://blog.sina.com.cn/s/blog_14d64daa10102xhx8.html 

OpenFOAM——气泡上升

接着我们设置controlDict文件 

OpenFOAM——气泡上升

为了加快计算速度,我们进行并行计算,回到算例文件夹下,打开终端,输入decomposePar进行分块

OpenFOAM——气泡上升
OpenFOAM——气泡上升    输入mpirun -np 4 interFoam -parallel开始计算 

OpenFOAM——气泡上升

等到计算结束

将计算结果导入paraview里面进行处理 

OpenFOAM——气泡上升

posted @ 2017-10-10 23:27  硫酸亚铜  阅读(1534)  评论(1编辑  收藏  举报