带副翼三维机翼绕流的Euler方程解

提出了一种求解带副翼偏转三维机翼绕的Ｅｕｌｅｒ方程计算方法。采用对接分区网格民分区求解算法的结合,有铲地求解绕引外形的复杂流动,提出了一种满足通量守恒的内边界耦吉选用ＶａｎＬｅｅｒ分裂格式离散玩粘通量项,并构造了一种Ｌｉｍｉｔｅｒ函数以保证ＴＶＤ性质。数值算例表明本方法是求解带操纵纵面偏转的机翼绕流的有效方法。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

大钝头短体跨音速大迎角不同流型的计算比较

 采用不同流型 (无粘,层流,湍流 )的数值模拟方法计算了绕大钝头短体外形的跨音速大迎角时的流态,并进行了比较分析。数值方法中选用 van Leer分裂格式离散无粘通量项,并构造了一种 Limiter函数以保证 TVD性质;采用中心差分格式来离散粘性通量项。湍流计算中采用了 B-L湍流模型。算例表明不同流型的流态共同特点是三维分离的非定常性,但不同流型计算所得流态的差异较大。     

微型飞行器小展弦比机翼的低雷诺数气动力特性分析

以固定机翼式微型飞行器的气动设计为研究背景,对小展弦比薄机翼的低雷诺数气动力特性进行了数值分析.不可压粘性绕流的求解采用了人工压缩性方法,其中对流项的离散应用了三阶迎风格式.湍流模型为Baldwin-Barth一方程模型.一无弯度薄翼的数值计算结果表明本文的数值方法是成功的.通过数值分析,揭示了弯度和展弦比对气动力特性的影响,也给出了不同弯度机翼流动结构的差别.     

一种边界非协调数值方法

研究了一种边界非协调的数值方法。该方法的提出是为了能够方便的求解复杂的运动边界绕流问题。空间上采用Galerkin有限体积空间离散格式。数值离散中,网格点分为内部计算点、IB点和外点三种类型。其中临近物面的IB点在嵌入物面边界条件时起到了作用,其流动变量值则是通过流场内解的线性插值获得。时间推进上采用结合人工压缩项的双时间步推进格式。静止和振荡圆柱不可压绕流的数值实验初步验证了本文方法的可靠性和准确性。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

Jameson/TVD格式应用旋翼翼型绕流分析

以一种常见的旋翼翼型为例 ,探讨了两种旋翼桨叶翼型绕流数值模拟方法 ,即 Jameson有限体积中心格式和有限体积通量修正 Jameson/ TVD混合格式。这两种格式的方法都采用有限体积法进行离散 ,并用龙格—库塔时间推进法求解。通过算例分析 ,比较了这两种方法的特点 ,为进一步有效地依靠求解湍流 Navier-Stokes方程来模拟完整的旋翼流场、旋翼阻力与功率奠定了基础     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

不可压缩流动的高精度非标准格子玻尔兹曼方法

首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法（FRLBM）。采用两种求解方法：一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程（直接法）;另一种是先执行碰撞步,再求解纯对流方程（分步法）。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。     

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

基于WENO-分段线性格式的直升机流场数值模拟

发展了一套适合于格心格式求解器的基于加权本质无振荡(WENO)-分段线性格式的旋翼/机身气动干扰高精度CFD计算方法。应用该方法对多种旋翼/机身组合模型前飞算例进行了数值模拟,计算得到的机身表面压力系数分布与实验结果吻合良好,说明了该方法对旋翼/机身气动干扰研究的有效性。之后进一步将该方法应用到悬停状态的X3构型复合式高速直升机旋翼/机翼/螺旋桨组合模型的复杂流场模拟中,并与悬停状态孤立旋翼和旋翼/机翼组合模型的流场进行了对比。研究发现:机翼对旋翼下洗流起阻滞作用,引起机翼下方不规则流动,且螺旋桨滑流与旋翼下洗流会相互干扰产生一定的偏折,旋翼下洗流速度更大,螺旋桨滑流会产生明显的向下偏折。      

表面涂漆对火箭尾翼热结构的影响

建立了运载火箭的气动加热工程计算方法,计算区域横跨连续流、稀薄过渡流和自由分子流,包括层流、转捩和湍流等各流态。重点研究了尾翼与芯级间的干扰加热,数值求解了自然正交曲线坐标系下尾翼内部二维单层和双层材料热传导方程,给出了典型部位的内外表面温度计算结果,表明表面涂漆对尾翼前缘并不能起到防热隔热效果。     

超声三速角翼背风区旋涡运动的数值模拟

本文采用杂交通量分裂的ＮＮＤ格式模拟了Ｍ∞＝１．９５、Ｒｅ＝９．５×１０５，α＝１０°、２０°的三角翼绕流流场，结果揭示了超声速旋涡沿其自身轴线的发展规律，当旋涡轴向速度为超声速且处于顺压区时，涡轴附近的横截面流线向外转，而由机翼前缘尖点处发出的截面流线向内卷，它们之间存在极限环。数值结果与张涵信的拓扑分析结果完全一致。     

三维自适应非结构网格的Euler方程解

 将 Ausm+ 迎风格式应用于三维非结构网格中求解 Euler方程。对单元变量进行重构以获得空间高阶精度,对时间域采用多步龙格库塔法推进,并采用了当地时间步长和隐式残差光顺技术来加速收敛。采用多点择优推进阵面法生成复杂曲面的三角形网格,利用推进阵面法生成四面体网格。采用网格自适应技术对网格进行局部加密,以减少总体网格数目,从而提高计算效率。最后给出了绕 ONERA M6机翼的跨音速流动及绕麻雀 导弹的超音速流动算例,结果表明了本方法的有效性。     

磁场对高超声速弱电离气体流动的影响

对偶极子磁场作用下的三维钝头体高超声速黏性绕流的化学非平衡流动进行了数值模拟。电导率利用组分公式计算,化学模型为7组元、6反应模型。应用诱导磁场方法将磁场与流场耦合,控制方程的空间离散采用有限体积法,扩散项用中心差分格式,对流项采用二阶迎风格式,时间推进采用3步Runge-Kutta法。计算结果表明,外加偶极子磁场使激波脱体距离增加,壁面摩擦系数和表面热流密度增大。与无磁场作用时相比,在0.153T外加磁场作用下,冻结流中的激波脱体距离增加约3倍,局部壁面摩擦系数最大增加63%,局部表面热流密度最大增加61%;而化学非平衡流中的激波脱体距离增加约0.5倍,局部壁面摩擦系数最大增加47%,局部表面热流密度最大增加31%,且化学非平衡流中激波层内温度的最大值约为冻结流中的64%。     

三角翼无尾布局全动翼尖的操纵性能研究

基于变前掠翼(VFSW)布局,采用Navier-Stokes控制方程的有限体积法离散格式,选取剪切应力输运(SST)湍流模型,对VFSW中三角翼飞行器全动翼尖(AMT)的流场进行数值分析。首先,通过未带机翼前缘延伸的三角翼试验模型验证了数值模拟算法的精度;其次,研究了三角翼无尾布局在超声速时AMT的操纵性能;最后,采用可视化方法分析了AMT的流场和作用机理。AMT计算结果表明:迎角对AMT偏航特性影响轻微,超声速时最大设计舵偏量的偏航力矩系数约为0.02,但偏航力矩和滚转力矩具有耦合性;耦合滚转力矩在局部大迎角时易反向,而舵面失升是滚转反向的根本原因;AMT的偏航作用线性较好,作动效率较高,消除不利滚转后是变前掠翼布局一种极具潜力的航向操纵面。     

大迎角振动三角翼跨声速粘性绕流数值分析

本文通过数值求解三维可压缩雷诺平均ＮＳ方程，研究发析了振动三角翼的三维复杂非定常流场，其中对流场项计算格式的空间离散在平行物物理面方向与物面法线方向的空间导数分别采用二阶ＮＮＤ格式和由ＮＮＤ格式高阶插值而得到的三阶迎风格式，时间推进为隐式ＬＵ分解方法，并使用了修正的代数ＢＬ湍流模型模拟高雷诺数的湍流流动。