谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.     

有限振幅扰动下高超声速钝楔绕流感受性

为了研究自由流脉冲波作用下高超音速流动感受性特征,采用高阶有限差分方法模拟钝楔超高音速绕流中来流慢声脉冲波作用下的非定常流场,并应用傅里叶频谱分析了边界层扰动模态的演化.结果显示,自由流脉冲扰动下,边界层不同模态波沿流线呈现不同变化规律.球头区边界层感受到的主要是基频波,在球头向非球头区过渡段,主导模态迅速向高频迁移;在非球头区,主导模态频率缓慢向高频迁移.总的来说,沿流向总扰动模态中的低频模态成分变小,高频,特别是二阶谐频附近模态成分显著地增大.沿流线发展,边界层不同扰动模态之间存在伴随能量迁移的模态竞争物理现象,主导模态的发展对其他模态有很强的抑制作用.      

跨音压气机叶栅流场的附面层修正

本文发展了二维压气机叶栅亚音、跨音流中粘性——无粘相互作用程序。无粘计算在弦向大扰动、弦法向小扰动的假设条件下用跨音松弛法求解;粘性附面层用Nash方法计算。考虑了叶栅尾迹的影响。最后,连同通道区和尾迹区一起,求得了包括冲波影响的叶栅的总损失系数。计算结果表明,在跨音情况下,粘性对冲波强度、位置以及全流场的速度分布都有不可忽视的影响。     

共轴刚性旋翼直升机着舰飞行特性

建立一套能够耦合非定常舰船流场的共轴刚性旋翼直升机飞行动力学模型,以研究该旋翼飞行器在着舰过程中的飞行特性。在舰艉流场模拟方面,采用分离涡方法以获得高精度非定常流场数据;在飞行力学建模方面,引入上下旋翼干扰因子建立共轴刚性旋翼诱导速度模型,并采用等效挥舞运动概念建立其挥舞运动方程;基于“单向耦合”思想构建了计算流体动力学（CFD）向飞行力学模型的数据传递策略,并分别以XH-59A直升机和UH-60A/SFS2组合为算例验证了飞行动力学模型和数据传递策略的正确性。以SFS2舰船模型和XH-59A共轴刚性旋翼直升机组合为研究对象,从直升机操纵余量和非定常载荷水平两方面分析了着舰过程中舰艉流场对共轴刚性旋翼直升机的扰动特征。时均研究结果显示：由于共轴上下旋翼受时均流场的扰动存在差异,飞行员在增大总距以维持高度稳定的同时,还需要减小差动总距以保持机头朝向的稳定。非定常水平研究结果表明：对于共轴构型直升机,舰船流场对拉力及俯仰力矩的非定常扰动,是引起飞行员工作载荷增大的主要因素。     

叶轮机械中转子角加速度对非定常流场的影响

本文推导了存在转子角加速度时,描述叶片排三维非定常流场的无粘基本气动方程组。并给出一些简化及变换形式。     

火焰筒头部流场研究之二——二元轴对称等温廻流区计算

本文将火焰筒头部的等温流场,用二元轴对称不可压定常流来处理并用数值方法解出初步结果。基本方法是用涡量流函数法,双方程粘度模型,采用交替方向隐含格式差分求解,在边界及迭代初值提出自己的看法,初步计算结果与[10]之实验定性一致。     

浮空气囊定常及非定常气动特性计算

对于0~20 km范围内不同外形的浮空气囊,采用虚拟压缩方法求解不可压缩非定常N-S(Navier-Stokes)方程,数值模拟了5个高度下气囊的定常和非定常流场.定常状态时,随高度及气囊外形的变化,模拟了气囊气动力及背风区复杂流动分离.结果表明,当气囊厚度率逐渐增大时,其背风区存在复杂的空间流动分离.运用双时间步方法及动网格技术,成功实现了气囊平移或俯仰运动状态下非定常流场的数值模拟.考察了不同的平移方向及速度对气囊非定常气动特性的影响,得到了气囊平移过程中气动力及流场的变化情况.对于气囊的俯仰运动,获得了不同俯仰速度及不同最大俯仰角时气囊的气动力及流场.计算结果表明气囊的非定常运动对其流场和气动力影响显著.     

一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。      

用于叶轮机械流场计算的交错网格有限体积法

在对各种交错网格离散格式的奇偶不耦合数值误差进行了详细的傅里叶分析和比较后,采用一种紧致的交错网格有限体积法将纳维-斯托克斯方程和k-ε湍流模型方程紧密的耦合在一起求解,发展了一套三维定常流程序.通过对NASA Rotor 67设计和非设计工况内部流场的数值模拟结果与实验的比较、分析,初步验证了程序的准确性和可靠性.      

新型主动液压滤波器的设计与分析

针对目前主动液压滤波器体积大、结构复杂等问题,提出了一种占用空间小、结构简单的新型主动液压滤波器. 采用计算流体动力学在时域内对滤波器所载流体域进行了压力分布计算,分别对三组可控参数激励位移幅值D_m及频率f_m作用下的流体域进行分析,通过比较分析发现,D_m和f_m的选择决定着该液压滤波器性能的优劣,如果二者与原始非定常流幅频特性匹配较好,则该液压滤波器滤波效果良好,反之,会使系统压力脉动变的极不稳定,这与波的干涉原理相互吻合. 算例计算结果表明,当激励频率与初始非定常流频率重合达525 Hz,且激励幅值达0.05 mm时,系统的压力脉动峰峰值可从2.26 MPa降到0.25 MPa,即脉动从±5.38%降至±0.6%,说明了该主动液压滤波器的有效性.     

绕丁坝流动结构实验研究

利用流动显示和片光源技术,对绕丁坝的流动进行了实验研究,以便了解丁坝参数对流态的影响.实验是在水槽中进行的.发现在丁坝的上游会形成马蹄涡,在丁坝的背风面会形成卡门涡.按雷诺Re数的大小可以把绕丁坝的流态分成4种类型:定常单个马蹄涡、非定常单个马蹄涡、非定常双马蹄涡Ⅰ及非定常双马蹄涡Ⅱ.实验表明丁坝的交角及丁坝的个数对丁坝的流态有重要影响.     

激波碰撞干扰流动非定常效应的数值研究

将AUSM+上风格式应用于数值模拟半球绕流激波碰撞干扰的流动中,揭示了该流动在定常来流条件下的非定常本质.分析弓形激波后超音速"喷流"的流动结构及其非定常效应的形成机制,提出一种新的激波碰撞干扰流动的非定常形成机理.      

轴流压气机内导叶/转子干涉对转子流场的影响

叶轮机械内的转静干涉现象多年来一直是一个研究的热点,然而现有的对于转 静干涉的研究多局限于对流动现象的描述,缺乏理论指导,因此无法在叶轮机械的设计体系 中包含这一非常重要的非定常流动现象,限制了叶轮机械性能的进一步提升.利用流动稳定 性及感受性理论来指导对这一问题的研究,能够找到利用这一现象提升叶轮机性能的规律. 对于一台单级低速轴流压气机内的导叶/转子干涉现象进行了实验研究.利用二维粒子图像 速度仪测量得到了导叶/转子干涉的定量流动细节,及其对于转子流场的影响.结果表明, 转子流场在定常来流条件下也表现出显著的非定常性.在引入导叶/转子干涉之后,转子流 场结构发生显著变化.动叶尾迹强度明显减弱,叶背流动分离得到抑制,叶尖泄漏涡也得到 了周期性的抑制.除了导叶尾迹对动叶流场的干涉作用以及气流角的周期性变化以外,动叶 非定常流动对于不同的导叶/动叶干涉激励频率的不同响应是引起压气机性能变化的主要原 因.      

基于小扰动理论的桨叶叶素气动载荷计算方法

针对直升机飞行仿真应用,采用小扰动理论建立了桨叶叶素气动载荷计算方法。该计算方法集成了弹性桨叶挥舞-摆振-扭转模态,耦合了旋翼动态入流、非定常动态失速气动模型和旋翼配平模型。整个算法采用标准的状态空间格式。为验证算法,以UH-60A直升机为例,在低速和高速2种定常飞行情况下,研究了偏航流效应对旋翼气功性能的影响,结果表明高速飞行时偏航流效应对旋翼气功性能影响较大。模拟了桨盘平面诱导速度分布以及桨叶叶素气动载荷,将仿真结果与飞行试验结果进行对比验证,本文建立的方法可以准确地求解定常飞行时桨叶的非定常气动载荷,捕捉其沿方位角的变化特征,对于大速度前飞仍能适用。      

对旋轴流风机非定常流场数值模拟

利用全隐格式的非稳态压力修正求解方法,分析和研究了对旋叶轮级间流场非定常流动情况,上游尾迹与下游叶片位势作用对流场的影响,并与动态实验结果进行了有益的对比.结果表明:非定常流场的计算更能反映叶轮机实际运行性能,如果能很好地组织对旋叶轮的非定常流动,还可以进一步提高其压比和效率.     

轴流压气机尾流撞击效应机理的数值模拟

通过精确给定转子进口边界条件,以模拟上游静子尾迹对转子流场的影响,这样利用单排的计算结果模拟转/静干涉引起的非定常流动对流场气动参数的影响,因为相邻叶片排的叶片数不同,从而避免传统方法中常用的叶片数约化存在的难以正确模拟非定常激励频率的问题,对计算时间和计算能力的要求也可减少.针对一个低速转子的数值模拟结果显示,如果能合理地组织来流尾迹的扫过频率和尾迹亏损的幅度,对转子气动性能的改善是很显著的,从而为风扇/压气机中控制和利用非定常流动提供了一个可行的方向.     

环量控制翼型非定常气动力建模

针对目前环量控制技术中射流参数与迎角对翼型气动特性的影响高度耦合，对应非定常气动力模型精度较差的研究现状，基于环量控制翼型强迫俯仰振动数值模拟数据，借助Kriging模型实现环量控制翼型的定常气动力插值，借助微分方程模型完成了适用于环量控制翼型的线性微分方程建模，采用两步线性回归参数辨识方法辨识线性微分方程模型中特征时间常数等参数，对高动量系数大振幅流动状态下的非线性影响进行修正。研究结果表明:基于Kriging模型实现的环量控制翼型定常气动力插值精度较传统气动导数模型高，建立的环量控制翼型非定常气动力模型能够精确预测不同流动状态下的气动力和力矩系数变化情况。      

后掠翼身干扰区流动特性及改善措施研究

利用流动显示及表面压力测量方法研究了后掠翼身干扰区的流动特性,并研究了用小边条等措施改善干扰区的流动特性的效果.结果表明,随着不同机翼后掠角、不同迎角及不同Re数对干扰区流动特性的影响,流态可以从一涡系变成多涡系,由定常变成非定常,而且在一定的Re数以后涡系会紊流化;翼身干扰区上游的的逆压梯度是导致边界层分离的物理原因,利用面积很小的边条可以降低干扰区局部的逆压梯度,可以导致干扰区的旋涡很弱,甚至不出现,这是很有实际意义的.