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

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

加权梯度限制器在节点DG格式中的应用

在使用节点间断Galerkin (DG)方法求解存在激波的流动问题时,为抑制Gibbs效应导致的非物理振荡,对一种加权梯度限制器进行研究,在二维非结构网格上对定常欧拉方程求解.限制器以原始变量作为限制对象,对其梯度进行限制,以加权形式进行重构.对翼型的跨音速流场、超燃冲压发动机内流场和前台阶模型流场进行模拟,结果表明限制器能很好抑制间断解附近的非物理振荡,有效避免计算中断,体现了DG格式精度高、对间断捕捉敏感的特点.     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

非结构网格通量重构算法下三种紧致WENO限制器对比研究

在通量重构算法框架下对比研究了3种最新提出的用于间断伽辽金算法的紧致型WENO(weighted essentially non-oscillatory)限制器，即简单WENO限制器、p阶加权WENO限制器和多精度WENO限制器。这3种WENO限制器能够较高精度地模拟流场并捕捉激波，其紧致性在于问题网格单元中的重构只涉及问题网格单元及其相邻网格单元。同时，使用保精度的保正限制器用以避免可能出现的非物理的密度和压力负值。最后，采用非结构四边形网格在双马赫反射、激波与涡相互作用、激波反射、黏性弓形激波和激波与层流边界层相互作用等多个二维算例中，对这3种WENO限制器进行分析比较。结果显示：多精度WENO限制器与p阶加权WENO限制器能够高精度模拟流场并捕捉激波，同时一定程度抑制通量重构方法本身在激波附近的数值伪振荡；p阶加权WENO限制器与多精度WENO限制器相比，其稳态收敛性相对更好；简单WENO限制器则性能较差。根据以上研究，提出亟须发展能够使高阶WENO限制器在稳态问题中收敛的间断探测器。     

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

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

高精度IPDG湍流模拟及通量格式数值特性北大核心CSCD

为更高精度数值求解复杂湍流问题,在内罚间断伽辽金(internal penalty discontinuous Galerkin,IPDG)方法框架内发展了基于SST k-ω模型的湍流模拟方法,通过对亚/跨/超声速工况下流场的数值计算与湍流特征捕捉,验证了方法的适用性,进而系统分析了AUSM、Lax-F、HLL和Roe 4种通量格式在IPDG湍流模拟中的数值特性。结果表明:AUSM格式在超声速工况下脱体激波面“褶皱”现象明显;Lax-F格式和HLL格式数值耗散大,在激波解析方面精度较低,且Lax-F格式在激波脚后诱导产生流动分离;Roe格式具备宽速域适用性,计算精度较高且能精确解析激波结构。     

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

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

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

一种求解激波问题的中心差分-WENO混合方法研究

为提高求解包含激波问题的计算精度和效率,发展了一种低耗散、高效率的中心差分-WENO(Weighted Essentially Non-Oscillatory)混合格式,Navier-Stokes方程的无黏项在光滑流场区域采用六阶中心差分格式离散,而间断附近采用五阶WENO格式求解;基于密度设计了一种新型格式开关实现两种格式在光滑区域和间断之间自动切换,确保数值解在间断附近基本无振荡。针对一维激波熵波作用问题验证了混合格式的低耗散特性;对二维Riemann问题的研究表明发展的基于密度的格式开关更为合理,混合格式的计算效率较WENO格式明显提高;将其应用到激波诱导燃烧问题中,混合格式能很好地捕捉定常流场中激波和化学反应锋面的位置以及非定常流场的振荡频率。     

WENO格式的一种新型光滑因子及其应用

光滑因子是构造WENO格式的关键，决定了WENO格式能否达到最优收敛精度以及在间断附近能否保持本质无振荡特性。针对广泛应用的五阶WENO格式，采用算子函数近似导数关系的方法，设计了一种新型光滑因子。在间断区域，与传统的光滑因子相比，新型光滑因子对间断的识别更准确。在光滑区域，新型光滑因子使得新型WENO格式的非线性权更接近线性权。理论分析和数值验证表明新型WENO格式即使在一阶临界点也能保持一致五阶精度。一维典型激波问题的数值结果表明，与现有WENO格式相比，新型WENO格式提高了短波的分辨率，降低了格式的耗散，同时能够准确识别间断。对于典型包含激波和剪切层的流动问题，新型WENO格式不仅能够准确地捕捉强激波，而且能够精细模拟剪切层和声波等流场结构。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

任意单元间断Galerkin有限元计算方法研究

基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。     

高精度两步TVD格式的构造及数值检验

基于通量分裂的思想,利用Taylor级数理论和第二步限制器函数的控制,构造了在光滑区域空间为三阶、时间为一阶精度的高分辨率全变差递减(TVD)格式.该格式在激波过渡点降为一阶迎风格式.并从单个线形方程推广到非线性方程及方程组情形.通过几个典型算例的计算,并与二阶TVD格式作了比较,表明该方法对流场中的激波有较高的分辨率,且是无波动的.     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

角度-空间间断有限元的热辐射限制器

由于航空发动机内高温热端部件是非规则形状的，其中的热辐射过程存在明显的间断遮蔽效应，如何实现该过程的精准模拟是一个挑战。采用间断有限元对辐射传递方程中角度和空间计算域进行离散，并基于分层限制的思想在角度和空间单元上引入巴斯-叶斯帕森限制器，对辐射强度数值解在角度和空间上的非物理振荡进行抑制。通过与文献中有限体积法和蒙特卡洛辐射模型的结果对比，证明了方法的有效性。此外，通过与辐射强度在角度上的解析解对比，验证了该限制器可以抑制辐射间断效应造成的数值振荡，消除辐射强度的非物理解。实现了任意空间位置上辐射强度在角度方向上的三维刻画。      

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

非结构网格划分的喷管流场计算中单调性约束条件的研究

提出了一种限制网格单元内梯度值大小的单调性约束条件,以减少间断附近数值解的振荡。这一约束条件在算法上具有易于实现的优点,因此能够很方便地用于非结构网格划分的流场计算中。