TY - GEN A1 - Ehrig, Rainald A1 - Nowak, Ulrich A1 - Deuflhard, Peter T1 - Massively Parallel Linearly-Implicit Extrapolation Algorithms as a Powerful Tool in Process Simulation. N2 - We study the parallelization of linearly--implicit extrapolation codes for the solution of large scale PDE systems and differential algebraic equations on distributed memory machines. The main advantage of these algorithms is that they enable adapativity both in time and space. Additive Krylov--Schwarz methods yield high parallel perfomance for such extrapolation methods. Our approach combines a slightly overlapping domain decomposition together with a polynomial block Neumann preconditioner and a reduced system technique. Furthermore we get important advantages through the explicit computation of the matrix--products of the preconditioner and the matrix of the linear system. The parallel algorithms exhibit scalability up to 64 processors already for medium--sized test problems. We show that the codes are really efficient in large application systems for chemical engineering problems. T3 - ZIB-Report - SC-97-43 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3127 ER - TY - GEN A1 - Lipnikov, Konstantin A1 - Deuflhard, Peter T1 - Domain Decomposition with Subdomain CCG for Material Jump Elliptic Problems. N2 - A combination of the cascadic conjugate gradient (CCG) method for homogeneous problems with a non-overlapping domain decomposition (DD) method is studied. Mortar finite elements on interfaces are applied to permit non-matching grids in neighboring subdomains. For material jump problems, the method is designed as an alternative to the cascadic methods. T3 - ZIB-Report - SC-97-36 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3057 ER - TY - GEN A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - GMERR - an Error Minimizing Variant of GMRES N2 - The paper analyzes a recently proposed iterative error minimizing method for the solution of linear systems. Sufficient and necessary conditions for convergence are studied, which show that the method essentially requires normal matrices. An efficient implementation similar to GMRES has been worked out in detail. Numerical tests on general non--normal matrices, of course, indicate that this approach is not competitive with GMRES. Summarizing, if error minimizing is important, one should rather choose CGNE. A computational niche for GMERR might be problems, where normal but non--symmetric matrices occur, like dissipative quantum mechanics. T3 - ZIB-Report - SC-97-63 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3323 ER - TY - GEN A1 - Beck, Rudolf A1 - Deuflhard, Peter A1 - Hiptmair, Ralf A1 - Hoppe, Ronald H. W. A1 - Wohlmuth, Barbara T1 - Adaptive Multilevel Methods for Edge Element Discretizations of Maxwell's Equations N2 - The focus of this paper is on the efficient solution of boundary value problems involving the double-- curl operator. Those arise in the computation of electromagnetic fields in various settings, for instance when solving the electric or magnetic wave equation with implicit timestepping, when tackling time--harmonic problems or in the context of eddy--current computations. Their discretization is based on on N\'ed\'elec's {\bf H(curl}; $\Omega$)--conforming edge elements on unstructured grids. In order to capture local effects and to guarantee a prescribed accuracy of the approximate solution adaptive refinement of the grid controlled by a posteriori error estimators is employed. The hierarchy of meshes created through adaptive refinement forms the foundation for the fast iterative solution of the resulting linear systems by a multigrid method. The guiding principle underlying the design of both the error estimators and the multigrid method is the separate treatment of the kernel of the curl--operator and its orthogonal complement. Only on the latter we have proper ellipticity of the problem. Yet, exploiting the existence of computationally available discrete potentials for edge element spaces, we can switch to an elliptic problem in potential space to deal with nullspace of curl. Thus both cases become amenable to strategies of error estimation and multigrid solution developed for second order elliptic problems. The efficacy of the approach is confirmed by numerical experiments which cover several model problems and an application to waveguide simulation. T3 - ZIB-Report - SC-97-66 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3356 ER - TY - GEN A1 - Deuflhard, Peter A1 - Seebass, Martin T1 - Adaptive Multilevel FEM as Decisive Tools in the Clinical Cancer Therapy Hyperthermia N2 - 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). T3 - ZIB-Report - SC-98-30 KW - hyperthermia KW - Maxwell's equations KW - nonlinear heat transfer KW - finite elements Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3739 ER - TY - GEN A1 - Deuflhard, Peter A1 - Weiser, Martin A1 - Seebass, Martin T1 - A New Nonlinear Elliptic Multilevel FEM Applied to Regional Hyperthermia N2 - In the clinical cancer therapy of regional hyperthermia nonlinear perfusion effects inside and outside the tumor seem to play a not negligible role. A stationary model of such effects leads to a nonlinear Helmholtz term within an elliptic boundary value problem. The present paper reports about the application of a recently designed adaptive multilevel FEM to this problem. For several 3D virtual patients, nonlinear versus linear model is studied. Moreover, the numerical efficiency of the new algorithm is compared with a former application of an adaptive FEM to the corresponding instationary model PDE. T3 - ZIB-Report - SC-98-35 KW - hyperthermia KW - nonlinear elliptic KW - multilevel FEM Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3785 ER - TY - GEN A1 - Weiser, Martin A1 - Schiela, Anton A1 - Deuflhard, Peter T1 - Asymptotic Mesh Independence of Newton's Method Revisited N2 - The paper presents a new affine invariant theory on asymptotic mesh independence of Newton's method in nonlinear PDEs. Compared to earlier attempts, the new approach is both much simpler and more natural from the algorithmic point of view. The theory is exemplified at collocation methods for ODE boundary value problems and at finite element methods for elliptic PDE problems. T3 - ZIB-Report - 03-13 KW - mesh independence KW - nonlinear partial differential equations KW - Newton method KW - finite element method KW - collocation method Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7352 ER - TY - GEN A1 - Deuflhard, Peter T1 - A Comparison of Related Concepts in Computational Chemistry and Mathematics N2 - This article studies the relation of the two scientific languages Chemistry and Mathematics via three selected comparisons: (a) QSSA versus dynamic ILDM in reaction kinetics, (b) lumping versus discrete Galerkin methods in polymer chemistry, and (c) geometrical conformations versus metastable conformations in drug design. The common clear message from these comparisons is that chemical intuition may pave the way for mathematical concepts just as chemical concepts may gain from mathematical precisioning. Along this line, significant improvements in chemical research and engineering have already been possible -- and can be further expected in the future from the dialogue between the two scientific languages. T3 - ZIB-Report - 02-46 KW - dimension reduction KW - model reduction KW - QSSA KW - ILDM KW - reaction kinetics KW - polyreaction kinetics KW - Hamilton systems KW - conformations KW - metastability KW - d Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7138 ER - TY - GEN A1 - Deuflhard, Peter T1 - Adaptive Pseudo-transient Continuation for Nonlinear Steady State Problems N2 - Pseudo--transient continuation methods are quite popular for the numerical solution of steady state problems, typically in PDEs. They are based on an embedding into a time dependent initial value problem. In the presence of dynamical invariants the Jacobian matrix of the nonlinear equation system is bound to be singular. The paper presents a convergence analysis which takes this property into account -- in contrast to known approaches. On the basis of the new analysis adaptive algorithms are suggested in detail. These include a variant with Jacobian approximations as well as inexact pseudo--transient continuation, both of which play an important role in discretized PDEs. Numerical experiments are left to future work. T3 - ZIB-Report - 02-14 KW - pseudo--transient continuation KW - linearly implicit KW - Euler discretization KW - stiff integration KW - contractivity of ordinary differential equations KW - lar Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-6814 ER - TY - GEN A1 - Deuflhard, Peter T1 - From Molecular Dynamics to Conformational Dynamics in Drug Design N2 - Computational drug design studies molecular recognition in the {\em virtual lab}. The arising Hamiltonian dynamics is known to be chaotic and ill-conditioned already after picoseconds, whereas times are $msec$ up to $min$. Classical molecular dynamics with long term trajectory computation gives, at best, information about time and statistical ensemble averages. The present paper surveys a recent new modeling approach called {\em conformational dynamics}, which is due to the author and Ch. Schütte. This approach achieves information about the dy time scales by telescoping a short term deterministic model with a statistical model. Examples of small biomolecules are included. T3 - ZIB-Report - 02-20 KW - molecular dynamics KW - conformational dynamics KW - drug design KW - Hamiltonian dynamics KW - almost invariant sets KW - transition operator KW - Markov chain KW - Perro Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-6878 ER - TY - GEN A1 - Hochmuth, Reinhard A1 - Deuflhard, Peter T1 - Multiscale Analysis for the Bio-Heat Transfer Equation N2 - The bio-heat transfer equation is a macroscopic model for describing the heat transfer in microvascular tissue. In [{\sl Deuflhard, Hochmuth 2002}] the authors applied homogenization techniques to derive the bio-heat transfer equation as asymptotic result of boundary value problems which provide a microscopic description for microvascular tissue. Here those results are generalized to a geometrical setting where the regions of blood are allowed to be connected. Moreover, asymptotic corrector results are derived. T3 - ZIB-Report - 03-08 KW - bio-heat equation KW - hyperthermia KW - homogenization KW - correctors KW - heat transfer KW - Robin boundary conditions Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7302 ER - TY - GEN A1 - Deuflhard, Peter A1 - Hochmuth, Reinhard T1 - Multiscale Analysis of Thermoregulation in the Human Microvascular System N2 - The bio-heat transfer equation is a macroscopic model for describing the heat transfer in microvascular tissue. So far the deduction of the Helmholtz term in the bio-heat transfer equation is not co role. In view of a future numerical application of this new mathematical model to treatment planning in hyperthermia we derive asymptotic estimates for first and second order correctors. T3 - ZIB-Report - 02-31 KW - bio-heat equation KW - hyperthermia KW - homogenization Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-6985 ER - TY - GEN A1 - Deuflhard, Peter A1 - Nowak, Ulrich A1 - Weiser, Martin T1 - Affine Invariant Adaptive Newton Codes for Discretized PDEs N2 - The paper deals with three different Newton algorithms that have recently been worked out in the general frame of affine invariance. Of particular interest is their performance in the numerical solution of discretized boundary value problems (BVPs) for nonlinear partial differential equations (PDEs). Exact Newton methods, where the arising linear systems are solved by direct elimination, and inexact Newton methods, where an inner iteration is used instead, are synoptically presented, both in affine invariant convergence theory and in numerical experiments. The three types of algorithms are: (a) affine covariant (formerly just called affine invariant) Newton algorithms, oriented toward the iterative errors, (b) affine contravariant Newton algorithms, based on iterative residual norms, and (c) affine conjugate Newton algorithms for convex optimization problems and discrete nonlinear elliptic PDEs. T3 - ZIB-Report - 02-33 KW - Affine invariant Newton methods KW - global Newton methods KW - inexact Newton methods KW - adaptive trust region methods KW - nonlinear partial differential equa Y1 - 2002 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7005 ER - TY - GEN A1 - Deuflhard, Peter A1 - Weber, Marcus T1 - Robust Perron Cluster Analysis in Conformation Dynamics N2 - The key to molecular conformation dynamics is the direct identification of metastable conformations, which are almost invariant sets of molecular dynamical systems. Once some reversible Markov operator has been discretized, a generalized symmetric stochastic matrix arises. This matrix can be treated by Perron cluster analysis, a rather recent method involving a Perron cluster eigenproblem. The paper presents an improved Perron cluster analysis algorithm, which is more robust than earlier suggestions. Numerical examples are included. T3 - ZIB-Report - 03-19 KW - Markov chains KW - cluster algorithms KW - Perron cluster analysis KW - conformation dynamics Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7415 ER - TY - GEN A1 - Deuflhard, Peter A1 - Schütte, Christof T1 - Molecular Conformation Dynamics and Computational Drug Design N2 - The paper surveys recent progress in the mathematical modelling and simulation of essential molecular dynamics. Particular emphasis is put on computational drug design wherein time scales of $msec$ up to $min$ play the dominant role. Classical long-term molecular dynamics computations, however, would run into ill-conditioned initial value problems already after time spans of only $psec=10^{-12} sec$. Therefore, in order to obtain results for times of pharmaceutical interest, a combined deterministic-stochastic model is needed. The concept advocated in this paper is the direct identification of metastable conformations together with their life times and their transition patterns. It can be interpreted as a {\em transfer operator} approach corresponding to some underlying hybrid Monte Carlo process, wherein short-term trajectories enter. Once this operator has been discretized, which is a hard problem of its own, a stochastic matrix arises. This matrix is then treated by {\em Perron cluster analysis}, a recently developed cluster analysis method involving the numerical solution of an eigenproblem for a Perron cluster of eigenvalues. In order to avoid the 'curse of dimension', the construction of appropriate boxes for the spatial discretization of the Markov operator requires careful consideration. As a biomolecular example we present a rather recent SARS protease inhibitor. T3 - ZIB-Report - 03-20 KW - conformation dynamics KW - Monte Carlo methods KW - transfer operations KW - Hamiltonian dynamics KW - Smoluchowski dynamics KW - metastable sets KW - Perron cluster an Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7427 ER - TY - GEN A1 - Deuflhard, Peter A1 - Nowak, Ulrich A1 - Pöhle, Uwe A1 - Schmidt, B. Ch. A1 - Weyer, J. T1 - Die Ausbreitung von HIV/AIDS in Ballungsgebieten. N2 - The rapidly increasing number of AIDS cases requires a realistic estimation of the future development of the HIV/AIDS disease. For that purpose we develop a large system of coupled nonlinear differential equations describing simultaneously the dynamics of the development of the disease, the population size, the gender and age structure. A set of 1650 coupled equations are linked by balanced parameters. The balancing procedure is described by a set of (formally) 2,178,000 additional algebraic conditions. As the considered system is stiff, it requires new special extrapolation methods combined with techniquees of dynamical sparsing for the solution of sparsely filled systems. According to our simulations we expect 19,0,000 deaths caused by AIDS in the Federal Republic of Germany (former territories) in the year 2000. Such an epidemical spread would tie up about 4-7 percent of the actual health care workers. T3 - ZIB-Report - TR-91-09 Y1 - 1991 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4844 ER - TY - GEN A1 - Ehrig, Rainald A1 - Nowak, Ulrich A1 - Deuflhard, Peter T1 - Highly scalable parallel linearly-implicit extrapolation algorithms N2 - We present parallel formulations of the well established extrapolation algorithms EULSIM and LIMEX and its implementation on a distributed memory architecture. The discretization of partial differential equations by the method of lines yields large banded systems, which can be efficiently solved in parallel only by iterative methods. Polynomial preconditioning with a Neumann series expansion combined with an overlapping domain decomposition appears as a very efficient, robust and highly scalable preconditioner for different iterative solvers. A further advantage of this preconditioner is that all computation can be restricted to the overlap region as long as the subdomain problems are solved exactly. With this approach the iterative algorithms operate on very short vectors, the length of the vectors depends only on the number of gridpoints in the overlap region and the number of processors, but not on the size of the linear system. As the most reliable and fast iterative methods based on this preconditioning scheme appeared GMRES or FOM and BICGSTAB. To further reduce the number of iterations in GMRES or FOM we can reuse the Krylov-spaces constructed in preceeding extrapolation steps. The implementation of the method within the program LIMEX results in a highly parallel and scalable program for solving differential algebraic problems getting an almost linear speedup up to 64 processors even for medium size problems. Results are presented for a difficult application from chemical engineering simulating the formation of aerosols in industrial gas exhaust purification. T3 - ZIB-Report - TR-96-11 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5405 ER - TY - GEN A1 - Deuflhard, Peter A1 - Bölling, Reinhard A1 - Bott, Raoul H. T1 - Materialien zur Euler-Vorlesung in Sanssouci. N2 - Aus dem Inhalt: Vorwort; Leonhard Euler - aus der Zeit seines Wirkens in Berlin; Introductory Remarks. T3 - ZIB-Report - TR-93-10 Y1 - 1993 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5029 ER - TY - GEN A1 - Deuflhard, Peter A1 - Hege, Hans-Christian A1 - Seebass, Martin T1 - Progress Towards a Combined MRI/Hyperthermia System N2 - Regional hyperthermia, a clinical cancer therapy, is the main topic of the Sonderforschungsbereich Hyperthermia: Scientific Methods and Clinical Applications'' at Berlin. In recent years, technological improvements towards a better concentration of heat to the desired target region have been achieved. These include a rather sophisticated integrated software environment for therapy planning and a new hyperthermia applicator. In a next step, a detailed closed loop monitoring of the actual treatment is to be developed. For this purpose the hyperthermia applicator is combined with an MRI system, which will allow to check the positioning of the patients and to measure individual blood perfusion as well as the 3D temperature distribution. The measurements will then be employed for an on-line control of the whole treatment. In this intended setting, new fast feedback control algorithms will come into play. T3 - ZIB-Report - 00-07 KW - hyperthermia KW - medical therapy planning KW - applicator design KW - interventional MRI Y1 - 2000 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5755 ER - TY - GEN A1 - Deuflhard, Peter A1 - Reich, Sebastian T1 - 2nd International Symposium Algorithms for Macromolecular Modelling N2 - The report contains the program and the collected abstracts of the 2nd International Symposium ``Algorithms for Macromolecular Modelling''. T3 - ZIB-Report - TR-97-03 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5455 ER - TY - GEN A1 - Deuflhard, Peter A1 - Louis, Alfred A1 - Schlegel, W. A1 - Seebass, Martin T1 - Workshop "Scientific Computing in der Medizin" SCMED '97 N2 - Die Arbeitstagung \glqq Scientific Computing in der Medizin\grqq , kurz SCMED \grq 97, findet vom 22. - 23. September 1997 im neu errichteten Gebäude des Konrad-Zuse-Zentrums auf dem Dahlemer naturwissenschaftlichen Campus der Freien Universtät Berlin statt. T3 - ZIB-Report - TR-97-05 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5474 ER -