绕HERMES外形高超声速无粘流动的数值模拟

本文采用时间相关法和隐式推进迭代方法求解了绕ＨＥＲＭＥＳ外形高超声速无粘流动。将抛物型方程和代数方程相结合生成计算网格，既节省了网格生成的时间，又能得到较合理的网格分布。在时间相关法计算中，采用了显式ＮＮＤ格式；在隐式推进迭代计算中，采用了推广的隐式ＮＮＤ格式。     

p型多重网格间断Galekin有限元方法研究

在二维非结构网格上,使用p型多重网格间断Galerkin方法求解定常可压缩欧拉方程。p型多重网格方法主要特征是通过对不同阶次多项式的近似解进行递归迭代求解。文中高阶近似(p0)上使用显式格式,在最低阶近似(p=0)上选用隐式格式,而非显式格式,从而在保证精度和占用较小内存的情况下加速收敛到定常解。运用该方法对NACA0012跨音速无粘流动进行数值模拟,数值结果表明:p型多重网格方法同单重显式Runge-Kutta方法相比,收敛速度能够提高6倍左右,并且精度保持不变。     

高阶DG格式多重网格方法构造中的不同隐式策略影响研究

DG方法是一种非常具有潜力的高精度方法,但其在对复杂外形的数值模拟方面仍存在内存需求量大、计算量巨大等不足.为了进一步提高DG方法求解Euler方程的效率,在传统p型多重网格的基础上,结合LU-SGS和GMRES两种隐式迭代方法,研究其整体加速性能.p型多重网格方法通过对不同阶次多项式近似解进行递归迭代求解,来达到加速收敛的目的.高阶近似(p＞0)使用显式龙格库塔格式,最低阶近似(p=0)使用隐式格式.对NACA0012翼型和ONERA M6机翼跨音速无粘流动进行数值模拟,结果表明:与显式TVD-RKDG时间格式相比,DG(p0)层上采用LU-SGS和GMRES的p型多重网格方法收敛速度均得到明显提高,且GMRES迭代法性能最佳,LU-SGS迭代法次之.     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

用高效矢通量分裂格式求解二维叶栅的跨音流动

提出了一种高效隐式矢通量分裂格式。该格式通过在显格式上添加合适的高阶小量, 避免了近似因子分解和系数矩阵运算, 大大地减少了计算量。同时该格式在收敛性及激波捕捉精度等方面均保持了矢通量分裂格式所具有的优势。将该格式用到了二维叶栅的无粘跨音流场的数值模拟中去, 结果令人满意, 证实了无矩阵运算、无近似因子分解的隐格式构造有一定的推广意义。      

二维非结构网格上的高阶间断有限元方法研究

对二维非结构网格上的求解欧拉方程的高阶间断有限元方法进行了研究.运用间断有限元基本理论,采用多项式基函数对解进行近似描述,欧拉方程的数值通量项运用高斯积分公式和Roe格式求取,时间推进采用传统四步Runge-Kutta方法.运用发展的间断有限元方法对NACA0012翼型在亚跨音速情况下的无粘流动和圆柱的低速无粘流动进行了数值计算,结果表明:间断有限元方法具有良好的收敛特性、很小的数值耗散和很好的激波捕捉能力.     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

解对流扩散方程的有限元逆风因子法

 本文提出了一种解对流扩散方程的有限元逆风因子法。这种格式的数值解的稳定性很强,计算收敛速度很快,所引入的虚假扩散较少。一系列的数值试验证实了上述结论。     

隐式通量分裂线松驰迭代法求解平面叶栅无粘绕流问题

本文工作是在 Mac Cormack通量分裂格式 [2 ] 基础上发展的一种有限面积通量分裂隐式格式 ,其特点是在求解隐式离散化方程时 ,采用往返扫描一次的 Gauss- Seidel线松驰迭代方法 ,避免了对时间步长的任何限制。为提高定常解精度 ,格式的显式右端项采用二阶精度的离散。在数值求解跨音速涡轮平面叶栅问题中 ,对壁面边界作了较仔细的隐式处理。数值计算表明本方法保持了 Mac Cormack格式具有的高收敛速率 ,(约 30个时间步即可达到定常解 )而每一时间步计算量约减少一半。数值结果与实验结果符合得很好。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

高超声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析，研究了Steger-Warming通矢量分裂、Van Leer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

控制体有限元方法中的一种高阶插值格式

针对控制体有限元方法对流动方向不平行于坐标轴方向的回流问题数值模拟采用各种迎风格式都容易产生假扩散的问题,本文给出了一种由插值单元及其迎风相邻单元所确定的高阶插值格式。通过对Re=1000时三种不同形状的二维四边形空穴流模拟,表明该格式与局部斜迎风格式的适当组合可有效地减小控制体有限元方法对回流问题数值模拟的假扩散。      

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

侧向多喷流流场数值研究

为研究侧向喷流对高超声速拦截弹喷流直接力控制性能的影响，发展了在有限体积离散方法框架下的三维非定常N—S方程，采用LU—SSOR隐式时间推进格式和残值项空间迎风型格式求解技术，对大气层内高超声速拦截弹多个姿态控制喷流发动机同时工作时的干扰流场进行了数值求解。计算给出了自由流与侧喷流、侧喷流与侧喷流之间复杂干扰流场结构图(等压云图、等马赫云图)以及喷流附近的压强分布，其结果与单一侧向喷流流场情况进行了比较和分析。结果显示，喷流前的压强场基本不受后喷流的影响，处于前喷流尾流区后的喷流的干扰效率将受到较大影响。     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。     

Numerical Simulation of Rotor Flow Field Based on Overset Grids and Several Spatial and Temporal Discretization Schemes

A numerical method based on solutions of Euler/Navier-Stokes (N-S) equations is developed for calculating the flow field over a rotor in hover. Jameson central scheme, van Leer flux-vector splitting scheme, advection upwind splitting method (AUSM) scheme, upwind AUSM/van Leer scheme, AUSM+ scheme and AUSMDV scheme are implemented for spatial discretization, and van Albada limiter is also applied. For temporal discretization, both explicit Runge-Kutta method and implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method are attempted. Simultaneously, overset grid technique is adopted. In detail, hole-map method is utilized to identify intergrid boundary points (IGBPs). Furthermore, aimed at identification issue of donor elements, inverse-map method is implemented. Eventually, blade surface pressure distributions derived from numerical simulation are validated compared with experimental data, showing that all the schemes mentioned above have the capability to predict the rotor flow field accurately. At the same time, vorticity contours are illustrated for analysis, and other characteristics are also analyzed.