Refine
Year of publication
Document Type
- ZIB-Report (29)
- In Proceedings (20)
- Article (18)
- Book chapter (1)
- In Collection (1)
Language
- English (69)
Is part of the Bibliography
- no (69)
Keywords
Institute
- Numerical Mathematics (34)
- Computational Medicine (33)
- ZIB Allgemein (18)
- Therapy Planning (8)
- Visual Data Analysis (8)
- Visual and Data-centric Computing (8)
In this paper we present a self--adaptive finite element method to solve flame propagation problems in 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. The proposed method is applied to an unsteady thermo--diffusive combustion model to demonstrate its potential for the solution of complicated problems.
An error controlled finite elemente method (FEM) for solving stationary Schrödinger equations in three space dimensions is proposed. The method is based on an adaptive space discretization into tetrahedra and local polynomial basis functions of order $p=1$--$5$ defined on these tetrahedra. According to a local error estimator the triangulation is automatically adapted to the solution. Numerical results for standard problems appearing in vibrational motion and molecular structure calculations are presented and discussed. Relative precisions better than 1e-8 are obtained. For equilateral H$_3^{++}$ the adaptive FEM turns out to be superior to global basis set expansions in the literature. Our precise FEM results exclude in a definite manner the stability or metastability of equilateral H$_3^{++}$ in its groundstate.
{\def\enorm {\mathop{\mbox{\boldmath{$|\!|$}}}\nolimits} Let $u \in H$ be the exact solution of a given self--adjoint elliptic boundary value problem, which is approximated by some $\tilde{u} \in {\cal S}$, $\cal S$ being a suitable finite element space. Efficient and reliable a posteriori estimates of the error $\enorm u - \tilde{u}\enorm $, measuring the (local) quality of $\tilde{u}$, play a crucial role in termination criteria and in the adaptive refinement of the underlying mesh. A well--known class of error estimates can be derived systematically by localizing the discretized defect problem using domain decomposition techniques. In the present paper, we provide a guideline for the theoretical analysis of such error estimates. We further clarify the relation to other concepts. Our analysis leads to new error estimates, which are specially suited to three space dimensions. The theoretical results are illustrated by numerical computations.}
We consider the discretization of obstacle problems for second order elliptic differential operators in three space dimensions by piecewise linear finite elements. Linearizing the discrete problems by suitable active set strategies, the resulting linear sub--problems are solved iteratively by preconditioned cg--iterations. We propose a variant of the BPX preconditioner and prove an $O(j)$ estimate for the resulting condition number. To allow for local mesh refinement we derive semi--local and local a posteriori error estimates. The theoretical results are illustrated by numerical computations.
We consider the approximate solution of selfadjoint elliptic problems in three space dimensions by piecewise linear finite elements with respect to a highly non-uniform tetrahedral mesh which is generated adaptively. The arising linear systems are solved iteratively by the conjugate gradient method provided with a multilevel preconditioner. Here, the accuracy of the iterative solution is coupled with the discretization error. as the performance of hierarchical bases preconditioners deteriorate in three space dimensions, the BPX preconditioner is used, taking special care of an efficient implementation. Reliable a-posteriori estimates for the discretization error are derived from a local comparison with the approximation resulting from piecewise quadratic elements. To illustrate the theoretical results, we consider a familiar model problem involving reentrant corners and a real-life problem arising from hyperthermia, a recent clinical method for cancer therapy.
Adaptive Solutions of Nonlinear Parabolic Equations with Application to Hyperthermia Treatments
(1997)
We present a self-adaptive finite element method to solve nonlinear evolution problems in 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. The proposed method is applied to hyperthermia treatments to demonstrate its potential for the solving of complicated problems.
The transient heating in an unsymmetrical coated hot--strip sensor was simulated with a self--adaptive finite element method. The first tests of this model show that it can determine with a small error the thermal conductivity of liquids, from the transient temperature rise in the hot--strip, deposited in a substrate and coated by an alumina spray.
We present a self--adaptive finite element method to solve combustion problems in 1D, 2D, and 3D. An implicit time integrator of Rosenbrock type is coupled with a multilevel approach in space. A posteriori error estimates are obtained by constructing locally higher order solutions involving all variables of the problem. Adaptive strategies such as step size control, spatial refinement and coarsening allow us to get economically an accurate solution. Various examples are presented to demonstrate practical applications of the proposed method.
KARDOS solves nonlinear evolution problems in 1, 2, and 3D. An adaptive multilevel finite element algorithm is used to solve the spatial problems arising from linearly implicit discretization methods in time. Local refinement and derefinement techniques are used to handle the development of the mesh over time. The software engineering techniques used to implement the modules of the KASKADE toolbox are reviewed and their application to the extended problem class is described. A notification system and dynamic construction of records are discussed and their values for the implementation of a mesh transfer operation are shown. The need for low-level and high--level interface elements of a module is discussed for the assembling procedure of KARDOS. At the end we will summarize our experiences.
KARDOS - User"s Guide
(2002)
The adaptive finite element code {\sc Kardos} solves nonlinear parabolic systems of partial differential equations. It is applied to a wide range of problems from physics, chemistry, and engineering in one, two, or three space dimensions. The implementation is based on the programming language C. Adaptive finite element techniques are employed to provide solvers of optimal complexity. This implies a posteriori error estimation, local mesh refinement, and preconditioning of linear systems. Linearely implicit time integrators of {\em Rosenbrock} type allow for controlling the time steps adaptively and for solving nonlinear problems without using {\em Newton's} iterations. The program has proved to be robust and reliable. The user's guide explains all details a user of {\sc Kardos} has to consider: the description of the partial differential equations with their boundary and initial conditions, the triangulation of the domain, and the setting of parameters controlling the numerical algorithm. A couple of examples makes familiar to problems which were treated with {\sc Kardos}. We are extending this guide continuously. The latest version is available by network: {\begin{rawhtml} <A href="http://www.zib.de/Numerik/software/kardos/"> <i> Downloads.</i></a> \end{rawhtml}}