High-accuracy schemes of the finite element method for systems of degenerate elliptic equations on an interval

DOI10.3103/S1066369X09070032MaRDI QIDQ735986zbMATH OpenFDO

Authors A. D. Lyashko, Sh. I. Tayupov, M. R. Timerbaev

Publication date 26 October 2009

Published in Russian Mathematics (Search for Journal in Brave)

Full work available at URL https://doi.org/10.3103/S1066369X09070032

continuation function convergence degenerate elliptic equation Dirichlet problem finite element method multiplicative singularity extraction system weighted Sobolev spaces

Mathematics Subject Classification ID

Boundary value problems for second-order elliptic equations (35J25) Second-order elliptic systems (35J47) Degenerate elliptic equations (35J70) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Error bounds for boundary value problems involving PDEs (65N15) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)

Review text

This paper deals with a finite element scheme for a system of degenerate second order elliptic equations in one dimension. It is known that the exact solution in such problems have unbounded derivatives near degeneration points. Therefore, the discretization of such problems by the standard finite element method leads to loss of convergence of approximate solution to the exact one near degeneration points of coefficients of the differential operator of the problems. The finite element scheme suggested by the authors is based on a representation of the exact solution in terms of some known functions and unknown but smooth function. In the case of the Dirichlet problem with homogeneous boundary conditions, the authors prove that the exact solution is a product of two functions. The first function is very simple and describes the behavior of the exact solution near the the degeneration points of coefficients of the differential equations, and, the second function is treated as a new desired smooth function. In the case of the Dirichlet problem with non homogeneous boundary conditions, the authors use the so called continuation function (of boundary values into the domain) to reduce the problem to a new problem with homogeneous Dirichlet boundary conditions. A construction for the continuation function is provided. By this way, the authors prove that the exact solution of the original problem can be represented as the sum of two terms. The first term is known and the second one is the product of a known function and a unknown smooth function. In both cases of homogeneous and non homogeneous boundary conditions and thanks to the above representation for the exact solution, the authors apply the usual finite element method to approximate the unknown smooth function which is the solution of a new system of smooth second order elliptic equations in one dimension. They remark that this new unknown function is smooth and can be approximated by a piecewise polynomial function with a higher convergence order. This last finite element approximation allows the authors, thanks to the above stated representation, to obtain an approximation for the exact solution of the original system, with optimal order in weighted Sobolev norms.

Recommendations

Cites work

Cited in

(8)

This page was built for publication: High-accuracy schemes of the finite element method for systems of degenerate elliptic equations on an interval

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q735986)