基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

基于混合网格的直升机旋翼/机身非定常干扰流场数值模拟方法

基于混合网格建立了一个适合于直升机旋翼/机身非定常干扰流场的数值模拟方法。其中,采用混合网格系统生成旋翼/机身干扰流场网格以降低网格生成难度;将二阶MUSCL格式和Roe格式相结合进行空间离散以提高流场的计算精度;同时引入隐式LU-SGS时间格式和OpenMP并行模式以提高干扰流场的计算效率。应用建立的计算方法,对Georgia和Robin旋翼/机身模型进行了数值模拟和分析,与试验结果良好的相关性表明了该方法的有效性。     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

基于嵌套网格技术和多种空间离散格式的旋翼流场数值模拟

给出了三维Euler/N-S方程数值模拟悬停状态旋翼流场的方法和模型.在空间离散方法上分别采用Jameson中心差分格式、Van Leer矢通量分裂格式、AUSM格式3种方法,同时加入了Van Albada限制器,并应用了嵌套网格方法.文中采用Hole-Map方法来确定洞边界,采用Inverse-Map方法来搜寻贡献单元.最后给出了旋翼桨叶表面压力分布的计算结果,并与实验数据进行了对比,证明以上3种空间离散格式在旋翼流场计算中的准确性.     

超声速化学反应流场迎风格式数值模拟及并行计算

基于有限体积迎风格式对超声速燃烧流场进行了的数值模拟.由于超声速燃烧流场绕流的复杂性,要求对多组分Euler/N-S方程求解的数值模拟方法应具有较高的计算精度及效率.本文引用辅助点方法建立了具有空间二阶精度的van Leer迎风矢通量分裂格式,并应用于超声速燃烧流场绕流的数值模拟.化学反应为氢气/空气十反应模型,采用考虑了化学反应特征时间的当地时间步长显式Runge-Kutta时间推进格式.对钝头体模型爆轰现象、后向台阶氢气喷射及二维内外流超声速燃烧流场模型进行了区域分裂技术的并行计算.计算结果与参考文献作了对比,得到了满意的结果.     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

弹射座椅系统外高速流场的数值模拟

对绕弹射座椅系统的高速流场进行了数值模拟.控制方程为雷诺平均的Navier-Stokes方程,空间离散为Osher逆风格式和五步龙格库塔格式,湍流模拟采用了基于M-SST两方程紊流模型的DES方法,计算中采用了区域分裂的并行技术.对来流马赫数分别为0.6,1.2的流场进行了计算,计算结果表明该方法可以较好地模拟弹射座椅系统外的高速流场.     

矩形喷管内外喷流流场的数值模拟

基于 Ｆｒａｖｅ 平均的 Ｎ Ｓ方程和 Ｂ／ Ｌ 湍流模型，采用了 Ｊａｍ ｅｓｏｎ 格式和 Ｍ Ｕ Ｓ Ｃ Ｌ 格式，对外流在亚声、超声和零 Ｍ ａｃｈ 数三种状态时矩形喷管的喷流流场进行了数值模拟，探讨了内外喷流及喷管内的流动特征，并与实验数据进行了比较，计算结果与实验结果吻合较好，验证了计算所使用的 Ｎ Ａ Ｐ Ａ 软件的可靠性，从而为用计算流体力学的方法模拟具有两股来流的三维湍流流场提供了一种有效的手段。     

地效翼风洞试验支架干扰数值分析

对地效翼移动地面风洞试验研究中的支架干扰进行了数值模拟和分析。分别对独立地效翼,带支架的地效翼及独立支架进行了数值模拟,计算采用可实现的κ-ε模型,通过求解定常不可压N-S方程,得出地效翼及支架周围流场分布情况。对几组计算结果比较分析了支架和地效翼的空气动力及由于干扰引起的空气动力,发现支架与地效翼之间的相互干扰随着地效翼迎角的增大而增强,如果忽略流动干扰造成的空气动力变化,地效翼升力误差很小,但阻力误差相对较大。同时对有干扰下和没有干扰下的流场进行了对比,分析了支架对翼尖涡流动及绕机翼流动的干扰。翼尖涡在地效翼翼尖附近的发展在0.5犮范围内基本不受支架的干扰;除支架对流场产生干扰外,移动带区域以外的固定地面附近粘性流动也对绕地效翼流动有一定的影响。本研究分析了风洞试验结果的可靠性,为地效翼风洞试验优化设计和地面效应风洞试验研究提供了参考。     

旋转直管流场的数值模拟

根据旋转坐标系下的 Favre平均 NS方程以及旋转条件下的流动特征 ,修改了 B/ L湍流模型 ,并提出了一种旋转条件下出口边界压力分布的给定方法 ,对已有的 NAPA通用湍流计算软件进行了改进和发展 ,并计算了旋转直管中的流动。通过与实验数据的比较 ,文中的计算结果与实验结果吻合较好 ,正确反映了流动特征和物理现象 ,验证了修改后的 NAPA软件的可靠性。     

壁压法在NH-2风洞中的应用

本文介绍了一种由Hackett在1981年提出的利用测量风洞壁面静压来进行大迎角、大堵塞洞壁干扰修正的方法,并用该方法编写了Fortran程序。对四个几何相似不同尺寸的模型,在NH-2风洞中进行了试验,试验和计算的结果表明,堵塞度从1.56％到16.7％的四个模型,其升力系数和阻力系数经洞壁干扰修正后非常接近。但必须指出,对于带平衡缝的闭口回流风洞,从壁面测量所得到的静压必须加以修正,才能成功地应用这种方法。     

低速开缝壁风洞洞壁干扰若干性能探讨

为适应大攻角、高升力等试验的需要,许多国家正在研究或开始使用低速开缝壁风洞。为论证在我所建造低速开缝壁风洞的可行性,用涡格法对三种机翼模型进行了计算,并用所得结论指导了试验,其结果与实验结果比较符合。经分析得出:低速开缝壁风洞升力洞壁干扰与开闭比之间的关系接近对数变化规律。最佳开闭比随缝的加深而增大,每一座低速开缝壁风洞都有一个工作范围,在此范围内洞壁干扰量及下洗量均比相应的低速闭口和开口风洞小得多,因此证明了低速开缝壁风洞的性能确实优于低速闭口和开口风洞。又由于它制造和使用不太复杂,因而正以第二代低速风洞的姿态出现在人们面前。     

侧滑机翼洞壁干扰修正量的计算方法

本文采用数值计算方法,计算了实验段截面为扁八角形的低、亚音速风洞中的侧滑薄机翼的气动特性及洞壁干扰修正量。 基于流谱观察实验,对带有后退侧缘的翼段用两种极端的气动模型来描述。计算结果表明,就自由流态的气动特性而言,两者的平均值与文[4]中的实验结果相符,而相应的洞壁干扰修正量,也与本文的实验结果相近。     

飞翼布局低速支撑干扰及局部外形畸变影响数值模拟

飞翼布局具有气动效率高、隐身性能好的优点，是未来军民用飞机的重要发展方向。该类布局模型风洞试验尾撑、腹撑及背撑等支撑形式的干扰量及由此带来的模型局部外形畸变影响较为复杂，目前还没有通用的支撑方案和试验修正方法。采用CFD数值模拟方法，分别对某小展弦比飞翼布局标模低速尾撑支杆干扰和尾部外形畸变影响进行了研究，结果表明:在常用角度范围内，所使用的数值模拟方法是可靠的，可用于风洞试验支撑方案的评估及支撑干扰的修正；对纵向特性，尾撑支杆干扰量和尾部外形畸变影响量相对全量较小；对横航向特性，尾撑支杆干扰量基本可忽略，尾部外形畸变影响量与全量相当。     

B/L湍流模型在强压力梯度流场计算中的应用

对在湍流计算中已得到广泛应用的Ｂａｌｄｗｉｎ／Ｌｏｍａｘ（Ｂ／Ｌ）代数湍流模型进行了修正。修正后的Ｂ／Ｌ模型考虑了压力梯度及物面引射速度对内层ＶａｎＤｒｉｅｓｔ衰减因子的影响。与此同时，文中还用修正的Ｂ／Ｌ湍流模型对亚声平板附面层流动、拉伐尔喷管内的正激波／湍流附面层干扰流动进行了计算。计算与实验结果比较表明，修正后的Ｂ／Ｌ模型不仅可以更准确地计算具有强压力梯度的激波／附面层干扰区内流动压力分布，而且还可以准确地预报激波／附面层干扰引起的流动分离，从而为Ｂ／Ｌ模型的工程应用开拓了领域。     

二维翼型实验的侧壁影响研究

本采用油流显示技术，测量模型弦向和展向压力分布的方法，对二维翼型风洞中的侧壁附面层影响进行了研究。为能更详细了解侧壁对模型的影响，采用不同弦长的模型和两个不同风洞的实验结果比较的方法，由油流显示揭示侧壁干扰流动图画，利用侧壁抽气和干扰修正方法，来消除侧壁对模型的影响，实验结果表明，选取合理的抽气方式可得到有效的二维实验结果。     

大飞机布局模型跨声速风洞实验尾支撑干扰研究

对某大飞机布局风洞实验尾支撑干扰开展了数值模拟和实验研究,发展的数值方法计算结果与风洞实验结果有很好的一致性.对于类似构型的飞机,在迎角-2°~6°范围,可认为尾支撑干扰量随迎角呈线性变化,采用前位叶片支撑作为辅助支撑带来的二次干扰量可以忽略,新型双天平辅助支撑系统试验进一步验证了这一结果;尾支撑对机身、尾翼、机翼等部件的绕流都有影响,干扰量随构型而变,对阻力、力矩影响较大,且随Ma数变化,因此不同构型实验数据需要单独修正.所发展的带风洞支撵系统的数值模拟软件能够满足工程应用要求,可用于支撑干扰修正研究以及风洞实验支撑系统优化设计.