New Approximate Method for Solving the Stokes Problem in a Domain with Corner Singularity

In this paper we introduce the notion of an Rν-generalized solution to the Stokes problem with singularity in a two-dimensional non-convex polygonal domain with one reentrant corner on its boundary in special weight sets. We construct a new approximate solution of the problem produced by weighted nite element method. An iterative process for solving the resulting system of linear algebraic equations with a block preconditioning of its matrix is proposed on the basis of the incomplete Uzawa algorithm and the generalized minimal residual method. Results of numerical experiments have shown that the convergence rate of the approximate Rν-generalized solution to an exact one is independent of the size of the reentrant corner on the boundary of the domain and equals to the rst degree of the grid size h in the norm of the weight space W1 2,ν(Ω) for the velocity eld components in contrast to the approximate solution produced by classical nite element or nite dierence schemes convergence to a generalized one no faster than at an O(h α) rate in the norm of the space W^1_2 (Ω) for the velocity eld components, where α < 1 and α depends on the size of the reentrant corner.В статье определено понятие R_u-обобщенного решения задачи Стокса с сингулярностью в двумерной невыпуклой многоугольной области с одним входящим углом на границе области в специальных весовых множествах. Построено новое приближенное решение задачи с помощью весового метода конечных элементов. Предложен итерационный процесс решения полученной системы линейных алгебраических уравнений с блочным переобуславливанием ее матрицы на основе неполного алгоритма Удзавы и обобщенного метода минимальных невязок. Результаты численных экспериментов показали, что скорость сходимости приближенного R_u-обобщенного решения к точному решению задачи не зависит от величины входящего угла на границе области и равна первой степени по шагу сетки h в норме весового пространства W^1_{2,u}(Omega) для компонент вектора скоростей, в отличие от стандартных конечно-элементных и конечно-разностных схем, приближенное решение которых сходится к точному решению задачи не быстрее чем со скоростью O(h^alpha) в норме пространства W^1_2(Omega) для компонент вектора скоростей, где alpha