TY - GEN A1 - Lehmann, Felix T1 - Inexaktheit in Newton-Lagrange-Verfahren für Optimierungsprobleme mit Elliptischen PDGL-Nebenbedingungen N2 - Bei der numerischen Lösung von Optimalsteuerungsproblemen mit elliptischen partiellen Differentialgleichungen als Nebenbedingung treten unvermeidlich Diskretisierungs- und Iterationsfehler auf. Man ist aus Aufwandsgründen daran interessiert die dabei entstehenden Fehler nicht sehr klein wählen zu müssen. In der Folge werden die linearisierten Nebenbedingungen in einem Composite-Step-Verfahren nicht exakt erfüllt. In dieser Arbeit wird der Einfluss dieser Ungenauigkeit auf das Konvergenzverhalten von Newton-Lagrange-Verfahren untersucht. Dabei sollen mehrere einschlägige lokale Konvergenzresultate diskutiert werden. Anschließend wird ein konkretes Composite-Step-Verfahren formuliert, in dem die Genauigkeit der inneren Iterationsverfahren adaptiv gesteuert werden kann. Am Ende der Arbeit wird an zwei Musterproblemen die hohe Übereinstimmung der analytischen Voraussagen und der tatsächlichen Performanz der dargestellten Methoden demonstriert. T3 - ZIB-Report - 13-37 KW - optimal control KW - PDE KW - inexact Newton Lagrange KW - adaptive Newton Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-41972 SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - Minion, Michael L. T1 - An Efficient Parallel-in-Time Method for Optimization with Parabolic PDEs N2 - To solve optimization problems with parabolic PDE constraints, often methods working on the reduced objective functional are used. They are computationally expensive due to the necessity of solving both the state equation and a backward-in-time adjoint equation to evaluate the reduced gradient in each iteration of the optimization method. In this study, we investigate the use of the parallel-in-time method PFASST in the setting of PDE constrained optimization. In order to develop an efficient fully time-parallel algorithm we discuss different options for applying PFASST to adjoint gradient computation, including the possibility of doing PFASST iterations on both the state and adjoint equations simultaneously. We also explore the additional gains in efficiency from reusing information from previous optimization iterations when solving each equation. Numerical results for both a linear and a non-linear reaction-diffusion optimal control problem demonstrate the parallel speedup and efficiency of different approaches. T3 - ZIB-Report - 19-04 KW - PDE-constrained optimization KW - parallel-in-time methods KW - PFASST Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-71490 SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - von Tycowicz, Christoph A1 - Polthier, Konrad A1 - Weiser, Martin T1 - Reducing Memory Requirements in Scientific Computing and Optimal Control N2 - In high accuracy numerical simulations and optimal control of time-dependent processes, often both many time steps and fine spatial discretizations are needed. Adjoint gradient computation, or post-processing of simulation results, requires the storage of the solution trajectories over the whole time, if necessary together with the adaptively refined spatial grids. In this paper we discuss various techniques to reduce the memory requirements, focusing first on the storage of the solution data, which typically are double precision floating point values. We highlight advantages and disadvantages of the different approaches. Moreover, we present an algorithm for the efficient storage of adaptively refined, hierarchic grids, and the integration with the compressed storage of solution data. T3 - ZIB-Report - 13-64 KW - optimal control KW - trajectory storage KW - mesh compression KW - compression Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42695 SN - 1438-0064 ER - TY - JOUR A1 - Powell, Gary A1 - Weiser, Martin T1 - Views, A new form of container adaptors JF - C/C++ Users Journal N2 - The C++ standard template library has many useful containers for data. The standard library includes two adpators, queue, and stack. The authors have extended this model along the lines of relational database semantics. Sometimes the analogy is striking, and we will point it out occasionally. An adaptor allows the standard algorithms to be used on a subset or modification of the data without having to copy the data elements into a new container. The authors provide many useful adaptors which can be used together to produce interesting views of data in a container. Y1 - 2000 VL - 18 IS - 4 SP - 40 EP - 51 ER - TY - GEN A1 - Sali, Adrian T1 - Coupling of Monodomain and Eikonal Models for Cardiac Electrophysiology N2 - The primary goal of this paper is to study the coupling of monodomain and eikonal models for the numerical simulation of cardiac electrophysiology. Eikonal models are nonlinear elliptic equations describing the excitation time of the cardiac tissue. They are often used as very fast approximations for monodomain or bidomain models - parabolic reaction-diffusion systems describing the excitation wavefront in terms of ionic currents. The excitation front is a thin region with high gradients, whereas excitation times vary over larger domains. Hence, eikonal equations can be solved on much coarser grids than monodomain equations. Moreover, as eikonal models are not time-dependent, no time integration is needed. Eikonal models are derived from monodomain models making additional assumptions and using certain approximations. While generally the approximation is rather good, several specific situations are not well captured by eikonal models. We consider coupling the two models, i.e. using the monodomain model in regions where more accurate results or the shape of the wavefront are needed, and the eikonal model in the remaining parts of the domain, where the excitation time is sufficient. Restricting the monodomain simulation to a small subdomain reduces the computational effort considerably. Numerical methods for the simulation of the individual models are presented, with the finite element method as the main ingredient. Coupling conditions as well as algorithms for implementing the coupling are explained. The approximation quality and efficiency of the coupled model is illustrated on simple geometries using an Aliev-Panfilov membrane model. T3 - ZIB-Report - 16-50 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60508 SN - 1438-0064 ER - TY - THES A1 - Sali, Adrian T1 - Coupling of Monodomain and Eikonal Models for Cardiac Electrophysiology N2 - The primary goal of this paper is to study the coupling of monodomain and eikonal models for the numerical simulation of cardiac electrophysiology. Eikonal models are nonlinear elliptic equations describing the excitation time of the cardiac tissue. They are often used as very fast approximations for monodomain or bidomain models - parabolic reaction-diffusion systems describing the excitation wavefront in terms of ionic currents. The excitation front is a thin region with high gradients, whereas excitation times vary over larger domains. Hence, eikonal equations can be solved on much coarser grids than monodomain equations. Moreover, as eikonal models are not time-dependent, no time integration is needed. Eikonal models are derived from monodomain models making additional assumptions and using certain approximations. While generally the approximation is rather good, several specific situations are not well captured by eikonal models. We consider coupling the two models, i.e. using the monodomain model in regions where more accurate results or the shape of the wavefront are needed, and the eikonal model in the remaining parts of the domain, where the excitation time is sufficient. Restricting the monodomain simulation to a small subdomain reduces the computational effort considerably. Numerical methods for the simulation of the individual models are presented, with the finite element method as the main ingredient. Coupling conditions as well as algorithms for implementing the coupling are explained. The approximation quality and efficiency of the coupled model is illustrated on simple geometries using an Aliev-Panfilov membrane model. Y1 - 2016 ER - TY - JOUR A1 - Weiser, Martin A1 - Ghosh, Sunayana T1 - Theoretically optimal inexact SDC methods JF - Communications in Applied Mathematics and Computational Science N2 - In several inital value problems with particularly expensive right hand side evaluation or implicit step computation, there is a trade-off between accuracy and computational effort. We consider inexact spectral deferred correction (SDC) methods for solving such initial value problems. SDC methods are interpreted as fixed point iterations and, due to their corrective iterative nature, allow to exploit the accuracy-work-tradeoff for a reduction of the total computational effort. On one hand we derive error models bounding the total error in terms of the evaluation errors. On the other hand, we define work models describing the computational effort in terms of the evaluation accuracy. Combining both, a theoretically optimal local tolerance selection is worked out by minimizing the total work subject to achieving the requested tolerance. The properties of optimal local tolerances and the predicted efficiency gain compared to simpler heuristics, and a reasonable practical performance, are illustrated on simple numerical examples. Y1 - 2018 U6 - https://doi.org/10.2140/camcos.2018.13.53 IS - 13-1 SP - 53 EP - 86 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Maierhofer, Christiane A1 - Richter, Regina T1 - Data Enhancement for Active Thermography N2 - Pulse thermography is a non-destructive testing method based on infrared imaging of transient thermal patterns. Heating the surface of the structure under test for a short period of time generates a non-stationary temperature distribution and thus a thermal contrast between the defect and the sound material. Due to measurement noise, preprocessing of the experimental data is necessary, before reconstruction algorithms can be applied. We propose a decomposition of the measured temperature into Green's function solutions to eliminate noise. T3 - ZIB-Report - 12-20 KW - data enhancement KW - active thermography KW - quantitative reconstruction KW - Green's functions Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-15243 UR - http://qirt.gel.ulaval.ca/archives/qirt2012/papers/QIRT-2012-167.pdf SN - 1438-0064 ER - TY - JOUR A1 - Hochmuth, Reinhard A1 - Deuflhard, Peter T1 - Multiscale Analysis for the Bio-Heat-Transfer Equation - The Nonisolated Case JF - Mathematical Models and Methods in Applied Sciences (M$^3$AS) Y1 - 2004 VL - 14 IS - 11 SP - 1621 EP - 1634 ER - TY - JOUR A1 - Heller, Markus O. A1 - Duda, Georg A1 - Ehrig, Rainald A1 - Schell, H. A1 - Seebeck, P. A1 - Taylor, William R. T1 - Muskuloskeletale Belastungen im Schafshinterlauf JF - Materialwissenschaft und Werkstofftechnik Y1 - 2005 IS - 36 SP - 775 EP - 780 ER - TY - CHAP A1 - Heller, Markus O. A1 - Duda, G. A1 - König, Christian A1 - Taylor, William R. A1 - Ehrig, Rainald A1 - Hinterwimmer, S. A1 - Graichen, H. T1 - Validation of a New Model of Human Knee Kinematics against in Vivo Data under Physiologic Neuromuscular Activation Patterns T2 - Proceedings of the 51th Annual Meeting of the Orthopedic Society Y1 - 2005 CY - Washington ER - TY - JOUR A1 - Hellmich, Christian A1 - Kober, Cornelia A1 - Erdmann, Bodo T1 - Micromechanics-based conversion of CT data into anisotropic elasticity tensors, applied to FE simulations of a mandible JF - Annals of Biomedical Engineering Y1 - 2008 U6 - https://doi.org/10.1007/s10439-007-9393-8 VL - 36 IS - 1 SP - 108 EP - 122 ER - TY - JOUR A1 - Heller, Markus O. A1 - Kratzenstein, Stefan A1 - Ehrig, Rainald A1 - Wassilew, Georgi A1 - Duda, Georg A1 - Taylor, William R. T1 - The weighted optimal common shape technique improves identification of the hip joint centre of rotation in vivo JF - Journal of Orthopedic Research Y1 - 2011 VL - 29 SP - 1470 EP - 5 ER - TY - JOUR A1 - Gellermann, Johanna A1 - Weihrauch, Mirko A1 - Cho, C. A1 - Wlodarczyk, Waldemar A1 - Fähling, Horst A1 - Felix, Roland A1 - Budach, Volker A1 - Weiser, Martin A1 - Nadobny, Johanna A1 - Wust, Peter T1 - Comparison of MR-thermography and planning calculations in phantoms JF - Medical Physics Y1 - 2006 VL - 33 SP - 3912 EP - 3920 ER - TY - JOUR A1 - Heller, Markus O. A1 - König, Christian A1 - Graichen, F. A1 - Hinterwimmer, S. A1 - Ehrig, Rainald A1 - Duda, G. A1 - Taylor, William R. T1 - A new model to predict in vivo human knee kinematics under physiological-like muscle activation JF - Journal of Biomechanics Y1 - 2007 SP - 45 EP - 53 ER - TY - RPRT A1 - Gladilin, Evgeny A1 - Peckar, W. A1 - Rohr, K. A1 - Stiehl, H. T1 - A Comparsion Between BEM and FEM for Elastic Registration of Medical Images Y1 - 1999 IS - FBI-HH-M-287/99 PB - FB Informatik CY - Hamburg, Germany ER - TY - JOUR A1 - Gänzler, Tobias A1 - Volkwein, S. A1 - Weiser, Martin T1 - SQP methods for parameter identification problems arising in hyperthermia JF - Optim. Methods Softw. Y1 - 2006 VL - 21 IS - 6 SP - 869 EP - 887 ER - TY - JOUR A1 - Ehrig, Rainald A1 - Heller, Markus O. A1 - Kratzenstein, Stefan A1 - Duda, Georg A1 - Trepczynski, Adam A1 - Taylor, William R. T1 - The SCoRE residual: A quality index to assess the accuracy of joint estimations JF - Journal of Biomechanics Y1 - 2011 U6 - https://doi.org/10.1016/j.jbiomech.2010.12.009 VL - 44 IS - 7 SP - 1400 EP - 1404 ER - TY - CHAP A1 - Erdmann, Bodo A1 - Kober, C. A1 - Lang, Jens A1 - Deuflhard, Peter A1 - Zeilhofer, J. A1 - Sader, Robert T1 - Efficient and Reliable Finite Element Methods for Simulation of the Human Mandible T2 - proceedings of 9th Workshop on The Finite Element Method in Biomedical Engineering, Biomechanics and Related Fields Y1 - 2002 CY - Ulm, Germany ER - TY - JOUR A1 - Erdmann, Bodo A1 - Ljubijankic, M. A1 - Nytsch-Geusen, C. T1 - Entwicklung einer Toolbox zur zwei- und dreidimensionalen instationären thermischen Bauteilsimulation mit KARDOS JF - Bauphsik Y1 - 2007 VL - 29 IS - 1 SP - 33 EP - 39 ER -