大水滴撞击壁面的动态特性数值模拟

采用 VOF 方法建立了大水滴撞击壁面的计算模型,模拟了大水滴以不同直径、不同速度撞击光滑壁面的动态撞击过程和撞击特性。计算结果表明：在大韦伯数情况下,水滴撞击光滑壁面会在铺展过程中发生边缘水滴飞溅;在水滴撞击壁面的收缩阶段,随着水滴直径的减小和水滴速度的提高,会发生液膜缓慢收缩、边缘液环和液膜分离、中心处部分液膜和边缘液膜分离、液膜完全破裂等不同结果;当水滴直径和撞击速度增大时,同一时刻水滴的铺展半径、最大铺展半径、最大铺展系数增大;水滴在壁面上达到最大铺展系数所用的无量纲时间随水滴直径增大而增大,同一直径水滴不同初始速度下达到最大铺展系数所用的无量纲时间变化较小。     

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

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

探头式结冰探测器安装位置分析

通过求解雷诺时均Navier-Stokes(N-S)方程获得空气流场,采用欧拉法求解水滴运动场获得水滴运动轨迹.通过计算来流液态水含量等值线距离部件表面的垂直距离,获得部件不同位置上的广义水滴遮蔽高度.研究了不同飞行条件和气象条件下的广义水滴遮蔽高度大小.根据探头式结冰探测器安装位置要求确定探头式结冰探测器的合理安装区域.通过分析不同的状态参数对广义水滴遮蔽高度的影响,为结冰探测器安装位置计算状态参数的选取提供了数据支持.      

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

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

运输机翼型结冰的计算和实验

采用数值计算和结冰风洞实验相结合的手段,对某运输机机翼剖面缩比翼型的结冰特性进行了研究.结冰实验在气动中心0.3m×0.2m结冰风洞中进行,实验模型为弦长0.18m的层流翼型.翼型的绕流流场通过求解低速黏流的时均Navier-Stokes(N-S)方程得到,采用拉格朗日法计算水滴撞击特性,在此基础上,求解结冰热力学模型...     

水滴撞击特性的重力影响分析

以NACA0012翼型为研究对象,采用拉格朗日法分析了重力对水滴撞击特性的影响,获得了过冷大水滴（SLD）的撞击特性。计算结果表明：当水滴直径超过50μm时,重力对水滴轨迹以及极限轨迹的初始位置影响很大;SLD的水收集系数和收集区域都较直径为10～40μm的水滴（美国联邦航空条例FAR25部附录C规定的防冰系统设计用水滴直径范围）有较大幅度的提高,该结论可为飞机防冰系统的设计以及冰风洞中喷嘴安装位置的确定提供参考。     

Roe格式在水滴欧拉方程数值求解中的应用

目前普遍采用欧拉法求解水滴运动控制方程以模拟结冰表面水滴撞击特性。水滴欧拉方程近似没有压力项的空气欧拉方程,但已应用于空气欧拉方程的空间离散格式多数不能直接应用于水滴欧拉方程,有研究者将经典Roe格式应用于水滴欧拉方程时出现了奇点问题。通过对水滴欧拉方程进行形式改造,基于中心格式验证改造方法的有效性,并应用Roe格式结合线性重构与限制器函数实现对水滴欧拉方程二阶精度的求解。将NACA23012翼型计算结果与FENSAP-ICE软件进行对比,两者水滴收集系数分布基本一致。     

过冷大水滴相继撞壁对结冰影响的实验研究

目前过冷大水滴（SLD）撞击结冰机理未能清楚解释SLD间的相互干扰对撞击结冰过程的影响。实验用高速相机拍摄双SLD偏移一定位置相继撞击壁面结冰的物理过程，研究双SLD落点间的偏移、相继撞击的时间间隔对结冰形态及结冰速度的影响。实验表明:双SLD间的相互干扰会抑制撞击后的回缩行为进而影响结冰形态，降低结冰速度，延长双SLD完全冻结时间:双SLD落点位置偏移量的大小对结冰形态和结冰速度的影响均很大，相继撞击时间间隔的长短对结冰形态影响不大，主要影响结冰速度。     

基于分段常值推力的水滴悬停构型控制策略

针对航天器编队飞行任务对相对运动控制的要求,研究了在分段常值推力控制下航天器受迫绕飞构型的设计与控制问题。首先,基于脉冲控制下的水滴悬停构型,提出了多段常值推力控制实现水滴悬停构型的打靶方程;将打靶方程转化为求解极值问题,采用最小二乘法来求解;分析了一段常值推力可行性。然后,以连续常值小推力控制方程为基础,推导了小邻域定理,分析了近距离相对运动条件下两段常值推力控制的可行性;针对可能出现求解精度差的问题,提出了小推力增量方程来修正精度,并证明在靠近理想解的情况下多次迭代可以趋近于理想解。最后,通过数值仿真实现常值小推力控制下的水滴悬停相对运动。数值仿真结果表明常值小推力控制策略可行,研究成果完善了航天器受迫绕飞构型设计与控制的相关理论,为工程应用提供参考。      

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

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