Isogeometric Analysis of PDEs and numerical implementation in the Finite Element Library LifeV
The aim of the project is double: to understand the flexibility of the Isogeometric Analysis tools through the solution of some PDEs problems; to test the improvement in the computational time given by a partial loops vectorization at compile-time of the LifeV IGA code. Three different appli- cations have been selected: the potential flow problem around an airfoil profile, the heat equation problem in a bent cylinder and the Laplace problem in a multi-patches geometry representing a blood vessels bifurcation. The geometries used are built through the NURBS package available with the software GeoPDEs. The numerical analysis of the first application is performed with both GeoPDEs and LifeV IGA code. The comparison between different implementations shows that the degrees of freedom loop vectorization at compile-time is able to reduce the matrix assembling time of around 20%. The automatic vectorization at compile-time of the loop on the elements requires too much computational effort without having a reasonable improvement in the running time performances. Unsteady problems and multi-patches geometry have not been tested with LifeV IGA code, but GeoPDEs results show the expected solutions.