求解双曲守恒律方程的高分辨率熵相容格式

为提高熵相容格式的精度,利用限制器机制构造高分辨率格式,将构造的通量限制器插入熵相容格式,得到一类高分辨率熵相容格式.构造Euler方程高分辨率熵相容格式时,对熵相容格式中的几个参数做简单调整,提高了接触间断处的分辨率.将所得格式的数值结果与熵相容格式的数值结果比较表明,构造的高分辨率熵相容格式具有稳健和基本无振荡等特性.     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…     

采用Power限制器的PNND和PWNND格式

为高精度捕捉激波等流场结构,引入一种Power限制器,对NND格式和WNND格式进行改进,分别得到二阶PNND(Power NND)格式和三阶PWNND(Power WNND)格式.该Power类型格式通过Power限制器对相邻待选模板上的一阶导数进行限制,改善了NND格式和WNND格式在间断附近的耗散效应.对各种格式的分析表明,在间断附近采用Power限制器的格式比原格式的表现要好,耗散小且捕捉间断精度高,其中PNND格式虽然只有二阶精度,但在所有算例中与三阶WNND格式的计算结果比较接近,在个别算例中甚至优于WNND格式.最后将PWNND格式应用到二维NACA0012翼型的强迫俯仰振动的数值模拟,计算结果与实验值、参考计算值吻合较好.     

耗散紧致格式的频谱特性研究与应用

采用Fourier分析方法,给出线性耗散紧致格式和基于m级Runge-Kutta时间积分方法的全离散格式的频谱特性,并应用五阶耗散紧致格式模拟典型高频波传播和超声速平面Couette流动的特征值问题及其稳定性边值问题,展示耗散紧致差分格式良好的频谱特性.     

求解可压Navier-Stokes方程的对称超紧致差分格式及其并行算法

考查了超紧致差分方法,并将其精度同传统差分格式和紧致差分格式做了比较,结果显示超紧致方法具有良好求解效率.用分块流水线方法设计了超紧致差分格式的并行算法,进行数值实验及并行性能分析.     

迎风紧致格式的混淆误差分析及其同谱方法的比较

对运用迎风紧致格式求解非线性方程时混淆误差产生的机理进行了研究,通过算例对五阶迎风紧致格式与谱方法进行了比较,发现在混淆误差的处理上迎风紧致格式优于谱方法.     

求解对流扩散方程的紧致修正方法

提出了求解对流扩散方程的紧致修正方法,该方法是在低阶离散格式的源项中,引入紧致修正项,从而构造高阶紧致修正格式,并进行求解.采用紧致修正方法对典型的对流扩散方程进行计算.结果表明,紧致修正方法虽然与二阶经典差分方法建立在相同的结点数上,但紧致修正方法的精度与紧致方法的精度相同,均具有四阶精度.所以紧致修正方法可以在少网...     

一种基于AUSM思想的通量分裂方法

根据对流迎风分裂(AUSM)思想提出一种通量分裂方法,称为K-CUSP格式.它与传统H-CUSP和E-CUSP格式的最大差异在于总能量的分裂:K-CUSP格式将无粘守恒通量中所有的运动学量分裂到对流项,所有的热力学量分裂到压力项,即总能量被分裂成动能和静焓.对于压力项的数值通量,采用一种新的界面构造方法.数值测试表明:①K-CUSP格式继承了FVS格式的简单性和稳健性.在激波后不易出现压力过冲,在膨胀区域没有振荡,优于AUSM和WPS格式;②K-CUSP格式继承了FDS格式的分辨率.激波间断的分辨率和H-CUSP、E-CUSP格式基本相同,接触间断的分辨率高于FVS格式,低于Roe、AUSM和WPS格式.AUSM和WPS格式在计算运动接触间断时,速度存在很大振荡,而新格式不存在振荡.     

求解Navier-Stokes方程组的组合紧致迎风格式

给出一种新的至少有四阶精度的组合紧致迎风(CCU)格式,该格式有较高的逼近解率,利用该组合迎风格式,提出一种新的适合于在交错网格系统下求解Navier-Stokes方程组的高精度紧致差分投影算法.用组合紧致迎风格式离散对流项,粘性项、压力梯度项以及压力Poisson方程均采用四阶对称型紧致差分格式逼近,算法的整体精度不低于四阶.通过对Taylor涡列、对流占优扩散问题和双周期双剪切层流动问题的计算表明,该算法适合于对复杂流体流动问题的数值模拟.     

应用于知通量分裂格式的一类限制器函数

在对绕大纯头体三维有粘跨音速流场进行计算中，ｖａｎＬｅｅｒ知通量分裂格式加通常的顺函数，比如ｖａｎ Ａｌｂａｄａ，ｖａｎＬｅｅｒ和ｍｉｎｍｏｄ型的函数总是引起数值解的“漂移”现象。针对此问题设计了一类新限制器函数，从而数值解稳定收敛到定态解。证明了此类函数可保证格式的ＴＶＤ性质，而且除了在解的极值点附近格式具有二阶精度。     

非结构网格下非定常流场双时间步长的加权ENO-强紧致杂交高分辨率格式

针对三维非定常、可压缩流场的Navier-Stokes方程组,本文提出一种新的双时间步长高精度快速迭代格式。该格式在时间上具有二阶精度,在空间离散上不低于三阶。在对流项与粘性项的处理上,本格式分别采用了加权ENO-强紧致格式与紧致四阶精度格式的思想。几个典型算例的实践表明:计算结果与相关实验数据比较吻合,初步表明了该算法可以在非结构网格下具有高效率与高分辨率的特征。     

A new family of high-order compact upwind difference schemes with good spectral resolution

This paper presents a new family of high-order compact upwind difference schemes. Unknowns included in the proposed schemes are not only the values of the function but also those of its first and higher derivatives. Derivative terms in the schemes appear only on the upwind side of the stencil. One can calculate all the first derivatives exactly as one solves explicit schemes when the boundary conditions of the problem are non-periodic. When the proposed schemes are applied to periodic problems, only periodic bi-diagonal matrix inversions or periodic block-bi-diagonal matrix inversions are required. Resolution optimization is used to enhance the spectral representation of the first derivative, and this produces a scheme with the highest spectral accuracy among all known compact schemes. For non-periodic boundary conditions, boundary schemes constructed in virtue of the assistant scheme make the schemes not only possess stability for any selective length scale on every point in the computational domain but also satisfy the principle of optimal resolution. Also, an improved shock-capturing method is developed. Finally, both the effectiveness of the new hybrid method and the accuracy of the proposed schemes are verified by executing four benchmark test cases.     

基于气动声学问题的高阶非线性紧致格式

文章基于线性中心紧致差分格式, 通过非线性加权插值的方法来求解网格中心处的函数值.这类格式保持了原有中心紧致差分格式的高阶精度和低耗散特性, 同时其分辨率也非常高, 由于其非线性插值的机制, 使得这类格式能够捕捉强激波, 所以这类新的高阶非线性紧致格式是一种较好的模拟湍流和气动声学等多尺度问题的方法.      

迎风紧致格式求解Hamilton-Jacobi方程

基于Hamilton-Jacobi(H-J)方程和双曲型守恒律之间的关系,将三阶和五阶迎风紧致格式推广应用于求解H-J方程,建立了高精度的H-J方程求解方法.给出了一维和二维典型数值算例的计算结果,其中包括一个平面激波作用下的Richtmyer Meshkov界面不稳定性问题.数值试验表明,在解的光滑区域该方法具有高精度,而在导数不连续的不光滑区域也获得了比较好的分辨效果.相比于同阶精度的WENO格式,本方法具有更小的数值耗散,从而有利于多尺度复杂流动的模拟中H-J方程的求解.     

Generalized finite compact difference scheme for shock/complex flowfield interaction

A class of generalized high order finite compact difference schemes is proposed for shock/vortex, shock/boundary layer interaction problems. The finite compact difference scheme takes the region between two shocks as a compact stencil. The high order WENO fluxes on shock stencils are used as the internal boundary fluxes for the compact scheme. A lemma based on the property of smoothness estimators on a 5-points stencil is given to detect the shock position. There is no free parameter introduced to switch the compact scheme and the WENO scheme. Some numerical experiments are given and they demonstrate that the present scheme has low dissipation due to the compact central differencing scheme used in the smooth regions.     

Robust HLLC Riemann solver with weighted average flux scheme for strong shock

Many researchers have reported failures of the approximate Riemann solvers in the presence of strong shock. This is believed to be due to perturbation transfer in the transverse direction of shock waves. We propose a simple and clear method to prevent such problems for the Harten–Lax–van Leer contact (HLLC) scheme. By defining a sensing function in the transverse direction of strong shock, the HLLC flux is switched to the Harten–Lax–van Leer (HLL) flux in that direction locally, and the magnitude of the additional dissipation is automatically determined using the HLL scheme. We combine the HLLC and HLL schemes in a single framework using a switching function. High-order accuracy is achieved using a weighted average flux (WAF) scheme, and a method for v-shear treatment is presented. The modified HLLC scheme is named HLLC–HLL. It is tested against a steady normal shock instability problem and Quirk’s test problems, and spurious solutions in the strong shock regions are successfully controlled.     

High-order kinetic flux vector splitting schemes in general coordinates for ideal quantum gas dynamics

A class of high-order kinetic flux vector splitting schemes are presented for solving ideal quantum gas dynamics based on quantum statistical mechanics. The collisionless quantum Boltzmann equation approach is adopted and both Bose–Einstein and Fermi–Dirac gases are considered. The formulas for the split flux vectors are derived based on the general three-dimensional distribution function in velocity space and formulas for lower dimensions can be directly deduced. General curvilinear coordinates are introduced to treat practical problems with general geometry. High-order accurate schemes using weighted essentially non-oscillatory methods are implemented. The resulting high resolution kinetic flux splitting schemes are tested for 1D shock tube flows and shock wave diffraction by a 2D wedge and by a circular cylinder in ideal quantum gases. Excellent results have been obtained for all examples computed.