Párování metod hraničních a konečných prvků v prostředí Netgen/NgSolve

Abstract

Párování metod konečných a hraničních prvků patří k efektivním numerickým metodám pro řešení okrajových eliptických úloh. V praxi se často setkáme s problémy, které potřebujeme řešit na neomezených oblastech. Tyto úlohy nelze přímo řešit pomocí metody konečných prvků (FEM), ale lze je efektivně řešit metodou hraničních prvků (BEM), kterou si v této práci podrobněji popíšeme. BEM je však efektivní jen dokud se v oblasti hledaného řešení jedná o Laplaceovu rovnici. Pokud je však úloha na této oblasti Poissonova typu, tak samotná metoda hraničních prvků ztrácí svou efektivitu a je nutné ji párovat s metodou konečných prvků. V této práci ukážeme využití naší BEM knihovny v prostředí Netgen/NGSolve a demonstrujeme ji i pro Johnsonovo-Nédélecovo FEM-BEM párování pro úlohu transmise. Funkčnost a efektivitu naší implementace ově-\break říme na různých problémech pomocí numerických experimentů.

Coupling finite and boundary element methods is one of the most efficient numerical methods for solving elliptic boundary value problems. In practice, we often encounter problems that need to be solved on unbounded domains. These problems cannot be directly solved using the finite element method (FEM), but can be efficiently solved using the boundary element method (BEM), which we will describe in detail in this thesis. However, the BEM is only effective as long as the problem on the region of interest is described by a Laplace equation. But if the problem on this region is of Poisson type, then the boundary element method alone loses its effectiveness and must be paired with the finite element method. In this thesis, we show the use of our BEM library in the Netgen/NGSolve environment and we also demonstrate the Johnson-Nedellec FEM-BEM coupling for the transmission problem. We verify the functionality and efficiency of our implementation on various problems through numerical experiments.