考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程，选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果，采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据，对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

悬停状态下倾转旋翼机向下载荷被动减缓措施研究

基于SA一方程紊流模型的DES方法,发展了一种翼型在-90°来流迎角下非定常大分离流动的数值模拟方法,系统地研究了后缘襟翼、克鲁格襟翼、前缘下垂以及扰流板四种被动控制技术在悬停状态下倾转旋翼机向下载荷减缓中的应用,并分析了各影响参数对阻力系数的影响。计算结果表明:当后缘襟翼偏角60°时,垂直阻力最小;克鲁格襟翼偏角85°时,阻力系数最小;前缘下垂偏角45°时,阻力系数最小;当扰流板高度h/c=0.05时,阻力系数最小。最优组合构型为后缘襟翼偏角60°,克鲁格襟翼偏角85°,其向下载荷相对于原始翼型减小了55.4%。     

阵风响应及减缓的非定常数值模拟

基于网格速度和动网格方法,通过求解非定常N-S方程计算翼型副翼偏转时的阵风减缓响应,研究开环控制律参数和副翼缝隙对阵风减缓效果的影响。对于单个1-cos阵风,随着偏转角度的增大,减缓效果变大;由于粘性的影响,CFD计算得到的升力峰值减小量明显呈非线性的趋势,减缓效果低于准定常线性理论。对于连续sin阵风,通过调整副翼偏转的频率和最大偏转角,升力能达到很好的减缓效果,此时,副翼偏转频率应接近阵风频率,最大偏转角存在一个最优值。在阵风环境下,副翼间隙对升力系数影响不大,而力矩系数增大一倍,间隙会降低副翼偏转对升力的减缓效果。     

基于多点约束的大展弦比机翼静气动弹性计算

基于雷诺平均N-S方程和多块结构网格技术计算大展弦比机翼的气动力.机翼采用梁模型结构,利用有限元方法计算结构变形.建立基于多点约束(MPC,Multi-Point Constrain)的气动、结构数据的双向传递方法.采用MPC数据插值方法快捷地实现自主研发的计算流体力学(CFD,Computational Fluid Dynamics)程序和NASTRAN计算结构力学(CSD,Computational Structural Dynamics)软件的耦合计算,发展出可以考虑结构非线性的静气动弹性CFD/CSD耦合计算方法,该方法比只考虑单方向变形的柔度法更精确.开展了某大展弦比机翼的静气动弹性数值计算分析,结果显示只考虑机翼纵向变形而忽略展向变形对计算结果会有一定影响,而该算例几何非线性影响较小.     

基于线性模型的风电场发电量计算与分析

针对相对平坦地形,采用线性化流体方程描述流场,并根据大气边界层特性进一步简化计算方法,即利用所谓地形影响系数快速求解感兴趣点(一般为风机位置)的流动速度;针对尾流影响,采用工程上广泛使用的修正Park模型来求解,尾流线性叠加并忽略较弱的尾流影响。文中介绍了一个已建典型海滨风电场的装机规模、地表粗糙度以及测风情况,并以此风电场为算例采用不同测风数据进行发电量计算,预测结果与实际发电量较好地吻合。最后探索了算例风电场各风机发电量不均衡原因,提出地形图大小对风场模拟准确性的影响,并从工程设计角度提出具体意见。     

分区对接网格下多段翼型N-S方程计算

对多段翼型的复杂粘性流场，各翼段物面及尾迹区生成C型网格．外场采用O型网格，中间用H型附加区过渡．自动生成高效的分区对接C-O-H型网格。采用时间精确有限体积高精度格式解可压缩N—S方程。层流、紊流的N—S方程求解结果，说明了本方法的有效性。