A CCD–ADI method for unsteady convection–diffusion equations

With a combined compact difference scheme for the spatial discretization and the Crank–Nicolson scheme for the temporal discretization, respectively, a high-order alternating direction implicit method (ADI) is proposed for solving unsteady two dimensional convection–diffusion equations. The method is sixth-order accurate in space and second-order accurate in time. The resulting matrix at each ADI computation step corresponds to a triple-tridiagonal system which can be effectively solved with a considerable saving in computing time. In practice, Richardson extrapolation is exploited to increase the temporal accuracy. The unconditional stability is proved by means of Fourier analysis for two dimensional convection–diffusion problems with periodic boundary conditions. Numerical experiments are conducted to demonstrate the efficiency of the proposed method. Moreover, the present method preserves the higher order accuracy for convection-dominated problems.     

High-order compact ADI (HOC-ADI) method for solving unsteady 2D Schrödinger equation

In this paper, a high-order compact (HOC) alternating direction implicit (ADI) method is proposed for the solution of the unsteady two-dimensional Schrödinger equation. The present method uses the fourth-order Padé compact difference approximation for the spatial discretization and the Crank-Nicolson scheme for the temporal discretization. The proposed HOC-ADI method has fourth-order accuracy in space and second-order accuracy in time. The resulting scheme in each ADI computation step corresponds to a tridiagonal system which can be solved by using the one-dimensional tridiagonal algorithm with a considerable saving in computing time. Numerical experiments are conducted to demonstrate its efficiency and accuracy and to compare it with analytic solutions and numerical results established by some other methods in the literature. The results show that the present HOC-ADI scheme gives highly accurate results with much better computational efficiency.     

A 15-point high-order compact scheme with multigrid computation for solving 3D convection diffusion equations

We propose a method with sixth-order accuracy to solve the three-dimensional (3D) convection diffusion equation. We first use a 15-point fourth-order compact discretization scheme to obtain fourth-order solutions on both fine and coarse grids using the multigrid method. Then an iterative mesh refinement technique combined with Richardson extrapolation is used to approximate the sixth-order accurate solution on the fine grid. Numerical results are presented for a variety of test cases to demonstrate the efficiency and accuracy of the proposed method, compared with the standard fourth-order compact scheme.     

High-order discontinuous Galerkin computation of axisymmetric transonic flows in safety relief valves

This paper presents a discontinuous Galerkin (DG) discretization of the compressible RANS and k–ω turbulence model equations for two-dimensional axisymmetric flows. The developed code has been applied to investigate the transonic flow in safety relief valves.This new DG implementation has evolved from the DG method presented in [1]. An “exact” Riemann solver is used to compute the interface numerical inviscid flux while the viscous flux discterization relies on the BRMPS scheme [2] and [3]. Control of oscillations of high-order solutions around shocks is obtained by means of a shock-capturing technique developed and assessed within the EU ADIGMA project [4].The code has been applied to compute the flow in a spring loaded safety valve at several back pressures and different disk lifts. The predicted device flow capacity and the pressure inside its bonnet have been checked against experimental data. The CFD simulations allow to clarify the complex flow patterns occurring and to explain the measured trends.     

Fourth-order central compact scheme for the numerical solution of incompressible Navier–Stokes equations

This paper provides an implicit central compact scheme for the numerical solution of incompressible Navier–Stokes equations. The solution procedure is based on the artificial compressibility method that transforms the governing equations into a hyperbolic-parabolic form. A fourth-order central compact scheme with a sixth-order numerical filtering is used for the discretization of convective terms and fourth-order central compact scheme for the viscous terms. Dual-time stepping approach is applied to time discretization with backward Euler difference scheme to the pseudo-time derivative, and three point second-order backward difference scheme to the physical time derivative. An approximate factorization-based alternating direction implicit scheme is used to solve the resulting block tridiagonal system of equations. The accuracy and efficiency of the proposed numerical method is verified by simulating several two-dimensional steady and unsteady benchmark problems.     

A spatially second order alternating direction implicit (ADI) method for solving three dimensional parabolic interface problems

A new matched alternating direction implicit (ADI) method is proposed in this paper for solving three-dimensional (3D) parabolic interface problems with discontinuous jumps and complex interfaces. This scheme inherits the merits of its ancestor for two-dimensional problems, while possesses several novel features, such as a non-orthogonal local coordinate system for decoupling the jump conditions, two-side estimation of tangential derivatives at an interface point, and a new Douglas–Rachford ADI formulation that minimizes the number of perturbation terms, to attack more challenging 3D problems. In time discretization, this new ADI method is found to be of first order and stable in various experiments. In space discretization, the matched ADI method achieves the second order accuracy based on simple Cartesian grids for various irregularly-shaped surfaces and spatial–temporal dependent jumps. Computationally, the matched ADI method is as efficient as the fastest implicit scheme based on the geometrical multigrid for solving 3D parabolic equations, in the sense that its complexity in each time step scales linearly with respect to the spatial degree of freedom $N$, i.e., $O\left(N\right)$. Furthermore, unlike iterative methods, the ADI method is an exact or non-iterative algebraic solver which guarantees to stop after a certain number of computations for a fixed $N$. Therefore, the proposed matched ADI method provides a very promising tool for solving 3D parabolic interface problems.     

High-order compact boundary value method for the solution of unsteady convection–diffusion problems

In this paper, we propose a new class of high-order accurate methods for solving the two-dimensional unsteady convection–diffusion equation. These techniques are based on the method of lines approach. We apply a compact finite difference approximation of fourth order for discretizing spatial derivatives and a boundary value method of fourth order for the time integration of the resulted linear system of ordinary differential equations. The proposed method has fourth-order accuracy in both space and time variables. Also this method is unconditionally stable due to the favorable stability property of boundary value methods. Numerical results obtained from solving several problems include problems encounter in many transport phenomena, problems with Gaussian pulse initial condition and problems with sharp discontinuity near the boundary, show that the compact finite difference approximation of fourth order and a boundary value method of fourth order give an efficient algorithm for solving such problems.     

Implicit methods for fluid dynamics

This paper represents the contributions to the development of implicit procedures for solving the equations of fluid dynamics made by Briley and McDonald (1975)[1], Beam and Warming (1976, 1978) [2] and [3] and Lombard et al. (1983) [4]. The contributions of Briley and McDonald and Beam and Warming are well known, but Lombard has not been fully recognized for his innovative contributions to flux vector splitting and use of the DDADI (Diagonally Dominate Alternating Direction Implicit) algorithm. Their contributions are presented herein.Fully implicit algorithms are applied to two complex flow problems of current interest, (1) hypersonic non-equilibrium flow about a blunt nosed body and (2) flow within an MFD (magneto-fluid dynamics) accelerator. These two applications would be exceeding costly without the use of fully implicit methods.     

Cooperative strategies of wireless access technologies: A game-theoretic analysis

The wireless Internet has to overcome the problem of spectrum scarcity as the number of mobile equipments could increase even by an order of magnitude in the next decade; the cooperation of mobile devices is foreseeable as a feasible solution to the problems. There exists a large body of literature on opportunistic ad hoc networking including Pelusi et al. (2006) [25], Chen et al. (2006) [26], Hui et al. (2005) [27]; however, the impact of the location of the devices on their access method selection is not yet appropriately dealt with. In this paper, we address this issue based on game-theoretic analyses. The key contribution of our work is threefold. First, we model the access method selection of mobile devices by extending the classical forwarding game with position, mobility, and availability of the devices. Second, we apply the model in game-theoretic analyses to better understand the optimal cooperation strategies in the presence of heterogeneous wireless technologies. We further extend our framework to include uncertainty. Finally, we present the applicability of the model in a cognitive radio scenario where complex structures of parameters are included.     

基于非等距网格高阶紧致差分格式的多重网格算法研究

本文结合非等距网格高精度紧致差分格式的优越性与多重网格方法的快速收敛性,求解二维对流扩散方程。研究结果表明,对于处理物理量在不同的空间方向呈现不同的性态特征或不同变化规律的物理问题时,用非等距网格离散的四阶紧致格式的多重网格算法和二阶中心差分格式的多重网格算法都比等距网格离散得高效。同时,在非等距网格下下,部分半粗化多重网格算法比完全粗化多重网格算法具有更高的计算效率。针对不同的松弛算子对误差残量的磨光效果比较研究表明,线松弛算子是最高效的。而且,非等距网格离散的高精度紧致格式的多重网格算法对于对流扩散问题中大网格雷诺数情形也是收敛的。     

一维对流扩散反应方程的高精度数值方法

通过将原方程变换为对流扩散方程,将所得方程的对流项采用四阶组合紧致迎风格式离散,扩散项采用四阶对称紧致格式离散之后,对得到的空间半离散格式采用四阶龙格库塔方法进行时间推进,得到了一种求解非定常对流扩散反应问题的高精度方法,其收敛阶为O(h4+τ4).经数值实验并与文献结果进行对比,表明该格式适用于对流占优问题的数值模拟,验证了格式的良好性能.     

A new approach to Lipschitz continuity in state constrained optimal control

For a linear-quadratic state constrained optimal control problem, it is proved in [11] that under an independence condition for the active constraints, the optimal control is Lipschitz continuous. We now give a new proof of this result based on an analysis of the Euler discretization given in [9]. There we exploit the Lipschitz continuity of the control to estimate the error in the Euler discretization. Here we show that the theory developed for the Euler discretization can be used to derive the Lipschitz continuity of the optimal control.     

An Extrapolation Cascadic Multigrid Method Combined with a Fourth-Order Compact Scheme for 3D Poisson Equation

Extrapolation cascadic multigrid (EXCMG) method is an efficient multigrid method which has mainly been used for solving the two-dimensional elliptic boundary value problems with linear finite element discretization in the existing literature. In this paper, we develop an EXCMG method to solve the three-dimensional Poisson equation on rectangular domains by using the compact finite difference (FD) method with unequal meshsizes in different coordinate directions. The resulting linear system from compact FD discretization is solved by the conjugate gradient (CG) method with a relative residual stopping criterion. By combining the Richardson extrapolation and tri-quartic Lagrange interpolation for the numerical solutions from two-level of grids (current and previous grids), we are able to produce an extremely accurate approximation of the actual numerical solution on the next finer grid, which can greatly reduce the number of relaxation sweeps needed. Additionally, a simple method based on the midpoint extrapolation formula is used for the fourth-order FD solutions on two-level of grids to achieve sixth-order accuracy on the entire fine grid cheaply and directly. The gradient of the numerical solution can also be easily obtained through solving a series of tridiagonal linear systems resulting from the fourth-order compact FD discretizations. Numerical results show that our EXCMG method is much more efficient than the classical V-cycle and W-cycle multigrid methods. Moreover, only few CG iterations are required on the finest grid to achieve full fourth-order accuracy in both the \(L^2\)-norm and \(L^{\infty }\)-norm for the solution and its gradient when the exact solution belongs to \(C^6\). Finally, numerical result shows that our EXCMG method is still effective when the exact solution has a lower regularity, which widens the scope of applicability of our EXCMG method.     

A high-order compact finite difference method and its extrapolation for fractional mobile/immobile convection–diffusion equations

This paper is concerned with high-order numerical methods for a class of fractional mobile/immobile convection–diffusion equations. The convection coefficient of the equation may be spatially variable. In order to overcome the difficulty caused by variable coefficient problems, we first transform the original equation into a special and equivalent form, which is then discretized by a fourth-order compact finite difference approximation for the spatial derivative and a second-order difference approximation for the time first derivative and the Caputo time fractional derivative. The local truncation error and the solvability of the resulting scheme are discussed in detail. The (almost) unconditional stability and convergence of the method are proved using a discrete energy analysis method. A Richardson extrapolation algorithm is presented to enhance the temporal accuracy of the computed solution from the second-order to the third-order. Applications using two model problems give numerical results that demonstrate the accuracy of the new method and the high efficiency of the Richardson extrapolation algorithm.     

Unconditionally stable ADI scheme of higher-order for linear hyperbolic equations

An unconditionally stable alternating direction implicit (ADI) method of higher-order in space is proposed for solving two- and three-dimensional linear hyperbolic equations. The method is fourth-order in space and second-order in time. The solution procedure consists of a multiple use of one-dimensional matrix solver which produces a computational cost effective solver. Numerical experiments are conducted to compare the new scheme with the existing scheme based on second-order spatial discretization. The effectiveness of the new scheme is exhibited from the numerical results.     

Confronting intractability via parameters

One approach to confronting computational hardness is to try to understand the contribution of various parameters to the running time of algorithms and the complexity of computational tasks. Almost no computational tasks in real life are specified by their size alone. It is not hard to imagine that some parameters contribute more intractability than others and it seems reasonable to develop a theory of computational complexity which seeks to exploit this fact. Such a theory should be able to address the needs of practitioners in algorithmics. The last twenty years have seen the development of such a theory. This theory has a large number of successes in terms of a rich collection of algorithmic techniques, both practical and theoretical, and a fine-grained intractability theory. Whilst the theory has been widely used in a number of areas of applications including computational biology, linguistics, VLSI design, learning theory and many others, knowledge of the area is highly varied. We hope that this article will show the basic theory and point at the wide array of techniques available. Naturally the treatment is condensed, and the reader who wants more should go to the texts of Downey and Fellows (1999) [2], Flum and Grohe (2006) [59], Niedermeier (2006) [28], and the upcoming undergraduate text (Downey and Fellows 2012) [278].     

A higher-order compact ADI method with monotone iterative procedure for systems of reaction-diffusion equations

This paper is concerned with an existing compact finite difference ADI method, published in the paper by Liao et al. (2002) [3], for solving systems of two-dimensional reaction-diffusion equations with nonlinear reaction terms. This method has an accuracy of fourth-order in space and second-order in time. The existence and uniqueness of its solution are investigated by the method of upper and lower solutions, without any monotone requirement on the nonlinear reaction terms. The convergence of the finite difference solution to the continuous solution is proved. An efficient monotone iterative algorithm is presented for solving the resulting discrete system, and some techniques for the construction of upper and lower solutions are discussed. An application using a model problem gives numerical results that demonstrate the high efficiency and advantages of the method.     

Numerical simulation of the three dimensional Allen–Cahn equation by the high-order compact ADI method

In this paper, a new linearized high-order compact difference method is presented for numerical simulation of three dimensional (3D) Allen–Cahn equation with three kinds of boundary conditions. The method, which is based on the Crank–Nicholson/Adams–Bashforth scheme combined with the Douglas–Gunn ADI method, is second order accurate in time and fourth order accurate in space and energy degradation. The main advantages of this method is that the nonlinear penalty term f(u)$f\left(u\right)$ is linear and an extra stabilizing term is added to alleviate the stability constraint while maintaining accuracy and simplicity. Numerical experiments are given to demonstrate the validity and applicability of the new method.     

Fast algorithms for high-order numerical methods for space-fractional diffusion equations

In this paper, fast numerical methods for solving space-fractional diffusion equations are studied in two stages. Firstly, a fast direct solver for an implicit finite difference scheme proposed by Hao et al. [A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805], which is fourth-order accurate in space and second-order accurate in time, is developed based on a circulant-and-skew-circulant (CS) representation of Toeplitz matrix inversion. Secondly, boundary value method with spatial discretization of Hao et al. [A fourth-order approximation of fractional derivatives with its applications, J. Comput. Phys. 281 (2015), pp. 787–805] is adopted to produce a numerical solution with higher order accuracy in time. Particularly, a method with fourth-order accuracy in both space and time can be achieved. GMRES method is employed for solving the discretized linear system with two preconditioners. Based on the CS representation of Toeplitz matrix inversion, the two preconditioners can be applied efficiently, and the convergence rate of the preconditioned GMRES method is proven to be fast. Numerical examples are given to support the theoretical analysis.     

High-resolution compact numerical method for the system of 2D quasi-linear elliptic boundary value problems and the solution of normal derivatives on an irrational domain with engineering applications

In this paper, we present a novel approach to attain fourth-order approximate solution of 2D quasi-linear elliptic partial differential equation on an irrational domain. In this approach, we use nine grid points with dissimilar mesh in a single compact cell. We also discuss appropriate fourth-order numerical methods for the solution of the normal derivatives on a dissimilar mesh. The method has been protracted for solving system of quasi-linear elliptic equations. The convergence analysis is discussed to authenticate the proposed numerical approximation. On engineering applications, we solve various test problems, such as linear convection–diffusion equation, Burgers’equation, Poisson equation in singular form, NS equations, bi- and tri-harmonic equations and quasi-linear elliptic equations to show the efficiency and accuracy of the proposed methods. A comprehensive comparative computational experiment shows the accuracy, reliability and credibility of the proposed computational approach.