CFDl论文.docVIP

计算流体动力学（CFD）结课论文 论文题目：冷、热水混合器内部流场数值模拟分析混合器内部流场数值模拟分析通过使用FLUENT软件的标准k-ε湍流模型对冷热水混合器进行维数值模拟,分析其内部流场变化情况。通过模拟,能真实反映混合器内部的复杂流动,工程热水恒温混合器,是为适应中央热水工程向大型化、自动化和人性化发展的技术要求而研发的,是为太阳能热水工程和各种生活热水供水系统专门配套的一种全自动洗浴水恒温控制设备。广泛适用于宾馆、饭店、学校、医院、厂矿、机关及洗浴中心、游泳池等大中小型生活热水系—斯托克斯（Navier-Stokes）方程的求解模块。用压力校正法作为低速不可压流动的计算方法，包括SIMPLE、SIMPLEC、PISO三种算法，采用有限体积法离散方程，其计算精度和稳定性都要优于传统编程中使用的有限差分法。而对于可压流动采用耦合法，即将连续性方程、动量方程以及能量联立求解。FLUENT软件主要由前处理、求解器以及后处理3大模块组成。采用自行研发的GAMBIT前处理软件来建立几何形状及生成网格，然后由FLUENT进行求解。 2．控制方程和数值模拟 2.1 控制方程与标准k-湍流模型 本文主要利用FLUENT-2D的非耦合、隐式求解器，针对在冷、热水混合器内的定常流动进行求解。在求解过程中，还会利用FLUENT的网格优化功能对网格进行优化，使所得到的解更加可信。分析冷水和热水分别在混合器的上侧和右下侧流入，在容器混合后经过左下侧渐缩管道流入等径的出流管，然后流入大气。下面是本文中给出实例的相关参数以及冷、热水混合器的结构图： 问题描述：一个冷、热水混合器的内部二维流动和热量交换问题。温度为350K的热水自上部的热水小管嘴流入，与自下部右侧小管嘴流入的温度为290K的冷水在混合器内进行热量与动量的交换后，自下部左侧的小管嘴流出。混合器结构如右图所示: 不可压缩流体的连续方程在空间直角坐标系中的表达式为： （1） 对于不可压缩黏性流体的N-S方程，在空间直角坐标系中的表达式为： （2） 上述方程（2）再加上连续性方程（1）原则上就可以求得不可压缩黏性流体流场的解。但由于N-S方程中出现了速度的二阶导数，它的普遍解在数学上还有困难，只有某些特殊情况才能使方程得到充分简化，求出近似解。 由于此流场处于湍流状态，因此采用标准k-方程模型。标准k-方程模型的湍动能k和耗散率方程如下： （3） （4） 式中：为由于平均速度梯度引起的湍动能，为由于浮力引起的湍动能，为可压缩湍流脉动膨胀对中耗散率的影响。 湍流黏性系数： （5） 在FLUENT中，作为默认值常数，=1.44，=1.92，=0.09，湍动能k与耗散率的湍流普朗特数分别为=1.0，=1.3。 2.2数值模拟步骤： 2.2.1利用GAMBIT建立混合器计算模型 利用坐标网格建立节点 在两个节点之间创建直线； 利用圆心和端点创建一段圆弧； 由边创建面； 对各条边定义网格节点的分布； 在面上创建网格； 定义边界类型； 利用FLUENT5/6输出网格； 2.2.2利用FLUENT-2D求解器进行求解 读入网格文件； 确定长度单位为cm； 确定流体材料及其物理属性； 确定边界条件； 计算初始化并设置监视器； 使用非耦合、隐式求解器进行求解； 利用图形显示方法查看流场和温度场分布； 使用能量方程的二阶差分格式重新计算，改善温度场的计算； 利用温度梯度改进网格，进一步改善温度场的计算。 2.2.3数值模拟的过程 a、前处理--利用GAMBIT建立网格模型 启动GAMBIT，本例所选用的求解器为FLUENT5/6，在菜单栏Display Grid XY平面内对X、Y轴分别设置网格节点，在GAMBIT工作区将显示20*20的网格，这样网格节点就创建好了。紧接着确定不同类型的边界交点及圆弧中心点。隐藏坐标网格，由点建立直线和圆弧边以及入流和出流管嘴（由点到线），分别创建管嘴和混合器的面。所建图形如下： 然后创建边界线的内部节点和结构化网格，并且设置边界类型，输出网格并保存文件，网格图如下图所示： b、利用FLUENT进行混合器内的流动与换热的仿真计算 启动FLUENT--2D求解器后，将在GAMBIT中建立的混合器网格文件读入，FLUENT在读入网格文件的同时，会在反馈信息框内显示如下信息： 紧接着进行网格检查（Grid--Check），在FLUENT