冲压发动机进气道压力振荡过程的数值研究

针对来流马赫数3.5,0°迎角飞行条件下冲压发动机的简化模型,数值模拟了进气道整流罩开启过程的非定常流动。计算采用AUSMPW有限体积格式,基于结构与非结构的混合网格技术,求解了非定常N-S方程。计算结果表明,在进气道整流罩打开以后,超声速气流会在燃烧室通道内形成振幅衰减较慢的低频剧烈振荡。在同样计算条件下采用Euler方程进行对比性研究,数值实验发现对这类大雷诺数、强非定常性的问题,Euler方程和N-S方程得到的结果相差不大。     

考虑螺旋桨滑流的机载探空仪分离安全性研究

为研究螺旋桨滑流效应下机载探空仪投放过程中气动干扰特性及探空仪运动规律,采用计算流体力学(CFD)方法对机载探空仪投放过程进行非定常数值模拟。通过动态嵌套网格技术模拟螺旋桨转动和探空仪运动,结合刚体六自由度运动方程求解三维非定常N-S(Navier-Stokes)方程,成功模拟了基于察/打一体无人机平台的探空仪自由投放过程,获得了包括探空仪运动姿态、运动轨迹等在内的数据信息。计算结果表明:探空仪初始弹射速度越大,其脱离螺旋桨危险区域的时间越短,对飞行越安全;螺旋桨滑流对机体方向气流有明显的加速效应,使探空仪更快脱离螺旋桨区域;三个常用的巡航飞行迎角对比发现,0°迎角下投放最优,2°迎角次之,-2°迎角不满足投放安全性要求。     

三维振荡叶栅内部非定常流动数值模拟研究

进一步发展了非定常雷诺平均N-S方程求解程序,并基于影响系数法对三种不同振型下的三维涡轮振荡叶片绕流问题进行了数值模拟研究.数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;振型对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,其中叶片运动导致叶栅各通道中流量的再分配以及出口气流角的变化等对叶片的负荷分布会产生较为显著的影响.     

跨声速开式空腔模型分离特性研究

采用非定常RANS方程求解器、非结构动网格技术及六自由度方程耦合的数值仿真技术,对跨声速条件下的开式空腔模型分离特性进行仿真研究,对比分析了高度、弹射力组合以及斜板措施对模型分离轨迹及姿态的影响。仿真结果表明,较高的飞行高度有利于改善模型的分离特性;弹射力对模型姿态影响较大,当弹射力组合提供一定的低头力矩时,模型分离下落较平稳;斜板措施可有效地改善模型下落姿态角,更利于安全分离。     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

单转子风扇的三维反问题气动设计

采用准三维流函数反问题设计方法与三维 N-S求解方法的相互迭代对单转子风扇进行气动设计。在叶片初步计算中得到叶栅进出口气流角沿径向分布 ,并将它作为本文的目标函数。采用准三维反问题求解方法 ,依次构造出各个 S1流面上的叶片几何形状和气流角分布。然后再采用 N-S方程的求解方法 ,对叶片进行全三维流场的数值计算。通过 N-S方程计算结果与目标函数的对比 ,重新修正叶片出口气流角分布 ,并作为下一次反问题设计的目标函数。经过反问题与 N-S方程求解的反复迭代 ,最终得到满足设计要求的叶型。     

进口气流角对三维振荡涡轮叶片非定常流动影响的数值模拟研究

进一步发展了非定常雷诺平均N—S方程求解程序,并基于影响系数法,在三种不同振型下就进口气流角对三维涡轮振荡叶片绕流的影响进行了数值模拟研究。数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;进口气流角对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,且在不同模态下,进口气流角对振荡叶片非定常流动的影响规律表现不同。     

卫星姿控发动机混合物羽流场分区耦合计算研究

研究求解喷管内流场N-S方程数值计算方法,发展基于N-S方程物面边界滑移流理论计算技术。提出求解羽流核心区轴对称DSMC模拟方法与远场三维DSMC仿真方案,发展多组元混合物羽流DSMC仿真方法。研究求解卫星姿控发动机内外近场、远场、倒流区和物面相互作用影响区多流域流场分区耦合计算技术,建立了一套用于求解混合物燃气羽流及对太阳电池帆板与卫星体表面撞击污染影响数值模拟方法。通过对分别安装于某在轨卫星不同位置两个典型姿控发动机燃气五组元混合物羽流计算研究及相关结果对比分析,证实本文数值方法可靠性。     

基于CFD/RBD方法的旋转弹轨迹仿真

旋转弹在飞行过程中,由于表面气动载荷分布不对称,可能会出现耦合共振、灾难偏航等现象,使导弹飞行失常偏离目标。针对此类问题,采用非定常N-S方程和刚体动力学方程耦合求解的方法,仿真旋转弹飞行轨迹。非定常气动力计算使用非结构混合网格N-S方程求解程序HUNS3D,刚体六自由度运动方程采用改进的四阶Adams预估校正法求解。将这种耦合求解方法应用于美国ARL旋转弹试验模型的轨迹仿真,并将计算得到的旋转弹姿态角、空间位置和气动力与相关文献进行对比分析,结果表明:本文方法能较精确地仿真旋转弹的飞行轨迹。     

翼型扑翼运动流场的数值模拟

基于非定常N-S方程,利用动网格方法和预处理方法对翼型在低速度下扑翼运动进行数值模拟。采用Jameson中心格式的有限体积法、双时间推进法求解非定常N-S方程;紊流模型采用基于SA紊流模型的DES模型。本文计算模拟翼型在低雷诺数下扑翼运动过程,计算所得推力系数结果与实验参考结果吻合较好,对扑翼飞行器设计有一定参考价值。     

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

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

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

运动边界非定常全速流场的数值模拟

将预处理引入到ALE格式的非定常可压缩N-S方程中,发展了一种能有效模拟含运动边界的三维非定常全速流场的数值方法.控制方程在非结构动态网格上用双时间步推进并用隐式、LU-SGS迭代求解;提出了一种特征变量的边界条件以适合预处理后的特征系统.几个典型算例的结果表明,预处理方法将可压缩流方法的计算能力拓展到了低马赫数范围,有效地提高了收敛速度和计算精度.     

大转角压气机静子叶栅附面层吹吸数值研究

针对大转角静子叶栅吸力面气流分离严重、总压损失较大的问题,根据二次流对主流进行引射的作用原理,采用附面层吹吸技术使主流自动贴近叶栅吸力面,并同时加速附面层内的低速气流。研究中使用FLUENT软件求解二维N-S方程,数值模拟采用吹吸气技术前后两种叶栅通道的流场。结果表明:采用附面层吹吸技术后叶栅的气动性能明显改善,气流转折能力增强,扩散因子增大,总压损失系数减小。     

迷宫通道内部流动和泄漏特性的数值分析

求解非定常N-S方程, 数值预测了深宽比H/W分别为0.1、0.3和0.5的单空腔直齿迷宫通道的内部流动, 计算结果展示了流动结构的振荡及出口泄漏量下降和回升的波动。3种结构的密封效果优劣与泄漏量波动幅度大小相对应, 其排列顺序是深宽比H/W分别为0.3、0.5和0.1的结构。还讨论了涡运动与涡合并、分离泡的迁移、射流的弯曲及对空腔壁面的冲击等复杂物理过程的内在联系, 分析了射流冲击所产生的多种效应对增加能耗减小泄漏的作用。      

多通道涡轮级的流场/温度场非定常数值模拟

通过求解二维 N-S方程对多通道涡轮级进口有温度畸变时的流场和温度场进行了非定常数值模拟。相对于简化成单通道涡轮级的数值模拟得到更精确的结果 ,尤其是对温度场的模拟。改变热斑数与导叶通道数比会对动叶时均温度分布产生较大影响 ,动叶压力面温度分布所受的影响更大。     

不同螺旋桨滑流数值模拟方法对比研究

采用动量激励盘理论,求解准定常RANS方程,模拟螺旋桨滑流对飞机气动特性的影响。分别基于螺旋桨+中心体、螺旋桨+短舱、螺旋桨+短舱+翼身组合体构型的非定常流场模拟得到的桨叶载荷,建立三种不同的激励盘载荷。将三种载荷的流场模拟结果与全机非定常模拟的时间平均结果进行对比,可知三种不同激励盘载荷均能准确地模拟螺旋桨滑流的影响,因此,可以通过最简单的螺旋桨+中心体构型建立激励盘模型,用于螺旋桨滑流工程预测。     

螺旋桨滑流与机翼气动干扰数值模拟研究

分别采用非定常方法和激励盘理论研究螺旋桨滑流对机翼的干扰作用,并对计算结果进行对比分析。非定常方法基于动态面搭接结构化网格,运用有限体积法求解非定常RANS方程。对于激励盘理论,首先采用点对点结构化网格,随单独螺旋桨旋转的方式进行非定常模拟得到螺旋桨桨叶压力分布,然后将桨叶压强分布转化为激励盘载荷分布,最后将载荷分布作为边界条件对干扰流场进行定常模拟。数值结果表明:除了在较大来流攻角时有所差异,两种方式得到的机翼气动力吻合良好;非定常方法可以得到螺旋桨与机翼的瞬时干扰流场,而激励盘理论则可以大大缩短数值模拟的时间并降低计算要求。     

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

预处理方法在混合笛卡尔网格中的应用研究

使用基于混合笛卡尔网格的预处理方法,对低马赫数下的定常/非定常流动问题进行了数值模拟研究。网格生成中,在物面附近生成贴体结构网格,其余计算区域使用笛卡尔网格进行填充,两套网格之间通过查找"贡献单元"实现流场数据间的传递。同时,使用基于密度的预处理方法,发展了一套可求解从低速(极低马赫数)到常规马赫数的N-S方程求解器,时间离散使用隐式双时间步LU-SGS格式,空间离散使用具有二阶精度的格心格式有限体积方法。非定常流动问题的数值模拟过程中,使用逆距离插值来进行新现网格单元上流场参数确定。通过对NACA0012翼型在低马赫数下的定常绕流和动态失速问题的数值模拟研究表明:使用预处理方法能够加快低马赫数下数值计算的收敛速度和提高计算的准确性,基于混合笛卡尔网格的N-S方程求解器能够模拟包含运动边界的低速不可压非定常流动问题。耦合预处理技术的混合笛卡尔网格方法给包含低速和常规马赫数的复杂运动边界问题的求解提供了新的思路。