高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

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

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

自适应的高精度中心﹣WE NO混合格式

为对可压缩流动进行高精度高分辨率的数值模拟,基于模板的光滑度量函数,构造了新型的加权因子,在此基础上将7阶精度的WENO格式与8阶中心格式进行加权组合,构造了一种自适应低耗散的中心-WENO混合格式(H-WENO7-CD8),并采用Fourier方法对离散格式数学特性进行了理论分析。该自适应的高精度中心-WENO混合格式相比于七阶的WENO格式具有更小的耗散误差。通过对激波-密度干扰问题、Rayleigh-Taylor不稳定性问题和双Mach反射问题的数值模拟,并与WENO-JS格式的计算结果对比,结果表明:该格式结合了WENO格式和中心格式的特点,能更好地捕捉激波,对物理脉动也具有较高的分辨率,适合于复杂流体流动的数值模拟。     

AUSM+格式的改进

为了提高AUSM 格式对流场计算中激波和剪切流的分辨率,减小数值振荡和激波后数值过冲问题,本文在原有AUSM 格式的基础上经过数值研究分析,对AUSM 格式进行了声速、特征马赫数的处理方式,对网格界面马赫数等进行了改进,并用改进后的AUSM 格式求解了含有脱体激波、斜激波、剪切流、膨胀扇、及湍流边界层等流动特征的复杂流动.经典型算例检验证实,改进后的AUSM 格式对激波和剪切流的分辨率得到了加强,求解的准确性和收敛性有一定的提高.改进后的AUSM 格式继承了原AUSM 格式优点,计算处理问题能力增强,可以有效地应用于亚声、跨声和高超声流动的工程实际问题中.     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

CPR方法与高分辨率二阶格式的混合算法

重构修正方法（correction procedure via reconstruction,CPR）具有紧致高效的优点,但对较强激波的捕捉能力还相对较弱,而加权紧致非线性格式（weighted compact nonlinear scheme,WCNS）具有很强的激波捕捉能力。将基于高阶WCNS插值的二阶格式引入到高阶CPR方法中,构造了一种高效高分辨率的混合激波捕捉格式。首先,基于非线性权偏离线性权的程度的激波侦测器侦测出问题单元,并在问题单元附近引入缓冲单元,其余单元则标记为光滑单元。然后,针对问题单元和缓冲单元采用二阶格式计算,光滑单元采用CPR方法计算,构造混合格式。通过对等熵涡问题、含激波的问题以及激波旋涡干扰问题的数值模拟,测试了混合格式的精度、激波捕捉能力和计算效率。数值模拟结果充分说明了该混合格式具有很强的激波捕捉能力,同时在光滑区具有高分辨特性,可以应用于高超声速流动问题的高效数值模拟中。相比于基于高阶WCNS插值的二阶格式,此格式具有更高的计算效率和更高的分辨率。     

先进上风格式跨声速CFD计算研究

介绍了空气动力学跨音速计算中应用广泛的几种上风格式,同时引入了JAMSON中心差分格式作为比较,针对激波分辨率,格式精度,计算效率等方面的问题对ONERAM6机翼进行了对比性计算,讨论了ROE的FDS格式,AUSM+格式和AUSMDV格式的计算性能,得出了很有意义的结论。     

有关有限差分高精度格式两个应用问题的讨论

激波装配法结合有限差分高精度格式是发展全场一致高精度算法的一种途径。在对流场内的间断全部进行装配后,对高精度格式的应用研究成为发展本算法的主要研究问题。本文将有限差分的高精度格式应用于贴体坐标系时发现,对均匀流场,高精度格式因不满足几何守恒律而产生的数值误差比一阶迎风格式大,初步分析认为是由于高精度格式所用的模板比一阶格式更宽,涉及的网格点数更多,从而引入了更多的误差。而作者提出的基于离散等价方程的相容性算法可消除这一误差。此外,本文在利用激波捕捉法求解正方形均匀网格上的正激波运动问题时发现因通量分裂格式的使用,在激波处会产生随着特征线传播的非物理波动,这一波动在激波与网格不完全匹配时表现为多维波动相互干扰的虚假"数值湍流"现象,高精度格式的高分辨率特性使得这一现象更加明显。这是因为激波捕捉法假设激波为空间连续函数,用于包含激波的流场时必然得到数值解表示的过渡区,导致可信度评估困难,使用激波装配法可以避免这一问题。     

高精度有限差分格式的色散优化及耗散控制

本文综述了我们在高精度有限差分格式的色散优化和耗散控制方面的研究进展。首先,我们提出了半离散有限差分格式色散和耗散相互独立的充分条件,实现了差分格式色散和耗散特性的独立调节。在此基础上,提出了色散最小、耗散可控的高精度差分格式,称为MDCD格式。MDCD格式已经得到了广泛应用,取得了很好的计算效果,但其主要缺点是耗散的调节依赖于经验。为了解决这一问题,我们进一步提出耗散的自适应调节方法。具有自适应耗散特性的高精度有限差分格式的基本特征是,差分格式的耗散能够随解的局部尺度自适应调节。为了构造这类格式,我们提出了一种新型的尺度识别器,它能够以等效无量纲波数的形式来定量衡量数值解的局部长度尺度。在此基础上,设计差分格式耗散参数与尺度识别器得到的等效无量纲波数之间的关系,从而构造了一类色散最小、耗散自适应的差分格式,称为MDAD格式。为了计算含有间断的问题,同时保持在光滑区的良好耗散特性,我们利用尺度识别器对一种经典的激波探测器进行改进,提出了一种新的激波探测器,并将自适应耗散格式与对应的WENO格式相混合,得到自适应耗散混合格式。近似色散关系显示该混合格式兼具高分辨率和鲁棒性。多个含间断流场的标准算例测试结果显示,自适应耗散混合格式具有良好的分辨率和激波捕捉能力。     

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

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

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

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

基于边界变差最小的高精度有限差分格式构造

从加权紧致非线性（WCNS）高精度有限差分格式出发,在其重构形式的基础上,根据边界变差最小（BVD）原理,遵循单元边界两侧重构物理量值之差最小的准则,在每个单元内通过两步空间重构,构造了一种新的高精度有限差分格式。一般对WCNS等加权非线性格式的改进都是基于改善色散耗散特性、优化非线性权、提高分辨率等单一途径,将它们作为重构候补函数进行结合,既保持了各自优势所在,又能控制格式整体黏性,所得格式具有丰富的应用场景。通过数值实验,将结果与单一格式进行对比,新格式既能在流场光滑区保证设计的精度,对激波等间断附近的振荡也有很好的抑制作用,提高了对高波数区的分辨率,而且在长时间计算后也有较为精确的结果。面向广泛发展的数值格式,还可以构造出其他新方法,对包含强间断和多尺度的流动问题可以获得更好的结果。     

适用于超声速的一种通量限制型紧致格式

紧致格式因其结构简单、在相同的网格点上能达到比非紧致格式更高的精度以及与谱方法相近的分辨率等优点,日益受到人们的重视。用紧致格式模拟超声速流场的主要问题之一是如何保证高阶紧致格式能光滑地捕捉到流场的各种间断。本文借鉴NND格式的思想,构造出一种总体上具有三阶精度的通量限制型紧致(LFC)格式,并成功地应用于含有激波、滑移面等复杂流动现象的数值模拟。计算结果表明这种格式不仅具有较高的精度和分辨率,而且还保证了在间断附近基本无虚假波动。     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

一种新的双重门限证实数字签名方案

在已有的(t,n)门限数字签名方案及证实数字签名的基础上,提出了新的双重门限证实数字签名方案(t1,n1;t2,n2)。新方案首次在签名证实阶段引入了门限方案,由一群符合一定条件的t2个人组成证实者(confirm er)。因为是由不同的各方参与进来,方案具有了更加广泛的可信任性,从而使得该方案适用于某些安全程度要求较高的情况。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。     

Scheme design of the CHANG’E-5T1 extended mission

Flight schemes for the CHANG’E-5T1 extended mission are investigated in this paper. In the flight scheme and trajectory design, the remaining propellant of the CHANG’E-5T1 mission is utilized. The CHANG’E-5T1 mission is firstly introduced with feasible flight goals derived based on the terminal trajectory and satellite status. The flight schemes are designed to include a lunar return and the libration points in the Sun-Earth/Moon and Earth-Moon systems, with an emphasis on the Earth-Moon triangle libration point thus far unexplored. Secondly, three schemes are proposed for the CHANG’E-5T1 extended mission with different flight goals. The direct libration point orbit transfer and injection method is adopted to solve the issue in the transfer trajectory design. Furthermore, an innovative concept is proposed to transfer from the Earth-Moon collinear libration point to the triangle point using the Sun-Earth/Moon libration point. Finally, the merits and drawbacks of the three schemes are discussed in terms of flight time, control energy and frequency, flight distance, and goal value. As a result, the scheme including a lunar return and the Earth-Moon L2 libration point is selected for the CHANG’E-5T1 extended mission. A flight to the Earth-Moon libration point is achieved, replicating the achievement of the ARTEMIS mission.     

气动声学直接数值模拟的空间差分格式分析

流场内声波的直接数值模拟对数值格式的耗散、色散特性提出了严格要求,首先介绍了两种适于声学模拟空间差分格式——紧致差分格式与色散关系保持(DRP)格式的基本思想,然后分析了格式的色散及其与四阶龙格-库塔法结合时的耗散特性。根据Fourier误差分析和数值模拟发现:无论从格式可有效求解的波数范围还是从与龙格-库塔时间推进结合后的耗散特性来看六阶紧致格式的可有效求解范围要大于DRP格式;不过,数值计算表明紧致格式的计算量远远超过DRP格式,对计算机资源提出了较高的要求。      

基于BVD原理的高保真空间重构方法

简要综述了一类基于单元边界变差最小化(Boundary Variation Diminishing,BVD)原理,设计双曲守恒律高保真数值格式的空间重构方法。BVD原理要求尽量减少通过重构得到的网格边界两侧物理量之间的差,从而能够有效地控制黎曼求解器中的数值黏性。BVD方法针对数值解的空间分布特征,选择多个函数作为空间重构的候补函数,并根据BVD判定准则从候补函数中选取最合适的函数进行空间重构。BVD判据不需要根据求解对象进行经验参数(阈值)的调整。选用适当的候补函数和BVD准则,可以完全避免现有算法中为抑制数值振荡而必须采用的非线性限制。BVD格式能在抑制数值振荡的同时,有效地控制数值耗散,可以对光滑解与间断解都获得高保真的计算结果。本文概述了BVD方法的基本思想、设计相关格式的基本思路,以及一些具有很强实用价值的BVD格式。并通过单相和两相可压缩流动的一些典型算例验证BVD格式的特点和优势。