| 查看: 1439 | 回复: 2 | ||
[求助]
圆柱形流化床fluent模拟问题
|
|
我的床体是圆柱筒内径100mm，外径108mm，高度1000mm，为了减少计算，只取400mm的高度作为计算区域的高度；床内球型颗粒直径0.52mm，颗粒密度1712kg/m3，床层初始填充率0.584，初始床层高度100mm，床层颗粒质量0.785kg；流化介质为标况下的空气，进口表观气体流速0.2m/s，均匀布风。我用ICEM CFD建立了三维几何模型和划分了O（Ogrid）型网格（图1），网格最小尺度1mm、最大尺度2mm（边界层网格依据k-ω（SST）湍流模型的y+值取0.5，第一层高度1mm），然后生成非结构网格导入fluent中进行4核CPU并行计算，并进行网格质量检查，如图2。 在fluent中采用Euler多相流DDPM（DEM）模型，同时采用k-ω（SST）湍流模型。颗粒的初始化采用surface injection方式，采用老师视频中的区域表面发射（底部的fluid-1区域），从0s开始，1e-9s结束，total flow rate为7.85e+8 kg/s，particle time step size为1e-7s。 边界条件：底部采用速度进口0.2m/s，出口采用压力出口，侧面是圆柱壁面条件。 求解：method里都选用了默认的一阶迎风格式，选了SIMPLEC求解器，并勾选了下角的N-phase volume equation；control里把动量的松弛因子改成了0.2，viscosity的松弛因子改为了0.8，其他采用默认。 初始化：把k值改成0.01，其他采取默认。 计算：选择了固定的time step size：1e-6,时间步为1，作为调试。 在调试过程中，发现连续方程的残差始终达不到1e-3以下（即使把每个时间步长内的迭代步增至500，图8-10），但是在report里flux的进口与出口的质量流率差值很小（图7），符合大家说的做为判断收敛的标准。不过，在report里，体积积分的mass中，fluid-1区域颗粒相显示的质量只有0.490637kg（图5），其mixture的质量加权平均mass-average（图3）为1025.17kg/m3，体积加权平均为1001.34kg/m3（图6）。 然后，我进行了2个时间补偿的调试计算（图11-12），连续方程的残差突然增大到1e+4量级，而且降幅也很小，并出现了particle parcels的incomplete不为零，计算的速度越来越慢。 我不知道为何会出现这种现象和错误，问题出在哪里，该如何解决？ 此外，我想问一下injection这种颗粒加入方式，采用surface里的一个区域在短时间内大量发射从而完成初始密集床层的堆积初始化，这样操作与实际完全可以符合吗？ 请大神们帮忙分析一下。谢谢！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 圆柱形流化床fluent模拟问题图片.docx
2019-10-22 13:01:02, 2.54 M
» 猜你喜欢
中国科学院深圳先进技术研究院-光纤传感课题组招生-中国科学院大学、深圳理工大学联培
已经有4人回复
-
调剂推荐
已经有4人回复
-
总分322求生物学/生化与分子/生物信息学相关调剂
已经有5人回复
-
0703化学338求调剂！
已经有3人回复
-
284求调剂
已经有7人回复
-
271求调剂
已经有6人回复
-
312求调剂
已经有9人回复
-
352求调剂
已经有4人回复
-
085602化学工程求调剂。
已经有3人回复
-
085600材料与化工306
已经有6人回复
蓝云思雨
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7370.6
- 散金: 75
- 红花: 1
- 帖子: 355
- 在线: 35.3小时
- 虫号: 2001929
- 注册: 2012-09-14
- 性别: MM
- 专业: 能源化工
2楼2021-03-06 21:40:10
蓝云思雨
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 7370.6
- 散金: 75
- 红花: 1
- 帖子: 355
- 在线: 35.3小时
- 虫号: 2001929
- 注册: 2012-09-14
- 性别: MM
- 专业: 能源化工
3楼2021-03-06 21:51:48
