@misc{BornemannYserentant, author = {Bornemann, Folkmar A. and Yserentant, Harry}, title = {A Basic Norm Equivalence for the Theory of Multilevel Methods.}, doi = {10.1007/BF01388699}, number = {SC-92-01}, abstract = {Subspace decompositions of finite element spaces based on \$L2\$-like orthogonal projections play an important role for the construction and analysis of multigrid like iterative methods. Recently several authors proved the equivalence of the associated discrete norms with the \$H^1\$-norm. The present report gives an elementary, self-contained derivation of this result which is based on the use of \$ K\$-functionals known from the theory of interpolation spaces. {\bf Keywords:} multilevel methods, nonuniform meshes, optimal convergence rates. {\bf AMS(MOS) Subject classifications:} 65N55, 65N30, 65N50.}, language = {en} } @misc{OevermannKlein, author = {Oevermann, Michael and Klein, Rupert}, title = {A cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-8996}, number = {06-05}, abstract = {We present a finite volume method for the solution of the two-dimensional Poisson equation \$ \nabla\cdot( \beta( {\mbox{\boldmath \$x\$}}) \nabla u({\mbox{\boldmath \$x\$}})) = f(\mbox{\boldmath \$x\$}) \$ with variable, discontinuous coefficients and solution discontinuities on irregular domains. The method uses bilinear ansatz functions on Cartesian grids for the solution \$u({\mbox{\boldmath \$x\$})\$ resulting in a compact nine-point stencil. The resulting linear problem has been solved with a standard multigrid solver. Singularities associated with vanishing partial volumes of intersected grid cells or the dual bilinear ansatz itself are removed by a two-step asymptotic approach. The method achieves second order of accuracy in the \$L^\infty\$ and \$L^2\$ norm.}, language = {en} } @phdthesis{Schmidt, author = {Schmidt, Frank}, title = {A New Approach to Coupled Interior-Exterior Helmholtz-Type Problems: Theory and Algorithms}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7128}, number = {02-45}, abstract = {The work presents a new approach to the numerical solution of time-harmonic and time-dependent scattering problems. We replace Sommerfeld's radiation condition valid for the Helmholtz equation by a more general concept called pole condition. The pole condition is based on the Laplace transform of the exterior solution and allows a characterization of outgoing waves. Both new insight into the analysis of scattering problems as well as new numerical algorithms are obtained.}, language = {en} } @misc{ZschiedrichKloseSchaedleetal., author = {Zschiedrich, Lin and Klose, Roland and Sch{\"a}dle, Achim and Schmidt, Frank}, title = {A new Finite Element realization of the Perfectly Matched Layer Method for Helmholtz scattering problems on polygonal domains in 2D}, doi = {10.1016/j.cam.2005.03.047}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7662}, number = {03-44}, abstract = {In this paper we propose a new finite element realization of the Perfectly Matched Layer method (PML-method). Our approach allows to deal with arbitrary shaped polygonal domains and with certain types of inhomogeneous exterior domains. Among the covered inhomogeneities are open waveguide structures playing an essential role in integrated optics. We give a detailed insight to implementation aspects. Numerical examples show exponential convergence behavior to the exact solution with the thickness of the PML sponge layer.}, language = {en} } @misc{OevermannScharfenbergKlein, author = {Oevermann, Michael and Scharfenberg, Carsten and Klein, Rupert}, title = {A sharp interface finite volume method for elliptic equations on Cartesian grids}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-10900}, number = {08-38}, abstract = {We present a second order sharp interface finite volume method for the solution of the three-dimensional poisson equation with variable coefficients on Cartesian grids. In particular, we focus on interface problems with discontinuities in the coefficient, the source term, the solution, and the fluxes across the interface. The method uses standard piecewiese trilinear finite elements for normal cells and a double piecewise trilinear ansatz for the solution on cells intersected by the interface resulting always in a compact 27-point stencil. Singularities associated with vanishing partial volumes of intersected grid cells are removed by a two-term asymptotic approach. In contrast to the 2D method presented by two of the authors in [M.~Oevermann, R.~Klein: A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, J.~Comp.~Phys.~219 (2006)] we use a minimization technique to determine the unknown coefficients of the double trilinear ansatz. This simplifies the treatment of the different cut-cell types and avoids additional special operations for degenerated interface topologies. The resulting set of linear equations has been solved with a BiCGSTAB solver preconditioned with an algebraic multigrid. In various testcases -- including large coefficient ratios and non-smooth interfaces -- the method achieves second order of accuracy in the L_inf and L_2 norm.}, language = {en} } @misc{Bornemann, author = {Bornemann, Folkmar A.}, title = {A Sharpened Condition Number Estimate for the BPX Preconditioner of Elliptic Finite Element Problems on Highly Nonuniform Triangulations.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-596}, number = {SC-91-09}, abstract = {In this paper it is shown that for highly nonuniformly refined triangulations the condition number of the BPX preconditioner for elliptic finite element problems grows at most linearly in the depth of refinement. This is achieved by viewing the computational available version of the BPX preconditioner as an abstract additive Schwarz method with exact solvers. {\bf AMS CLASSIFICATION:} 65F10, 65F35, 65N20, 65N30.}, language = {en} } @misc{BraessDeuflhardLipnikov, author = {Braess, Dietrich and Deuflhard, Peter and Lipnikov, Konstantin}, title = {A Subspace Cascadic Multigrid Method for Mortar Elements}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3954}, number = {SC-99-07}, abstract = {A cascadic multigrid (CMG) method for elliptic problems with strong material jumps is proposed and analyzed. Non--matching grids at interfaces between subdomains are allowed and treated by mortar elements. The arising saddle point problems are solved by a subspace confined conjugate gradient method as smoother for the CMG. Details of algorithmic realization including adaptivity are elaborated. Numerical results illustrate the efficiency of this CMG algorithm.}, language = {en} } @misc{SchmidtFrieseZschiedrichetal., author = {Schmidt, Frank and Friese, Tilmann and Zschiedrich, Lin and Deuflhard, Peter}, title = {Adaptive Multigrid Methods for the Vectorial Maxwell Eigenvalue Problem for Optical Waveguide Design}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6228}, number = {00-54}, abstract = {This paper has been motivated by the need for a fast robust adaptive multigrid method to solve the vectorial Maxwell eigenvalue problem arising from the design of optical chips. Our nonlinear multigrid methods are based on a previous method for the scalar Helmholtz equation, which must be modified to cope with the null space of the Maxwell operator due to the divergence condition. We present two different approaches. First, we present a multigrid algorithm based on an edge element discretization of time-harmonic Maxwell's equations, including the divergence condition. Second, an explicit elimination of longitudinal magnetic components leads to a nodal discretization known to avoid discrete \emph{spurious modes} also and a vectorial eigenvalue problem, for which we present a multigrid solver. Numerical examples show that the edge element discretization clearly outperforms the nodal element approach.}, language = {en} } @misc{DeuflhardSeebass, author = {Deuflhard, Peter and Seebass, Martin}, title = {Adaptive Multilevel FEM as Decisive Tools in the Clinical Cancer Therapy Hyperthermia}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3739}, number = {SC-98-30}, abstract = {The paper surveys recent progress in a joint mathematical-medical project on cancer therapy planning. Within so-called regional hyperthermia the computational task is to tune a set of coupled radiofrequency antennas such that a carefully measured tumor is locally heated, but any outside hot spots are avoided. A mathematical model of the whole clinical system -- air, applicator with antennas, water bolus, individual patient body -- involves Maxwell's equations in inhomogeneous media and a parabolic bioheat transfer equation, which represents a simplified model of heat transfer in the human body (ignoring strong blood vessel heat transport). Both PDEs need to be computed fast and to medical reliability (!) on a workstation within a clinical environment. This requirement triggered a series of new algorithmic developments to be reported here, among which is an adaptive multilevel FEM for Maxwell's equations, which dominates the numerical simulation time. In total, however, the main bulk of computation time (see Table 3 in Section 4 below) still goes into segmentation -- a necessary preprocessing step in the construction a 3D virtual patient from the input of a stack of 2D computed tomograms (left out here).}, language = {en} } @misc{WeiserDeuflhardErdmann, author = {Weiser, Martin and Deuflhard, Peter and Erdmann, Bodo}, title = {Affine conjugate adaptive Newton methods for nonlinear elastomechanics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7768}, number = {04-01}, abstract = {The paper extends affine conjugate Newton methods from convex to nonconvex minimization, with particular emphasis on PDE problems originating from compressible hyperelasticity. Based on well-known schemes from finite dimensional nonlinear optimization, three different algorithmic variants are worked out in a function space setting, which permits an adaptive multilevel finite element implementation. These algorithms are tested on two well-known 3D test problems and a real-life example from surgical operation planning.}, language = {en} }