基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

用规范型直接法研究立方非线性机翼的颤振

研究了在不可压缩流中具有立方非线性俯仰刚度的二元机翼的颤振问题。首先，对系统平衡点进行了特征值分析，得到在Hopf分岔点附近的颤振参数方程。然后，应用规范形直接法的Maple推导程序，计算得到Hopf分岔的规范形．并数值模拟验证了超、亚临界类型分岔情形下响应对初值的依赖性。最后．分析了线性和非线性刚度系数对颤振响应拓扑结构的影响。     

一种前后缘连续变弯度翼型气动特性分析

提出了一种基于厚度不变的翼型前后缘连续偏转变形规律,并以NACA0015翼型为例,实现了翼型变弯度的参数化。以柔性伸缩蒙皮支撑结构和机械结构实现了可连续光滑偏转后缘的变弯度翼型,验证了变形规律。以前后缘偏转角度为参数,计算并分析了各个变形状态下翼型扰流流场和气动特性,讨论了变形参数对气动特性的影响,研究了气动特性的产生机理。研究结果表明:在大迎角下,前缘偏转角对翼型失速有一定的抑制作用;在中小迎角范围内,翼型升阻比随着后缘偏转角的增大而增大,且不论迎角如何变化,总可以通过前后缘偏转来获得较高的升阻比。     

基于CFD研究等离子体环量控制流场特性

为了研究等离子体气动激励射流在翼型钝尾缘产生Coanda效应时对流场特性的影响,基于雷诺平均N-S方程中添加体积力源项的唯象学仿真方法,建立并优化了等离子体环量控制的数值计算模型,研究了基于分离点的不同激励器位置对流场的影响规律。结果表明,优化激励器位置后环量控制的效费比ΔCL／Cμ最高可以达到113．49。分析发现,当等离子体激励产生的激励射流刚好可以将分离点推移至压力面时,机翼有效弯度增加最为明显,激励效果最好。     

耦合梯度与分级Kriging模型的高效气动优化方法

Kriging代理模型中引入梯度信息能够提高模型的预测精度,但常规耦合梯度的方法都有不足之处。本文结合分级Kriging模型,提出了一种变可信度分级Kriging模型耦合梯度（GEHK）信息的新方法。首先利用梯度信息,选取扰动步长得到初始样本点附近的派生点,以派生点拟合出能够预测目标函数趋势的低精度Kriging模型。然后以初始样本点修正该模型得到高精度的Kriging模型。翼型减阻优化设计算例表明,与常规耦合梯度的Kriging模型相比,基于分级Kriging的梯度耦合方法对于扰动步长引入的误差不敏感,明显提高了模型预测精度,优化效率因此提升并使得目标函数值下降得更加迅速。相比Euler解作为低精度数据的常规分级Kriging模型,由梯度信息得到的派生点为模型提供了更准确的全局趋势预测,取得了更好的优化结果。本文方法成功应用于翼型多点减阻优化问题,说明该方法对复杂设计问题有很好的适应性。基于分级Kriging模型的耦合梯度方法克服了常规方法的缺点,提高了模型全局拟合精度,是一种优化效率更高的Kriging方法。     

旋翼非定常动态失速的三维效应分析

采用建立的高精度计算流体动力学(CFD)方法,针对旋翼非定常动态失速的三维（3D）效应特性进行研究。以Helishape 7AD旋翼为基准,开展三维效应对旋翼非定常气动特性的影响分析。研究了来流马赫数对旋翼翼型动态失速特性影响。在此基础上,针对三维情形旋翼动态失速非定常涡流动特性及诱导分离特征开展了数值分析,通过与二维情形对比表明:受旋翼旋转、轴向诱导速度等三维效应影响,旋翼桨叶剖面动态失速涡的产生、对流和脱落明显滞后于二维翼型情形,并且涡强度也更弱。越靠近桨叶内段,桨叶剖面非定常动态失速特性与二维旋翼翼型情形的差距越为明显。      

低雷诺数下协同射流关键参数对翼型气动性能的影响

协同射流是一种高效的新型主动流动控制技术,至今缺乏关键参数对翼型气动性能影响规律的系统研究。通过引入螺旋桨激励盘模型,发展了一种新的协同射流翼型流动模拟方法,使得射流反作用力计算更符合实际。在低雷诺数条件下,以NACA6415为基准翼型开展了射流动量系数、开口尺寸和位置等关键参数对翼型气动性能的影响规律研究,并探讨了相应的物理机制。结果表明：大迎角分离流状态下,射流动量系数对翼型气动性能的影响规律比小迎角附着流状态更复杂;随着吹气口尺寸增加,气泵功率系数先减后增,有效升阻比先增后减;随着吸气口尺寸增加,气泵功率系数逐渐减小,有效升阻比先增加后趋于平稳;吹/吸气口位置对翼型气动性能和气泵功率系数的影响很小。     

基于高阶谐波平衡的跨声速颤振高效预测方法

跨声速流场激波及其诱导的附面层分离等非线性因素导致跨声速颤振边界很难被准确预测,尤其是目前工程常用的偶极子格网法,在跨声速时该方法的预测精度大幅下降。在雷诺平均Navier-Stokes方程流场求解器的框架内,利用结构模态建立广义结构运动方程,利用径向基函数建立模态振型的插值方法,结合径向基函数和无限插值两种网格变形方法的优点实现高效高鲁棒性网格变形方法,从而实现颤振时间推进分析流程,利用国际颤振标模AGARD445.6机翼验证程序在跨声速颤振边界预测中的可靠性。然而,时域方法在气动/结构反复迭代,需要耗费大量的计算资源和时间。为了提高颤振预测效率,基于高阶谐波平衡（HOHB）方法快速获得广义力影响系数矩阵,利用该矩阵建立频域模态位移和气动力之间的关系,实现高效颤振频域分析方法。通过二维翼型和三维机翼算例进行验证,结果表明：在不对计算精度产生明显影响的前提下,HOHB方法能够提高颤振预测效率约6倍。     

新的解跨音速Euler方程的隐式杂交方法

 本文提出一个新的杂交格式,它由Jameson-Turkel LU分解法和Harten TVD（total variation diminishing）格式的修正数值通量耦合而成。LU分解用于Euler方程隐式部分的离散化而Harten的数值通量控制所谓的稳态解残差。精心构造LU算子,可使矩阵算子L和U很容易逐点推进求逆;以Harten的数值通量计算杂交格式右边的通量,这样做在a=1/2的条件下,可保证收敛的稳态解具有二阶精度并能高分辨率地捕获激波。数值实验表明,本文格式对于稳态计算是有生命力的并能保证较高的激波分辨率。     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe