TY - GEN A1 - Weiser, Martin A1 - Götschel, Sebastian T1 - State Trajectory Compression for Optimal Control with Parabolic PDEs N2 - In optimal control problems with nonlinear time-dependent 3D PDEs, full 4D discretizations are usually prohibitive due to the storage requirement. For this reason gradient and quasi-Newton methods working on the reduced functional are often employed. The computation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, again requiring the storage of a full 4D data set. We propose a lossy compression algorithm using an inexact but cheap predictor for the state data, with additional entropy coding of prediction errors. As the data is used inside a discretized, iterative algorithm, lossy coding maintaining an error bound is sufficient. T3 - ZIB-Report - 10-05 KW - optimal control KW - adjoint gradient computation KW - trajectory storage Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11676 SN - 1438-0064 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 - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Schiela, Anton T1 - Solving Optimal Control Problems with the Kaskade 7 Finite Element Toolbox N2 - This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study. T3 - ZIB-Report - 10-25 KW - partial differential equations KW - optimal control KW - finite elements KW - generic programming KW - adaptive methods Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11909 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 - GEN A1 - Götschel, Sebastian A1 - Minion, Michael L. T1 - Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST N2 - In gradient-based methods for parabolic optimal control problems, it is necessary to solve both the state equation and a backward-in-time adjoint equation in each iteration of the optimization method. In order to facilitate fully parallel gradient-type and nonlinear conjugate gradient methods for the solution of such optimal control problems, we discuss the application of the parallel-in-time method PFASST to adjoint gradient computation. In addition to enabling time parallelism, PFASST provides high flexibility for handling nonlinear equations, as well as potential extra computational savings from reusing previous solutions in the optimization loop. The approach is demonstrated here for a model reaction-diffusion optimal control problem. T3 - ZIB-Report - 17-51 Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-64989 SN - 1438-0064 ER - TY - JOUR A1 - Fischer, Lisa A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Lossy data compression reduces communication time in hybrid time-parallel integrators JF - Comput. Vis. Sci. N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1007/s00791-018-0293-2 VL - 19 IS - 1 SP - 19 EP - 30 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Maierhofer, Christiane A1 - Richter, Regina ED - Cardone, Gennaro T1 - Data Enhancement for Active Thermography T2 - E-book Proceedings, 11th International Conference on Quantitative Infrared Thermography, Naples 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. Y1 - 2012 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Schiela, Anton ED - Dedner, A. ED - Flemisch, B. ED - Klöfkorn, R. T1 - Solving Optimal Control Problems with the Kaskade 7 Finite Element Toolbox T2 - Advances in DUNE N2 - This paper presents concepts and implementation of the finite element toolbox Kaskade 7, a flexible C++ code for solving elliptic and parabolic PDE systems. Issues such as problem formulation, assembly and adaptivity are discussed at the example of optimal control problems. Trajectory compression for parabolic optimization problems is considered as a case study. Y1 - 2012 SP - 101 EP - 112 PB - Springer ER - TY - JOUR A1 - Weiser, Martin A1 - Götschel, Sebastian T1 - State Trajectory Compression for Optimal Control with Parabolic PDEs JF - SIAM J. Sci. Comput. N2 - In optimal control problems with nonlinear time-dependent 3D PDEs, full 4D discretizations are usually prohibitive due to the storage requirement. For this reason gradient and quasi-Newton methods working on the reduced functional are often employed. The computation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, again requiring the storage of a full 4D data set. We propose a lossy compression algorithm using an inexact but cheap predictor for the state data, with additional entropy coding of prediction errors. As the data is used inside a discretized, iterative algorithm, lossy coding maintaining an error bound is sufficient. Y1 - 2012 U6 - https://doi.org/10.1137/11082172X VL - 34 IS - 1 SP - A161 EP - A184 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - State Trajectory Compression in Optimal Control JF - PAMM N2 - In optimal control problems with nonlinear time-dependent 3D PDEs, the computation of the reduced gradient by adjoint methods requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. Since the state enters into the adjoint equation, the storage of a 4D discretization is necessary. We propose a lossy compression algorithm using a cheap predictor for the state data, with additional entropy coding of prediction errors. Analytical and numerical results indicate that compression factors around 30 can be obtained without exceeding the FE discretization error. Y1 - 2010 U6 - https://doi.org/10.1002/pamm.201010282 VL - 10 IS - 1 SP - 579 EP - 580 ER - TY - GEN A1 - Fischer, Lisa A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Lossy data compression reduces communication time in hybrid time-parallel integrators N2 - 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. T3 - ZIB-Report - 17-25 KW - parallel-in-time integration KW - hybrid parareal KW - convergence KW - lossy compression KW - MPI Y1 - 2017 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-63961 SN - 1438-0064 ER - TY - CHAP A1 - Müller, Jan A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Maierhofer, Christiane T1 - Thermografie mit optimierter Anregung für die quantitative Untersuchung von Delaminationen in kohlenstofffaserverstärkten Kunststoffen T2 - NDT.net Proc. DGZfP 2017 N2 - Da kohlenstofffaserverstärkte Kunststoffe (CFK) in anspruchsvollen sicherheitsrelevanten Einsatzgebieten wie im Automobilbau und in der Luftfahrt eingesetzt werden, besteht ein zunehmender Bedarf an zerstörungsfreien Prüfmethoden. Ziel ist die Gewährleistung der Sicherheit und Zuverlässigkeit der eingesetzten Bauteile. Aktive Thermografieverfahren ermöglichen die effiziente Prüfung großer Flächen mit hoher Auflösung in wenigen Arbeitsschritten. Ein wichtiges Teilgebiet der Prüfungen ist die Ortung und Charakterisierung von Delaminationen, die sowohl bereits in der Fertigung als auch während der Nutzung eines Bauteils auftreten können, und dessen strukturelle Integrität schwächen. ;In diesem Beitrag werden CFK-Strukturen mit künstlichen und natü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ältnis (CNR) der detektierten Fehlstellen optimiert. Verglichen werden anschließend die Nachweisempfindlichkeit, das CNR und die Ortsauflösung der zu charakterisierenden Delaminationen für die unterschiedlichen Anregungs- und Auswertungstechniken. Ergänzt werden die Experimente durch numerische Simulationen des dreidimensionalen Wärmetransportes. Y1 - 2017 ER - TY - GEN A1 - Götschel, Sebastian A1 - Minion, Michael L. T1 - Parallel-in-Time for Parabolic Optimal Control Problems Using PFASST BT - DD 2017. T2 - Domain Decomposition Methods in Science and Engineering XXIV N2 - In gradient-based methods for parabolic optimal control problems, it is necessary to solve both the state equation and a backward-in-time adjoint equation in each iteration of the optimization method. In order to facilitate fully parallel gradient-type and nonlinear conjugate gradient methods for the solution of such optimal control problems, we discuss the application of the parallel-in-time method PFASST to adjoint gradient computation. In addition to enabling time parallelism, PFASST provides high flexibility for handling nonlinear equations, as well as potential extra computational savings from reusing previous solutions in the optimization loop. The approach is demonstrated here for a model reaction-diffusion optimal control problem. Y1 - 2018 U6 - https://doi.org/10.1007/978-3-319-93873-8_34 SP - 363 EP - 371 PB - Springer CY - Cham ER - TY - CHAP A1 - Götschel, Sebastian A1 - Höhne, Christian A1 - Kolkoori, Sanjeevareddy A1 - Mitzscherling, Steffen A1 - Prager, Jens A1 - Weiser, Martin T1 - Ray Tracing Boundary Value Problems: Simulation and SAFT Reconstruction for Ultrasonic Testing T2 - Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016) Y1 - 2016 UR - https://www.wcndt2016.com/portals/wcndt/bb/Fr1H4.pdf ER - TY - CHAP A1 - Götschel, Sebastian A1 - Maierhofer, Christiane A1 - Müller, Jan A1 - Rothbart, Nick A1 - Weiser, Martin T1 - Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites T2 - Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016) Y1 - 2016 UR - https://www.wcndt2016.com/portals/wcndt/bb/Th4C4.pdf ER - TY - CHAP A1 - Müller, Jan A1 - Götschel, Sebastian A1 - Maierhofer, Christiane A1 - Weiser, Martin T1 - Determining the material parameters for the reconstruction of defects in carbon fiber reinforced polymers from data measured by flash thermography T2 - AIP Conference Proceedings Y1 - 2017 U6 - https://doi.org/10.1063/1.4974671 VL - 1806 ER - TY - GEN A1 - Götschel, Sebastian A1 - Nagaiah, Chamakuri A1 - Kunisch, Karl A1 - Weiser, Martin T1 - Lossy Compression in Optimal Control of Cardiac Defibrillation N2 - This paper presents efficient computational techniques for solving an optimization problem in cardiac defibrillation governed by the monodomain equations. Time-dependent electrical currents injected at different spatial positions act as the control. Inexact Newton-CG methods are used, with reduced gradient computation by adjoint solves. In order to reduce the computational complexity, adaptive mesh refinement for state and adjoint equations is performed. To reduce the high storage and bandwidth demand imposed by adjoint gradient and Hessian-vector evaluations, a lossy compression technique for storing trajectory data is applied. An adaptive choice of quantization tolerance based on error estimates is developed in order to ensure convergence. The efficiency of the proposed approach is demonstrated on numerical examples. T3 - ZIB-Report - 13-26 KW - monodomain model KW - defibrillation KW - optimal control KW - Newton-CG KW - trajectory storage KW - compression Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-18566 SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Lossy Compression for PDE-constrained Optimization: Adaptive Error Control N2 - For the solution of optimal control problems governed by nonlinear parabolic PDEs, methods working on the reduced objective functional are often employed to avoid a full spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the ad-joint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories have to be stored. To get numerical results which are accurate enough, in many case very fine discretizations in time and space are necessary, which leads to a significant amount of data to be stored and transmitted to mass storage. Lossy compression methods were developed to overcome the storage problem by reducing the accuracy of the stored trajectories. The inexact data induces errors in the reduced gradient and reduced Hessian. In this paper, we analyze the influence of such a lossy trajectory compression method on Newton-CG methods for optimal control of parabolic PDEs and design an adaptive strategy for choosing appropriate quantization tolerances. T3 - ZIB-Report - 13-27 KW - optimal control KW - semilinear parabolic PDEs KW - Newton-CG KW - trajectory storage KW - lossy compression Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-18575 SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - Maierhofer, Christiane A1 - Müller, Jan P. A1 - Rothbart, Nick A1 - Weiser, Martin T1 - Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites N2 - Carbon-fiber reinforced composites are becoming more and more important in the production of light-weight structures, e.g., in the automotive and aerospace industry. Thermography is often used for non-destructive testing of these products, especially to detect delaminations between different layers of the composite. In this presentation, we aim at methods for defect reconstruction from thermographic measurements of such carbon-fiber reinforced composites. The reconstruction results shall not only allow to locate defects, but also give a quantitative characterization of the defect properties. We discuss the simulation of the measurement process using finite element methods, as well as the experimental validation on flat bottom holes. Especially in pulse thermography, thin boundary layers with steep temperature gradients occurring at the heated surface need to be resolved. Here we use the combination of a 1D analytical solution combined with numerical solution of the remaining defect equation. We use the simulations to identify material parameters from the measurements. Finally, fast heuristics for reconstructing defect geometries are applied to the acquired data, and compared for their accuracy and utility in detecting different defects like back surface defects or delaminations. T3 - ZIB-Report - 16-13 KW - thermography KW - inverse problems KW - parameter identification KW - geometry identification KW - nondestructive testing Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-58374 UR - https://www.wcndt2016.com/portals/wcndt/bb/Th4C4.pdf SN - 1438-0064 ER - TY - GEN A1 - Götschel, Sebastian A1 - Höhne, Christian A1 - Kolkoori, Sanjeevareddy A1 - Mitzscherling, Steffen A1 - Prager, Jens A1 - Weiser, Martin T1 - Ray Tracing Boundary Value Problems: Simulation and SAFT Reconstruction for Ultrasonic Testing N2 - The application of advanced imaging techniques for the ultrasonic inspection of inhomogeneous anisotropic materials like austenitic and dissimilar welds requires information about acoustic wave propagation through the material, in particular travel times between two points in the material. Forward ray tracing is a popular approach to determine traveling paths and arrival times but is ill suited for inverse problems since a large number of rays have to be computed in order to arrive at prescribed end points. In this contribution we discuss boundary value problems for acoustic rays, where the ray path between two given points is determined by solving the eikonal equation. The implementation of such a two point boundary value ray tracer for sound field simulations through an austenitic weld is described and its efficiency as well as the obtained results are compared to those of a forward ray tracer. The results are validated by comparison with experimental results and commercially available UT simulation tools. As an application, we discuss an implementation of the method for SAFT (Synthetic Aperture Focusing Technique) reconstruction. The ray tracer calculates the required travel time through the anisotropic columnar grain structure of the austenitic weld. There, the formulation of ray tracing as a boundary value problem allows a straightforward derivation of the ray path from a given transducer position to any pixel in the reconstruction area and reduces the computational cost considerably. T3 - ZIB-Report - 16-14 KW - ray tracing KW - ODE boundary value problems KW - ultrasonics KW - nondestructive testing Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-58386 UR - https://www.wcndt2016.com/portals/wcndt/bb/Fr1H4.pdf SN - 1438-0064 ER - TY - CHAP A1 - Götschel, Sebastian A1 - Tycowicz, Christoph von A1 - Polthier, Konrad A1 - Weiser, Martin ED - Carraro, T. ED - Geiger, M. ED - Koerkel, S. ED - Rannacher, R. T1 - Reducing Memory Requirements in Scientific Computing and Optimal Control T2 - Multiple Shooting and Time Domain Decomposition Methods Y1 - 2015 SP - 263 EP - 287 PB - Springer ER - TY - JOUR A1 - Götschel, Sebastian A1 - Nagaiah, Chamakuri A1 - Kunisch, Karl A1 - Weiser, Martin T1 - Lossy Compression in Optimal Control of Cardiac Defibrillation JF - J. Sci. Comput. N2 - This paper presents efficient computational techniques for solving an optimization problem in cardiac defibrillation governed by the monodomain equations. Time-dependent electrical currents injected at different spatial positions act as the control. Inexact Newton-CG methods are used, with reduced gradient computation by adjoint solves. In order to reduce the computational complexity, adaptive mesh refinement for state and adjoint equations is performed. To reduce the high storage and bandwidth demand imposed by adjoint gradient and Hessian-vector evaluations, a lossy compression technique for storing trajectory data is applied. An adaptive choice of quantization tolerance based on error estimates is developed in order to ensure convergence. The efficiency of the proposed approach is demonstrated on numerical examples. Y1 - 2014 U6 - https://doi.org/10.1007/s10915-013-9785-x VL - 60 IS - 1 SP - 35 EP - 59 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Lossy Compression for PDE-constrained Optimization: Adaptive Error Control JF - Comput. Optim. Appl. N2 - For the solution of optimal control problems governed by nonlinear parabolic PDEs, methods working on the reduced objective functional are often employed to avoid a full spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the ad-joint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories have to be stored. To get numerical results which are accurate enough, in many case very fine discretizations in time and space are necessary, which leads to a significant amount of data to be stored and transmitted to mass storage. Lossy compression methods were developed to overcome the storage problem by reducing the accuracy of the stored trajectories. The inexact data induces errors in the reduced gradient and reduced Hessian. In this paper, we analyze the influence of such a lossy trajectory compression method on Newton-CG methods for optimal control of parabolic PDEs and design an adaptive strategy for choosing appropriate quantization tolerances. Y1 - 2015 VL - 62 IS - 1 SP - 131 EP - 155 PB - Springer ER - TY - THES A1 - Götschel, Sebastian T1 - Adaptive Lossy Trajectory Compression for Optimal Control of Parabolic PDEs N2 - Optimal control problems governed by nonlinear, time-dependent PDEs on three-dimensional spatial domains are an important tool in many fields, ranging from engineering applications to medicine. For the solution of such optimization problems, methods working on the reduced objective functional are often employed to avoid a full spatio-temporal discretization of the problem. The evaluation of the reduced gradient requires one solve of the state equation forward in time, and one backward solve of the adjoint equation. The state enters into the adjoint equation, requiring the storage of a full 4D data set. If Newton-CG methods are used, two additional trajectories have to be stored. To get numerical results that are accurate enough, in many cases very fine discretizations in time and space are necessary, leading to a significant amount of data to be stored and transmitted to mass storage. This thesis deals with the development and analysis of methods for lossy compression of such finite element solutions. The algorithms are based on a change of basis to reduce correlations in the data, combined with quantization. This is achieved by transforming the finite element coefficient vector from the nodal to the hierarchical basis, followed by rounding the coefficients to a prescribed precision. Due to the inexact reconstruction, and thus inexact data for the adjoint equation, the error induced in the reduced gradient, and reduced Hessian, has to be controlled, to not impede convergence of the optimization. Accuracy requirements of different optimization methods are analyzed, and computable error estimates for the influence of lossy trajectory storage are derived. These tools are used to adaptively control the accuracy of the compressed data. The efficiency of the algorithms is demonstrated on several numerical examples, ranging from a simple linear, scalar equation to a semi-linear system of reaction-diffusion equations. In all examples considerable reductions in storage space and bandwidth requirements are achieved, without significantly influencing the convergence behavior of the optimization methods. Finally, to go beyond pointwise error control, the hierarchical basis transform can be replaced by more sophisticated wavelet transforms. Numerical experiments indicate that choosing suitable norms for error control allows higher compression factors. KW - optimal control KW - semi-linear parabolic PDEs KW - adjoint gradient computation KW - trajectory storage KW - lossy compression Y1 - 2015 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000098552 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Weiser, Martin A1 - Maierhofer, Christiane A1 - Richter, Regina A1 - Röllig, Mathias T1 - Fast Defect Shape Reconstruction Based on the Travel Time in Pulse Thermography JF - Nondestructive Testing of Materials and Structures Y1 - 2013 VL - 6 SP - 83 EP - 89 ER - TY - JOUR A1 - Weber, Marie-Christin A1 - Fischer, Lisa A1 - Damerau, Alexandra A1 - Ponomarev, Igor A1 - Pfeiffenberger, Moritz A1 - Gaber, Timo A1 - Götschel, Sebastian A1 - Lang, Jens A1 - Röblitz, Susanna A1 - Buttgereit, Frank A1 - Ehrig, Rainald A1 - Lang, Annemarie T1 - In vitro and in silico modeling of cellular and matrix-related changes during the early phase of osteoarthritis JF - BioRxiv N2 - Understanding the pathophysiological processes of osteoarthritis (OA) require adequate model systems. Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific parts of the disease. This study aimed to combine in vitro and in silico modeling to describe cellular and matrix-related changes during the early phase of OA. We developed an in vitro OA model based on scaffold-free cartilage-like constructs (SFCCs), which was mathematically modeled using a partial differential equation (PDE) system to resemble the processes during the onset of OA. SFCCs were produced from mesenchymal stromal cells and analyzed weekly by histology and qPCR to characterize the cellular and matrix-related composition. To simulate the early phase of OA, SFCCs were treated with interleukin-1β (IL-1β), tumor necrosis factor α (TNFα) and examined after 3 weeks or cultivated another 3 weeks without inflammatory cytokines to validate the regeneration potential. Mathematical modeling was performed in parallel to the in vitro experiments. SFCCs expressed cartilage-specific markers, and after stimulation an increased expression of inflammatory markers, matrix degrading enzymes, a loss of collagen II (Col-2) and a reduced cell density was observed which could be partially reversed by retraction of stimulation. Based on the PDEs, the distribution processes within the SFCCs, including those of IL-1β, Col-2 degradation and cell number reduction was simulated. By combining in vitro and in silico methods, we aimed to develop a valid, efficient alternative approach to examine and predict disease progression and new therapeutic strategies. Y1 - 2019 U6 - https://doi.org/10.1101/725317 ER - TY - JOUR A1 - Weber, Marie-Christin A1 - Fischer, Lisa A1 - Damerau, Alexandra A1 - Ponomarev, Igor A1 - Pfeiffenberger, Moritz A1 - Gaber, Timo A1 - Götschel, Sebastian A1 - Lang, Jens A1 - Röblitz, Susanna A1 - Buttgereit, Frank A1 - Ehrig, Rainald A1 - Lang, Annemarie T1 - Macroscale mesenchymal condensation to study cytokine-driven cellular and matrix-related changes during cartilage degradation JF - Biofabrication N2 - Understanding the pathophysiological processes of cartilage degradation requires adequate model systems to develop therapeutic strategies towards osteoarthritis (OA). Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific disease aspects. This study aimed to combine in vitro and in silico modeling based on a tissue-engineering approach using mesenchymal condensation to mimic cytokine-induced cellular and matrix-related changes during cartilage degradation. Thus, scaffold-free cartilage-like constructs (SFCCs) were produced based on self-organization of mesenchymal stromal cells (mesenchymal condensation) and i) characterized regarding their cellular and matrix composition or secondly ii) treated with interleukin-1β (IL-1β) and tumor necrosis factor α (TNFα) for 3 weeks to simulate OA-related matrix degradation. In addition, an existing mathematical model based on partial differential equations was optimized and transferred to the underlying settings to simulate distribution of IL-1β, type II collagen degradation and cell number reduction. By combining in vitro and in silico methods, we aim to develop a valid, efficient alternative approach to examine and predict disease progression and effects of new therapeutics. Y1 - 2020 U6 - https://doi.org/10.1088/1758-5090/aba08f VL - 12 IS - 4 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Minion, Michael L. T1 - An Efficient Parallel-in-Time Method for Optimization with Parabolic PDEs JF - SIAM J. Sci. Comput. 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 the 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 nonlinear reaction-diffusion optimal control problem demonstrate the parallel speedup and efficiency of different approaches. Y1 - 2019 U6 - https://doi.org/10.1137/19M1239313 VL - 41 IS - 6 SP - C603 EP - C626 ER - TY - GEN A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Lossy Compression for Large Scale PDE Problems N2 - Solvers for partial differential equations (PDE) are one of the cornerstones of computational science. For large problems, they involve huge amounts of data that needs to be stored and transmitted on all levels of the memory hierarchy. Often, bandwidth is the limiting factor due to relatively small arithmetic intensity, and increasingly so due to the growing disparity between computing power and bandwidth. Consequently, data compression techniques have been investigated and tailored towards the specific requirements of PDE solvers during the last decades. This paper surveys data compression challenges and corresponding solution approaches for PDE problems, covering all levels of the memory hierarchy from mass storage up to main memory. Exemplarily, we illustrate concepts at particular methods, and give references to alternatives. T3 - ZIB-Report - 19-32 KW - partial differential equation KW - data compression KW - floating point compression KW - lossy compression Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-73817 SN - 1438-0064 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Weiser, Martin T1 - Compression Challenges in Large Scale Partial Differential Equation Solvers JF - Algorithms N2 - Solvers for partial differential equations (PDEs) are one of the cornerstones of computational science. For large problems, they involve huge amounts of data that need to be stored and transmitted on all levels of the memory hierarchy. Often, bandwidth is the limiting factor due to the relatively small arithmetic intensity, and increasingly due to the growing disparity between computing power and bandwidth. Consequently, data compression techniques have been investigated and tailored towards the specific requirements of PDE solvers over the recent decades. This paper surveys data compression challenges and discusses examples of corresponding solution approaches for PDE problems, covering all levels of the memory hierarchy from mass storage up to the main memory. We illustrate concepts for particular methods, with examples, and give references to alternatives. Y1 - 2019 U6 - https://doi.org/10.3390/a12090197 VL - 12 IS - 9 SP - 197 ER - TY - JOUR A1 - Götschel, Sebastian A1 - Schiela, Anton A1 - Weiser, Martin T1 - Kaskade 7 - a Flexible Finite Element Toolbox JF - Computers and Mathematics with Applications N2 - Kaskade 7 is a finite element toolbox for the solution of stationary or transient systems of partial differential equations, aimed at supporting application-oriented research in numerical analysis and scientific computing. The library is written in C++ and is based on the \textsc{Dune} interface. The code is independent of spatial dimension and works with different grid managers. An important feature is the mix-and-match approach to discretizing systems of PDEs with different ansatz and test spaces for all variables. We describe the mathematical concepts behind the library as well as its structure, illustrating its use at several examples on the way. Y1 - 2021 U6 - https://doi.org/10.1016/j.camwa.2020.02.011 VL - 81 SP - 444 EP - 458 ER - TY - GEN A1 - Götschel, Sebastian A1 - Schiela, Anton A1 - Weiser, Martin T1 - Kaskade 7 -- a Flexible Finite Element Toolbox N2 - Kaskade 7 is a finite element toolbox for the solution of stationary or transient systems of partial differential equations, aimed at supporting application-oriented research in numerical analysis and scientific computing. The library is written in C++ and is based on the Dune interface. The code is independent of spatial dimension and works with different grid managers. An important feature is the mix-and-match approach to discretizing systems of PDEs with different ansatz and test spaces for all variables. We describe the mathematical concepts behind the library as well as its structure, illustrating its use at several examples on the way. T3 - ZIB-Report - 19-48 KW - finite elements KW - generic programming KW - partial differential equations Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-74616 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 -