JACOBI PSEUDOSPECTRAL METHOD FOR FOURTH ORDER PROBLEMS

In this paper, we investigate Jacobi pseudospectral method for fourth order problems. We establish some basic results on the Jacobi-Gauss-type interpolations in non-uniformly weighted Sobolev spaces, which serve as important tools in analysis of numerical quadratures, and numerical methods of differential and integral equations. Then we propose Jacobi pseudospectral schemes for several singular problems and multiple-dimensional problems of fourth order. Numerical results demonstrate the spectral accuracy of these schemes, and coincide well with theoretical analysis.     

Singly diagonally implicit runge‐kutta method for time‐dependent reaction‐diffusion equation

In this article, an efficient fourth‐order accurate numerical method based on Padé approximation in space and singly diagonally implicit Runge‐Kutta method in time is proposed to solve the time‐dependent one‐dimensional reaction‐diffusion equation. In this scheme, we first approximate the spatial derivative using the second‐order central finite difference then improve it to fourth‐order by applying Padé approximation. A three stage fourth‐order singly diagonally implicit Runge‐Kutta method is then used to solve the resulting system of ordinary differential equations. It is also shown that the scheme is unconditionally stable, and is suitable for stiff problems. Several numerical examples are solved by the scheme and the efficiency and accuracy of the new scheme are compared with two widely used high‐order compact finite difference methods. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1423–1441, 2011     

Numerical experiments on a hybrid WENO5 filter for shock‐capturing

In the present paper, a hybrid filter is introduced for high accurate numerical simulation of shock‐containing flows. The fourth‐order compact finite difference scheme is used for the spatial discretization and the third‐order Runge–Kutta scheme is used for the time integration. After each time‐step, the hybrid filter is applied on the results. The filter is composed of a linear sixth‐order filter and the dissipative part of a fifth‐order weighted essentially nonoscillatory scheme (WENO5). The classic WENO5 scheme and the WENO5 scheme with adaptive order (WENO5‐AO) are used to form the hybrid filter. Using a shock‐detecting sensor, the hybrid filter reduces to the linear sixth‐order filter in smooth regions for damping high frequency waves and reduces to the WENO5 filter at shocks in order to eliminate unwanted oscillations produced by the nondissipative spatial discretization method. The filter performance and accuracy of the results are examined through several test cases including the advection, Euler and Navier–Stokes equations. The results are compared with that of a hybrid second‐order filter and also that of the WENO5 and WENO5‐AO schemes.     

An efficient spectral method and rigorous error analysis based on dimension reduction scheme for fourth order problems

In this paper, we propose an effective spectral method based on dimension reduction scheme for fourth order problems in polar geometric domains. First, the original problem is decomposed into a series of one‐dimensional fourth order problems by polar coordinate transformation and the orthogonal properties of Fourier basis function. Then the weak form and the corresponding discrete scheme of each one‐dimensional fourth order problem are derived by introducing polar conditions and appropriate weighted Sobolev spaces. In addition, we define the projection operators in the weighted Sobolev space and give its approximation properties, and further prove the error estimation of each one‐dimensional fourth order problem. Finally, we provide some numerical examples, and the numerical results show the effectiveness of our algorithm and the correctness of the theoretical results.     

High order compact solution of the one‐space‐dimensional linear hyperbolic equation

In this article, we introduce a high‐order accurate method for solving one‐space dimensional linear hyperbolic equation. We apply a compact finite difference approximation of fourth order for discretizing spatial derivative of linear hyperbolic equation and collocation method for the time component. The main property of this method additional to its high‐order accuracy due to the fourth order discretization of spatial derivative, is its unconditionally stability. In this technique the solution is approximated by a polynomial at each grid point that its coefficients are determined by solving a linear system of equations. Numerical results show that the compact finite difference approximation of fourth order and collocation method produce a very efficient method for solving the one‐space‐dimensional linear hyperbolic equation. We compare the numerical results of this paper with numerical results of (Mohanty, 3 .© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

An exponential high‐order compact ADI method for 3D unsteady convection–diffusion problems

In this article, we develop an exponential high order compact alternating direction implicit (EHOC ADI) method for solving three dimensional (3D) unsteady convection–diffusion equations. The method, which requires only a regular seven‐point 3D stencil similar to that in the standard second‐order methods, is second order accurate in time and fourth‐order accurate in space and unconditionally stable. The resulting EHOC ADI scheme in each alternating direction implicit (ADI) solution step corresponding to a strictly diagonally dominant matrix equation can be solved by the application of the one‐dimensional tridiagonal Thomas algorithm with a considerable saving in computing time. Numerical experiments for three test problems are carried out to demonstrate the performance of the present method and to compare it with the classical Douglas–Gunn ADI method and the Karaa's high‐order compact ADI method. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Legendre spectral method for solving mixed boundary value problems on rectangle

In this paper, we develop a spectral method for mixed inhomogeneous Dirichlet/Neumann/Robin boundary value problems defined on rectangle. Some results on two‐dimensional Legendre approximation in Jacobi‐weighted Sobolev space are established. As examples of applications, spectral schemes are provided for two model problems with mixed inhomogeneous boundary conditions. The spectral accuracy in space of proposed algorithms are proved. Efficient implementations are presented. Numerical results demonstrate their high accuracy and confirm the theoretical analysis well. Copyright © 2016 John Wiley & Sons, Ltd.     

High‐order finite difference schemes for numerical solutions of the generalized Burgers–Huxley equation

In this article, up to tenth‐order finite difference schemes are proposed to solve the generalized Burgers–Huxley equation. The schemes based on high‐order differences are presented using Taylor series expansion. To establish the numerical solutions of the corresponding equation, the high‐order schemes in space and a fourth‐order Runge‐Kutta scheme in time have been combined. Numerical experiments have been conducted to demonstrate the high‐order accuracy of the current algorithms with relatively minimal computational effort. The results showed that use of the present approaches in the simulation is very applicable for the solution of the generalized Burgers–Huxley equation. The current results are also seen to be more accurate than some results given in the literature. The proposed algorithms are seen to be very good alternatives to existing approaches for such physical applications. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1313‐1326, 2011     

A series of high‐order quasi‐compact schemes for space fractional diffusion equations based on the superconvergent approximations for fractional derivatives

Based on the superconvergent approximation at some point (depending on the fractional order α, but not belonging to the mesh points) for Grünwald discretization to fractional derivative, we develop a series of high‐order quasi‐compact schemes for space fractional diffusion equations. Because of the quasi‐compactness of the derived schemes, no points beyond the domain are used for all the high‐order schemes including second‐order, third‐order, fourth‐order, and even higher‐order schemes; moreover, the algebraic equations for all the high‐order schemes have the completely same matrix structure. The stability and convergence analysis for some typical schemes are made; the techniques of treating the fractional derivatives with nonhomogeneous boundaries are introduced; and extensive numerical experiments are performed to confirm the theoretical analysis or verify the convergence orders. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1345–1381, 2015     

Spectral problems with mixed Dirichlet–Neumann boundary conditions: Isospectrality and beyond

In this study, an implicit semi-discrete higher order compact (HOC) scheme, with an averaged time discretization, has been presented for the numerical solution of unsteady two-dimensional (2D) Schrödinger equation. The scheme is second order accurate in time and fourth order accurate in space. The results of numerical experiments are presented, and are compared with analytical solutions and well established numerical results of some other finite difference schemes. In all cases, the present scheme produces highly accurate results with much better computational efficiency.     

Fast direct solver for Poisson equation in a 2D elliptical domain

In this article, we extend our previous work 3 for developing some fast Poisson solvers on 2D polar and spherical geometries to an elliptical domain. Instead of solving the equation in an irregular Cartesian geometry, we formulate the equation in elliptical coordinates. The solver relies on representing the solution as a truncated Fourier series, then solving the differential equations of Fourier coefficients by finite difference discretizations. Using a grid by shifting half mesh away from the pole and incorporating the derived numerical boundary value, the difficulty of coordinate singularity can be elevated easily. Unlike the case of 2D disk domain, the present difference equation for each Fourier mode is coupled with its conjugate mode through the numerical boundary value near the pole; thus, those two modes are solved simultaneously. Both second‐ and fourth‐order accurate schemes for Dirichlet and Neumann problems are presented. In particular, the fourth‐order accuracy can be achieved by a three‐point compact stencil which is in contrast to a five‐point long stencil for the disk case. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 72–81, 2004     

A CLASS OF EXPLICIT RATIONAL SYMPLECTIC INTEGRATORS

In this paper, a class of rational explicit symplectic integrators for one-dimensional oscillatory Hamiltonian problems are presented. These methods are zero-dissipative, and of first algebraic order and high phase-lag order. By means of composition technique, we construct second and fourth order methods with high phase-lag order of this type. Based on our ideas, three applicable explicit symplectic schemes with algebraic order one, two and four are derived, respectively. We report some numerical results to illustrate the good performance of our methods.     

二维抛物型方程的高精度多重网格解法

提出了数值求解二维抛物型方程的一种新的高精度加权平均紧隐格式，利用Fourier分析方法证明了该格式是无条件稳定的，为了克服传统迭代法在求解隐格式是收敛速度慢的缺陷，利用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的精确性和可靠性。     

Minimal dissipation hybrid bicompact schemes for hyperbolic equations

New monotonicity-preserving hybrid schemes are proposed for multidimensional hyperbolic equations. They are convex combinations of high-order accurate central bicompact schemes and upwind schemes of first-order accuracy in time and space. The weighting coefficients in these combinations depend on the local difference between the solutions produced by the high- and low-order accurate schemes at the current space-time point. The bicompact schemes are third-order accurate in time, while having the fourth order of accuracy and the first difference order in space. At every time level, they can be solved by marching in each spatial variable without using spatial splitting. The upwind schemes have minimal dissipation among all monotone schemes constructed on a minimum space-time stencil. The hybrid schemes constructed has been successfully tested as applied to a number of two-dimensional gas dynamics benchmark problems.     

A class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes

In this paper, we developed a class of the fourth order accurate finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes based on the work (Computers & Fluids, 34: 642–663 (2005)) by Qiu and Shu, with Total Variation Diminishing Runge-Kutta time discretization method for the two-dimensional hyperbolic conservation laws. The key idea of HWENO is to evolve both with the solution and its derivative, which allows for using Hermite interpolation in the reconstruction phase, resulting in a more compact stencil at the expense of the additional work. The main difference between this work and the formal one is the procedure to reconstruct the derivative terms. Comparing with the original HWENO schemes of Qiu and Shu, one major advantage of new HWENO schemes is its robust in computation of problem with strong shocks. Extensive numerical experiments are performed to illustrate the capability of the method. Corresponding author This work was partially supported by the National Natural Science Foundation of China (Grant No. 10671097), the European project ADIGMA on the development of innovative solution algorithms for aerodynamic simulations, Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry and the Natural Science Foundation of Jiangsu Province (Grant No. BK2006511)     

Symmetrical weighted essentially non‐oscillatory‐flux limiter schemes for Hamilton–Jacobi equations

In this paper, we propose a new scheme that combines weighted essentially non‐oscillatory (WENO) procedures together with monotone upwind schemes to approximate the viscosity solution of the Hamilton–Jacobi equations. In one‐dimensional (1D) case, first, we obtain an optimum polynomial on a four‐point stencil. This optimum polynomial is third‐order accurate in regions of smoothness. Next, we modify a second‐order ENO polynomial by choosing an additional point inside the stencil in order to obtain the highest accuracy when combined with the Harten–Osher reconstruction‐evolution method limiter. Finally, the optimum polynomial is considered as a symmetric and convex combination of three polynomials with ideal weights. Following the methodology of the classic WENO procedure, then, we calculate the non‐oscillatory weights with the ideal weights. Numerical experiments in 1D and 2D are performed to compare the capability of the hybrid scheme to WENO schemes. Copyright © 2015 John Wiley & Sons, Ltd.     

High order Hybrid central—WENO finite difference scheme for conservation laws

In this article we present a high resolution hybrid central finite difference—WENO scheme for the solution of conservation laws, in particular, those related to shock–turbulence interaction problems. A sixth order central finite difference scheme is conjugated with a fifth order weighted essentially non-oscillatory WENO scheme in a grid-based adaptive way. High order multi-resolution analysis is used to detect the high gradients regions of the numerical solution in order to capture the shocks with the WENO scheme while the smooth regions are computed with the more efficient and accurate central finite difference scheme. The application of high order filtering to mitigate the dispersion error of central finite difference schemes is also discussed. Numerical experiments with the 1D compressible Euler equations are shown.     

Convergence analysis of some high–order time accurate schemes for a finite volume method for second order hyperbolic equations on general nonconforming multidimensional spatial meshes

The present work is an extension of our previous work (Bradji, Numer Methods Partial Differ Equations, to appear) which dealt with error analysis of a finite volume scheme of a first convergence order (both in time and space) for second‐order hyperbolic equations on general nonconforming multidimensional spatial meshes introduced recently in (Eymard et al. IMAJ Numer Anal 30(2010), 1009–1043). We aim in this article to get some higher‐order time accurate schemes for a finite volume method for second‐order hyperbolic equations using the same class of spatial generic meshes stated above. We derive a family of finite volume schemes approximating the wave equation, as a model for second‐order hyperbolic equations, in which the discretization in time is performed using a one‐parameter scheme of the Newmark's method. We prove that the error estimate of these finite volume schemes is of order two (or four) in time and it is of optimal order in space. These error estimates are analyzed in several norms which allow us to derive approximations for the exact solution and its first derivatives whose the convergence order is two (or four) in time and it is optimal in space. We prove in particular, when the discrete flux is calculated using a stabilized discrete gradient, that the convergence order is \begin{align*}k^2+h_\mathcal{D}\end{align*} or \begin{align*}k^4+h_\mathcal{D}\end{align*}, where \begin{align*}h_\mathcal{D}\end{align*} (resp. k) is the mesh size of the spatial (resp. time) discretization. These estimates are valid under the regularity assumption \begin{align*}u\in C^4(\lbrack 0,T\rbrack;C^2(\overline{\Omega}))\end{align*}, when the schemes are second‐order accurate in time, and \begin{align*}u\in C^6(\lbrack 0,T\rbrack;C^2(\overline{\Omega}))\end{align*}, when the schemes are four‐order accurate in time for the exact solution u. The proof of these error estimates is based essentially on a comparison between the finite volume approximate solution and an auxiliary finite volume approximation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A compact finite difference scheme for solving a three‐dimensional heat transport equation in a thin film

Heat transport at the microscale is of vital importance in microtechnology applications. The heat transport equation differs from the traditional heat diffusion equation in having a second‐order derivative of temperature with respect to time and a third‐order mixed derivative of temperature with respect to space and time. In this study, we develop a high‐order compact finite difference scheme for the heat transport equation at the microscale. It is shown by the discrete Fourier analysis method that the scheme is unconditionally stable. Numerical results show that the solution is accurate. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 441–458, 2000     

A high‐order B‐spline collocation scheme for solving a nonhomogeneous time‐fractional diffusion equation

A high‐accuracy numerical approach for a nonhomogeneous time‐fractional diffusion equation with Neumann and Dirichlet boundary conditions is described in this paper. The time‐fractional derivative is described in the sense of Riemann‐Liouville and discretized by the backward Euler scheme. A fourth‐order optimal cubic B‐spline collocation (OCBSC) method is used to discretize the space variable. The stability analysis with respect to time discretization is carried out, and it is shown that the method is unconditionally stable. Convergence analysis of the method is performed. Two numerical examples are considered to demonstrate the performance of the method and validate the theoretical results. It is shown that the proposed method is of order O(Δx⁴ + Δt^{2 ? α}) convergence, where α ∈ (0,1) . Moreover, the impact of fractional‐order derivative on the solution profile is investigated. Numerical results obtained by the present method are compared with those obtained by the method based on standard cubic B‐spline collocation method. The CPU time for present numerical method and the method based on cubic B‐spline collocation method are provided.