气相多组分反应流场的一种新算法.pdfVIP

中国工程热物理学会 984073 燃烧学术会议 气相多组分反应流场的～种新算法球 任安禄，彭杰，戴志潜 (浙江大学力学系，杭州，310027) 0．摘要 本文提出在气相多组分反应流场中用时间推进法积分组分和运动方程．把混合气连续方程同状态 方程相联系，在同一时间步把密度的调整量转化为压力的调整量，计算绝对压力场。来获得气相多组 分反应流动的数值解。我们首先计算了单～气体在均匀温度场下的流场和大温差下的气流场，以证明 压力汁算方法的可行性。再计算在大温差下的有化学反应的多组分炉内流场，均获得成功。 关键诃：多纽分扩散方程绝对压力场计算 数值模亍cf { 本文受国家自然基助 1·前言 作过～些炉膛流场的计算工作．提出用粒子强度源模型[6](和STP模型近似[2])计 算颗粒相，用投影法计算流场，用区域法计算幅射传热。并且在流场计算中沿流线加 平衡粘性．这种方法对于回旋流有较高的精度[7，8]。但是如何把流场计算方法扩展到 热态反应流场计算中去却是十分困难的，而把它与SI、IPLE方法连接起来又不可能。所 以发展一种非SDIPLE方法的热态反应流场计算方法是十分必要的。 在热态流场计算中必须满足状态方程，而用于不可压流场汁算的压力场是相对压力 场，它不满足状态方程，已失去了作为流场变量的依据。流场变量需为绝对压力，但 它的量级为105，相对于其它变量，此量级容易引起计算的发散。为了解决这一困难， 我们采取了两种措施：一是无量纲化方程以缩小压力量与速度量等的量级差别；=是 通过状态方程把密度的调整量转化为压力的调整量，直接计算绝对压力场。用这种方 法我们成功地计算了常温下的流场．结果和以前用投影法的计算结果吻合很好，然后 我们又成功地计算了炉内大温差下单一气体的无反应流场以及有反应的热态流场(由 于只是原理性的，我们假设炉内只有五种组分，它们是：C02，CH。，0，，H：0，N：其 中Ⅳ，不参加反应)。因此．我们相信这种方法可以发展为一种不同于SIMPLE算法的 炉内流场、反应流场新计算方法。 从实际应用的角度来说．我们计算的结构形状比较简单的切圆燃烧疗型炉膛[6]并 且组分的选择仅限于原理性．但是，我们认为把这种方法推向实际运用是不困难的， 考虑污染物产：生也是完全可能的。 2-控制方程 ” ’ a iz ＼vax)0v＼’av)az＼vez) ≯分别代表】，Ufv，w，k，占．r，cpT，其中IJ，v，w分别为三个坐标方 向的速度分量，k为湍动能．s为湍动能耗散率，r是组分的质量百分比，C。为等压 气体比热容，T代表温度，只+S。分别为下列各项： O+Sp。 Ⅳ一176 一瓦dp+丢(比到Yx)+毒(儿期+象(“期q。 一 一OP_Pdy+丢(以期+号(，t期+尝(“期螺+瓦l以冽+万l以刮+瓦l“冽娟一 一笔+去(以期+号(雎期+差(肛期+‰ ∥ 肛 ，，．．．，。L 盟钙 手一矛 ¨习盟鹤 如4筹剖筹乜爿 咝+S。 一w；·dH。七S。 S，为煤粉燃烧产生的质量源项．S，。为气体产生或消耗的源项。其中，挥发份的 产生来自于傅维标的煤粉热解通用模型[9-11]。 Rl、／M 状态方程：P=p p为气体混合物密度，R为气体常数，T为温度，M为气体分子量。 通用热解模型： ，”，鲁=Ldaf--lllv!bexp(一