@article{HammerschmidtWeiserSantiagoetal., author = {Hammerschmidt, Martin and Weiser, Martin and Santiago, Xavier Garcia and Zschiedrich, Lin and Bodermann, Bernd and Burger, Sven}, title = {Quantifying parameter uncertainties in optical scatterometry using Bayesian inversion}, series = {Proc. SPIE}, volume = {10330}, journal = {Proc. SPIE}, doi = {10.1117/12.2270596}, pages = {1033004}, language = {en} } @misc{HammerschmidtWeiserSantiagoetal., author = {Hammerschmidt, Martin and Weiser, Martin and Santiago, Xavier Garcia and Zschiedrich, Lin and Bodermann, Bernd and Burger, Sven}, title = {Quantifying parameter uncertainties in optical scatterometry using Bayesian inversion}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64704}, abstract = {We present a Newton-like method to solve inverse problems and to quantify parameter uncertainties. We apply the method to parameter reconstruction in optical scatterometry, where we take into account a priori information and measurement uncertainties using a Bayesian approach. Further, we discuss the influence of numerical accuracy on the reconstruction result.}, language = {en} } @article{MoldenhauerWeiserZachow, author = {Moldenhauer, Marian and Weiser, Martin and Zachow, Stefan}, title = {Adaptive Algorithms for Optimal Hip Implant Positioning}, series = {PAMM}, volume = {17}, journal = {PAMM}, number = {1}, doi = {10.1002/pamm.201710071}, pages = {203 -- 204}, abstract = {In an aging society where the number of joint replacements rises, it is important to also increase the longevity of implants. In particular hip implants have a lifetime of at most 15 years. This derives primarily from pain due to implant migration, wear, inflammation, and dislocation, which is affected by the positioning of the implant during the surgery. Current joint replacement practice uses 2D software tools and relies on the experience of surgeons. Especially the 2D tools fail to take the patients' natural range of motion as well as stress distribution in the 3D joint induced by different daily motions into account. Optimizing the hip joint implant position for all possible parametrized motions under the constraint of a contact problem is prohibitively expensive as there are too many motions and every position change demands a recalculation of the contact problem. For the reduction of the computational effort, we use adaptive refinement on the parameter domain coupled with the interpolation method of Kriging. A coarse initial grid is to be locally refined using goal-oriented error estimation, reducing locally high variances. This approach will be combined with multi-grid optimization such that numerical errors are reduced.}, language = {en} } @article{LubkollSchielaWeiser, author = {Lubkoll, Lars and Schiela, Anton and Weiser, Martin}, title = {An affine covariant composite step method for optimization with PDEs as equality constraints}, series = {Optimization Methods and Software}, volume = {32}, journal = {Optimization Methods and Software}, number = {5}, doi = {10.1080/10556788.2016.1241783}, pages = {1132 -- 1161}, abstract = {We propose a composite step method, designed for equality constrained optimization with partial differential equations. Focus is laid on the construction of a globalization scheme, which is based on cubic regularization of the objective and an affine covariant damped Newton method for feasibility. We show finite termination of the inner loop and fast local convergence of the algorithm. We discuss preconditioning strategies for the iterative solution of the arising linear systems with projected conjugate gradient. Numerical results are shown for optimal control problems subject to a nonlinear heat equation and subject to nonlinear elastic equations arising from an implant design problem in craniofacial surgery.}, language = {en} } @misc{FischerGoetschelWeiser, author = {Fischer, Lisa and G{\"o}tschel, Sebastian and Weiser, Martin}, title = {Lossy data compression reduces communication time in hybrid time-parallel integrators}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63961}, abstract = {Parallel in time methods for solving initial value problems are a means to increase the parallelism of numerical simulations. Hybrid parareal schemes interleaving the parallel in time iteration with an iterative solution of the individual time steps are among the most efficient methods for general nonlinear problems. Despite the hiding of communication time behind computation, communication has in certain situations a significant impact on the total runtime. Here we present strict, yet no sharp, error bounds for hybrid parareal methods with inexact communication due to lossy data compression, and derive theoretical estimates of the impact of compression on parallel efficiency of the algorithms. These and some computational experiments suggest that compression is a viable method to make hybrid parareal schemes robust with respect to low bandwidth setups.}, 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} } @inproceedings{MuellerGoetschelWeiseretal., author = {M{\"u}ller, Jan and G{\"o}tschel, Sebastian and Weiser, Martin and Maierhofer, Christiane}, title = {Thermografie mit optimierter Anregung f{\"u}r die quantitative Untersuchung von Delaminationen in kohlenstofffaserverst{\"a}rkten Kunststoffen}, series = {NDT.net Proc. DGZfP 2017}, booktitle = {NDT.net Proc. DGZfP 2017}, abstract = {Da kohlenstofffaserverst{\"a}rkte Kunststoffe (CFK) in anspruchsvollen sicherheitsrelevanten Einsatzgebieten wie im Automobilbau und in der Luftfahrt eingesetzt werden, besteht ein zunehmender Bedarf an zerst{\"o}rungsfreien Pr{\"u}fmethoden. Ziel ist die Gew{\"a}hrleistung der Sicherheit und Zuverl{\"a}ssigkeit der eingesetzten Bauteile. Aktive Thermografieverfahren erm{\"o}glichen die effiziente Pr{\"u}fung großer Fl{\"a}chen mit hoher Aufl{\"o}sung in wenigen Arbeitsschritten. Ein wichtiges Teilgebiet der Pr{\"u}fungen ist die Ortung und Charakterisierung von Delaminationen, die sowohl bereits in der Fertigung als auch w{\"a}hrend der Nutzung eines Bauteils auftreten k{\"o}nnen, und dessen strukturelle Integrit{\"a}t schw{\"a}chen. ;In diesem Beitrag werden CFK-Strukturen mit k{\"u}nstlichen und nat{\"u}rlichen Delaminationen mit Hilfe unterschiedlich zeitlich modulierter Strahlungsquellen experimentell untersucht. Verwendet werden dabei Anregungen mit Blitzlampen und mit frequenzmodulierten Halogenlampen. Mittels Filterfunktionen im Zeit- und Frequenzbereich wird das Kontrast-zu-Rausch-Verh{\"a}ltnis (CNR) der detektierten Fehlstellen optimiert. Verglichen werden anschließend die Nachweisempfindlichkeit, das CNR und die Ortsaufl{\"o}sung der zu charakterisierenden Delaminationen f{\"u}r die unterschiedlichen Anregungs- und Auswertungstechniken. Erg{\"a}nzt werden die Experimente durch numerische Simulationen des dreidimensionalen W{\"a}rmetransportes.}, language = {de} } @inproceedings{MuellerGoetschelMaierhoferetal., author = {M{\"u}ller, Jan and G{\"o}tschel, Sebastian and Maierhofer, Christiane and Weiser, Martin}, title = {Determining the material parameters for the reconstruction of defects in carbon fiber reinforced polymers from data measured by flash thermography}, series = {AIP Conference Proceedings}, volume = {1806}, booktitle = {AIP Conference Proceedings}, doi = {10.1063/1.4974671}, language = {en} } @article{SchenklMuggenthalerHubigetal.2017, author = {Schenkl, Sebastian and Muggenthaler, Holger and Hubig, Michael and Erdmann, Bodo and Weiser, Martin and Zachow, Stefan and Heinrich, Andreas and G{\"u}ttler, Felix Victor and Teichgr{\"a}ber, Ulf and Mall, Gita}, title = {Automatic CT-based finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis}, series = {International Journal of Legal Medicine}, volume = {131}, journal = {International Journal of Legal Medicine}, number = {3}, doi = {doi:10.1007/s00414-016-1523-0}, pages = {699 -- 712}, year = {2017}, abstract = {Temperature based death time estimation is based either on simple phenomenological models of corpse cooling or on detailed physical heat transfer models. The latter are much more complex, but allow a higher accuracy of death time estimation as in principle all relevant cooling mechanisms can be taken into account. Here, a complete work flow for finite element based cooling simulation models is presented. The following steps are demonstrated on CT-phantoms: • CT-scan • Segmentation of the CT images for thermodynamically relevant features of individual geometries • Conversion of the segmentation result into a Finite Element (FE) simulation model • Computation of the model cooling curve • Calculation of the cooling time For the first time in FE-based cooling time estimation the steps from the CT image over segmentation to FE model generation are semi-automatically performed. The cooling time calculation results are compared to cooling measurements performed on the phantoms under controlled conditions. In this context, the method is validated using different CTphantoms. Some of the CT phantoms thermodynamic material parameters had to be experimentally determined via independent experiments. Moreover the impact of geometry and material parameter uncertainties on the estimated cooling time is investigated by a sensitivity analysis.}, language = {en} } @inproceedings{WeiserScacchi, author = {Weiser, Martin and Scacchi, Simone}, title = {Spectral Deferred Correction methods for adaptive electro-mechanical coupling in cardiac simulation}, series = {G. Russo et al.(eds.) Progress in Industrial Mathematics at ECMI 2014}, booktitle = {G. Russo et al.(eds.) Progress in Industrial Mathematics at ECMI 2014}, publisher = {Springer}, doi = {10.1007/978-3-319-23413-7_42}, pages = {321 -- 328}, abstract = {We investigate spectral deferred correction (SDC) methods for time stepping and their interplay with spatio-temporal adaptivity, applied to the solution of the cardiac electro-mechanical coupling model. This model consists of the Monodomain equations, a reaction-diffusion system modeling the cardiac bioelectrical activity, coupled with a quasi-static mechanical model describing the contraction and relaxation of the cardiac muscle. The numerical approximation of the cardiac electro-mechanical coupling is a challenging multiphysics problem, because it exhibits very different spatial and temporal scales. Therefore, spatio-temporal adaptivity is a promising approach to reduce the computational complexity. SDC methods are simple iterative methods for solving collocation systems. We exploit their flexibility for combining them in various ways with spatio-temporal adaptivity. The accuracy and computational complexity of the resulting methods are studied on some numerical examples.}, language = {en} }