满足几何守恒律的WENO格式及其应用

对几何守恒律的来源进行了分析,发展了一种满足几何守恒律的WENO格式,并应用于翼型层流分离现象的数值模拟中。为消除网格质量影响,采用守恒型方法计算网格导数,并将标准的WENO格式分解为中心差分部分和数值耗散部分。算例计算结果表明,几何守恒律对高精度有限差分WENO格式计算结果具有重要影响,本文方法能够消除网格导数计算误差,保证来流保持性。将本文方法应用于SD7003翼型层流分离现象的数值模拟中,计算结果与文献中计算及试验数据吻合较好,同时能够精细捕捉小尺度流场结构,准确模拟翼型层流分离现象中的复杂流动过程。     

求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。     

双同守恒律方程的加权本质无振荡格式新进展

近几年,在计算流体力学中,高精度、高分辨率的加权本质无振荡(weighted essentially non-oscillatory , WENO)格式得到很大的发展.WENO格式的主要思想是通过低阶的数值流通量的凸组合重构得到高阶的逼近,并且在间断附近具有本质无振荡的性质.本文综合介绍了双曲守恒律方程的有限差分和有限体积迎风型WENO,中心WENO,紧致中心WENO以及优化的WENO格式等,讨论了负权的处理和多维问题的解决方法.最后,通过一些算例证明WENO格式的高精度,本质无振荡的性质.图6参40      

基于简单WENO-RKDG的溃坝问题数值模拟

对溃坝问题水流间断面的高精度、高分辨率数值模拟是水动力学的重要内容。简单加权本质无振荡(WENO)限制器由"问题单元"及其相邻单元的解重构"问题单元"的解,从而抑制数值解的非物理振荡,能够很好地模拟间断问题。本文详细介绍了简单WENO限制器的基本原理和过程。将简单WENO限制器-Runge-Kutta间断Galerkin方法应用于二维浅水控制方程的求解中,对二维矩形明渠中大坝瞬间全溃、局部溃塌所致的水流运动进行了数值模拟,并将数值计算结果与理论分析进行了比较。计算结果表明,方法能够清晰地捕捉到溃坝全过程中的间断,没有非物理的振荡现象发生,简单WENO限制器-RKDG方法能够很好地模拟溃坝波的演进过程。     

改进WENO格式及其在爆轰计算中的应用

对已有的一种改进型WENO-M格式进行了理论上的修正,得到了新的WENO权函数。与原WENO-M格式相比,新得到的WENO-M-P格式减少了约9%的CPU计算时间,同时也保证了在一阶极值点处不降低精度（经典的WENO格式在一阶极值点处精度下降）,仍然保持5阶精度。为了验证修正后的格式,采用二步反应模型数值模拟了几组一维和二维爆轰：在一维爆轰模拟中,对比了WENO-M格式和WENO-M-P格式在一定温度下的最低起爆压强以及一定压强下的最低起爆温度;在二维旋转爆轰的模拟中,对比了WENO-M格式和WENO-M-P格式的计算效果。结果表明：在同样的起爆温度下,WENO-M 格式和 WENO-M-P 格式的最低起爆压强均高于 WENO格式,但WENO-M-P格式比WENO-M格式的最低起爆压强低;三种格式在不同压强下的最低起爆温度相同;WENO-M-P格式具有和WENO-M格式一样的计算效果,优势是节省了计算时间。     

无粘可压缩流动的改进高精度方法

为更准确捕捉复杂流场的流动细节,通过对WENO格式的光滑因子进行改进,发展了一种新的五阶WENO格式。对三阶ENO格式进行加权可以得到五阶WENO格式,但是不同的加权处理,WENO格式在极值处保持加权基本无振荡的效果不同,本文构造了二阶精度的局部光滑因子,及不含一阶二阶导数的高阶全局光滑因子,从而实现WENO格式在极值处有五阶精度。基于改进五阶WENO格式,对一维对流方程、一维和二维可压缩无粘问题进行算例验证,并与传统WENO-JS格式和WENO-Z格式进行比较。计算结果表明,改进五阶WENO格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

应用于计算气动声学的优化有限紧致格式

对于含间断的计算气动声学问题,数值计算的格式不仅要求低耗散低色散的设计,对短波具有较高的分辨率,还要求能捕捉激波.中心紧致格式具有高精度,具有无耗散和低色散特征,但不能捕捉间断和激波;WENO格式处理间断较为成功,而耗散和色散误差相对较大.有限紧致格式可以将紧致格式与WENO格式相结合构造成混合格式,利用光滑因子之间的关系对激波区域进行自动判断,将传统的全域求解的紧致格式划分为有限的局部紧致求解,间断点上的激波捕捉铜梁自动作为局部紧致求解的边界通量,在在光滑区域具有紧致格式的高精度低耗散性质,在激波附近不产生非物理振荡.本文利用有限紧致格式思想,构造了新的适合于气动声学问题的优化有限紧致格式,将其应用于计算气动声学一维标准测试问题,对相关格式的模拟性能进行了评估,显示该格式在宽频声波传播和含有间断的声波传播模拟方面具有优势.     

基于混合型WENO格式数值模拟四类Riemann问题

本文基于三阶WENO格式和三阶WENO-Z格式,利用有限体积法研究了同一格式在不同方向使用不同的求解器,以及在不同方向采用相同精度的不同格式数值模拟4类Riemann问题,分析各向异性对数值计算结果的影响．数值结果表明,无论是使用不同的格式,还是使用不同的求解器进行数值模拟,都会导致数值结果不同程度地失去图像对称性．     

一种改进的中心-WENO混合格式

基于中心差分与WENO格式混合可以改善WENO格式耗散特性的思想,在理论推导的基础上,给出了一种用于激波捕捉计算的守恒型中心-WENO混合格式,该混合格式可视为三阶WENO格式和二阶中心差分格式的加权平均。在数值研究现有加权函数的基础上,给出了适用于该混合格式的加权函数,使其能够自适应地调整数值耗散以捕捉激波间断。数值结果表明:与三阶WENO格式相比,混合格式HY3_4能够降低数值耗散,更陡峭地捕捉间断,对复杂流场结构具有较高的分辨率;混合格式HY3_5对于包含高压比激波间断流场结构,能给出无振荡、低耗散的结果。     

非稳态Hamilton-Jacobi方程的7阶加权紧致非线性格式

Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.