环流喷管的轴对称湍流流场计算

对火箭发动机环流喷管进行了二维数值仿真,用数值方法求解Navier－Stokes方程,数值格式为显格式时间推进格式,湍流模型选择可靠性较好的两层代数湍流模型。由于所研究的流场边界形状复杂,计算中采用了分区的方法。通过多个算例的对比,确定了使环形喷管推力效率较高的几何形状。画出了环流喷管内部的计算网格图、速度矢量图、等马赫线图、等压线图及喷管上壁面及轴线的压力分布图。     

数值模拟二维跨音速叶栅可压粘性流动

用有限体积法数值模拟了任意回转流面跨音速压气机叶栅中的粘性流动，四步龙格－库塔法用于N－S方程的时间推进。给出了一种抑制叶栅计算中“数值失速”的方法，并采用隐式残差光顺技术加速收敛。用所发展的方法计算了三种叶栅在宽广来流Ma数条件下的流动，计算结果与实验结果吻合良好，流场激波具有较高的分辨率。     

带尖针杆的钝体粘性绕流的数值模拟

用无波动、无自由参数耗散的差分（ＮＮＤ）格式数值求解了Ｎ－Ｓ方程,模拟了尖锥柱与球头组合体的高超音速粘性流动;清晰地给出了激波、膨胀波、剪切层和旋涡等流动图象。计算给出的定常流场与实验图象甚为一致。通过数值模拟研究了激波、旋涡的形成和发展过程以及相互作用情况。发现分离区中的旋涡存在分叉、合并及二次分离等演变历程。钝体主激波和分离激波从开始形成到相互撞击;发展到在与分离区相互干扰下向前、向后运动和变形;最后趋于定常状态。     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

微型飞行器的仿生力学 ——蝴蝶飞行的气动力特性

研究一种蝴蝶(Morpho peleides)前飞时的气动力特性.在运动重叠网格上数值求解Navier-Stokes方程,对蝴蝶前飞时左、右翅膀的拍动运动以及跟随身体一起的俯仰运动进行计算.结果表明:蝴蝶主要用"阻力原理"作拍动飞行,即平衡身体重量的举力和克服身体阻力的推力均主要由翅膀的阻力提供.蝴蝶翅在下拍中产生很大的瞬态阻力(平行于拍动运动的力),对流动结构分析表明,产生此力的机制如下:每次下拍中产生了一个由前缘涡, 翅端涡及起动涡构成的强"涡环",其包含一个沿拍动方向的射流,产生此射流的反作用力即翅膀的阻力.平衡身体重量的举力主要由翅膀下拍中产生的阻力提供.上拍时(由于身体上仰,上拍实际是向后和向上拍动的),翅也产生阻力,但较下拍时小的多.平衡身体阻力的推力主要由翅膀上拍中产生的阻力提供.      

基于IPDG方法的嵌套网格技术

在内罚间断伽辽金（IPDG）方法框架内应用嵌套网格技术求解了二维可压缩Navier-Stokes方程。通过把黏性通量作为辅助变量实现了Navier-Stokes方程的降阶,继而用间断有限元方法进行离散得到空间半离散方程。时间推进采用Newton-Krylov隐式方法。利用库埃特流动精确解验证了该方法的精度。在此基础上,对包括有黏NACA0012翼型绕流、定常和非定常的圆柱绕流在内的若干典型测试算例开展数值模拟,进一步验证了该方法的鲁棒性和可靠性。在嵌套网格中应用了h型网格自适应技术用于提高激波分辨率,对NACA0012翼型无黏跨声速绕流开展数值模拟,结果表明：自适应之后的单元数相比无自适应时只增加了8.4%,但是激波分辨率得到了显著提高,激波附近的计算结果也与实验值符合得更好,从而验证了该方法的有效性。     

以风轮气动性能为目标的风力机翼型优化设计

提出了一种优化风力机翼型气动特性的设计方法。在现有风力机翼型的基础上,通过改变其最大厚度位置至后缘部分的外形进行优化计算。采用N-S方程计算翼型的流场及气动特性,以提高风力机切向力同时减小其法向力为目的,用二维线性插值的方法进行多目标优化。以DU 93W 210翼型为例,针对风力机翼型气动性能本身的矛盾性,根据实际需求提出几种方案加以优化。结果表明,在不同的性能需求下,本文所提出的设计方法能取得较好的效果。     

基于混合动网格的二维非定常粘性流动数值模拟

利用基于Delaunay三角化的动网格方法,在混合网格上求解了带运动边界的二维非定常Navier-Stokes方程。在有限体积法格心格式的基础上,采用带人工耗散项的Jameson中心格式以及双时间步长推进方法对Navier-Stokes方程中的空间项和时间项分别进行了离散。湍流模型为Spalart-Allmaras方程模型。对在跨声速粘性绕流中进行俯仰振动以及后缘摆动的翼型进行了数值模拟,得到了令人满意的结果。     

基于高阶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格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

基于重叠网格技术和多重网格算法的悬停旋翼粘性绕流数值模拟

发展了一种基于重叠网格的多重网格算法,并应用于悬停旋翼粘性绕流的Navier-Stokes方程数值模拟.多重网格算法在重叠网格上实施主要面临如下两个问题:粗网格间几何关系的建立,即洞单元、边界单元的确定;洞边界存在的情况下如何进行粗细网格间流场信息的传递.本文针对以上问题给出了相应的解决方案.采用两层网格的V循环,应用有限体积法对一跨声速悬停旋翼粘性流场进行了数值计算.计算结果表明该多重网格方法可以显著加快悬停旋翼粘性流场数值解收敛速度,加速比约为3.