基于STAR―CCM+的某型车SCR系统流场分析.docVIP

基于STAR―CCM+的某型车SCR系统流场分析 　　一、引言 　　随着社会的发展，汽车的保有量快速增加，由此引发的空气、噪声污染也越来越引起社会的关注，汽车的SCR后处理系统在设计过程中对排气背压的要求十分严格，因为后处理系统的背压不仅与车辆的排气噪声有关，而且关乎到车辆的经济性和动力性，如果排气背压过大，会使排气不畅，导致气缸内气体无法正常排出，影响进气量，使发动机性能下降，进而影响整个车辆的油耗经济性。如果靠传统的试验方法，不仅耗时耗力，也难以同时达到很好的效果，往往是顾此失彼。如果借助数值模拟的方法，先从理论的角度分析SCR系统包括压力场、速度场等各种性能，可以大大缩短开发周期，提高开发效率，同时避免设计过程的性能缺陷。文章利用CFD方法，对某型车SCR系统总成进行了流场模拟计算，得到了SCR系统排气总背压，并分析了各个部件的背压及不同尾气流量下的背压，为该型车SCR系统结构布置和设计提供预测性的意见。 　　二、SCR系统流场分析控制方程 　　车辆SCR后处理系统内部结构一般比较复杂，内部流场也相应地不规律，所以属于紊流流场问题。文中选用STAR-CCM+中提供的Realizable k-ε湍流模型进行数值计算，湍流控制方程为三维不可压缩雷诺时均Navier-Stokes方程。 　　由于在对Navier-Stokes方程进行雷诺时均处理过程中，引入了与湍流脉动值有关的新的变量项- （雷诺应力项），为了使方程组封闭，必须对雷诺应力做出假设，即建立应力表达式或引入新的湍流模型方程。通过这些表达式或湍流模型，把湍流的脉动值与时均值联系起来。研究表明，标准 k -ε模型应用于时均应变率特别大的情况时，有可能导致负的正应力，为了使流动符合湍流流动的规律，需要对正应力进行某种数学约束，从而提出了Realizable k-ε模型中关于k和ε的运输方程。 　　三、模型建立及边界条件 　　1.模型建立 　　计算模型为某型车SCR后处理系统，采用CATIA建立其三维模型，如图1所示。车辆尾气进入排气管，再通过壳体内薄壁多孔管流出，经过薄壁隔板的两个孔流出后到达消声陶瓷蜂窝载体上端所在的腔体，从消声陶瓷蜂窝载体的上端流入下端，最后通过排气尾管排放到大气中。为了便于考察模型细节处的流场情况并保证其计算的准确性，尽量保证计算模型与实际三维模型的一致性，并采用1∶1实际模型进行数值计算。前处理在ANSA软件中进行处理，最终形成Triangle网格数为1 70222，SCR系统的面网格模型如图2所示。 　　在ANSA软件中将SCR系统的面网格转化为stl格式，导入STAR-CCM+中进行下一步体网格的生成。在STAR-CCM+中指定网格尺寸和边界层数及厚度，由于STAR-CCM+特有的多面体网格具备很多优势，既保证计算精度，提高计算的收敛性和稳定性，又能控制网格总体数量，节约计算时间，所以网格类型选择Polyhedral多面体网格，最终生成的体网格数约为284456，如图3所示。 　　2.消声陶瓷蜂窝载体 　　由于陶瓷载体管道太细、太多，如果这些管道生成体网格进行流体计算，计算机容量难以承受且没有必要。从工程角度来看，人们更多的是关心载体整体性能及对气流的影响。由于蜂窝载体的管道具有相同的几何形状，且管道的直径远小于载体的直径，载体的气流可以认为是定常不可压层流，可以把蜂窝载体用当量连续法进行模拟，即把管道内的气流在整个载体内看作连续分布的气流。载体内的流动假设为多孔介质内的流动，流体只沿轴向流动，不存在径向质量交换。在STAR-CCM+中，多孔介质的压降公式可表示为： 　　四、额定工况SCR系统计算结果分析 　　1.边界条件的设定 　　该型车发动机额定工况：发动机转速1 900r/min，扭矩1 267N?m，尾气流量1 465kg/h。SCR后处理系统在数值模拟过程中，排气管道来流方向的横截端面为入口边界，设定入口边界为流量入口，流量为 g=1 465kg/h，湍流强度为0.01，湍流长度为0.01m；出口边界为压力出口，出口压力为p=0Pa（相对于大气压）；消声陶瓷蜂窝载体的惯性阻力系数为3.842kg/m 4，粘性阻力系数为718.764kg/（m 3s）SCR后处理系统固壁面采用无滑移边界条件，求解过程中各松弛因子：速度项为0.7，压力项为0.3。 　　2.计算结果及分析 　　（1）表面压力分析。 　　通过STAR-CCM+软件求得该型车SCR后处理系统在额定工况下的排气背压为6.682kPa。图4为该模型整体表面压力分布云图。图5为SCR系统排气管的表面压力分布云图。 　　由图4、图5可以看出，SCR后处理系统在排气管有转弯外侧位置表面压力较大，这是由于尾气从发动机排出后高速进入排气管，遇到转