High Order of Accuracy for Poisson Equation Obtained by Grouping of Repeated Richardson Extrapolation with Fourth Order Schemes

In this article,we improve the order of precision of the two-dimensional Poisson equation by combining extrapolation techniques with high order schemes.The high order solutions obtained traditionally generate non-sparse matrices and the calculation time is very high.We can obtain sparse matrices by applying compact schemes.In this article,we compare compact and exponential finite difference schemes of fourth order.The numerical solutions are calculated in quadruple precision(Real*16 or extended precision)in FORTRAN language,and iteratively obtained until reaching the round-off error magnitude around 1.0E−32.This procedure is performed to ensure that there is no iteration error.The Repeated Richardson Extrapolation(RRE)method combines numerical solutions in different grids,determining higher orders of accuracy.The main contribution of this work is based on a process that initializes with fourth order solutions combining with RRE in order to find solutions of sixth,eighth,and tenth order of precision.The multigrid Full Approximation Scheme(FAS)is also applied to accelerate the convergence and obtain the numerical solutions on the fine grids.

On the Robustness of the xy-Zebra-Gauss-Seidel Smoother on an Anisotropic Diffusion Problem

Studies of problems involving physical anisotropy are applied in sciences and engineering,for instance,when the thermal conductivity depends on the direction.In this study,the multigrid method was used in order to accelerate the convergence of the iterative methods used to solve this type of problem.The asymptotic convergence factor of the multigrid was determined empirically(computer aided)and also by employing local Fourier analysis(LFA).The mathematical model studied was the 2D anisotropic diffusion equation,in whichε>0 was the coefficient of a nisotropy.The equation was discretized by the Finite Difference Method(FDM)and Central Differencing Scheme(CDS).Correction Scheme(CS),pointwise Gauss-Seidel smoothers(Lexicographic and Red-Black ordering),and line Gauss-Seidel smoothers(Lexicographic and Zebra ordering)in x and y directions were used for building the multigrid.The best asymptotic convergence factor was obtained by the Gauss-Seidel method in the direction x for 0<ε<<1 and in the direction y forε>>1.In this sense,an xy-zebra-GS smoother was proposed,which proved to be efficient and robust for the different anisotropy coefficients.Moreover,the convergence factors calculated empirically and by LFA are in agreement.