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.

In this paper, a uniformly elliptic second order boundary value problem in 2D is discretized by the p-version of the finite element method. An inexact Dirichlet-Dirichlet domain decomposition pre-conditioner for the system of linear algebraic equations is investigated. The solver for the problem in the sub-domains and a pre-conditioner for the Schur-complement are proposed as ingredients for the inexact DD-pre-conditioner. Finally, several numerical experiments are given.

On a large class of two-dimensional anisotropic meshes, the inf-sup condition (stability) is proved for the triangular and quadrilateral finite element pairs suggested by Bernardi/Raugel and Fortin. As a consequence the pairs ${\cal P}_2-{\cal P}_0$, ${\cal Q}_2-{\cal P}_0$, and ${\cal Q}_2^\prime-{\cal P}_0$ turn out to be stable independent of the aspect ratio of the elements.

The paper is concerned with the Nitsche mortaring in the framework of domain decomposition where non-matching meshes and weak continuity of the finite element approximation at the interface are admitted. The approach is applied to singularly perturbed reaction-diffusion problems in 2D. Non-matching meshes of triangles being anisotropic in the boundary layers are applied. Some properties as well as error estimates of the Nitsche mortar finite element schemes are proved. In particular, using a suitable degree of anisotropy of triangles in the boundary layers of a rectangle, we derive convergence rates as known for the conforming finite element method in presence of regular solutions. Numerical examples illustrate the approach and the results.

The preprint delivers an efficient solution technique for the numerical simulation of crack propagation of 2D linear elastic formulations based on finite elements together with the conjugate gradient method in order to solve the corresponding linear equation systems. The developed iterative numerical approach using hierarchical preconditioners comprehends the interesting feature that the hierarchical data structure will not be destroyed during crack propagation. Thus, one gets the possibility to simulate crack advance in a very effective numerical manner including adaptive mesh refinement and mesh coarsening. Test examples are presented to illustrate the efficiency of the given approach. Numerical simulations of crack propagation are compared with experimental data.