网格划分对柴油喷雾油束形状多维数值模拟影响地研究.pdfVIP

网格划分对柴油喷雾油束 形状多维数值模拟影响的研究 赵昌普苏万华张晓宇 天津大学内燃机燃烧学国家重点实验室 摘要用CFD多维数值分析软件对柴油喷雾的油束特性进行了研究，计算结果 表明计算网格的划分直接影响到计算结果的正确性，影响到油束的几何形扶。在建 立计算网格时，一定要注意网格的均匀一致，避免网格密度的突变；特殊几何形状 壁面附近的局部贴体坐标网格的范围不宜太大，否则会引起计算结果的畸变。 关键词 柴油机 喷雾特性 计算网格 数值模拟 内燃机缸内气流运动、喷雾及燃烧过程的改善是降低排放、改善柴油栅陛能的重要方面， 币j数仿模拟是揭示其详细物理过程的重要手段。因数值摸拟具有成本低、速度快、资料完备 以及模拟真实条件和理想条件的能力，其应用日益广泛cI。】。但计算结果取决于计算网格及 数学模型，如果所采用的计算网格或数学模型不合适，那么，即便是一个完善而又令人满意 自I数值方法，也只能产生毫无价值的结果。本文旨在探索计算网格划分方式对柴油喷雾模拟 计算结果的影响，以便合理地建立柴油喷雾模拟计算计算网格、准确预报柴油机喷雾混合过 1计算方法 柴油机的喷雾、燃烧及流动过程属于瞬态问题，故采用具有两个压力校正步、对压力场 校正精度高．且比较容易得到收敛解的PISO算法。PISO算法的基本思想和结构都与SIMPLE Solution 非常棚似。其名称由PressureImplicit bySplitOperator的字头组成，主要用于非稳态 汁算，有时也用于相关变量间存在强耦合的场合(例如浮力流)。PISO的计算步骤如下：(1) 估汁压力场p聿，并估计一个迭代初始速度场；(2)求魑动量方程，得到速度场u+，V8，w8； 7和p2， 7方程，得到p1。，用Pl’校正速度值；(4)求解p27方程，得到p2，1用pt (3)求解Pl 校j}：压力，即p：p++p17+p27；(5)求解那些通过源项、流体物性等影响流场的其它变量(如 温度、浓度等)的离散化方程： (6)将新求得的速度场及新的物性、源项等代人动量方程， 用新得到的压力P当作一个新的试探压力场p+，返回第二步，重复全部过程，直到得到收敛 斛； 【7)得到收敛的速度场后，求解剩下的其它待求的变量。其计算流程图参阅文献【6】。 041094 2网格划分对计算结果的影响 计算网格的划分是否合理直接影响到计算结果的正确性‘”，在计算过程中必须高度重视。 下面用图l(a)、(b)、(c)三种不同的网格对喷雾撞壁过程进行了模拟计算。表1为模 拟计算主要参数，图2为计算结果。计算域的初始网格数为30×70×30，网格边长为lmm。 为了提高计算精度，揭示限流沿(Bump)【8’附近及附面层的详细物理过程，对Bump附近及 附面层的网格进行了细化，壁面垂直方向网格最小尺寸达0．125ntm。喷雾模型使用Huh模型， Ultra 液滴破碎模型使用Reitz／Diwakar模型，撞壁模型使用Bai模型9】。在Sun80工作站 (UltraSPARC--11450双CPU：4M高速缓存；双18GSCSI磁盘；2G内存；中文操作界面， Solaris 8操作系统)上使用多维数值分析软件STAR，CD进行计算，时间步长及计算步数分 8小时。 别设为0．0lms和250步，三种计算网格分别需CPU时间大约19、2l及1 三维透视圈 二维局部图 三维透视图 二维局部图 三维透视图 二维局部图 a)原计算网格 (b)Bump周围细化后网格 (c)重划局部贴体网格 囤l 三种不同的计算两格 由图2(a)可以看出，当采用图1(a)的正方体网