Dealing with periodic boundary conditions for 1D, 2D and 3D isogeometric finite element method

Marcin Łoś¹ , Maciej Paszyński¹, Lisandro Dalcin², Victor Calo²

¹AGH University of Science and Technology, al. Mickiewicza 30, 30-059 Krakow, Poland.
²King Abdullah University of Science and Technology, Thuwal, Saudi Arabia.

DOI:

https://doi.org/10.7494/cmms.2015.1.0525

Abstract:

In this paper we analyze the problem of implementing periodic boundary conditions in the isogeometric finite element method (ISO-FEM). The ISO-FEM method uses the B-spline-based basis functions, which facilitates usage of the same basis functions for approximation of the geometry as well as for the numerical solution of the modeled physical phenomena. The usage of the B-spline based basis functions results in C^(p-1) global continuity of the solution. The drawback is a difficulty in implementing the periodic boundary conditions, and special dedicated methods are necessary. In this paper we present two algorithms implementing the periodic boundary conditions. The first one is an iterative algorithm that utilizes widely available block-diagonal LAPACK solver. The second one is a modification of the multi-frontal solver algorithm itself, and it requires a dedicated solver with its source code modified accordingly. The presented methods can be applied in one, two or three-dimensional isogeometric finite element method.

Cite as:

Łoś, M., Paszyński, M., Dalcin, L., & Calo, V. (2015). Dealing with periodic boundary conditions for 1D, 2D and 3D isogeometric finite element method. Computer Methods in Materials Science, 15(1), 213-218. https://doi.org/10.7494/cmms.2015.1.0525

Article (PDF):