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

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

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

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

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

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

三维高超声速无粘定常绕流的数值模拟

本文采用一种简单有效的通量分裂结合一种二阶TVD格式的数值通量的方法,提出一种隐式的迎风有限体积格式,并利用这种格式,从气体动力学非定常Euler方程组出发,数值模拟了三维不对称物体的高超声速无粘定常绕流。数值结果表明此格式具有分辨率较高和收敛速度较快的优点。     

多车种LWR交通流模型的半离散中心迎风格式

对多车种LWR交通流模型,给出一种半离散中心迎风格式,该格式以五阶WENO-Z重构和半离散中心迎风数值通量为基础.WENO-Z重构方法的引入提高了格式的精度,并保证格式具有基本无振荡的性质.时间的离散采用保持强稳定性的Runge-Kutta方法.通过数值算例验证了格式的有效性.     

粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

随机扰动对拟小波方法求解对流扩散方程的影响

引进拟小波方法数值求解对流扩散方程，研究结果表明，计算带宽W有一个极值，当计算带宽W取该极值时，该方程的拟小波解的精度最高，且好于迎风格式。当边界发生随机不等幅扰动时，对于积分时间较长的情况，拟小波格式的效果要稍逊于迎风格式；当边界发生随机等幅扰动时，若计算带宽W取大于等于20的整数时，方程拟小波解的精度与迎风格式相同；当参数受到随机扰动时，W取10时的拟小波解的均方根误差要小于迎风格式；在初值发生随机扰动且计算带宽W取10时，方程的拟小波解的精度最高，好于迎风格式。     

非结构网格上的迎风有限元格式

在非结构网格上提出一种基于修正积分区域的迎风有限元格式,它与一阶迎风差分格式相当,可应用于构造各种不同的数值格式。     

摄动有限体积法重构近似高精度的意义

研讨有限体积(FV)方法重构近似高精度的作用问题.FV方法中积分近似采用中点规则为二阶精度时,重构近似高精度(精度高于二阶)的意义和作用是一个有争议的问题.利用数值摄动技术[1,2]构造了标量输运方程的积分近似为二阶精度、重构近似为任意阶精度的迎风型和中心型摄动有限体积(PFV)格式.迎风PFV格式无条件满足对流有界准则(CBC),中心型PFV格式为正型格式,两者均不会产生数值振荡解.利用PFV格式求解模型方程的数值结果表明:与一阶迎风和二阶中心格式相比,PFV格式精度高、对解的间断分辨率高、稳定性好、雷诺数的适用范围大,数值地"证实"重构近似高精度和PFV格式的实际意义和好处.     

一类迎风紧致差分格式全离散特性分析

采用Ｆｏｕｒｉｅｒ分析方法，通过显式多步Ｒｕｎｇｅ－Ｋｕｔｔａ时间离攻一维线性对流方程，导出了一类高精度迎风紧致格式全离散色散关系式，详细分析了不同ＣＦＬ数下所研究差分格式的耗散、色散及相应的相速度、群速度等特性。以数值实验显示了格式较高的计算精度和分辨率。     

New embedded pairs of explicit Runge-Kutta methods with FSAL properties adapted to the numerical integration of oscillatory problems

A higher order Runge-Kutta (pair) method specially adapted to the numerical integration of IVPs with oscillatory solutions is presented. This method is based on the adapted methods proposed by Franco (see Ref. [J.M. Franco, Appl. Numer. Math. 50 (2004) 427]). We give explicit method (up to order 5) as well as pairs of embedded Runge-Kutta methods of order 5 and 4 designed using the FSAL properties. The stability of the new methods is analyzed. The numerical experiments are carried out to show the efficiency and robustness of our methods in comparison with some efficient methods.     

A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems

We present a fully second order implicit/explicit time integration technique for solving hydrodynamics coupled with nonlinear heat conduction problems. The idea is to hybridize an implicit and an explicit discretization in such a way to achieve second order time convergent calculations. In this scope, the hydrodynamics equations are discretized explicitly making use of the capability of well-understood explicit schemes. On the other hand, the nonlinear heat conduction is solved implicitly. Such methods are often referred to as IMEX methods [2], [1], [3]. The Jacobian-Free Newton Krylov (JFNK) method (e.g. [10], [9]) is applied to the problem in such a way as to render a nonlinearly iterated IMEX method. We solve three test problems in order to validate the numerical order of the scheme. For each test, we established second order time convergence. We support these numerical results with a modified equation analysis (MEA) [21], [20]. The set of equations studied here constitute a base model for radiation hydrodynamics.     

Explicit multi-symplectic methods for Klein–Gordon–Schrödinger equations

In this paper, we propose explicit multi-symplectic schemes for Klein–Gordon–Schrödinger equation by concatenating suitable symplectic Runge–Kutta-type methods and symplectic Runge–Kutta–Nyström-type methods for discretizing every partial derivative in each sub-equation. It is further shown that methods constructed in this way are multi-symplectic and preserve exactly the discrete charge conservation law provided appropriate boundary conditions. In the aim of the commonly practical applications, a novel 2-order one-parameter family of explicit multi-symplectic schemes through such concatenation is constructed, and the numerous numerical experiments and comparisons are presented to show the efficiency and some advantages of the our newly derived methods. Furthermore, some high-order explicit multi-symplectic schemes of such category are given as well, good performances and efficiencies and some significant advantages for preserving the important invariants are investigated by means of numerical experiments.     

A NUMERICAL COMPARISON OF ALTERNATIVE GALERKIN METHODS FOR EIGENVALUE ESTIMATION

This paper provides a concise and unified comparison of four distinct variations of Galerkin's method. The four Galerkin variations are the explicit (traditional), implicit, quadratic implicit, and diagonalized-implicit Galerkin methods. Results indicate that the explicit Galerkin method is superior to all implicit formulations. Among the implicit methods, the quadratic implicit and diagonalized-implicit perform equivalently and are significantly better than the standard implicit method. The explicit method is recommended in all the cases except those in which the number of trial functions is so large that numerical conditioning affects the eigenvalue estimate. In this case, the diagonalized-implicit method is recommended.     

A second order self-consistent IMEX method for radiation hydrodynamics

We present a second order self-consistent implicit/explicit (methods that use the combination of implicit and explicit discretizations are often referred to as IMEX (implicit/explicit) methods ,  and ) time integration technique for solving radiation hydrodynamics problems. The operators of the radiation hydrodynamics are splitted as such that the hydrodynamics equations are solved explicitly making use of the capability of well-understood explicit schemes. On the other hand, the radiation diffusion part is solved implicitly. The idea of the self-consistent IMEX method is to hybridize the implicit and explicit time discretizations in a nonlinearly consistent way to achieve second order time convergent calculations. In our self-consistent IMEX method, we solve the hydrodynamics equations inside the implicit block as part of the nonlinear function evaluation making use of the Jacobian-free Newton Krylov (JFNK) method ,  and . This is done to avoid order reductions in time convergence due to the operator splitting. We present results from several test calculations in order to validate the numerical order of our scheme. For each test, we have established second order time convergence.     

Energy stable numerical methods for hyperbolic partial differential equations using overlapping domain decomposition

Overlapping domain decomposition methods, otherwise known as overset grid or chimera methods, are useful for simplifying the discretization of partial differential equations in or around complex geometries. Though in wide use, such methods are prone to numerical instability unless numerical diffusion or some other form of regularization is used, especially for higher-order methods. To address this shortcoming, high-order, provably energy stable, overlapping domain decomposition methods are derived for hyperbolic initial boundary value problems. The overlap is treated by splitting the domain into pieces and using generalized summation-by-parts derivative operators and polynomial interpolation. New implicit and explicit operators are derived that do not require regularization for stability in the linear limit. Applications to linear and nonlinear problems in one and two dimensions are presented, where it is found the explicit operators are preferred to the implicit ones.     

Upwind monotonic interpolation methods for the solution of the time dependent radiative transfer equation

We describe the use of upwind monotonic interpolation methods for the solution of the time-dependent radiative transfer equation in both optically thin and thick media. These methods, originally developed to solve Eulerian advection problems in hydrodynamics, have the ability to propagate sharp features in the flow with very little numerical diffusion. We consider the implementation of both explicit and implicit versions of the method. The explicit version is able to keep radiation fronts resolved to only a few zones wide when higher order interpolation methods are used. Although traditional implementations of the implicit version suffer from large numerical diffusion, we describe an implicit method which considerably reduces this diffusion.     

Numerical stability for laser radiation computations

The numerical stability for solving the coupled differential equations for laser radiation intensity and gain is investigated. Simplified lasing conditions are considered under which a closed-form solution is obtained and numerical stability can be determined. The analytical solution illustrates the effect of the important physical parameters, and it provides a basis for evaluation of numerical methods. It is shown that the explicit numerical procedure is inherently unstable, and that stability requires an implicit approach. Partially-implicit methods, which are conditionally stable, are also discussed.     

Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge–Kutta–Nyström (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge–Kutta method, which is equivalent to the well-known Störmer–Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine–Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation.     

A discrete action principle for electrodynamics and the construction of explicit symplectic integrators for linear,non-dispersive media

In this work, we derive a discrete action principle for electrodynamics that can be used to construct explicit symplectic integrators for Maxwell’s equations. Different integrators are constructed depending on the choice of discrete Lagrangian used to approximate the action. By combining discrete Lagrangians in an explicit symplectic partitioned Runge–Kutta method, an integrator capable of achieving any order of accuracy is obtained. Using the von Neumann stability analysis, we show that the integrators greatly increase the numerical stability and reduce the numerical dispersion compared to other methods. For practical purposes, we demonstrate how to implement the integrators using many features of the finite-difference time-domain method. However, our approach is also applicable to other spatial discretizations, such as those used in finite element methods. Using this implementation, numerical examples are presented that demonstrate the ability of the integrators to efficiently reduce and maintain a minimal amount of numerical dispersion, particularly when the time-step is less than the stability limit. The integrators are therefore advantageous for modeling large, inhomogeneous computational domains.