旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

翼型表面明冰的数值模拟

提出一种翼型表面明冰的数值模拟方法.基于欧拉两相流理论,建立空气、水滴两相流控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击;采用考虑了粗糙度影响的边界层积分法求解结冰表面的对流换热特性,基于经典Messinger模型提出明冰数值模拟的热力学模型;利用冰层时间推进法模拟积冰过程,采用积冰法向生长假设生成积冰外形.对NACA 0012翼型在不同环境温度下的明冰冰形进行了预测,并与实验结果进行了对比.      

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.      

多段翼构型的积冰数值模拟

基于非结构网格提出一种针对多段翼构型的积冰数值模拟方法.通过求解气体欧拉方程得到多段翼的绕流流场,基于空气绕流流场的气流速度分布,采用迭代法对水滴控制方程进行数值求解,得到多段翼各翼面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对某典型三段翼霜冰和明冰的结冰情况进行了预测研究,根据多段翼的流场特点提出一种明冰冻结模型的求解思路,并分析了积冰对多段翼性能的影响.      

求解可压缩流动的低耗散E-CUSP迎风格式性能分析

采用对比分析的方法,对求解可压缩流动的低耗散E-CUSP(LDE)迎风格式的性能进行测评,目的是进一步更好地将其用于分离涡数值模拟。对比格式采用广泛使用的Roe格式。控制方程为Euler和雷诺平均Navier-Stokes方程耦合Spalart-Allmaras一方程湍流模型,采用全隐式无分裂Gauss-Seidel线松驰时间推进方法。结果表明,LDE格式具有强的捕捉激波能力,计算结果与Roe格式一致,和实验数据吻合良好。相对于Roe格式,LDE格式具有较快的收敛速度。     

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程;四阶龙格一库塔法求解水滴运动轨迹方程;逐次二分法寻找水滴撞击极限,然后计算局部水滴收集率;分析翼型表面的传质、传热项,建立相应质量、能量方程并迭代求解,得到翼型表面结冰质量;根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好,结冰过程中的特征参数与冰型及相应冰形特征吻合。     

结冰面水滴收集率欧拉计算方法研究及应用

发展了一套用欧拉法计算飞机结冰表面水滴收集率的算法。先用CFD方法计算空气场,在获得空气流场分布的基础上,求解水滴相控制方程,进而得到水滴收集率。空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制。计算了典型条件下圆柱和NACA0012翼型的水滴收集率,得到了与实验或文献一致的结果,初步表明本文发展的有限体积算法以及相应的计算代码是有效的。最后给出了本文方法在三维复杂外形水滴收集率计算中的应用,对三段翼、T型尾翼的水滴收集率进行了计算,获得了相应的分布规律。     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

冰风洞试验中水滴的传热传质计算研究

进行冰风洞试验时,出口过冷水滴的参数往往使用风洞来流空气参数进行计算,忽略了真实情况下水滴与空气的传热传质过程,致使试验数据产生偏差,并影响测试结果的可靠性.通过对冰风洞试验中水滴的运动过程进行研究,考察了水滴与空气主流间的传热与传质现象,建立了水滴运动过程中参数变化的控制方程.基于该方程编制了冰风洞水滴粒径温度变化分析软件,计算分析了过冷水滴在行进过程中相关参数的变化曲线,比较了不同初始条件及各参数对水滴温度、直径和速度的影响.结果表明:不同环境下水滴在运动过程中温度、尺寸和速度变化受来流温度、速度、相对湿度和水滴初始温度的影响,呈现出不同的特点,不能忽略水滴状态参数在冰风洞内的变化.     

跳点有限体积法及其在跨音速流场计算中的应用

[1]和[2]曾用几种差分格式求解一维N-S方程和欧拉方程,计算表明,Hopscotch格式(跳点格式)计算速度最快。为探讨这种格式在多维流场特别是跨音速流场计算中应用的效果,本文将此差分格式应用于有限体上,成为跳点有限体法,并用该方法求解二维欧拉方程组。对若干个喷管和叶栅所做的计算得到了满意的结果,计算时间较目前常用的差分格式少得多。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

过冷大水滴动力学特性对结冰影响数值研究

对比常规水滴,过冷大水滴结冰对飞机气动性能和飞行安全具有更大的危害性。在常规水滴结冰过程数值模拟基础上,针对过冷大水滴条件下各种动力学特性及其对成冰过程的影响进行了数值研究。采用结构化网格和中心有限体积法求解 N-S 方程获得空气流场,用拉格朗日法求解水滴流场,基于改进的 Messinger 热力学模型完成结冰过程模拟;对于过冷大水滴,采用泰勒类比理论和水滴碰撞模型,对比分析了水滴动力学特性及其对结冰过程的影响;通过数值模拟获得了水滴 Langmuir D 分布的撞击特性及结冰计算结果,研究了水滴多尺度分布对成冰过程的影响规律。通过一些典型结冰算例的数值模拟,并与参考文献结果及实验数据的对比分析,说明了本文计算方法及计算模型模拟过冷大水滴结冰是可行和正确的。     

过冷大水滴飞溅特性数值分析

针对过冷大水滴（SLD）撞击结冰表面的飞溅现象,开发了基于欧拉法的二维算法程序,该程序将LEWICE飞溅模型和FENSAP飞溅模型耦合到水滴流场控制方程中进行求解,实现了对SLD撞击特性的准定常计算。根据水滴飞溅模型中对于飞溅生成子水滴质量和速度等信息的定义,通过定义控制方程中的源项,将其耦合到控制方程的求解中,进而研究水滴飞溅效应对于撞击特性的影响。研究发现,由于水滴的飞溅作用,导致翼型壁面附近的液态水含量（LWC）分布发生变化;通过对上述两种典型的SLD飞溅模型进行计算分析,发现了SLD撞击特性的一些共性和不同,并揭示了SLD撞击效应的一般特点。