New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

A sixth-order dual preserving algorithm for the Camassa-Holm equation

The paper presents a sixth-order numerical algorithm for studying the completely integrable Camassa-Holm (CH) equation. The proposed sixth-order accurate method preserves both the dispersion relation and the Hamiltonians of the CH equation. The CH equation in this study is written as an evolution equation, involving only the first-order spatial derivatives, coupled with the Helmholtz equation. We propose a two-step method that first solves the evolution equation by a sixth-order symplectic Runge-Kutta method and then solves the Helmholtz equation using a three-point sixth-order compact scheme. The first-order derivative terms in the first step are approximated by a sixth-order dispersion-relation-preserving scheme that preserves the physically inherent dispersive nature. The compact Helmholtz solver, on the other hand, allows us to use relatively few nodal points in a stencil, while achieving a higher-order accuracy. The sixth-order symplectic Runge-Kutta time integrator is preferable for an equation that possesses a Hamiltonian structure. We illustrate the ability of the proposed scheme by examining examples involving peakon or peakon-like solutions. We compare the computed solutions with exact solutions or asymptotic predictions. We also demonstrate the ability of the symplectic time integrator to preserve the Hamiltonians. Finally, via a smooth travelling wave problem, we compare the accuracy, elapsed computing time, and rate of convergence among the proposed method, a second-order two-step algorithm, and a completely integrable particle method.     

A new sixth-order scheme for nonlinear equations

In this paper we present a new efficient sixth-order scheme for nonlinear equations. The method is compared to several members of the family of methods developed by Neta (1979) [B. Neta, A sixth-order family of methods for nonlinear equations, Int. J. Comput. Math. 7 (1979) 157-161]. It is shown that the new method is an improvement over this well known scheme.     

Convergence analysis of nonic-spline solutions for special nonlinear sixth-order boundary value problems

In this paper, we use nonic-spline polynomial method for the numerical solution of special nonlinear sixth-order two-point boundary value problems. The main idea is to use the conditions of continuity as discretization equations for the sixth-order boundary value problem. The end conditions are derived for defined spline. A new approach for convergence analysis of the presented method discussed. Some examples are solved to illustrate the applications of method, and to compare the computed results with other existing known methods.     

Fast and Robust Sixth Order Multigrid Computation for 3D Convection Diffusion Equation

We present a sixth-order explicit compact finite difference scheme to solve the three-dimensional (3D) convection-diffusion equation. We first use a multiscale multigrid method to solve the linear systems arising from a 19-point fourth-order discretization scheme to compute the fourth-order solutions on both a coarse grid and a fine grid. Then an operator-based interpolation scheme combined with an extrapolation technique is used to approximate the sixth-order accurate solution on the fine grid. Since the multigrid method using a standard point relaxation smoother may fail to achieve the optimal grid-independent convergence rate for solving convection-diffusion equations with a high Reynolds number, we implement the plane relaxation smoother in the multigrid solver to achieve better grid independency. Supporting numerical results are presented to demonstrate the efficiency and accuracy of the sixth-order compact (SOC) scheme, compared with the previously published fourth-order compact (FOC) scheme.     

New variants of Jarratt’s method with sixth-order convergence

In this paper, by using the two-variable Taylor expansion formula, we introduce some new variants of Jarratt’s method with sixth-order convergence for solving univariate nonlinear equations. The proposed methods contain some recent improvements of Jarratt’s method. Furthermore, a new variant of Jarratt’s method with sixth-order convergence for solving systems of nonlinear equations is proposed only with an additional evaluation for the involved function, and not requiring the computation of new inverse. Numerical comparisons are made to show the performance of the presented methods.     

三维泊松方程基于Richardson外推法的高阶紧致差分方法

基于Richardson外推法提出了数值求解三维泊松方程的高阶紧致差分方法.方法通过利用四阶和六阶紧致差分格式,分别在细网格和粗网格上求解,然后利用Richardson外推技术和算子插值方法,得到三维泊松方程在细网格上的六阶和八阶精度的数值解.数值实验结果验证了该方法的精确性和有效性.     

Septic spline solutions of sixth-order boundary value problems

Septic spline is used for the numerical solution of the sixth-order linear, special case boundary value problem. End conditions for the definition of septic spline are derived, consistent with the sixth-order boundary value problem. The algorithm developed approximates the solution and their higher-order derivatives. The method has also been proved to be second-order convergent. Three examples are considered for the numerical illustrations of the method developed. The method developed in this paper is also compared with that developed in [M. El-Gamel, J.R. Cannon, J. Latour, A.I. Zayed, Sinc-Galerkin method for solving linear sixth order boundary-value problems, Mathematics of Computation 73, 247 (2003) 1325–1343], as well and is observed to be better.     

Two new families of sixth-order methods for solving non-linear equations

In this paper, we developed two new families of sixth-order methods for solving simple roots of non-linear equations. Per iteration these methods require two evaluations of the function and two evaluations of the first-order derivatives, which implies that the efficiency indexes of our methods are 1.565. These methods have more advantages than Newton’s method and other methods with the same convergence order, as shown in the illustration examples. Finally, using the developing methodology described in this paper, two new families of improvements of Jarratt method with sixth-order convergence are derived in a straightforward manner. Notice that Kou’s method in [Jisheng Kou, Yitian Li, An improvement of the Jarratt method, Appl. Math. Comput. 189 (2007) 1816-1821] and Wang’s method in [Xiuhua Wang, Jisheng Kou, Yitian Li, A variant of Jarratt method with sixth-order convergence, Appl. Math. Comput. 204 (2008) 14-19] are the special cases of the new improvements.     

Higher order difference methods for second order two-point boundary-value problems

Using methods of numerical integration, difference schemes of sixth order with off-step points have been obtained and applied to the second order differential equation, with and without mixed boundary conditions. Numerical results obtained by these methods have been compared with those obtained by using h⁴-extrapolation of the Numerov method. It is found that the sixth-order method based upon 4-point Lobatto quadrature with approximation 2 is more economical and accurate from the computational viewpoint than the existing sixth-order methods.     

A Sixth-order Tridiagonal Finite Difference Method for General Non-linear Two-point Boundary Value Problems

We present a sixth-order finite difference method for the generalsecond-order non-linear differential equation Y"=f(x, y, y')subject to the boundary conditions y(a) = A, y(b) = B. In thecase of linear differential equations, our finite differencescheme leads to tridiagonal linear systems. We establish, underappropriate conditions, O(h⁶)-convergence of the finite differencescheme. Numerical examples are given to illustrate the methodand its sixth-order convergence.     

Nested Runge-Kutta methods for the differential equation y”=F(x,y)

Nested Runge-Kutta methods of fourth, fifth, and sixth orders of accuracy are constructed for numerical solution of the Cauchy problem for the ordinary differential equation y=f(x, y). The first two stages of the fifth- and sixth-order methods determine the fourth-order method, and the ftrst three stages of the sixth-order method determine the fifthorder method. Numerical values of the parameters for these methods are determined and two test problems are solved numerically.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 57, pp. 42–48, 1985.     

A sixth order diagonally implicit symmetric and symplectic Runge-Kutta method for solving Hamiltonian systems

The paper is concerned with construction of symmetric and symplectic Runge-Kutta methods for Hamiltonian systems. Based on the symplectic and symmetrical properties, a sixth-order diagonally implicit symmetric and symplectic Runge-Kutta method with seven stages is presented, the proposed method proved to be P-stable. Numerical experiments with some Hamiltonian oscillatory problems are presented to show the proposed method is as competitive as the existing Runge-Kutta methods in scientic literature.     

Modified Householder iterative method for nonlinear equations

In this paper, we suggest and analyze a new two-step predictor–corrector type iterative method for solving nonlinear equations of the type f(x)=0. This new method includes the two-step Newton method as a special case. We show that this new two-step method is a sixth-order convergent method. Several examples are given to illustrate the efficiency of this new method and its comparison with other sixth-order methods. This method can be considered as a significant improvement of the Newton method and its variant forms.     

Blow-up and extinction for a sixth-order parabolic equation

In this paper, we consider an initial-boundary value problem for a sixth-order parabolic equation. We use the modified method of potential wells to study the relationship which the equation solutions existence, blow-up and the asymptotic behavior with initial conditions.     

N-soliton solutions for the integrable bidirectional sixth-order Sawada-Kotera equation

In this work, the integrable bidirectional sixth-order Sawada-Kotera equation is examined. The equation considered is a KdV6 equation that was derived from the fifth order Sawada-Kotera equation. Multiple soliton solutions and multiple singular soliton solutions are formally derived for this equation. The Cole-Hopf transformation method combined with the Hirota’s bilinear method are used to determine the two sets of solutions, where each set has a distinct structure.     

New exact solutions for a generalization of the Korteweg-de Vries equation (KdV6)

The generalized tanh-coth method is used to construct periodic and soliton solutions for a new integrable system, which has been derived from an integrable sixth-order nonlinear wave equation (KdV6). The system is formed by two equations. One of the equations may be considered as a Korteweg-de Vries equation with a source and the second equation is a third-order linear differential equation.     

Semilocal convergence of a sixth-order Jarratt method in Banach spaces

In this paper, we study the semilocal convergence for a sixth-order variant of the Jarratt method for solving nonlinear equations in Banach spaces. The semilocal convergence of this method is established by using recurrence relations. We derive the recurrence relations for the method, and then prove an existence-uniqueness theorem, along with a priori error bounds which demonstrates the R-order of the method. Finally, we give some numerical applications to demonstrate our approach.     

Sinc-Galerkin method for solving linear sixth-order boundary-value problems

There are few techniques available to numerically solve sixth-order boundary-value problems with two-point boundary conditions. In this paper we show that the Sinc-Galerkin method is a very effective tool in numerically solving such problems. The method is then tested on examples with homogeneous and nonhomogeneous boundary conditions and a comparison with the modified decomposition method is made. It is shown that the Sinc-Galerkin method yields better results.

     

Multiple soliton solutions for the sixth-order Ramani equation and a coupled Ramani equation

In this work, the completely integrable sixth-order nonlinear Ramani equation and a coupled Ramani equation are studied. Multiple soliton solutions and multiple singular soliton solutions are formally derived for these two equations. The Hirota’s bilinear method is used to determine the two distinct structures of solutions. The resonance relations for the three cases are investigated.