基于线性弹性体的网格变形方法研究

网格变形方法在计算流体力学的各个方面都发挥着极其重要的作用,是各种应用问题解决方案的基本组成部分.研究和实现了基于线性弹性体的网格变形方法,对基于弹性体法的网格变形方法进行了综述,从弹性方程理论、高效求解方法等方面进行了介绍,给出了方法的最优参数选择,并从基本变形能力、多段翼型设计、非定常动态特性模拟等方面,展示基于线性弹性体网格变形方法的研究和应用进展.结果表明:实现的线弹性体网格变形方法特别善于处理几何大变形,尤其是网格扭转比较大的应用,并且变形后网格品质高、鲁棒性强.     

非线性气动恢复力与飞行器动稳定性

复杂流动中,飞行器非定常气动数值计算后的参数辨识、飞行器模型的风洞定常测量、风洞动态(风洞自由飞)测量均显示了飞行器非线性气动刚度项(气动恢复力项、如俯仰力矩 攻角曲线Cm θ)将出现局部静不稳定( Cm/ θ>0)。　　飞行器气动恢复力对运动稳定性影响以往讨论较少,本文讨论了非线性气动恢复力局部出现静不稳定条件下的动稳定特性,其特点是突跃性、单边性与不确定性。飞行器飞行遥测数据分析支持了上述动稳定特性的描述。     

带控制舵弹体气动特性分析

本文用数值模拟方法研究了带控制舵弹体的气动特性。该问题外形复杂，在控制舵和弹身之间会出现明显的缝隙。为了解决此类外形的气动力计算问题，我们发展了结构网格和非结构网格的混合网格技术。这样一方面充分发挥了贴体结构网格质量好、效率高的优点，另一方面充分发挥了非结构网格几何灵活性好、易于模拟复杂外形的优点。数值计算与风洞实验相当吻合，说明本方法具有良好的精度。最后对两种外形进行大量计算并进行气动特性分析，给出了静稳定特性和配平特性的研究结果。     

低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。      

仿鸟柔性扑翼气动特性与能耗的数值研究

建立了适当的三维仿鸟柔性扑翼模型,并以配平重力和平衡阻力为条件,数值计算了它的低雷诺数非定常流场.研究揭示了翼面初始扭转角度、动态俯仰幅度等重要设计参数与飞行性能的关系,表明扑翼平面的初始扭转程度、扑翼柔性材料的选择以及两者之间的合理搭配对扑翼机的成功飞行至关重要.研究分析了仿鸟扑翼的流场涡结构、升力推力产生原理,下扑过程附着上翼面的前缘涡是升力产生的重要机制.对扑翼气动功率的比较分析也发现,人造扑翼机需要的气动功率明显高出同等大小的鸟类,在效率方面尚不及扑翼飞行生物.     

非对称激励对叶栅流场时空结构及特性的影响

基于非定常耦合流动的研究采用导叶非对称周向分布时,可以通过非对称尾流撞击效应实现比对称周向分布更好的气动性能.为了探索非对称激励下的非定常耦合流动特点,首先以单个槽道的时间非均匀射流激励形成的间歇耦合流动研究为基础,分析了时间非均匀激励对二维叶栅流场时空结构的影响.然后通过对采用非对称导叶的某低速轴流压气机的叶中截面二维流场进行1/4周的数值模拟,分析了周向非均匀分布导叶对下游转子流场时空结构的影响,最后通过实验对前面所得到的结论进行验证.认为通过周向非均匀导叶实现流场的耦合流动时,激励区的长度决定了激励效果的好坏,而缓冲区内加入离散导叶可以提升转子流场时空结构品质改善压气机性能.      

缝隙-腔体密封结构在高速气流冲击下的整体流动、传热特性分析

 根据高马赫数流场特征参数变化快、固态场特征参数变化慢的特性,编制开发了针对高速气流的"半解耦"显式流固耦合近似计算程序,并通过数值方法分析了该方法的计算误差,分析结果显示其能够较准确地模拟高速流场与固态场长时间非稳态耦合问题.在此基础上,运用"半解耦"流固耦合方法数值模拟了高速气流横掠缝隙-腔体典型密封结构的非稳态过程,并与相关实验测试数据进行了对比,验证了程序可靠性.随后,进一步分析了气流侵入密封结构的主要特性,总结了密封结构内、外流场中气流温度、压力和速度的分布特征以及其随时间的变化规律,研究了密封结构中加热板气动热流随时间的变化规律,探讨了密封结构中固体温度场分布特征及其随时间积累的变化规律等.最终,计算结果说明了密封体的结构布局对其内部热状况的决定性影响.     

纵列式双旋翼直升机气动力和流场的数值模拟

纵列式双旋翼直升机的旋翼尾迹受前飞自由流、机身和旋翼之间的干扰影响,旋翼尾迹不仅影响机身的空气动力,也影响旋翼自身的空气动力特性.应用欧拉方程与计算流体力学(CFD)方法,采用自适应非结构直角网格技术生成计算域网格,数值模拟纵列式双旋翼直升机全机流场,给出了考虑旋翼影响的机身表面压力分布以及全机流场特性.     

离心叶轮与扩压器相互干扰数值模拟

采用定常与非定常数值模拟方法,对某小型离心压气机进行数值计算,首先比较了定常结果与非定常时均结果的特性;然后结合离心叶轮出口和径向扩压器进口典型位置非定常压力脉动的频域特性,研究了离心叶轮与扩压器之间的转静干涉作用;最后着重分析了离心叶轮带与不带分流叶片对径向扩压器进口压力脉动频率特性的影响。结果表明:在相同流量下,非定常时均结果的压比偏高约3.7%,效率偏高0.1%;在带有分流叶片的离心压气机中,离心叶轮叶片出口与扩压器叶片进口的脉动除了叶片通过频率的非定常脉动外,分流叶片尾缘会产生与之强度相当的2倍叶片通过频率的非定常压力脉动。同时,分流叶片会对径向扩压器进口产生更强的0.5倍叶片通过频率非定常压力脉动,这些频率的非定常压力脉动很可能会对离心压气机叶片的疲劳强度造成影响。