Unified multilevel adaptive finite element methods for elliptic problems
Many elliptic partial differential equations can be solved numerically with near optimal efficiency through the uses of adaptive refinement and multigrid solution techniques. It is our goal to develop a more unified approach to the combined process of adaptive refinement and multigrid solution which can be with high order finite elements. The basic step of the refinement process is the bisection of a pair of triangles, which corresponds to the addition of one or more basis functions to the approximation space. An approximation of the resulting change in the solution is used as an error indicator to determine which triangles to divide. The multigrid iteration uses a red-black Gauss--Seidel relaxation in which the black relaxations are used only locally. The grid transfers use the change between the nodal and hierarchical bases. This multigrid iteration requires only O(/Nu/) operations, even for highly nonuniform grids, and is defined for any finite element space. The full multigrid method is an optimal blending of the processes of adaptive refinement and multigrid iteration. So as to minimize the number of operations required, the duration of the refinement phase is based on increasing the dimension of the approximation space by some fixed factor which is determined to be the largest possible for the given error-reducing power of the multigrid iteration. The result is an algorithm which (i) uses only O(/Nu/) operations with a reasonable constant of proportionality, (ii) solves the discrete system to the accuracy of the discretization error, (iii) is able to achieve the optimal order of convergence of the discretization error in the presence of singularities. Numerical experiments confirm this for linear, quadratic and cubic elements. It is believed that the method can also be applied to more practical problems involving systems of PDE'S, time dependence, and three spatial dimensions. 39 refs.
- Research Organization:
- Illinois Univ., Urbana (USA). Dept. of Computer Science
- DOE Contract Number:
- FG02-87ER25026
- OSTI ID:
- 7230871
- Report Number(s):
- DOE/ER/25026-21; UIUCDCS-R-88-1436; UILU-ENG-88-1741; ON: DE88011938
- Resource Relation:
- Other Information: Thesis (Ph.D.)
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multilevel adaptive solution procedure for material nonlinear problems in visual programming environment
An efficient flexible-order model for 3D nonlinear water waves
Related Subjects
PARTIAL DIFFERENTIAL EQUATIONS
FINITE ELEMENT METHOD
ALGORITHMS
ELLIPTICAL CONFIGURATION
NUMERICAL SOLUTION
CONFIGURATION
DIFFERENTIAL EQUATIONS
EQUATIONS
MATHEMATICAL LOGIC
990230* - Mathematics & Mathematical Models- (1987-1989)