一种具有较小数值耗散的对称TVD格式

为了减小TVD格式的数值粘性,提高TVD格式的分辨度,提出了一种具有较小数值耗散的Harten型TVD格式.数值模拟结果表明,这种Harten型TVD格式比原来的Harten格式^[1]及其改进形式Harten-Yee^[2]的TVD格式能更好地模拟出边界层及涡流.     

2阶精度时间后差隐式格式TVD充分条件

非定常流动问题计算中常用到含三个时间层的二阶精度时间后差隐式格式,并且希望构成TVD格式,然而理论上的问题多年一直没有解决。本文找到了解决办法,构造了这种类型的隐式TVD格式,证明了其为TVD的充分条件。理论结果为计算所验证,并表明通常未采取本文对时间差分处理方法的格式尚不具备TVD性质。     

非定常隐式TVD格式的基本理论和验证

用于非定常问题的时间二阶精度的隐式TVD格式的理论多年来一直是个空白。本文按新的途径构造了一种时间二阶精度的隐式TVD格式,以发展了的手段证明了其为无条件(任意时间步长)TVD的,且限量因子Φ不受传统上界的限制,具体构造了这样的新的Φ。所得理论结果均为计算所验证,格式显示了良好的稳定性和分辨率。     

具有TVD性质的三阶精度GODUNOV格式在粘性流场计算中的应用

本文发展了一种具有TVD性质的三阶精度的Godunov格式。隐式部分采用迎风对角线形隐式近似因式分解法。并引入了粘性通量的简化算法。显式部分采用三阶精度TVD格式。为进一步增强格式的稳定性及对间断的捕捉能力,在单元边界上构造Riemann问题。应用上述方法对某型涡轮高压级及压气机进行了数值模拟。     

用TVD格式求解一维滑坡问题

本文用二阶精度显式TVD格式计算了一维滑坡问题,得到了和实际值基本一致的结果。     

湍流模式对转捩的初步研究

采用二阶时间精度的LU-SGS-τTS目的进行时间推进求解Navier-Stokes方程组,用基于Roe的六阶对称TVD格式离散对流项,并采用Baldwin-Lomax和Spalart-Allmaras模式对平板边界层、Aerospatial A-profile翼型低速大迎角及RAE-2822翼型跨声速激波/边界层干扰流动进行数值分析,自动捕捉流动的转捩过程,获取了有意义的结果.     

涡轮叶栅气热耦合数值方法与验证

应用三阶精度TVD格式和自由型曲面阿格技术,对涡轮内冷叶片进行气热耦合数值模拟.结合商用软件采用多种湍流模型计算方法的结果进行实验验证和对比分析,结果表明,商用软件采用传统湍流模型对温度边界层的模拟是不合理的,采用考虑转捩的湍流模型,调整转捩起始雷诺数,可以较好地模拟温度边界层的热传导.本文对非转捩区边界层温度的计算结果与实验结果吻合较好.     

非定常喷流有限差分方法数值试验

以欠膨胀自由喷流初期流动为例,采用Euler方程和Beam Warming TVD有限差分格式,对比分析了Beam WarmingAF方法,隐式亚迭代方法和简化Runge Kutta五步格式的非定常流场描述能力,结果表明:(1)隐式近似因式分解方法基本上可以描述非定常流动现象;(2)隐式亚迭代一阶时间精度格式会导致流场结构的变化,其精度可能是不足的;(3)隐式亚迭代二阶时间精度与简化Runge Kutta五步格式的计算结果一致,可以认为是计算非定常问题的适当方法.     

隐式格式求解拟压缩性非定常不可压Navier-Stokes方程

采用Rogers发展的双时间步拟压缩方法,数值求解不可压非定常问题.数值通量分别采用三阶精度Roe格式和二阶精度Harten-Yee的TVD格式离散.为了加快收敛,提高求解效率,试验了几种隐式格式(ADI-LU,LGS,LU-SGS).针对经典的低雷诺数(Re=200)圆柱绕流问题,比较了不同隐式方法的计算结果和求解效率,以及两种数值离散格式计算结果的异同.最后采用Roe格式数值求解了两种典型的低速非定常流动问题:绕转动圆柱(ω=1)低雷诺数流动;NACA0015翼型等速拉起数值模拟.     

涡流数值模拟中的计算格式粘性分辨率探讨

采用理论分析和三角翼数值模拟实验相结合的方法,研究了目前常用的CFD计算格式——Jameson中心格式、RoeFDS格式、VanLeerFVS格式、迎风TVD格式等的粘性分辨率,论述了Jameson中心格式和RoeFDS格式在粘性分辨率上具有的优势,以及VanLeerFVS格式和迎风TVD格式在此方面的不足     

Multi-Symplectic Method for the Zakharov-Kuznetsov Equation

A new scheme for the Zakharov-Kuznetsov (ZK) equation with the accuracy order of $\mathcal{O}(∆t^2+∆x+∆y^2)$ is proposed. The multi-symplectic conservation property of the new scheme is proved. The backward error analysis of the new multi-symplectic scheme is also implemented. The solitary wave evolution behaviors of the Zakharov-Kunetsov equation are investigated by the new multi-symplectic scheme. The accuracy of the scheme is analyzed.     

径向对称非线性Schrdinger方程的守恒差分格式

对径向对称的非线性Schrdinger方程提出了一个新的守恒差分格式,这是一个三层格式,它不需迭代求解,因此提高了计算速度,同时也较好地保持了方程的两个守恒律。文中证明了格式的收敛性与稳定性,数值计算结果表明,该差分格式是有效的。     

An improved monotone finite volume scheme for diffusion equation on polygonal meshes

We construct a new nonlinear monotone finite volume scheme for diffusion equation on polygonal meshes. The new scheme uses the cell-edge unknowns instead of cell-vertex unknowns as the auxiliary unknowns in order to improve the accuracy of monotone scheme. Our scheme is locally conservative and has only cell-centered unknowns. Numerical results are presented to show how our scheme works for preserving positivity on various distorted meshes. Specially, numerical results show that the new scheme is robust, and more accurate than the existing monotone scheme on some kinds of meshes.     

一种新的LU型隐式格式及其应用

发展了一种新的ＬＵ型隐式格式。新格式除包含ＬＵ－ＳＧＳ格式全部特色外，尽管采用了相似变换而不是近似处理来精确地构造通量矢量的迎风Ｊａｃｏｂｉａｎ矩阵，仍可避免块对角矩阵求逆．且左端项仍可保持矢量化处理。新格式显示出更快的收敛性和更高的稳定性且没有多余的数值耗散和自由参数，尤其适用于求解三维非定常流动问题。     

Third-order Energy Stable WENO scheme

A new third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme for scalar and vector hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be linearly stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. We also present new weight functions which drastically improve the accuracy of the third-order ESWENO scheme. Based on a truncation error analysis, we show that the ESWENO scheme is design-order accurate for smooth solutions with any number of vanishing derivatives, if its tuning parameters satisfy certain constraints. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO scheme of Jiang and Shu in terms of accuracy, while providing essentially non-oscillatory solutions near strong discontinuities.     

阻性管道消声器消声量方程的一个新算法

将阻性管道消声器消声量特征方程（组）转化为常微分方程（组），由此给出消声量计算的一个新方法。该方法毋需迭代因而不存在计算稳定性问题，计算量小，初值的迭取也很简单，适合于扁矩形、任意矩形及圆形等各种截面管道消声器，尤其对多变量、考虑平均流影响的情形更为有效。     

A 5th order monotonicity-preserving upwind compact difference scheme

Based on an upwind compact difference scheme and the idea of monotonicity-preserving, a 5th order monotonicity-preserving upwind compact difference scheme (m-UCD5) is proposed. The new difference scheme not only retains the advantage of good resolution of high wave number but also avoids the Gibbs phenomenon of the original upwind compact difference scheme. Compared with the classical 5th order WENO difference scheme, the new difference scheme is simpler and small in diffusion and computation load. By emplo...     

一种处理漫射因子的新方法

提出了一种新的处理漫射因子的拟合方案,通过与前人提出的方案相比,具有较高的拟合精度;然后,将该拟合方案和前人提出的方案同时放入辐射传输模式中来计算大气加热率,结果表明:新方案的结果总体上优于前人的结果,尤其是对45 km以上的中间层大气. 关键词： 漫射因子 指数积分 加热率     

An implicit FDTD scheme for the propagation of VLF–LF radio waves in the Earth–ionosphere waveguide

A new finite-difference time-domain scheme is presented for the propagation of VLF–LF radio waves in the Earth–ionosphere waveguide. The new scheme relies on the implicit solution of the auxiliary equation that governs the current density in the ionosphere. The advantages and drawbacks of the new scheme are discussed. Its main advantage is its stability condition, which is the same as that of the FDTD method in a vacuum. This permits the time step of the calculation to be increased and then the overall computational time to be reduced. Numerical experiments demonstrate the accuracy of the new scheme and the reduction of the computational time.     

A conservative numerical scheme for solving an autonomous Hamiltonian system

A new numerical scheme is proposed for solving Hamilton’s equations that possesses the properties of symplecticity. Just as in all symplectic schemes known to date, in this scheme the conservation laws of momentum and angular momentum are satisfied exactly. A property that distinguishes this scheme from known schemes is proved: in the new scheme, the energy conservation law is satisfied for a system of linear oscillators. The new numerical scheme is implicit and has the third order of accuracy with respect to the integration step. An algorithm is presented by which the accuracy of the scheme can be increased up to the fifth and higher orders. Exact and numerical solutions to the two-body problem, calculated by known schemes and by the scheme proposed here, are compared.