发动机进口支板过冷水滴撞击特性的数值模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083513 发动机进口支板过冷水滴撞击特性的 数值模拟 范训 吉洪湖 胡娅萍 曹广州 （南京航空航天大学能源与动力学院205 教研室 210016） 联系电话：0252205，Email: fanxunking-666@163.com 摘要：本文采用数值模拟的方法对航空发动机进口支板的过冷水滴撞击特性进行了数值模拟。首先对 空气－过冷水滴两相流动的空气流场进行了求解，在此基础上，将过冷水滴以离散相模型的形式加入 流场，与空气相进行耦合求解，跟踪其运动轨迹，得到水滴的撞击特性。计算结果表明不同的过冷水 滴直径和水滴速度对于水滴撞击极限和水收集系数有较大的影响。 关键字： 发动机进口支板 离散相模型 水滴撞击特性 数值模拟 1 引言 飞机在含有过冷水滴的云中飞行时，当达到结冰气象条件时，飞机的迎风部件如机翼以 及发动机的进气道、进气部件和动力装置均会发生结冰。发动机进气道及进气部件的结冰， 指的是进气道前缘、发动机压气机前的整流罩、支撑及第一级压气机前的导流叶片等部件的 结冰[1]。在涡扇发动机中，风扇出口的导向叶片结冰，会导致气流的阻力增加，使得风扇的 操作线达到滞止。压气机导向叶片结冰，会导致气体流量降低，流体的加速度减慢，当结冰 脱落时会破坏下游的机构，导致发动机熄火[2] 。因此有必要针对发动机部件的结冰进行研 究。在结冰数值模拟中，首先必须计算结冰区域和结冰量，这主要取决于水滴撞击特性的分 析研究。目前国内外，对于水滴撞击特性的数值计算方法主要有两种：拉格朗日法和欧拉法 [3] 。而针对机翼的水滴撞击特性的研究已经比较成熟。对于发动机部件的研究还比较少， 本文利用拉格朗日法针对发动机的进口支板进行了水滴撞击特性的研究。 2 控制方程 在用拉格朗法日建立控制方程时，做了如下假设： ① 过冷水滴与空气间有速度滑移； ② 把过冷水滴视为分散群，考虑其湍流扩散； ③ 过冷水滴按一定初始位置沿着各自的轨道运动，彼此互不干扰，可追踪其沿轨道的速度、 温度以及水滴直径的变化，在拉格朗日坐标中考虑水滴影响； ④ 把过冷水滴对气相的质量、动量相互作用当作连续分布的物质源、动量源加入到气相方 程中。 空气控制方程： 在笛卡儿坐标系中，定常的雷诺平均 N-S 方程： 连续性方程： ∂ ρu m (1) ( ) i l ∂x i 动量方程： ∂(ρu u ) ∂p ∂ ∂u ∂u i j i j F (2) =− + [μ( + )]+ j x x x x x ∂ ∂ ∂ ∂ ∂ j i j j i 式中，ml 为质量源项，Fj 为动量源项。 水滴运动轨迹，通过积分拉氏坐标系下的水滴作用力微分方程来求解得到。水滴的作用力平 衡方程在笛卡儿坐标系下的形式为： du g (ρ −ρ) p F u u x