Efficient iterative methods and solvers for FEM analysis

Tomáš Luber

Doctoral thesis

Abstract:

Hlavním tématem této práce je trojpolní formulace Biotova modelu poroelasticity a předpodmiňovače k řešení soustav vzniklých diskretizací modelu pomocí metody konečných prvků. Pro popis stavu porézního prostředí jsou použita pole posunutí, toku a tlaku kapaliny. Biotův model je představen nejprve v klasické formulaci této trojpolní formulace a následně je přestavena variační formulace problému pro krok diskretizace v čase. Hlavní náplní práce jsou předpodmiňovače pro řešení systému rovnic pro jednotlivé časové kroky. Je představeno několik blokově diagonálních předpodmiňovačů založených na Schurových doplňcích, které jsou analyzovány na úrovni slabé formulace a následně na maticové úrovni. Stejný přístup je poté aplikovám pro předpodmiňovače pro Biotův-Barenblattův model, který je zobecněním Biotova modelu. Výsledky analýzy čísel podmíněnosti jednotlivých předpodmiňovačů jsou podpořeny numerickými experimenty. Přínos práce je založen na publikovaných výsledcích autora, které popisují konstrukci a analýzu předpodmínění založeného na Schurových doplňcích. V práci jsou navíc tyto výsledky zobecněné tím, že většina analýzy se přesouvá do slabá formulace. Prezentace se tak stává přehlednější a umožňuje jednodušší odvození výsledků nezávisejících na konkretní diskretizaci.

Abstract:

The main topic of this thesis is the three-field formulation of Biot's model of poroelasticity and preconditioners to solve the equations of the model after finite element discretization. The three field-formulation that is used describes the state of the porous media using displacement, fluid flux and pressure. Using this formulation we introduce a classical formulation of the model, use implicit time discretization and introduce the weak formulation for the timestep problem. The main focus is on the preconditioners for the timestep problem. Several block diagonal preconditioners based on Schur complements are presented and analysed on the level of the weak formulation and on the matrix level. The same approach is then applied to Biot-Barenblatt model that generalizes Biot's model. The results on condition number estimates are illustrated by numerical experiments. The author's contributions are based on the coauthored articles that cover the construction of Schur complement preconditioners. This thesis generalizes them by taking most of the analysis to the weak formulation which provides clearer presentation and discretization-independent results.

Keywords

Numerická matematika, poroelasticita, smíšené konečné prvky, předpodmínění, Schurův doplněk

Keywords

Numerical mathematics, poroelasticity, mixed finite elements, preconditioning, Schur complement

Language used: English

Date on which the thesis was submitted / produced: 28. 2. 2022

Thesis defence

• Date of defence: 29. 6. 2022
• Supervisor: Stanislav Sysala
• Reader: Dalibor Lukáš, Václav Kučera, Jan Stebel

VSB – Technical University of Ostrava

Fakulta elektrotechniky a informatiky

Doctoral programme / field:
Informatika, komunikační technologie a aplikovaná matematika / Výpočetní a aplikovaná matematika