fluent汽车模型cfd模拟..docVIP

基于FLUENT的汽车模型CFD模拟 本次作业是利用CFD对流经结构外表面流动进行研究，模拟汽车运动中流经汽车表面的气体动力学现象。 目的： 创建对单一小车模型的CFD模拟； 进行两车模型CFD模拟并研究车间距对模拟的影响； 了解湍流流动模型以及对y+； 创建对单一小车加尾翼模型CFD模拟。 操作过程 单车模型 （1）创建独立性网格 小车形状以及计算域如图1，小车长度为1。为了不影响计算，计算域设定的很大，长27，宽为5。 图1 小车形状及计算域 由于小车以100km/h的速度运动，相当于小车不动，地面以27.778m/s的速度向x方向运动，此时空气流动为湍流，因此在小车的边界上要设置边界层网格。设置第一层为0.0029，成长因子为1.09，层数为6，如图2。 图2 边界层网格 因为主要研究小车附近的空气流动特性，从小车到计算域边界应用sizing function。小车地面到地面的sizing function参数设为start size：0.035，growth rate：1.15，max size：0.1；小车其他边到计算域边界参数设为start size：0.035，growth rate：1.2，max size：20，如图3。 图3 sizing function设定 创建面网格，非结构网格，选择三角形网格，为了创建独立性网格，将网格密度逐次加倍得到种网格如图4。 图4 网格 对网格进行检查，发现密度最大的第四种网格的aspect ratio大于6的有13个，edge ratio大于6的有8个如图5，其他网格质量较好。 图5 网格检查 定义边界： Sky: Symmetry Inlet: Velocity inlet Outlet: Pressure outlet Ground: Wall Car_front: Wall Car_top: Wall Car_back: Wall Car_bottom: Wall （2）应用fluent对小车网格进行计算 分别对四种不同密度的网格进行计算。读入网格，检查完毕后需要对网格进行光滑处理，如图6。 图6 smooth/swap 因为空气速度比较大，是湍流流动，所以要选用湍流模型，这里选用k-e模型，各参数选用默认值，如图7。 图7 湍流模型 设定边界条件，inlet为速度入口，速度为27.778，outlet为压力出口，ground为移动壁面相对速度为27.778，其他参数如图8。 图8 定义边界条件 从入口初始化，并计算，监视残差与阻力系数，如图9。 图9 残差与阻力系数 Continuity都是随着迭代次数增加一直减小，其他残差随迭代次数增加先增大后减小。阻力系数是从计算开始立刻减小，然后增加，再迭代10次左右以后保持不变。 对近壁面Y+值作图。四种网格小车边界层的y+值如图。 图10 y+值 四种网格小车近壁面处的y+值大部分都在75-175之间，符合要求。 （3）分析结果 分别计算四种网格小车的阻力系数和升力系数如图。 图11 阻力系数 图12 升力系数 第一种到第四种网格阻力系数分别为0.4898、0.4666、0.45396和0.45382，升力系数分别为-0.143、-0.211、-0.236、-0.193。阻力为负值说明小车受到下压力。后两个阻力系数已经很接近，而且第四种网格中已经存在质量不好的网格所以认为第三种网格为独立性网格。车前端面所受到的阻力最大，占阻力的绝大部分，顶部的阻力几乎为零，并且是负值，说明顶部受力方向为车的行进方向。选择第三中网格，其压力、速度、湍流动能以及速度矢量图如下，同时顶部也提供了小车的绝大部分下压力。 图13 压力图 图14 速度图 图15 湍流动能 图16 速度矢量 图17 速度矢量（部分） 压力图中，空气压缩在车头处，压力明显高于其他区域，车前部与顶部交汇处存在明显的压力梯度，在车顶部最高处压力最低，压力能转化为动能，车底部压力也较低，车后部压力较均匀，同时在车后部形成了一个高压区。速度图中，车头区域速度较低，车顶部和底部的空气流动速度很高。紧挨车后部空气流动速度很低，车后部在很大范围内速度比较小。车头与顶部的交汇处湍流动能最大，车顶部的前端较小。同时可以看到车后部产生了漩涡，车头处气流分别流向顶部和底部所以有一个低速区。 双车模型 两小车的间距分别为0.5、1、2个车间距，模型如图。 图18 双车模型 应用FLUENT分别计算。Y+值如图。 图19 0.5倍车距前车、后车Y+值 图20 1倍车间距前车、后车Y+值 图2