Article information
2018 , Volume 23, ¹ 2, p.102-116
Rukavishnikov A.V.
Stabilized numerical method for solving the Oseen type problem with a singularity
Purpose. To construct modified approximation approach using the finite element method and to perform numerical analysis for a two dimensional problem on the flow of a viscous inhomogeneous fluids - the Oseen type problem, that is obtained by sampling in time and linearizing the incompressible Navier-Stokes equations. To consider the convection dominated flow case. Methodology. Based on the domain decomposition method with a smooth curvilinear boundary between subdomains, a stabilization nonconformal finite element method is constructed that satisfies the inf-sup-stability condition. To solve the resulting system of linear algebraic equations, an iterative process is considered that uses the decomposition of the vector in the Krylov subspace with minimal inviscidity, with a block preconditioning of the matrix. Findings. The results of the numerical experiments demonstrate the robustness of the considered method for different (even small) discontinuous values of viscosity. The differences between finite element and exact solutions for the velocity field and pressure in the norms of the grid spaces decrease as 𝒪( h ) for each curvilinear interface, that agrees with the theoretical estimate which is proved by the author. Originality/value. Using the stabilized finite element method is often good and better then the classical finite element method for convection dominated flows with discontinuous viscosity and density.
[full text] Keywords: domain decomposition method, Oseen type problem, stabilized finite element method, discontinuous coefficients
doi: 10.25743/ICT.2018.23.12803
Author(s): Rukavishnikov Alexey Victorovich PhD. , Associate Professor Position: Leading research officer Office: Computing Center of the Far-Eastern Branch Russian Academy of Sciences Address: 680000, Russia, Khabarovsk, 65, Kim Yu Chen Str.
Phone Office: (4212) 70-43-42 SPIN-code: 7680-1450 References:
[1] Braess, D., Dahmen, W., Wieners, C. A multigrid algorithm for the mortar finite element methods. SIAM Journal on Numerical Analysis. 1999; 37(1):48-69.
[2] Bernardi, C., Maday, Y., Patera, A. A New nonconforming approach to domain decomposition: The mortar element method. Nonlinear Partial Differential Equations and their Applications. Eds H. Brezis and J.L. Lions. Paris: Pitman; 1994: 13-51.
[3] Wohlmuth, B. Hierarchical a posteriori error estimators for mortar finite element methods with Lagrange multipliers. SIAM Journal on Numerical Analysis. 1999; 36(5):1636-1658.
[4] Flemisch, B., Melenk, J.M., Wohlmuth, B. Mortar methods with curved interfaces. Applied Numerical Mathematics. 2005; (54):339–361.
[5] Huang, J., Zou, J. A mortar element method for elliptic problems with discontinuous coefficients. IMA Journal of Numerical Analysis. 2002; (22):549-576.
[6] Rukavishnikov, A.V., Rukavishnikov, V.A. On the nonconformal finite element method for the Stokes problem with a discontinuous coefficient. Journal of Applied and Industrial Mathematics. 2007;. Vol. 10. 𝒩 4, P. 104-117. (In Russ.)
[7] Ben Belgacem, F. The mixed mortar finite element method for the incompressible Stokes problem: convergence analysis. SIAM Journal on Numerical Analysis. 2000; 37(4):1085-1100.
[8] Rukavishnikov, A.V. On construction of a numerical method for the Stokes problem with a discontinuous coefficient of viscosity. Computational Technologies. 2009; 14(2):110-123. (In Russ.)
[9] Rukavishnikov, A.V. Nonconformal finite element method for a fluid dynamics problem with a curved interface. Zhurnal Vychislitel'noi Matematiki i Matematicheskoi Fiziki. 2012; 52(6):1072-1094. (In Russ.)
[10] Rukavishnikov, A.V. Domain decomposition method and numerical analysis of a fluid dynamics problem. Computational Mathematics and Mathematical Physics. 2014; 54(9):1459–1480.
[11] Brezzi, F., Fortin, M. Mixed and hybrid finite element methods. New-York: Springer-Verlag; 1991: 362.
[12] Ciarlet, P. The Finite element method for elliptic problems. Amsterdam: North-Holand; 1978: 530.
[13] Kuznetsov, Yu. A. Efficient iterative solvers for elliptic finite element problems on nonmatching grids. Russian Journal of Numerical Analysis and Mathematical Modelling. 1995; (10):187–211.
[14] Olshanskii, M.A., Reusken, A. Grad-div stabilization for Stokes equations. Mathematics of Computation. 2004; (73):1699–1718.
[15] Grisvard, P. Elliptic problems in nonsmooth domains. Boston: Pitman; 1985: 410.
[16] Bramble, J.H., Pasciak, J.E., Vassilev, A.T. Analysis of the inexact Uzawa algorithm for saddle point problems. SIAM Journal on Numerical Analysis. 1997; (34):1072–1092.
[17] Saad, Y. Iterative methods for sparse linear systems. New Jersey: PWS Publishing Company; 1996: 417.
[18] Little, L., Saad, Y., Smoch, L. Block, LU preconditioners for symmetric and nonsymmetric saddle point problems. SIAM Journal on Scientific Computing. 2003; (25):729–748.
[19] Olshanskii, M.A. Reusken A. Analysis of a Stokes interface problem. Numerische Mathematik. 2006; (103):129-149.
Bibliography link: Rukavishnikov A.V. Stabilized numerical method for solving the Oseen type problem with a singularity // Computational technologies. 2018. V. 23. ¹ 2. P. 102-116
|