Introduction In recent years, domain decomposition methods have been used extensively to efficiently solve boundary value problems for partial differential equations in complex{form domains. On the other hand, multilevel techniques on hierarchical data structures also have developed into an effective tool for the construction and analysis of fast solvers. But direct realization of multilevel techniques on a parallel computer system for the global problem in the original domain involves difficult communication problems. I this paper, we present and analyze a combination of these two approaches: domain decomposition and multilevel decomposition on hierarchical structures to design optimal preconditioning operators.

Using material models on the basis of the flow theory of plasticity the asymptotic behaviour of solid mechanics solutions in crack tips, interface corners etc. strongly depends on the local realized load trajectory. For incrementally proportional load paths the equations determining the asymptotic fields are very simple ones. The paper considers two-dimensional statements in the neighbourhood of an interface corner consisting of two material ranges. At a distance from the corner the finite element nodes of a regular net are established in a polar co-ordinate system together with the displacement degrees of freedom. The main idea of the presented singular and non-singular stress and deformation field calculation at interface corners characterizes an replacement of the corner neighbourhood effect to the surrounding body by introducing stiffness actions which in usual manner can be assembled together with the other element stiffness matrices to the global stiffness matrix of the body. According to this there exists an in teresting invariant stiffness independence in corner and crack neighbourhoods. The applied technique allows extensions to non-proportional local load increments simplifying the mathematical calculations for the presentation of stress and strain fields in this general case. All computations are made on modern parallel computers. Concrete examples show the advantages of the presented approach.

We study two interacting particles in a random potential chain by a transfer matrix method which allows a correct handling of the symmetry of the two- particle wave function, but introduces an artificial ¨bag¨ interaction. The dependence of the two-particle localization length lambta 2on disorder, interaction strength and range is investigated. Our results demonstrate that the recently proposed enhancement of lambta 2 as compared to the results for single particles is vanishingly small for a Hubbard interaction. For longer-range interactions, we observe a small enhancement but with a different disorder dependence than proposed previously.

This work describes the creation of detailed FEM models of the selected area. Primarily is focused on the process of creating a detailed FEM model of the part of airplane using the static condensation. This work also contains a description of the process stress analysis of part of the fuselage of the airplane EV-55 Outback.

This paper deals with parallelization methods for time-dependent problems where the time steps are shared out among the processors. A Full Multigrid technique serves as solution algorithm, hence information of the preceding time step and of the coarser grid is necessary to compute the solution at each new grid level. Applying the usual extrapolation formula to process this information, the parallelization will not be very efficient. We developed another extrapolation technique which causes a much higher parallelization effect. Test examples show that no essential loss of exactness appears, such that the method presented here shall be well-applicable.

The purpose of this paper is to supply a summary of library subroutines and functions for parallel MIMD computers. The subroutines have been developed at the University of Chemnitz during a period of the last five years. In detail, they are concerned with vector operations, inter-processor communication and simple graphic output to workstations. One of the most valuable features is the machine-independence of the communication subroutines proposed in this paper for a hypercube topology of the parallel processors (excepting a kernel of only two primitive system-dependend operations). They were implemented and tested for different hardware and operating systems including transputer, nCube, KSR, PVM. The vector subroutines are optimized by the use of C language and enrolled loops (BLAS1-like). The paper includes hints for using the libraries with both Fortran and C programs.

The paper gives an overview over local inequalities for anisotropic simplicial Lagrangian finite elements. The main original contributions are the estimates for higher derivatives of the interpolation error, the formulation of the assumptions on admissible anisotropic finite elements in terms of geometrical conditions in the three-dimensional case, and an anisotropic variant of the inverse inequality. An application of anisotropic meshes in the context of a stabilized Galerkin method for a convection-diffusion problem is given.

Nonsteady Navier-Stokes equations represent a differential-algebraic system of strangeness index one after any spatial discretization. Since such systems are hard to treat in their original form, most approaches use some kind of index reduction. Processing this index reduction it is important to take care of the manifolds contained in the differential-algebraic equation (DAE). We investigate for several discretization schemes for the Navier-Stokes equations how the consideration of the manifolds is taken into account and propose a variant of solving these equations along the lines of the theoretically best index reduction. Applying this technique, the error of the time discretisation depends only on the method applied for solving the DAE.

The characterization of the local ill-posedness and the local degree of nonlinearity are of particular importance for the stable solution of nonlinear ill-posed problems. We present assertions concerning the interdependence between the ill-posedness of the nonlinear problem and its linearization. Moreover, we show that the concept of the degree of nonlinearity com bined with source conditions can be used to characterize the local ill-posedness and to derive a posteriori estimates for nonlinear ill-posed problems. A posteriori estimates are widely used in finite element and multigrid methods for the solution of nonlinear partial differential equations, but these techniques are in general not applicable to inverse an ill-posed problems. Additionally we show for the well-known Landweber method and the iteratively regularized Gauss-Newton method that they satisfy a posteriori estimates under source conditions; this can be used to prove convergence rates results.

A class of Lagrange-Newton-SQP methods is investigated for optimal control problems governed by semilinear parabolic initial- boundary value problems. Distributed and boundary controls are given, restricted by pointwise upper and lower bounds. The convergence of the method is discussed in appropriate Banach spaces. Based on a weak second order sufficient optimality condition for the reference solution, local quadratic convergence is proved. The proof is based on the theory of Newton methods for generalized equations in Banach spaces.