限制器对高超声速体表面热流数值模拟的影响

采用Roe的FDS（Flux Difference Splitting）差分格式结合min mod,van Leer和Osher-C限制器,对高超声速体外部流场进行数值模拟.研究了限制器对高超声速体气动力/热数值模拟的区别,对比分析了3种限制器对三维模型表面热流数值模拟结果的影响.研究表明,限制器对数值模拟精度有重要影响,且对气动热数值模拟的影响大于对气动力的影响;对简单或较复杂的气动外形,采用Roe的FDS格式结合min mod限制器获得的表面热特性比它结合van Leer和Osher-C限制器时要好,可为热防护设计提供可靠结果.      

适用于Mie-Grüneisen状态方程的一种通用激波限制器

为了求解含Mie-Grüneisen状态方程多介质可压缩流动问题,针对界面追踪方法,发展了一种通用的激波限制器。根据平均守恒误差及L1误差随参数的变化情况,确定了参数的范围。由于表达式不含有状态方程中的物性参数,所以它适用于任何一种形式的状态方程。通过具体算例分析了激波限制器的开关闭状态。数值结果表明文中构造的激波限制器可以用于计算含强激波的一般状态方程的多介质可压缩流动问题。     

通量限制器与LED性质、空间离散精度

本文从标量形式的ＵＳＬＩＰ格式出发,从理论上分析了通量限制器与ＬＥＤ性质和空间离散精度的关系。利用Ｒｏｅ平均矩阵的Ｕ特性和广义数值Ｒｉｅｍａｎｎ不变量的概念,将此类格式应用于二维翼型跨声流场的数值模拟,针对提供的来流条件,考查了不同的通量限制器对流场计算精度的影响。     

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin(DG)有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用,其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题,均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题,发展了r型网格自适应方法,实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值,并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值,实现网格沿激波方向各向异性自适应加密,并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明：基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波,提高模拟精度,具有良好的收敛性和鲁棒性。     

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。      

一种新的间断侦测器及其在DGM中的应用

根据单元交界面左右变量的差别,提出了一种新的间断侦测器构造方法。该间断侦测器的构造原理简单,编程实现容易。针对一维和二维Euler方程,我们将此间断侦测器用于间断Galerkin格式的数值计算中。数值实验表明本文构造的间断侦测器能够准确捕捉到激波的位置,从而只在间断区域引入限制器,在减少计算量、保证光滑区计算精度的同时,提高了激波等强间断的分辨率,能够明显地改善流动精细结构的模拟精度。     

间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明：如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。     

一种基于间断Galerkin格式的新型限制器设计

使用高阶间断Galerkin格式求解守恒律方程组时,激波附近的Gibbs效应容易导致非物理解的产生。为抑制这一现象,必须构造合理的限制器对数值解进行处理。目前间断Galerkin格式中的限制器多源于有限体积法,在非结构网格上只对低阶导数项进行限制,对高阶导数项则很难给出普适判据。文章对间断Galerkin解进行广义Fourier展开,实现不同频域范围内的谱分解;在新的模板坐标系下描述各阶方向导数的变化规律;结合当前单元和相邻单元的信息,分层限制各阶方向导数,实现对非物理解的抑制。通过求解Euler方程,对二维Riemann问题、翼型周围的亚、跨声速流动问题、前台阶问题以及超燃冲压发动机内流场激波反射问题进行数值模拟,检验了新型限制器的可靠性以及向高阶格式推广的可行性。