University, Diss., 2005--Kassel.

Techn. Hochsch., Diss., 2005--Aachen.

Techn. Universiẗat, Diss., 2006--Darmstadt.

Zugl.: Stuttgart, Univ., Diss., 2008

After some remarks on the parallel implementation of the Finite Element package FEAP, our realisation of the parallel CG-algorithm is sketched. From a technical point of view, a hierarchical preconditioner with and without additional global crosspoint preconditioning is presented. The numerical properties of this preconditioners are discussed and compared to a Schur-complement-preconditioning, using a wide range of data from computations on technical and academic examples from elasticity.

Es werden Fragen der effektiven Kommunikation in parallelen FEM-Systemen behandelt. Durch geschickte Partitionierung des betrachteten Gebietes und Verteilung auf die vorhandenen Prozessoren kann man versuchen, die Kommunikationslast zu minimieren, z. B. mit dem Programmsystem chaco. Ein möglichst paralleler Datenaustausch wird durch Anordnung der Kommunikationsschritte in so genannten Linkleveln versucht. In der vorliegenden Arbeit wird ausgehend von der MPI-basierten Koppelrandkommunikation ein Split-Algorithmus vorgestellt, der versucht, die Koppelranddaten großer Kommunikatoren auf die kleinerer Sub-Kommunikatoren abzubilden und damit die Anzahl der zu übertragenden Datenpakete zu minimieren.

The continuum theory of large deformation elasto-plasticity is summarized as far as it is necessary for the numerical treatment with the Finite-Element-Method. Using the calculus of modern differential geometry and functional analysis, the fundamental equations are derived and the proof of most of them is shortly outlined. It was not our aim to give a contribution to the development of the theory, rather to show the theoretical background and the assumptions to be made in state of the art elasto-plasticity.

Singularly perturbed reaction-diffusion problems exhibit in general solutions with anisotropic features, e.g. strong boundary and/or interior layers. This anisotropy is reflected in the discretization by using meshes with anisotropic elements. The quality of the numerical solution rests on the robustness of the a posteriori error estimator with respect to both the perturbation parameters of the problem and the anisotropy of the mesh. An estimator that has shown to be one of the most reliable for reaction-diffusion problem is the <i>equilibrated residual method</i> and its modification done by Ainsworth and Babuška for singularly perturbed problem. However, even the modified method is not robust in the case of anisotropic meshes. The present work modifies the equilibrated residual method for anisotropic meshes. The resulting error estimator is equivalent to the equilibrated residual method in the case of isotropic meshes and is proved to be robust on anisotropic meshes as well. A numerical example confirms the theory.

In this paper, we consider domain decomposition preconditioners for a system of linear algebraic equations arising from the <i>p</i>-version of the fem. We propose several multi-level preconditioners for the Dirichlet problems in the sub-domains in two and three dimensions. It is proved that the condition number of the preconditioned system is bounded by a constant independent of the polynomial degree. The proof uses interpretations of the <i>p</i>-version element stiffness matrix and mass matrix on [-1,1] as <i>h</i>-version stiffness matrix and weighted mass matrix. The analysis requires wavelet methods.

In the adaptive finite element method, the solution of a p.d.e. is approximated from finer and finer meshes, which are controlled by error estimators. So, starting from a given coarse mesh, some elements are subdivided a couple of times. We investigate the question of avoiding instabilities which limit this process from the fact that nodal coordinates of one element coincide in more and more leading digits. In a previous paper the stable calculation of the Jacobian matrices of the element mapping was given for straight line triangles, quadrilaterals and hexahedrons. Here, we generalize this ideas to linear and quadratic triangles on curved boundaries.