A cell-centered lagrangian scheme in two-dimensional cylindrical geometry

A new Lagrangian cell-centered scheme for two-dimensional compressible flows in planar geometry is proposed by Maire et al.The main new feature of the algorithm is that the vertex velocities and the numerical fluxes through the cell interfaces are all evaluated in a coherent manner contrary to standard approaches.In this paper the method introduced by Maire et al.is extended for the equations of Lagrangian gas dynamics in cylindrical symmetry.Two different schemes are proposed,whose difference is that one uses volume weighting and the other area weighting in the discretization of the momentum equation.In the both schemes the conservation of total energy is ensured,and the nodal solver is adopted which has the same formulation as that in Cartesian coordinates.The volume weighting scheme preserves the momentum conservation and the area-weighting scheme preserves spherical symmetry.The numerical examples demonstrate our theoretical considerations and the robustness of the new method.     

Iterative methods and high-order difference schemes for 2D elliptic problems with mixed derivative

We derive a fourth-order compact finite difference scheme for a two-dimensional elliptic problem with a mixed derivative and constant coefficients. We conduct experimental study on numerical solution of the problem discretized by the present compact scheme and the traditional second-order central difference scheme. We study the computed accuracy achieved by each scheme and the performance of the Gauss-Seidel iterative method, the preconditioned GMRES iterative method, and the multigrid method, for solving linear systems arising from the difference schemes.     

同位非结构和结构网格摄动有限体积法（PFV）求解二维Navier-Stokes方程组

根据NS方程组的一阶迎风和二阶中心有限体积（UFV和CFV）格式,导出NS方程组迎风和中心摄动有限体积(UPFV和CPFV)格式．该格式通过把控制体界面质量通量摄动展开成网格间距的幂级数,并由守恒方程本身求得幂级数系数而获得．迎风和中心摄动有限体积格式使用了与一阶迎风和二阶中心格式相同的基点数和相同的表达形式,宜于计算机编程．顶盖驱动方腔流和驻点流标量输运的数值实验证明,迎风PFV格式比一阶UFV、二阶CFV格式有更高的精度,更高的分辨率．尤其是良好的鲁棒特性．对顶盖驱动方腔流,在Re数从102到104范围内,亚松弛系数可在0.3～0.8任取,收敛性能良好．     

A formally fourth-order accurate compact scheme for 3D Poisson equation in cylindrical coordinates

In this paper, we extend our previous work (M.-C. Lai, A simple compact fourth-order Poisson solver on polar geometry, J. Comput. Phys. 182 (2002) 337–345) to 3D cases. More precisely, we present a spectral/finite difference scheme for Poisson equation in cylindrical coordinates. The scheme relies on the truncated Fourier series expansion, where the partial differential equations of Fourier coefficients are solved by a formally fourth-order accurate compact difference discretization. Here the formal fourth-order accuracy means that the scheme is exactly fourth-order accurate while the poles are excluded and is third-order accurate otherwise. Despite the degradation of one order of accuracy due to the presence of poles, the scheme handles the poles naturally; thus, no pole condition is needed. The resulting linear system is then solved by the Bi-CGSTAB method with the preconditioner arising from the second-order discretization which shows the scalability with the problem size.     

MUSTA: A multi-stage numerical flux

Numerical methods for conservation laws constructed in the framework of finite volume and discontinuous Galerkin finite elements require, as the building block, a monotone numerical flux. In this paper we present some preliminary results on the MUSTA approach [E.F. Toro, Multi-stage predictor–corrector fluxes for hyperbolic equations, Technical Report NI03037-NPA, Isaac Newton Institute for Mathematical Sciences, University of Cambridge, UK, 17th June, 2003] for constructing upwind numerical fluxes. The scheme may be interpreted as an un-conventional approximate Riemann solver that has simplicity and generality as its main features. When used in its first-order mode we observe that the scheme achieves the accuracy of the Godunov method used in conjunction with the exact Riemann solver, which is the reference first-order method for hyperbolic systems. At least for the scalar model hyperbolic equation, the Godunov scheme is the best of all first-order monote schemes, it has the smallest truncation error. Extensions of the scheme of this paper are realized in the framework of existing approaches. Here we present a second-order TVD (TVD for the scalar case) extension and show numerical results for the two-dimensional Euler equations on non-Cartesian geometries. The schemes find their best justification when solving very complex systems for which the solution of the Riemann problem, in the classical sense, is too complex, too costly or is simply unavailable.     

Flux-splitting schemes for parabolic problems

Splitting with respect to space variables can be used in solving boundary value problems for second-order parabolic equations. Classical alternating direction methods and locally one-dimensional schemes could be examples of this approach. For problems with rapidly varying coefficients, a convenient tool is the use of fluxes (directional derivatives) as independent variables. The original equation is written as a system in which not only the desired solution but also directional derivatives (fluxes) are unknowns. In this paper, locally one-dimensional additional schemes (splitting schemes) for second-order parabolic equations are examined. By writing the original equation in flux variables, certain two-level locally one-dimensional schemes are derived. The unconditional stability of locally one-dimensional flux schemes of the first and second approximation order with respect to time is proved.     

Implicit Solution of a Two-Dimensional Parabolic Inverse Problem with Temperature Overspecification

Three different implicit finite difference schemes for solving the two-dimensional parabolic inverse problem with temperature overspecification are considered. These schemes are developed for indentifying the control parameter which produces, at any given time, a desired temperature distribution at a given point in the spatial domain. The numerical methods discussed, are based on the second-order (5,1) Backward Time Centered Space (BTCS) implicit formula, and the second-order (5,5) Crank-Nicolson implicit finite difference formula and the fourth-order (9,9) implicit scheme. These finite difference schemes are unconditionally stable. The (9,9) implicit formula takes a huge amount of CPU time, but its fourth-order accuracy is significant. The results of a numerical experiment are presented, and the accuracy and central processor (CPU) times needed for each of the methods are discussed and compared. The implicit finite difference schemes use more central processor times than the explicit finite difference techniques, but they are stable for every diffusion number.     

High-order ADI orthogonal spline collocation method for a new 2D fractional integro-differential problem

We use the generalized L1 approximation for the Caputo fractional derivative, the second-order fractional quadrature rule approximation for the integral term, and a classical Crank-Nicolson alternating direction implicit (ADI) scheme for the time discretization of a new two-dimensional (2D) fractional integro-differential equation, in combination with a space discretization by an arbitrary-order orthogonal spline collocation (OSC) method. The stability of a Crank-Nicolson ADI OSC scheme is rigourously established, and error estimate is also derived. Finally, some numerical tests are given.     

Comparative study of 1D, 2D and 3D simplified gas kinetic schemes for simulation of inviscid compressible flows

With assumption that all the particles in the phase velocity space are concentrated on a circle and on a sphere, the circular function-based gas kinetic scheme and sphere function-based gas kinetic scheme have been developed by Shu and his coworkers [21], [22], [23]. These schemes are simpler than the Maxwellian function-based gas kinetic schemes. The simplicity is due to the fact that the integral domain of phase velocity of circular function and sphere function is a finite region while the integral domain of Maxwellian distribution function is infinite. In this work, the 1D delta function-based gas kinetic scheme is also developed to form a complete set of the simplified gas kinetic schemes. The 1D, 2D and 3D simplified gas kinetic schemes can be viewed as the truly 1D, 2D and 3D flux solvers since they are based on the multi-dimensional Boltzmann equation. On the other hand, to solve the 3D flow problem, the tangential velocities are needed to be approximated by some ways for the 1D and 2D simplified gas kinetic schemes, and to solve the 1D flow problem, the tangential velocities should be taken as zero for the 2D and 3D simplified gas kinetic schemes. The performances of these three schemes for simulation of inviscid compressible flows are investigated in this work by their application to solve the test problems from 1D to 3D cases. Numerical results showed that the efficiency of the delta function-based gas kinetic scheme is slightly superior to that of the circular function- and sphere function-based gas kinetic schemes, while its stability is inferior significantly to the latter. For simulation of the 3D hypersonic flows, the sphere function-based gas kinetic scheme could be the best choice.     

High-order methods for elliptic equations with variable coefficients

In this article, we give a simple method for developing finite difference schemes on a uniform square gird. We consider a general, two-dimensional, second-order, partial differential equation with variable coefficients. In the case of a nine-point scheme, we obtain the known results of Young and Dauwalder in a fairly elegant fashion. We show how this can be extended to obtain fourth-order schemes on thirteen points. We derive two such schemes which are attractive because they can be adapted quite easily bnto obtain formulas for gird points near the boundary. In addition to this, these formulas only require nine evaluations for the typical forcing function. Numerical examples are given to demonstrate the performance of one of the fourth-order schemes.     

Numerical modelling of free-surface shallow flows over irregular topography with complex geometry

A well-balanced Godunov-type finite volume algorithm is developed for modelling free-surface shallow flows over irregular topography with complex geometry. The algorithm is based on a new formulation of the classical shallow water equations in hyperbolic conservation form. Unstructured triangular grids are used to achieve the adaptability of the grid to the geometry of the problem and to facilitate localised refinement. The numerical fluxes are calculated using HLLC approximate Riemann solver, and the MUSCL-Hancock predictor–corrector scheme is adopted to achieve the second-order accuracy both in space and in time where the solutions are continuous, and to achieve high-resolution results where the solutions are discontinuous. The novelties of the algorithm include preserving well-balanced property without any additional correction terms and the wet/dry front treatments. The good performance of the algorithm is demonstrated by comparing numerical and theoretical results of several benchmark problems, including the preservation of still water over a two-dimensional hump, the idealised dam-break flow over a frictionless flat rectangular channel, the circular dam-break, and the shock wave from oblique wall. Besides, two laboratory dam-break cases are used for model validation. Furthermore, a practical application related to dam-break flood wave propagation over highly irregular topography with complex geometry is presented. The results show that the algorithm can correctly account for free-surface shallow flows with respect to its effectiveness and robustness thus has bright application prospects.     

保持球对称的相容流体算法的构造

重点研究如何在二维柱坐标系中保持流体运动的球对称问题．该问题需要对整个流体算法做细致的分析与设计,从而获得不同的方法．基于Caramana介绍的修正压强梯度算子构造了保持球对称的相容流体算法,该算法还包括了边人工粘性方法、子网格压力的计算方法,以及预估步与校正步相结合的策略．并且,考虑了总能量的守恒性,讨论了网格节点质量保持不变的必要性,以及基于相容格式证明了边人工粘性的耗散性．数值模拟结果与已知数据的比对验证了程序的正确性、稳定性.     

Mesh-centered finite differences from nodal finite elements for elliptic problems

After it is shown that the classical five-point mesh-centered finite difference scheme can be derived from a low-order nodal finite element scheme by using nonstandard quadrature formulae, higher-order block mesh-centered finite difference schemes for second-order elliptic problems are derived from higher-order nodal finite elements with nonstandard quadrature formulae as before, combined to a procedure known as “transverse integration.” Numerical experiments with uniform and nonuniform meshes and different types of boundary conditions confirm the theoretical predictions, in discrete as well as continuous norms. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 439–465, 1998     

Finite-element method for time-dependent euler equation

This paper presents finite element methods to approximate inviscid incompressible flow problems. First we emphasize the conservation properties of these problems, and we show that finite element methods appear as a very natural way to find conservative schemes such as Arakawa's scheme. We give convergence theorems and an error analysis of finite element discretization schemes. We turn then to the time differencing problem. We derive stability and convergence results for a second-order semi-implicit scheme and for the leap-frog scheme.     

Fast second-order accurate difference schemes for time distributed-order and Riesz space fractional diffusion equations

The aim of this paper is to develop fast second-order accurate difference schemes for solving one- and two-dimensional time distributed-order and Riesz space fractional diffusion equations. We adopt the same measures for one- and two-dimensional problems as follows: we first transform the time distributed-order fractional diffusion problem into the multi-term time-space fractional diffusion problem with the composite trapezoid formula. Then, we propose a second-order accurate difference scheme based on the interpolation approximation on a special point to solve the resultant problem. Meanwhile, the unconditional stability and convergence of the new difference scheme in $L_2$-norm are proved. Furthermore, we find that the discretizations lead to a series of Toeplitz systems which can be efficiently solved by Krylov subspace methods with suitable circulant preconditioners. Finally, numerical results are presented to show the effectiveness of the proposed difference methods and demonstrate the fast convergence of our preconditioned Krylov subspace methods.     

Development and Comparison of Numerical Fluxes for LWDG Methods

The discontinuous Galerkin （DO） or local discontinuous Galerkin （LDG） method is a spatial discretization procedure for convection-diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters. The Lax- Wendroff time discretization procedure is an altemative method for time discretization to the popular total variation diminishing （TVD） Runge-Kutta time discretizations. In this paper, we develop fluxes for the method of DG with Lax-Wendroff time discretization procedure （LWDG） based on different numerical fluxes for finite volume or finite difference schemes, including the first-order monotone fluxes such as the Lax-Friedfichs flux, Godunov flux, the Engquist-Osher flux etc. and the second-order TVD fluxes. We systematically investigate the performance of the LWDG methods based on these different numerical fluxes for convection terms with the objective of obtaining better performance by choosing suitable numerical fluxes. The detailed numerical study is mainly performed for the one-dimensional system case, addressing the issues of CPU cost, accuracy, non-oscillatory property, and resolution of discontinuities. Numerical tests are also performed for two dimensional systems.     

Quasi-acoustic scheme for the Euler equations of gas dynamics

We suggest an original scheme and an algorithm for the numerical solution of the Euler equations of gas dynamics. The construction of the scheme is based on the mass, momentum, and energy conservation laws. The flux computation is carried out by summation of elementary fluxes formed by small-amplitude running waves that satisfy the linearized equations of gas dynamics. The scheme contains no artificial regularizers, has second-order accuracy on smooth solutions, and is quasimonotone in a neighborhood of the discontinuities. Examples of one- and two-dimensional computations are given.     

High-order accurate implicit running schemes

An approach to the construction of high-order accurate implicit predictor-corrector schemes is proposed. The accuracy is improved by choosing a special time integration step for computing numerical fluxes through cell interfaces by using an unconditionally stable implicit scheme. For smooth solutions of advection equations with constant coefficients, the scheme is second-order accurate. Implicit difference schemes for multidimensional advection equations are constructed on the basis of Godunov’s method with splitting over spatial variables as applied to the computation of “large” values at an intermediate layer. The numerical solutions obtained for advection equations and the radiative transfer equations in a vacuum are compared with their exact solutions. The comparison results confirm that the approach is efficient and that the accuracy of the implicit predictor-corrector schemes is improved.     

Central ADER schemes for hyperbolic conservation laws

In this paper we first briefly review the very high order ADER methods for solving hyperbolic conservation laws. ADER methods use high order polynomial reconstruction of the solution and upwind fluxes as the building block. They use a first order upwind Godunov and the upwind second order weighted average (WAF) fluxes. As well known the upwind methods are more accurate than central schemes. However, the superior accuracy of the ADER upwind schemes comes at a cost, one must solve exactly or approximately the Riemann problems (RP). Conventional Riemann solvers are usually complex and are not available for many hyperbolic problems of practical interest. In this paper we propose to use two central fluxes, instead of upwind fluxes, as the building block in ADER scheme. These are the monotone first order Lax-Friedrich (LXF) and the third order TVD flux. The resulting schemes are called central ADER schemes. Accuracy of the new schemes is established. Numerical implementations of the new schemes are carried out on the scalar conservation laws with a linear flux, nonlinear convex flux and non-convex flux. The results demonstrate that the proposed scheme, with LXF flux, is comparable to those using first and second order upwind fluxes while the scheme, with third order TVD flux, is superior to those using upwind fluxes. When compared with the state of art ADER schemes, our central ADER schemes are faster, more accurate, Riemann solver free, very simple to implement and need less computer memory. A way to extend these schemes to general systems of nonlinear hyperbolic conservation laws in one and two dimensions is presented.     

Galerkin-Legendre spectral schemes for nonlinear space fractional Schrödinger equation

In the paper, we first propose a Crank-Nicolson Galerkin-Legendre (CN-GL) spectral scheme for the one-dimensional nonlinear space fractional Schrödinger equation. Convergence with spectral accuracy is proved for the spectral approximation. Further, a Crank-Nicolson ADI Galerkin-Legendre spectral method for the two-dimensional nonlinear space fractional Schrödinger equation is developed. The proposed schemes are shown to be efficient with second-order accuracy in time and spectral accuracy in space which are higher than some recently studied methods. Moreover, some numerical results are demonstrated to justify the theoretical analysis.