We describe an optimization process specially designed for regional hyperthermia of deep-seated tumors in order to achieve desired steady-state temperature distributions. A nonlinear three-dimensional heat transfer model based on temperature-dependent blood perfusion is applied to predict the temperature. Using linearly implicit methods in time and adaptive multilevel finite elements in space, we are able to integrate efficiently the instationary nonlinear heat equation with high accuracy. Optimal heating is obtained by minimizing an integral objective function which measures the distance between desired and model predicted temperatures. A sequence of minima is calculated from successively improved constant-rate perfusion models employing a damped Newton method in an inner iteration. We compare temperature distributions for two individual patients calculated on coarse and fine spatial grids and present numerical results of optimizations for a Sigma 60 Applicator of the BSD 2000 Hyperthermia System.

Abstract. Let u e H be the exact solution of a given selfadjoint elliptic boundary value problem, which is approximated by some S, S being a suitable finite-element space. Efficient and reliable a posteriod estimates of the error u II, measuring the (local) quality of , 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 by using domain decomposition techniques. In this 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.

Adaptive numerical methods in space and time are introduced and studied for multiscale cardiac reaction-diffusion models in three dimensions. The evolution of a complete heartbeat, from the excitation to the recovery phase, is simulated with both the anisotropic Bidomain and Monodomain models, coupled with either a variant of the simple FitzHugh-Nagumo model or the more complex phase-I Luo-Rudy ionic model. The simulations are performed with the kardos library, that employs adaptive finite elements in space and adaptive linearly implicit methods in time. The numerical results show that this adaptive method successfully solves these complex cardiac reaction-diffusion models on three-dimensional domains of moderate sizes. By automatically adapting the spatial meshes and time steps to the proper scales in each phase of the heartbeat, the method accurately resolves the evolution of the intra-and extra-cellular potentials, gating variables and ion concentrations during the excitation, plateau and recovery phases.Keywords: reaction-diffusion equations, cardiac Bidomain and Monodomain models, adaptive finite elements, adaptive time integration Recent advances in contemporary cardiac electrophysiology are progressively revealing the complex multiscale structure of the bioelectrical activity of the

SUMMARYWe consider the approximate solution of self-adjoint 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 deteriorates in three space dimensions, the BPX preconditioner is used, taking special care of an efficient implementation. Reliable o 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.

No abstract