Robust Solution of Singularly Perturbed Problems Using Multigrid
نویسندگان
چکیده
We consider the problem of solving linear systems of equations that arise in the numerical solution of singularly perturbed ordinary and partial differential equations of reaction-diffusion type. Standard discretization techniques are not suitable for such problems and, so, specially tailored methods are required, usually involving adapted or fitted meshes that resolve important features such as boundary and/or interior layers. In this study, we consider classical finite difference schemes on the layer adapted meshes of Shishkin and Bakhvalov. We show that standard direct solvers exhibit poor scaling behaviour, with respect to the perturbation parameter, when solving the resulting linear systems. We propose and prove optimality of a new block-structured preconditioning approach that is robust for small values of the perturbation parameter, and compares favourably with standard robust multigrid preconditioners for these linear systems. We also derive stopping criteria which ensure that the potential accuracy of the layer-resolving meshes is achieved.
منابع مشابه
Robust Solution of Singularly Perturbed Problems Using Multigrid Methods; Analysis and Numerical Results in One and Two Dimensions∗
We consider the problem of solving linear systems of equations that arise in the numerical solution of singularly perturbed ordinary and partial differential equations of reaction-diffusion type. Standard discretization techniques are not suitable for such problems and, so, specially tailored methods are required, usually involving adapted or fitted meshes that resolve important features such a...
متن کاملRobust Solution of Singularly Perturbed Problems Using Multigrid Methods
We consider the problem of solving linear systems of equations that arise in the numerical solution of singularly perturbed ordinary and partial differential equations of reaction-diffusion type. Standard discretization techniques are not suitable for such problems and, so, specially tailored methods are required, usually involving adapted or fitted meshes that resolve important features such a...
متن کاملA hybrid method for singularly perturbed delay boundary value problems exhibiting a right boundary layer
The aim of this paper is to present a numerical method for singularly perturbed convection-diffusion problems with a delay. The method is a combination of the asymptotic expansion technique and the reproducing kernel method (RKM). First an asymptotic expansion for the solution of the given singularly perturbed delayed boundary value problem is constructed. Then the reduced regular delayed diffe...
متن کاملNumerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type
In this paper, we have proposed a numerical method for singularly perturbed fourth order ordinary differential equations of convection-diffusion type. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference method. In order to get a numerical solution for the derivative of the solution, the given interval is divided in...
متن کاملA first-order system Petrov–Galerkin discretization for a reaction–diffusion problem on a fitted mesh
We consider the numerical solution, by a Petrov–Galerkin finite-element method, of a singularly perturbed reaction–diffusion differential equation posed on the unit square. In Lin & Stynes (2012, A balanced finite element method for singularly perturbed reaction-diffusion problems. SIAM J. Numer. Anal., 50, 2729–2743), it is argued that the natural energy norm, associated with a standard Galerk...
متن کامل