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文件，在此处设置空气和水的物性参数，内容如下： 

turbulenceProperties文件和g文件无需修改

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

p_rgh文件当中的内容如下： 

U文件当中的内容如下： 

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

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

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

接着我们设置controlDict文件 

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

    输入mpirun -np 4 interFoam -parallel开始计算 

等到计算结束

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