论文部分内容阅读
微型飞行器因尺寸小、重量轻、隐蔽性好等特点而具有良好的应用前景,其概念一经提出即迅速成为科学界当前研究的前沿和热点课题。微型飞行器不是常规飞行器的简单缩小,其尺寸的微型化带来许多亟待解决的难题,低雷诺数空气动力学即为其中最为重要和最为关键的问题之一。本文的主要内容是固定翼/扑翼低雷诺数气动特性的数值和试验研究,包括非定常低雷诺数流场的数值计算方法、动态网格生成、固定翼/扑翼低雷诺数流场机理等研究内容。传统可压缩流计算方法用于求解低速流场时会面临收敛慢、精度低的“刚性”问题,因此,本文将预处理技术拓展到非定常可压缩N-S方程的求解中,发展了一种能有效模拟含运动边界的三维非定常低雷诺数流场的数值方法。控制方程在非结构动态网格上用双时间步推进并用无矩阵计算的隐式LU-SGS迭代求解;由于预处理改变了控制方程的特征值和特征向量,也提出了一种特征变量的边界条件以适合预处理后的特征系统。众多典型算例的结果表明,预处理方法将可压缩流方法的计算能力拓展到了低马赫数范围,有效地提高了收敛速度和计算精度。对运动边界问题,动态网格生成是其非定常流场计算的一项重要工作。基于Delaunay图映射的动态网格变形方法无需迭代计算,效率高,稳定性好。但对大幅运动复杂构形的动边界问题,背景图容易交叉,重新生成背景图和重新定位网格点不仅费时而且会导致网格质量的严重下降。本文提出一种双重Delaunay图映射的动网格生成方法,通过在初始背景图中添加辅助点,生成虚实两种背景图和虚实两种映射关系。分别根据虚映射关系和虚背景图、实映射关系和实背景图,移动辅助点和网格点。复杂构形的扑翼算例表明,双重图映射方法多付出极少的内存代价即可避免背景图交叉引起的问题,提高了动网格生成的效率和质量,增强了处理大变形复杂动边界问题的强健性。本文也发展了一种非结构嵌套网格方法,并提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据壁面的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。数值模拟了小展弦比固定翼微型飞行器的低雷诺数流场,考察了小展弦比机翼在低雷诺数条件下的升阻力特性、流场流动机理,分析了翼尖涡和翼梢小翼对气动特性的影响,并对大攻角下流场的非定常现象进行了相关的讨论。数值研究的结果表明,低雷诺数小展弦比机翼的流场具有较复杂的三维效应,翼尖涡的存在对流场结构和升力、阻力的产生都有重要影响,涡的脱落和相互作用使得流场在大攻角下具有一定的非定常特征。建立了适当的三维仿鸟柔性扑翼模型,并以配平重力和平衡阻力为条件,数值计算了它的低雷诺数非定常流场,研究分析了仿鸟扑翼的流场涡结构、升力推力产生机制和气动功率。研究了二维柔性扑翼在不同斯德鲁哈尔数下的推力特性以及流场尾涡的分布特点。研究认为,扑翼平面的初始扭转程度、扑翼柔性材料的选择以及两者之间的合理搭配对扑翼机的成功飞行至关重要;下扑过程附着上翼面的前缘涡是升力产生的重要机制。数值研究了二维前后置双扑动翼的非定常流场,考察了前后翼扑动的相位差、水平距离、频率和攻角等参数对其气动力和气动效率的影响,并对前后翼之间存在的气动干扰进行了分析。蜻蜓能通过改变前后翼间的相位实现对不同飞行姿态的控制。本文数值模拟了三维蜻蜓模型在悬停(前进比J=0)和中等前飞速度(J=0.3)两种状态下的流场,每种状态各计算从0o至360o间隔30o的13个不同相位。考察气动力和气动功率随相位的变化关系以及前后翼间的气动干扰,也分析了流场涡结构和气动力的产生机制。结果发现,平均举力和平均气动功率随相位的变化呈U型分布,举力足够用于平衡重力,气动功率也符合蜻蜓真实飞行的统计数据。90o~270o的较宽相位范围内,两翼间的气动干扰较为强烈且比较稳定,该范围内能量消耗较小,功率最省,在一定程度上解释了蜻蜓翼具有特定相位的现象。对仿鸟扑翼微型飞行器进行了流场显示试验。在低速风洞中,采用高速摄影仪实时拍摄了扑翼机的烟线流场,获得了对仿鸟扑翼非定常流场的直观认识。