本文摘要(由AI生成):
本文介绍了使用CFX软件对某模型进行流体动力学计算的过程。首先,通过设置入口边界条件、出口边界条件和对称边界等,完成了模型的边界定义。然后,在求解控制中设置了相关参数,并启动计算。计算完成后,使用CFD-Post对结果进行了后处理,得到了压力降的计算结果,并与解析值进行了比较。文章最后提供了下载链接和版权声明。通过本案例的计算,验证了CFX软件在流体动力学计算方面的应用能力和精度。
本案例利用CFX计算管道中湍流流动压降。
本案例计算空气流经光滑的水平管道,管道长度2m,半径0.002 m,空气密度1.225 kg/m3,粘度1.7894e-5 kg/m-s,管道入口速度50 m/s,出口压力 0 Pa。计算管道的压降。
采用轴对称模型。入口采用速度边界,50m/s;出口为压力出口,静压 0Pa;采用稳态求解。湍流模型采用SST k-omega模型。
流经管道的雷诺数:
为湍流流动,采用布拉休斯公式计算摩擦系数(F. M. White. Fluid Mechanics(7nd ed). P366页)
则压力降:
3.1 启动CFX
启动CFX后,软件界面如下图所示,设置Working Directory为合适的英文路径。
点击按钮CFX-Pre启动CFX前处理界面
点击菜单File → New Case新建case
弹出对话框中选择General,点击OK按钮
3.2 读取网格
选择菜单File → Import → Mesh,弹出对话框中选择网格文件VM03.cas
注:对于2D轴对称模型,通常可以在Fluent中设置好轴对称后输出为cas文件再导入CFX,这样CFX会自动构建楔形体。否则若直接导入msh文件,则导入的模型为长方体。
导入的模型如图所示。
3.3 设置材料参数
鼠标双击模型树节点Materials > Air at 25 C打开材料属性编辑面板
如下图所示,材料设置面板中修改Density为1.225 kg/m3
设置Dynamic Viscosity为1.789E-5 Pa.s
点击OK按钮关闭面板
2.4 设置计算域
鼠标双击模型树节点Flow Analysis 1 > Default Domain打开区域设置
切换到Fluid Models标签页,设置Heat Transfer为IsoThermal,设置Turbulence为Shear Stress Transport,点击OK按钮关闭面板
注:本案例不考虑传热
切换至标签Initialization,如下图所示进行全局初始化,点击OK按钮关闭对话框
注:本案例为稳态计算,也可以采用自动初始化
2.5 边界条件
1、入口边界
右键选择模型树节点Default Domain,点击弹出菜单项Insert → Boundary插入边界
将边界命名为inlet
在Basic Settings标签页中设置Location为inlet
切换至Boundary Details标签页,设置Mass And Momentum指定方式为Cart. Vel. Components,设置U为50 m/s,其他两方向速度为0
设置Turbulence指定方式为Intensity and length Scale,设置Fractional Intensity为0.01,指定Eddy Length Scale为0.002 m
其他参数保持默认设置,点击OK按钮关闭设置面板
2、出口边界
右键选择模型树节点Default Domain,点击弹出菜单项Insert → Boundary插入边界outlet
在Basic Settings标签页中设置Location为outlet
切换至Boundary Details标签页,设置设置Mass And Momentum指定方式为Static Pressure,设置Relative Pressure为0 Pa
点击OK按钮关闭对话框
3、设置对称边界
右键选择模型树节点Default Domain,点击弹出菜单项Insert → Boundary插入边界symmetry
在Basic Settings标签页中设置Location为periodic.1,如下图所示
其他参数保持默认设置,点击OK按钮关闭设置面板
注:CFX读入Fluent的二维cas文件后会多出一些边界面,必须将其设置为对称面。剩下未定义的面默认为wall边界,本案例采用默认设置。
2.6 Solution Controls
鼠标双击模型树节点Solver Control,打开设置面板
如下图所示的设置面板中,设置Turbulence Numerics Option为High Resolution,设置Max. Iterations为400,设置Length Scale Option为Aggressive,设置Residual Target为1e-7
点击OK按钮关闭面板
点击菜单项File → Save Case保存文件
点击工具栏按钮Define Run写入求解器文件
点击Save按钮启动CFX求解管理器
2.7 开始计算
求解管理器对话框中设置Run Mode为Intel MPI Local Parallel,点击下侧的加号可增加CPU数量,如下图所示设置CPU数量为6个
点击按钮-Start Run开始计算
计算完毕后弹出提示对话框,激活选项Post-Process Results及Shut down CFX-Solver Manager,点击OK按钮启动CFD-Post
CFD-Post中左侧模型面板中切换至Expressions标签页,在任意空白位置点击鼠标右键,选择弹出菜单项New新建表达式
弹出对话框中为表达式命名为PressureDrop,点击OK按钮关闭对话框
输入表达式areaAve(Total Pressure)@inlet -areaAve(Total Pressure)@outlet,点击按钮Apply,则Value显示为计算结果。
如图所示,计算得到的压力降为20618.3 Pa,与解析值22366 Pa相比,误差7.81%,本案例网格较为粗糙,增加网格密度有助于提高计算精度。
下载链接:见附件
声明:原创文章,本文已授权,首发CFD之道公 众号,部分图片源自网络,如有不当请联系我们,欢迎分享,禁止私自转载,转载请联系我们。