Parallelization of the solution of elliptic boundary value problems using TFETI domain decomposition method

Abstract

Import 03/08/2012Práce se zabývá paralelní implementací TFETI metody rozložení oblasti, která je vyvíjena na Katedře aplikované matematiky VŠB - TU Ostrava pro řešení rozsáhlých inženýrských problémů, a dále její aplikací na řešení eliptických okrajových úloh. Po diskretizaci oblasti pomocí metody konečných prvků se oblast rozdělí na nepřekrývající se podoblasti tvořené elementy sítě. Toto vede na soustavu rovnic s blokovou maticí tuhosti, což usnadňuje paralelizaci jejího řešení. Pomocí Lagrangeových multiplikátorů jsou pak vynucovány nejen lepící podmínky, ale i Dirichletovy okrajové podmínky. Metoda je tak ve srovnání s metodou FETI efektivnější a snažší na implementaci. Je zde také provedena efektivní regularizace matic tuhosti podoblastí a jsou zavedeny ortogonální projektory, které zlepšují podmíněnost duální úlohy s Lagrangeovými multiplikátory. Na numerických experimentech je ověřena numerická i paralelní škálovatelnost metody.The bachelor thesis deals with parallel implementation of the TFETI domain decomposition method, which is developed at the Department of Applied Mathematics VŠB - TU Ostrava to solve large engineering problems. It deals also with application of TFETI on the solution of elliptic boundary value problems. After finite element space discretization we decompose the domain into nonoverlapping subdomains defined by mesh elements. This leads to a system of linear equations with block-diagonal system matrix which enables its effective parallel solution. The gluing conditions and even the Dirichlet boundary conditions are enforced by the Lagrange multipliers, so the method is more effective and simpler for implementation than FETI. The effective regularization of the stiffness matrices of the subdomains is performed and the condition number of the dual problem with Lagrange multipliers is improved by using orthogonal projectors. Both numerical and parallel scalability of TFETI are tested in the numerical experiments.470 - Katedra aplikované matematikyvýborn