基于迎风格式数值模拟三角翼涡流场

本文应用基于有限体积方法的高精度迎风格式模拟了三角翼在大迎角条件下的绕流流场。选用不同几何外形的三角翼，由计算结果可以清楚的看到前缘自由剪切层的发展，涡的形成及其复杂的相互作用。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

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

采用杂交通量分裂的ＮＮＤ格式模拟了跨，超声速三角翼背风区的旋涡运动结构。结果表明，旋涡的截面结构依赖于涡轴上的速度与密度的乘积沿轴向的导数λ＝「（１／ρ）（δρｕ１／δｌ）」，当它大于零时，涡心附近的截面流线向里转，而当它小于零时，涡心附近的截面流线向外转，λ由大于零变到小于零，涡心附近将出现一个稳定的极限环，多一次变号，将多出现一个极限环。     

三角翼大攻角绕流数值模拟研究

采用Harten-Yee的二阶精度隐式TVD格式对N-S方程进行数值模拟,研究了后掠角为65°和70°的三角翼绕流的流场特性。计算结果表明：TVD格式中限制器函数的选取对于三角翼绕流结构的计算结果具有较重要的影响。结合张涵信的旋涡定性分析理论对旋涡的发展过程进行研究,此外,还对二次涡破裂现象和网格分布影响等问题进行了探讨。     

基于迎风格式的双曲型网格生成方法

双曲型网格生成方法是一种十分适合于重叠网格方法应用的网格生成方法。为了避免双曲型方程将初始外形的不光滑性继续向外传播，本文采用了自动隐含数值耗散项的迎风格式进行数值离散。在法向推进方向，采用了隐式格式，使得网格推进步长的大小不受稳定性条件的限制。对实际应用中可能出现的边界条件进行了分类和处理。几个典型外形网格生成的质量是令人满意的，说明本文方法是成功的。     

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

跨声速下迎风格式的比较研究

研究了不同迎风格式的耗散机理,通过求解N-S方程,比较了其耗散性能.采用SST,EASM k-ω湍流模型,分析了湍流数值模拟中不同格式的计算精度,重点比较了迎风格式抑制壁面振荡的能力.选取跨声速收敛-扩张喷管内流场算例,从不同角度考察了格式的计算性能,数值计算结果和试验结果吻合,获得了不同格式耗散性能的比较结果.     

利用紧致格式捕捉间断的数值方法研究

将通量限制或加权的思想引入到紧致格式中,构造了两类传统方法与紧致格式混合组成的差分格式:通量限制与加权型差分格式。通过方波、组合波、定常激波、非定常Sod问题、Shu问题和Lax问题上的计算,以及与精确解的比较,结果表明这两种方法在间断的捕捉上具有高精度、高分辨率,而在计算方波、组合波或Shu问题和Lax问题上,加权格式具有更大的优势。     

绕三角翼旋涡流动的数值模拟

本文通过高清晰度的Ｎ－Ｓ方程数值模拟对绕大后掠三角机翼的旋涡流动结构及旋涡破裂进行了研究,给出了机翼背风面上复杂的旋涡流动结构,绘制了流动横截面上的截面流线。算例的计算采用了隐式近似分解格式。     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

三角翼动态气动特性低速实验研究

三角翼快速上仰和正弦振荡非定常气动力的实验研究是在西北工业大学Φ1．5m开口、低速风洞中进行的。60°三角翼迎角从0°到90°的快速上仰实验结果表明，随上仰速率的提高，最大升力系数和失速迎角随之增大。本文还对大振幅正弦振荡三角翼的减缩频率，平均迎角，振幅和俯仰轴位置对机翼瞬态载荷的影响进行了研究。     

上仰三角翼的流动与频率特性

在水洞和风洞中分别完成了上仰三角翼的流动显示和动态测压试验,目的是研究上仰-停止三角翼的流动和频率特性。水洞中的流动显示试验结果表明,在不同的上仰迎角范围和不同的简化频率k值下,三角翼前缘涡破裂点的位置变化是不相同的。通过分析在风洞中测得的三角翼做上仰运动时的压力信号的结果表明,在破裂涡区中,在相同的名义迎角下,无量纲螺旋波的主频比静态时的高,且k值越大无量纲螺旋波频率越高。不同后掠角的三角翼均能得到相似的结论。     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

基于DES方法的三角翼激波-涡干扰流场数值模拟

采用基于Spalart-Allmaras湍流模型的脱体涡模拟（DES）方法,数值求解Navier-Stokes方程,模拟绕尖前缘三角翼的跨音速流动,并对三角翼上翼面的复杂激波-旋涡干扰流场进行了分析。与NASA兰利研究中心的NTF风洞实验结果对比分析表明,DES方法能很好地模拟跨音速三角翼上的旋涡流动。随着攻角由中度攻角增加到大攻角,支架附近的激波越来越强,对主分离涡的干扰作用越来越大,直至出现激波干扰导致的涡破裂。激波的形状、位置及涡破裂位置均与实验结果吻合良好。     

预处理AUSMDV格式在微型扑翼绕流数值模拟中的应用

由于微型扑翼的低速低雷诺特性,对其粘性绕流的数值模拟存在较大困难。传统的中心空间离散格式不具有格式耗散而需要加入带有自由度的人工耗散,给计算带来不确定因素;此外,在低速低雷诺条件下,RANS方程组的条件数太大,导致计算收敛缓慢甚至不收敛,需要引入预处理方法加速收敛。本文使用具有低耗散特性的AUSMDV空间离散格式结合预处理方法对扑翼的非定常粘性绕流进行了数值模拟,非定常推进采用了含双时间迭代的LU-SGS隐式时间推进,湍流模型采用了BL模型。计算结果表明,加入预处理能够克服AUSMDV格式难以收敛的缺点,并且与中心格式相比,基于预处理的AUSMDV格式能够更准确地模拟扑翼的非定常粘性绕流,通过将本文计算结果与文献结果对比验证了本文方法的有效性。在此基础上研究了扑翼相关参数对气动特性的影响,对进一步了解扑翼气动力产生机理,指导扑翼的气动设计具有参考价值。     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

基于光纤传感网络的可变体机翼应变场数值模拟及实验验证

在不同飞行姿态及飞行条件下,可变体机翼智能结构及驱动器根据外部条件及控制指令,可以实现翼体自适应动态变形,表现出良好的空气动力学性能.由于翼体结构应力场与温度场的分布和作用形式复杂多样,因此需要对结构进行多物理场耦合分析,并准确获取翼体结构应变分布的有效信息.为此,采用COMSOL Multiphysic多物理场数值分...     

带配平翼钝体高超声速粘性绕流的数值模拟

采用ＮＮＤ差分格式，通过求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了带配平翼钝体的高超声速粘性绕流。文中首先以球头驻点压强和热流为参数，讨论了在差分方程左端采用ＬＵ－ＳＧＳ、ＬＵ－ＡＤＩ和对角化ＡＤＩ三种不同隐式算法收敛效率的异同。然后以带配平翼钝体的高超声速粘性绕流为模型，对所研制的程序进行了计算验证。在计算中采用了代数方法和求解椭圆型方程方法相结合的网格生成技术，针对配平翼外形给出了贴体性、