雾化模拟算例.doc

静电雾化模拟算例 问题描述本文利用FLUENT的DPM模型对带电液体的雾化情况进行研究。计算区域是一个直径100mm，高70mm的圆柱，简化为二维模型为100mm×70mm的平面。喷头支撑结构分为上下两段，毛细孔径为0.5mm，带电液体从毛细管喷出。 燃烧器启用DPM模型，先计算连续相，在利用UDF计算离散相 一 利用Gambit建立雾化模型 第一步：启动gambit并选定求解器（FLUENT/UNS) 第二步：创建雾化模型 操作：Operation→Tools→Coordinate System 在弹出的Display Grid 对话框中，输入X,Y的值，分别是100，70，点击Apply 图1 雾化区域计算图 第三步：建立喷嘴 喷嘴支撑结构分为上下两部分，上段尺寸为5mm×5mm，下段为3mm×3mm,喷头直径为0.5mm，长10mm。按照点、线、面的顺序逐步生成，如图2所示。 图2 喷嘴及支撑结构 第三步：划分网格 网格划分采用TGrid类型，喷头附近网格划分密集Intervai size为0.3，四周稀疏Intervai size为1，这样可以减少计算量。划分后的网格如图3所示。图3 网格划分图 第四步：设置边界类型 操作：ZONES →SPECIFY BOUNDARY TYPES 打开边界类型设置对话框如图4所示． 边界名称 边界类型 液体进口inlet2 VELOCITY-INLET 支撑结构及喷头Wall Wall 接收板Wall Wall 空气入口inlet1 VELOCITY-INLET 空气出口outlet PRESSURE-OUT 图4边界条件对话框 图5 边界条件设置 第五步：输出2D网格 操作：File→Export→Mesh．．．．．．输出3D网格，完成Gambit前处理． 二 利用FLUENT-2D求解器进行模拟计算 第一步：与网格相关的操作 　　1．读入网格文件 　操作：File→Read→Case．．． 　　在读网格文件后，将在FLUENT的console窗口中，报告网格和其他一些相关文件信息． 　　2．检查网格 　　操作：Grid→Check 　　网格检查列出网格的最小和最大的x与y值，并报告其他许多关于被检查网格的特征或错误，比如，网格体积必须不为负。 　　3．网格比例设置 　　FLUENT的缺省单位是m．若网格是以cm单位建立的，在Scale Grid面板中应选用相应的比例关系　　操作：Grid→Scale．．． （１）在Units Conversion中的下拉列表中选cm表示网格以厘米生成　　（２）点击Scale． 　　4．显示网格 　　操作：Display→Grid．．．　 图6 雾化模型的网格显示图 第二步：设置求解模型 　　1．定义计算域为2D，且保持缺省的求解器 　　操作：Define→Models→Solver．．．打开“Solver”对话框如图7所示． 图7求解器对话框 2．空气相选用层流模型 　　操作：Define→Models→Viscous．．． 　　打开“Viscous　Mode”对话框如图8所示． 图8 计算模型对话框 第三步：流体材料设置 操作：Define→Materials．．． 选择理想气体参数，点击Close关闭此面板。 　图9 材料对话框 第四步：边界条件设置 1．打开边界条件面板 操作：Define→Boundary Conditions．．． 打开“Boundary Conditions”对话框如图10所示． 图10 　边界条件对话框 设定空气进口inlet1的边界条件如图11所示，赋予如“Velocity Inlet”面板所示的进口边界条件． 图11　空气速度进口边界条件 4．设定空气出口pressureoutlet边界条件见图12所示． 图12 压力出口边界条件 5.设定支撑结构和喷嘴外壁的边界条件接收板Wall2的设定同上。 图13 壁面边界条件 第五步：初始化并求解 1．设定初场 操作：Sovle→Initialize→Initialize．．．见图14所示． 图14　求解初始化 （１）在Computer From下拉列表中选择inlet1． （2）点击Init设定变量初值，然后关闭面板． 2．设定松弛因子 操作：Solve→Controls→Solution．．．保持默认值。 图15　求解控制对话框 3．在计算期间打开残差的图形监视图 　 操作：Solve→Monitors→Residua