In this paper, we study the numerical solution of the Stokes system in deformed axisymmetric geometries. In the azimuthal direction the discretization is carried out by using truncated Fourier series, thus reducing the dimension of the problem. The resulting two-dimensional problems are discretized using the spectral element method which is based on the variational formulation in primitive variables. The meridian domain is subdivided into elements, in each of which the solution is approximated by truncated polynomial series. The results of numerical experiments for several geometries are presented.

This work is devoted to the study of a two-dimensional vectorPoisson equation with the normal component of the unknown andthe value of the divergence of the unknown prescribed simultaneouslyon the entire boundary.These two scalar boundary conditions appear prima faciealternative in a standard variational framework. An originalvariational formulation of this boundary value problemis proposed here. Furthermore, an uncoupled solution algorithm isintroduced together with its finite element approximation.The numerical scheme has been implemented and appliedto solve a simple test problem.