@misc{KoberErdmannSaderetal.2003, author = {Kober, Cornelia and Erdmann, Bodo and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Simulation of the Human Mandible: Comparison of Bone Mineral Density and Stress/Strain Profiles due to Masticatory Muscles' Traction}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7458}, number = {03-23}, year = {2003}, abstract = {The correlation of the inner architecture of bone and its functional loading was already stated by Wolff in 1892. Our objective is to demonstrate this interdependence in the case of the human mandible. For this purpose, stress/strain profiles occuring at a human lateral bite were simulated. Additionally, by a combination of computer graphics modules, a three--dimensional volumetric visualization of bone mineral density could be given. Qualitative correspondences between the density profile of the jaw and the simulated stress/strain profiles could be pointed out. In the long run, this might enable the use of the simulation for diagnosis and prognosis. The solution of the underlying partial differential equations describing linear elastic material behaviour was provided by an adaptive finite element method. Estimates of the discretization errors, local grid refinement, and multilevel techniques guaranteed the reliability and efficiency of the method.}, language = {en} } @misc{ErdmannLangRoitzsch1993, author = {Erdmann, Bodo and Lang, Jens and Roitzsch, Rainer}, title = {Kaskade Manual - Version 2.0.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4970}, number = {TR-93-05}, year = {1993}, abstract = {The KASKADE toolbox defines an interface to a set of C subroutines which can be used to implement adaptive multilevel Finite Element Methods solving systems of elliptic equations in two and three space dimensions. The manual contains the description of the data structures and subroutines. The main modules of the toolbox are a runtime environment, triangulation and node handling, assembling, direct and iterative solvers for the linear systems, error estimators, refinement strategies, and graphic utilities. Additionally, we included appendices on the basic command language interface, on file formats, and on the definition of the partial differential equations which can be solved. The software is available on the ZIB ftp--server {\tt elib} in the directory {\tt pub/kaskade}. TR 93--5 supersedes TR 89--4 and TR 89--05.}, language = {en} } @misc{BeckErdmannRoitzsch1995, author = {Beck, Rudolf and Erdmann, Bodo and Roitzsch, Rainer}, title = {KASKADE 3.0 User's Guide}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-5269}, number = {TR-95-11}, year = {1995}, abstract = {KASKADE 3.x was developed for the solution of partial differential equations in one, two, or three space dimensions. Its object-oriented implementation concept is based on the programming language C++\$\,\$.~Adaptive finite element techniques are employed to provide solution procedures of optimal computational complexity. This implies a posteriori error estimation, local mesh refinement and multilevel preconditioning. The program was designed both as a platform for further developments of adaptive multilevel codes and as a tool to tackle practical problems. Up to now we have implemented scalar problem types like stationary or transient heat conduction. The latter one is solved with the Rothe method, enabling adaptivity both in space and time. Some nonlinear phenomena like obstacle problems or two-phase Stefan problems are incorporated as well. Extensions to vector-valued functions and complex arithmetic are provided. This report helps to work with KASKADE Especially we \begin{itemize} \setlength{\parskip}{0ex} \item [{\bf --}] study a set of examples, \item [{\bf --}] explain how to define a user's problem and \item [{\bf --}] introduce a graphical user interface. \end{itemize} We are extending this guide continuously. The latest version is available by network.}, language = {en} } @misc{BeckErdmannRoitzsch1995, author = {Beck, Rudolf and Erdmann, Bodo and Roitzsch, Rainer}, title = {KASKADE 3.0 - An Object Oriented Adaptive Finite Element Code}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-5199}, number = {TR-95-04}, year = {1995}, abstract = {KASKADE 3.0 was developed for the solution of partial differential equations in one, two, or three space dimensions. Its object-oriented implementation concept is based on the programming language C++\$\,\$.~Adaptive finite element techniques are employed to provide solution procedures of optimal computational complexity. This implies a posteriori error estimation, local mesh refinement and multilevel preconditioning. The program was designed both as a platform for further developments of adaptive multilevel codes and as a tool to tackle practical problems. Up to now we have implemented scalar problem types like stationary or transient heat conduction. The latter one is solved with the Rothe method, enabling adaptivity both in space and time. Some nonlinear phenomena like obstacle problems or two-phase Stefan problems are incorporated as well. Extensions to vector-valued functions and complex arithmetic are provided. \%Such free boundary problems ... We have implemented several iterative solvers for both symmetric and unsymmetric systems together with multiplicative and additive multilevel preconditioners. Systems arising from the nonlinear problems can be solved with lately developed monotone multigrid methods. \%An object-oriented concept was chosen for KASKADE~3.0, based on the programming \%language C++\$\,\$. This should provide the desired extensibilty and clearly \%reflect the structure of the code. \%A direct sparse matrix solver (Harwell MA28) is included.}, language = {en} } @misc{BornemannErdmannRoitzsch1991, author = {Bornemann, Folkmar A. and Erdmann, Bodo and Roitzsch, Rainer}, title = {KASKADE - Numerical Experiments.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4762}, number = {TR-91-01}, year = {1991}, abstract = {The C-implementation of KASKADE, an adaptive solver for linear elliptic differential equations in 2D, is object of a set of numerical experiments to analyze the use of resources (time and memory) with respect to numerical accuracy. We study the dependency of the reliability, robustness, and efficiency of the program from the parameters controlling the algorithm.}, language = {en} } @misc{LangErdmannSeebass1998, author = {Lang, Jens and Erdmann, Bodo and Seebass, Martin}, title = {Impact of Nonlinear Heat Transfer on Temperature Control in Regional Hyperthermia}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3426}, number = {SC-97-73}, year = {1998}, abstract = {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 object 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.}, language = {en} } @misc{ErdmannLangSeebass1997, author = {Erdmann, Bodo and Lang, Jens and Seebass, Martin}, title = {Optimization of Temperature Distributions for Regional Hyperthermia Based on a Nonlinear Heat Transfer Model}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3283}, number = {SC-97-59}, year = {1997}, abstract = {We describe an optimization process specially designed for regional hyperthermia of deap 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. Optimal heating is obtained by minimizing an integral object function which measures the distance between desired and model predicted temperatures. Sequential minima are calculated from successively improved constant--rate perfusion models employing a damped Newton method in an inner iteration. Numerical results for a Sigma 60 applicator are presented. This work has been supported by Deutsche Forschungsgemeinschaft (DFG) within the Sonderforschungsbereich 273 \glqq Hyperthermie: Methodik und Klinik \grqq .}, language = {en} } @misc{KoberErdmannHellmichetal.2004, author = {Kober, Cornelia and Erdmann, Bodo and Hellmich, Christian and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Anisotropic Simulation of the Human Mandible}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7873}, number = {04-12}, year = {2004}, abstract = {We focus on the role of anisotropic elasticity in the simulation of the load distribution in a human mandible due to a lateral bite on the leftmost premolar. Based on experimental evidence, we adopt ``local''" orthotropy of the elastic properties of the bone tissue. Since the trajectories of anisotropic elasticity are not accessible from Computer Tomographic (CT) data, they will be reconstructed from (i) the organ's geometry and (ii) from coherent structures which can be recognized from the spatial distribution of the CT values. A sensitivity analysis comprising various 3D FE simulations reveals the relevance of elastic anisotropy for the load carrying behavior of a human mandible: Comparison of the load distributions in isotropic and anisotropic simulations indicates that anisotropy seems to ``spare''" the mandible from loading. Moreover, a maximum degree of anisotropy leads to kind of an load minimization of the mandible, expressed by a minimum of different norms of local strain, evaluated throughout the organ. Thus, we may suggest that anisotropy is not only relevant, but also in some sense ``optimal''.}, language = {en} } @misc{KoberErdmannLangetal.2004, author = {Kober, Cornelia and Erdmann, Bodo and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian}, title = {Adaptive Finite Element Simulation of the Human Mandible Using a New Physiological Model of the Masticatory Muscles}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7917}, number = {04-16}, year = {2004}, abstract = {Structural mechanics simulation of bony organs is of general medical and biomechanical interest, because of the interdependence of the inner architecture of bone and its functional loading already stated by Wolff in 1892. This work is part of a detailed research project concerning the human mandible. By adaptive finite element techniques, stress/strain profiles occurring in the bony structure under biting were simulated. Estimates of the discretization errors, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method. In general, our simulation requires a representation of the organ's geometry, an appropriate material description, and the load case due to teeth, muscle, or joint forces. In this paper, we want to focus on the influence of the masticatory system. Our goal is to capture the physiological situation as far as possible. By means of visualization techniques developed by the group, we are able to extract individual muscle fibres from computed tomography data. By a special algorithm, the fibres are expanded to fanlike (esp. for the musc. temporalis) coherent vector fields similar to the anatomical reality. The activity of the fibres can be adapted according to compartmentalisation of the muscles as measured by electromyological experiments. A refined sensitivity analysis proved remarkable impact of the presented approach on the simulation results.}, language = {en} } @article{AckermannErdmannRoitzsch1994, author = {Ackermann, J{\"o}rg and Erdmann, Bodo and Roitzsch, Rainer}, title = {A self-adaptive multilevel finite element method for the stationary Schr{\"o}dinger equation in three space dimensions}, volume = {101}, journal = {J. Chem. Phys.}, pages = {7643 -- 7650}, year = {1994}, language = {en} } @inproceedings{BeckErdmannRoitzsch1997, author = {Beck, Rudolf and Erdmann, Bodo and Roitzsch, Rainer}, title = {An Object-Oriented Adaptive Finite Element Code}, booktitle = {Modern Software Tools for Scientific Computing}, editor = {Arge, Erlend and Bruaset, Are Magnus and Langtangen, Hans Petter}, publisher = {Birkh{\"a}user, Boston}, pages = {105 -- 124}, year = {1997}, language = {en} } @article{ColliFranzoneDeuflhardErdmannetal.2006, author = {Colli Franzone, Piero and Deuflhard, Peter and Erdmann, Bodo and Lang, Jens and Pavarino, Luca Franco}, title = {Adaptivity in Space and Time for Reaction-Diffusion Systems in Electrocardiology}, volume = {28}, journal = {SIAM J. Sc. Comp.}, pages = {942 -- 962}, year = {2006}, language = {en} } @inproceedings{ZachowHierlErdmann2004, author = {Zachow, Stefan and Hierl, Thomas and Erdmann, Bodo}, title = {A quantitative evaluation of 3D soft tissue prediction in maxillofacial surgery planning}, booktitle = {Proc. 3. Jahrestagung der Deutschen Gesellschaft f{\"u}r Computer- und Roboter-assistierte Chirurgie e.V.}, address = {M{\"u}nchen}, year = {2004}, language = {en} } @inproceedings{ZachowHierlErdmann2004, author = {Zachow, Stefan and Hierl, Thomas and Erdmann, Bodo}, title = {On the Predictability of tissue changes after osteotomy planning in maxillofacial surgery}, booktitle = {Computer Assisted Radiology and Surgery (CARS)}, address = {Chicago, USA}, doi = {10.1016/j.ics.2004.03.043}, pages = {648 -- 653}, year = {2004}, language = {en} } @inproceedings{ZachowHierlErdmann2004, author = {Zachow, Stefan and Hierl, Thomas and Erdmann, Bodo}, title = {{\"U}ber die Genauigkeit einer 3D Weichgewebepr{\"a}diktion in der MKG-Cirurgie}, booktitle = {Workshop 'Bildverarbeitung f{\"u}r die Medizin' (BVM)}, address = {Berlin, Germany}, pages = {75 -- 79}, year = {2004}, language = {en} } @misc{WeiserErdmannSchenkletal.2017, author = {Weiser, Martin and Erdmann, Bodo and Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Mall, Gita and Zachow, Stefan}, title = {Uncertainty in Temperature-Based Determination of Time of Death}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63818}, year = {2017}, abstract = {Temperature-based estimation of time of death (ToD) can be per- formed either with the help of simple phenomenological models of corpse cooling or with detailed mechanistic (thermodynamic) heat transfer mod- els. The latter are much more complex, but allow a higher accuracy of ToD estimation as in principle all relevant cooling mechanisms can be taken into account. The potentially higher accuracy depends on the accuracy of tissue and environmental parameters as well as on the geometric resolution. We in- vestigate the impact of parameter variations and geometry representation on the estimated ToD based on a highly detailed 3D corpse model, that has been segmented and geometrically reconstructed from a computed to- mography (CT) data set, differentiating various organs and tissue types. From that we identify the most crucial parameters to measure or estimate, and obtain a local uncertainty quantifcation for the ToD.}, language = {en} } @misc{WeiserErdmannDeuflhard2010, author = {Weiser, Martin and Erdmann, Bodo and Deuflhard, Peter}, title = {On Efficiency and Accuracy in Cardioelectric Simulation}, journal = {Progress in Industrial Mathematics at ECMI 2008}, editor = {Wilson, E. and Fitt, A. and Ockendon, H. and Norbury, J.}, publisher = {Springer}, pages = {371 -- 376}, year = {2010}, abstract = {Reasons for the failure of adaptive methods to deliver improved efficiency when integrating monodomain models for myocardiac excitation are discussed. Two closely related techniques for reducing the computational complexity of linearly implicit integrators, deliberate sparsing and splitting, are investigated with respect to their impact on computing time and accuracy.}, language = {en} } @article{WeiserRoelligArndtetal.2010, author = {Weiser, Martin and R{\"o}llig, Mathias and Arndt, Ralf and Erdmann, Bodo}, title = {Development and test of a numerical model for pulse thermography in civil engineering}, volume = {46}, journal = {Heat and Mass Transfer}, number = {11-12}, pages = {1419 -- 1428}, year = {2010}, abstract = {Pulse thermography of concrete structures is used in civil engineering for detecting voids, honeycombing and delamination. The physical situation is readily modeled by Fourier's law. Despite the simplicity of the PDE structure, quantitatively realistic numerical 3D simulation faces two major obstacles. First, the short heating pulse induces a thin boundary layer at the heated surface which encapsulates all information and therefore has to be resolved faithfully. Even with adaptive mesh refinement techniques, obtaining useful accuracies requires an unsatisfactorily fine discretization. Second, bulk material parameters and boundary conditions are barely known exactly. We address both issues by a semi-analytic reformulation of the heat transport problem and by parameter identification. Numerical results are compared with measurements of test specimens.}, language = {en} } @misc{LangErdmann2000, author = {Lang, Jens and Erdmann, Bodo}, title = {Adaptive Linearly Implicit Methods for Heat and Mass Transfer Problems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-5892}, number = {00-21}, year = {2000}, abstract = {Dynamical process simulation of complex real-life problems often requires the use of modern algorithms, which automatically adapt both the time and space discretization in order to get error-controlled approximations of the solution. In this paper, a combination of linearly implicit time integrators of Rosenbrock type and adaptive multilevel finite elements based on a posteriori error estimates is presented. This approach has proven to work quite satisfactorily for a wide range of challenging practical problems. We show the performance of our adaptive method for two applications that arise in the study of flame balls and brine transport in porous media.}, language = {en} } @misc{ErdmannKoberLangetal.2001, author = {Erdmann, Bodo and Kober, Cornelia and Lang, Jens and Sader, Robert and Zeilhofer, Hans-Florian and Deuflhard, Peter}, title = {Efficient and Reliable Finite Element Methods for Simulation of the Human Mandible}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-6403}, number = {01-14}, year = {2001}, abstract = {By computed tomography data (CT), the individual geometry of the mandible is quite well reproduced, also the separation between cortical and trabecular bone. Using anatomical knowledge about the architecture and the functional potential of the masticatory muscles, realistic situations were approximated. The solution of the underlying partial differential equations describing linear elastic material behaviour is provided by an adaptive finite element method. Estimations of the discretization error, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method.}, language = {en} }