三维多相搅拌机械Fluent仿真实例
本帖最后由 sookzhou 于 2009-12-15 14:11 编辑做了好几个月的搅拌机械,过程中走了一些弯路,知道初学者的艰辛,今天抽时间把简单的实例简单总结一下,提供给有用的人.
1 建模
a.搅拌桨的生成
b.搅拌釜的生成
c.网格划分
1)这里对静止区域和2个转动区域进行网格划分
2）网格质量的检查
d.边界条件的设置
1）设置实体区域，这里要设置3个实体区域：静止区域（zone sta），转动区域上1（zone 1 ），转动区域下2（zone 2）；
2）设置边界条件：这里要设置之前编号号的12个面为interfac，搅拌轴其他的区域(axis)为wall,自由液面（water face）为symmetry,通气管入口（gasinlet）为速度进口（在gamit中没有velocity inlet类型。在fluent可以改），他们的设置同前面；
3）点击标题栏的solver-----fluent 5/6即选择运算器为fluent;
e.保存与输出
1）点击标题栏的file-----export—mesh，输入文件名名确定即可，mesh文件的保存路径极为开始时指定的工作目录；
2）点击标题栏的file-----save,然后点击标题栏的file-----exit即退出gambit
2.fluent操作
A .网格操作
1) 启动fluent后选择3位求解器即可,依次点击file----read---case,找到刚才gmbit工作目录下的mesh文件；网格读入完成；可以选择display—grid，选择相应得区域，按display即可；
2）依次点击grid—check，检查最小网格体积是否小于0；
3）依次点击grid---scale，fluent默认的长度单位为m,改为gambit默认的单位mm,依次点change length units,scale,close ;
4）依次点击frid---swmooth/swap grid,分别点击smooth,swap循环,直至swap信息中出现number faces swaps: 0为止,此时关闭对话框即可;这步操作是对网格光顺以及对等角倾斜度(skewness) 高的地方交换网格以便于后面的运算;
B .模型设置
1) 点击define—model---solver,因为后面所用的MRF为稳态处理法,因此选取默认的3维稳态求解器即可;
2) 点击define—model---multiphas,选择 Mixture多相流模型,相数为2,其余默认即可
3) 点击define—model---vicious model选择标准k-e模型即可
4) 点击define—materials,选择液态水(H2O) ,fluent默认的工作介质为空气;
5) 点define—phase,确定phase1为主相,设置空气(air) 为第二相,设置气泡直径为0.003mm;点mass选项卡,设置传输数目为1,从phase2 to phase 1,mechanism设为7.5e-5;这里的数据,别的保持默认,确定即可;
6) 点define---operatin conditions 设置大气压强及重力;
7)点define---grid interface,选择转动动区域与静止区域相交面(即在gambit里定义的interface边界) ,设置6组即可;这个设置是为了使两个区域的数据能进行数据交换;
C.边界条件设置
点define—boundary types
1)对zone sta，phase选mixture,点set，默认设置，点ok即可;
2)对zone 2(下面转动区域)，phase选mixture，点set,在运动选项卡，运动类型选MRF,设置转速为4rad/s,其余默认；这里即设置了转动区域与静止区域的处理运用MRF的方法；同理对zone 1 做同样的设置
3)对wall:（上层搅拌桨）, phase选mixture,点set,点运动选项卡，设为运动的墙，运动类型为转动，设置它与zone 1相对静止（相对速度为0），别的默认；对wall（下层搅拌桨）也运用相同的处理方法；
4)点 axis, phase选mixture，点set，设为运动的墙，绝对转动速度为4rad/s,由于axis的相对区域为zone sta 是静止的，因此用绝对速度设置较好；
5）点gaslet，phase选mixture，将gaslet的边界条件类型改为velocity inlet, 定即可；phase选phase 1,set,
将速度方式改为组分（components）,各项绝对速度为0；phase选phase 2,set,
，将速度方式改为组分（components），Z向绝对速度为0.01m/s,体积分数为1；这表示出口处置有气体，且气体沿Z轴正向以0.01m/s的绝对速度喷出；
6）其余的边界条件默认即可；
D.运算
1）点solve—controls---solutions，压力速度耦合采用SIMPLEC算法，其余不变；
2）点solve—innitialize--- initialize,,选择从所有区域开始计算（compute from all zones），依次点init,apply,close
3）点solve—moniter—residual设置，
4）点solve—iretate,出现的对话框，设置迭代次数，进行迭代即可。
5) 当计算收敛后，计算会自动停止 很想学习学习，可以把模型和CAS文件发我邮箱吗？jksmnh@gmail.com 有收获，谢谢！我也想试一下你做的这个算例能帮我也发一份吗？dcfeng85@sina.com en,很详细，用的滑移网格吗，帮顶一下 三维的 CAS 和 DAT文件都很大
发不了 楼主，你好，能留个联系方式么？
或者加我QQ：31996626
谢谢，期待向你请教 向您请教 caiyuan0418@126.com QQ：33220707 向您请教 感谢 谢谢分享....1017496054@qq.com 谢谢楼主，我也是做搅拌器的，也走了好多弯路 我的qq是421332493，期待和你交流一下 可以学习一下，楼主最近又做欧拉欧拉的没有？液固两相应该怎么弄？ 非常感谢你的分享，可以留下你的邮箱或者qq吗 非常感谢分享 :handshake 说的太详细了，谢谢分享！！！ 楼主，我做列车制动盘叶片的三维流场分析，能加一下我的qq号码447063684吗，有问题想请教一下，谢谢！ 楼主，会模拟混合时间吗我也是弄搅拌桨 混合时间弄了好久没弄出来 还是发cas上来吧。 同道中人，期待与首楼主交流 大家都是做搅拌的，该建个搅拌群，方便交流！
页:
[1]
2