临近空间旋翼气动性能分析研究崔钊.doc

第二十八届（2012）全国直升机年会论文 1 二维翼型数值计算 低雷诺数下绕翼型的边界层流动粘性效应显著，边界层厚度与特征尺寸相比常常达到了一个不容忽视的程度，流动主要趋于层流特征，稳定性差，易发生分离。计算网格采用双曲线方程生成的C型网格，采用此种方法生成网格在壁面区域的网格单元正交性好，能够较好的模拟边界层流动，计算外边界为翼型的14倍弦长。三种网格密度进行计算比较，以排除网格密度对计算结果的影响，300×119的网格密度进行计算考虑低雷诺数下流场的气流粘性影响，控制方程采用定常、不可压纳维－斯托克斯（N-S）方程。层流粘性系数通过萨特兰公式给出，湍流粘性系数由湍流模型给定。湍流模型采用Spalart-Allmaras模型，压力-速度耦合采用SIMPLE算法。动量方程、能量方程、湍流粘度均采用二阶迎风格式耗散。入口边界给定速度，出口边界给定大气压强，翼型壁面条件为无滑移边界条件。 图1 翼型流场网格 图2 0.08马赫下升力系数 图3 0.08马赫下阻力系数湍流模型采用Spalart-Allmaras模型，。动量方程、能量方程、湍流粘度均采用二阶迎风格式耗散。入口边界，翼型壁面为无滑移边界条件。 图4 0.55马赫下升力系数 图5 0.55马赫下阻力系数 图6 升力系数随马赫数变化 图7 阻力系数随马赫数变化 临近空间两万米处大气密度仅为海平面的7.2%，三万米处的大气密度仅为海平面的1.5%，如表1所示。极低的大气密度使得旋翼桨叶翼型所处的雷诺数明显低于海平面的水平，翼型升阻特性变差；而旋翼产生的拉力显著小于海平面处产生的拉力。 表1 设计环境参数 巡航高度 大气温度 大气密度 大气动黏度 大气压强 当地音速 0 288 1.225 101325 340 20000 217 0.088 5475 295 30000 227 0.018 1172 302 图8 桨尖0.5马赫时旋翼桨尖雷诺数随高度变化 图9 桨尖0.6马赫时旋翼桨尖雷诺数随高度变化 本别采用本模型与未考虑旋流的影响和未计入小角度假设影响的气动模型，计算临近空间旋翼的升阻特性，计算结果如下： 图10 旋翼升力系数比较 图11 旋翼扭矩系数比较 图12 不同桨叶半径下旋翼升力系数随实度变化 图13 不同总距下实度对功率系数的影响 图14 不同总距下旋翼实度对旋翼功率载荷的影响 图15 不同总距下旋翼实度对旋翼桨盘载荷的影响 CX10B_104Z）资助，特此表示衷心感谢。 参 考 文 献 Young L.A., Chen R.T.N., Aiken E.W., and Briggs.G.A.. Design Opportunities and Challenges in the Development of Vertical Lift Planetary Aerial Vehicles[C]. American Helicopter Society (AHS) Vertical Lift Aircraft Design Conference, San Francisco, CA, January 2000 Young. L.A., Briggs.G.A., Derby.M.R., and Aiken. E.W.. Use of Vertical Lift Planetary Aerial Vehicles for the Exploration of Mars[R]. NASA Headquarters and Lunar and Planetary Institute Workshop on Mars Exploration Concepts, LPI Contribution # 1062, Houston, TX, July 2000 Young. L.A.. Vertical Lift -- Not Just For Terrestrial Flight[C]. AHS / AIAA / SAE / RaeS International Powered Lift Conference, Arlington, VA, October 2000 Young. L.A., Johnson.J.L., Demblewski. R., and Andrews.J.. New Concepts and Perspec