本案例利用Eulerian多相流模型及组分输运模型模拟臭氧在流化床中的分解过程。 ![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_1_2023102304170772.png)
- 利用UDF定义修正的Syamlal-O’Brien阻力模型
1 问题描述该问题涉及流化床臭氧分解的瞬态启动。液相是臭氧和空气的混合物,而固相由直径为87.75微米的沙粒组成。流化床原理图如下图所示,该区域被建模为二维平面圆柱壳。![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_2_20231023041707822_wm.png)
2 Fluent设置- 以2D、Double Precision方式启动Fluent
- 读取网格文件2-D-FBed_Ozone_Cartesian.msh.gz
2.1 General设置![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_3_20231023041708541_wm.png)
2.2 Models设置![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_4_20231023041709838_wm.png)
- 激活Species Transport模型,启用Volumetric体积反应
![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_5_20231023041710932_wm.png)
2.3 Materials设置![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_6_20231023041711800_wm.png)
- 修改混合材料mixture-template的属性
![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_7_20231023041713400_wm.png)
![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_8_20231023041714697_wm.png)
![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_9_20231023041717322_wm.png)
![](http://image109.360doc.com/DownloadImg/2023/10/2316/274254993_10_20231023041718760_wm.png)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
注:用于流化床颗粒的固体颗粒,需要用fluid类型的材料 2.4 定义组分![](http://pubimage.360doc.com/wz/default.gif)
- 如下图所示设置phase-1为mixture-template
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
2.5 编译并加载UDF![](http://pubimage.360doc.com/wz/default.gif)
- 添加UDF文件rrate.c及bp_drag.c,编译并加载
![](http://pubimage.360doc.com/wz/default.gif)
- 点击按钮Function Hooks…打开设置对话框
![](http://pubimage.360doc.com/wz/default.gif)
- 如下图所示设置volume Reaction Rate为rrate::libudf
![](http://pubimage.360doc.com/wz/default.gif)
2.6 定义相![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
- 进入Collisions标签页,设置相间碰撞恢复系数为0.8
![](http://pubimage.360doc.com/wz/default.gif)
2.7 操作条件![](http://pubimage.360doc.com/wz/default.gif)
2.8 边界条件![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
- 打开边界leftwall的phase-2相设置对话框
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
2.9 Region- 右键选择Cell Registers,选择菜单New → Region…打开区域新建对话框
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
2.10 Controls- 右键选择节点Controls,点击弹出菜单项Equations…打开设置对话框
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
2.11 初始化计算![](http://pubimage.360doc.com/wz/default.gif)
- 点击Patch…按钮打开对话框,如下图所示Patch区域的Phase-2体积分数为0.52
![](http://pubimage.360doc.com/wz/default.gif)
- Patch完毕后查看phase-2体积分数分布,如下图所示
![](http://pubimage.360doc.com/wz/default.gif)
2.11 设置Autosave![](http://pubimage.360doc.com/wz/default.gif)
2.12 求解计算- 设置时间步数为10000,时间步长为0.001,设置max iteration/Time step为40
![](http://pubimage.360doc.com/wz/default.gif)
3 计算结果![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
![](http://pubimage.360doc.com/wz/default.gif)
|