AUTOMATIC NUMERICAL ELEMENT GENERATION BY BOUNDARY-FITTED CURVILINEAR COORDINATE SYSTEM

Abstract

A method for dividing a two-dimensional multi-connected region of a complex shape into a set of triangular elements is developed. A region of a complex shape in the physical plane is divided into some simply connected subregions, and each subregion is mapped onto a square region in the transformed plane. The inverse functions of the mapping are determined by the solution to elliptic partial differential equations with the Dirichlet boundary conditions. After the square region is divided into a set of finite elements, each element is inversely mapped onto the subregions by use of the functions. The finite element data for the global region are made of those for the divided subregions.