TY - JOUR A1 - Kleefeld, Andreas A1 - Pieronek, J. T1 - Elastic transmission eigenvalues and their computation via the method of fundamental solutions JF - Applicable Analysis N2 - A stabilized version of the fundamental solution method to catch ill-conditioning effects is investigated with focus on the computation of complex-valued elastic interior transmission eigenvalues in two dimensions for homogeneous and isotropic media. Its algorithm can be implemented very shortly and adopts to many similar partial differential equation-based eigenproblems as long as the underlying fundamental solution function can be easily generated. We develop a corroborative approximation analysis which also implicates new basic results for transmission eigenfunctions and present some numerical examples which together prove successful feasibility of our eigenvalue recovery approach. KW - elastic scattering KW - method of fundamental solutions KW - Interior transmission eigenvalues Y1 - 2020 U6 - https://doi.org/10.1080/00036811.2020.1721473 SN - 1563-504X VL - 100 IS - 16 SP - 3445 EP - 3462 PB - Taylore & Francis CY - London ER - TY - CHAP A1 - Abele, Daniel A1 - Kleefeld, Andreas ED - Constanda, Christian T1 - New Numerical Results for the Optimization of Neumann Eigenvalues T2 - Computational and Analytic Methods in Science and Engineering N2 - We present new numerical results for shape optimization problems of interior Neumann eigenvalues. This field is not well understood from a theoretical standpoint. The existence of shape maximizers is not proven beyond the first two eigenvalues, so we study the problem numerically. We describe a method to compute the eigenvalues for a given shape that combines the boundary element method with an algorithm for nonlinear eigenvalues. As numerical optimization requires many such evaluations, we put a focus on the efficiency of the method and the implemented routine. The method is well suited for parallelization. Using the resulting fast routines and a specialized parametrization of the shapes, we found improved maxima for several eigenvalues. Y1 - 2020 SN - 978-3-030-48185-8 (Print) SN - 978-3-030-48186-5 (Online) U6 - https://doi.org/10.1007/978-3-030-48186-5_1 SP - 1 EP - 20 PB - Birkhäuser CY - Cham ER - TY - GEN A1 - Burgeth, Bernhard A1 - Kleefeld, Andreas A1 - Naegel, Benoît A1 - Perret, Benjamin T1 - Editorial — Special Issue: ISMM 2019 T2 - Mathematical Morphology - Theory and Applications N2 - This editorial presents the Special Issue dedicated to the conference ISMM 2019 and summarizes the articles published in this Special Issue. Y1 - 2020 U6 - https://doi.org/10.1515/mathm-2020-0200 SN - 2353-3390 VL - 4 IS - 1 SP - 159 EP - 161 PB - De Gruyter CY - Warschau ER - TY - JOUR A1 - Martín-Vaquero, J. A1 - Kleefeld, Andreas T1 - Solving nonlinear parabolic PDEs in several dimensions: Parallelized ESERK codes JF - Journal of Computational Physics N2 - There is a very large number of very important situations which can be modeled with nonlinear parabolic partial differential equations (PDEs) in several dimensions. In general, these PDEs can be solved by discretizing in the spatial variables and transforming them into huge systems of ordinary differential equations (ODEs), which are very stiff. Therefore, standard explicit methods require a large number of iterations to solve stiff problems. But implicit schemes are computationally very expensive when solving huge systems of nonlinear ODEs. Several families of Extrapolated Stabilized Explicit Runge-Kutta schemes (ESERK) with different order of accuracy (3 to 6) are derived and analyzed in this work. They are explicit methods, with stability regions extended, along the negative real semi-axis, quadratically with respect to the number of stages s, hence they can be considered to solve stiff problems much faster than traditional explicit schemes. Additionally, they allow the adaptation of the step length easily with a very small cost. Two new families of ESERK schemes (ESERK3 and ESERK6) are derived, and analyzed, in this work. Each family has more than 50 new schemes, with up to 84.000 stages in the case of ESERK6. For the first time, we also parallelized all these new variable step length and variable number of stages algorithms (ESERK3, ESERK4, ESERK5, and ESERK6). These parallelized strategies allow to decrease times significantly, as it is discussed and also shown numerically in two problems. Thus, the new codes provide very good results compared to other well-known ODE solvers. Finally, a new strategy is proposed to increase the efficiency of these schemes, and it is discussed the idea of combining ESERK families in one code, because typically, stiff problems have different zones and according to them and the requested tolerance the optimum order of convergence is different. KW - Multi-dimensional partial differential equations KW - Higher-order codes KW - Nonlinear PDEs Y1 - 2020 U6 - https://doi.org/10.1016/j.jcp.2020.109771 SN - 0021-9991 IS - 423 PB - Elsevier CY - Amsterdam ER - TY - CHAP A1 - Kleefeld, Andreas ED - Constanda, Christian T1 - Numerical calculation of interior transmission eigenvalues with mixed boundary conditions T2 - Computational and Analytic Methods in Science and Engineering N2 - Interior transmission eigenvalue problems for the Helmholtz equation play an important role in inverse wave scattering. Some distribution properties of those eigenvalues in the complex plane are reviewed. Further, a new scattering model for the interior transmission eigenvalue problem with mixed boundary conditions is described and an efficient algorithm for computing the interior transmission eigenvalues is proposed. Finally, extensive numerical results for a variety of two-dimensional scatterers are presented to show the validity of the proposed scheme. Y1 - 2020 SN - 978-3-030-48185-8 (Hardcover) U6 - https://doi.org/10.1007/978-3-030-48186-5_9 SP - 173 EP - 195 PB - Birkhäuser CY - Cham ER - TY - JOUR A1 - Asante-Asamani, E.O. A1 - Kleefeld, Andreas A1 - Wade, B.A. T1 - A second-order exponential time differencing scheme for non-linear reaction-diffusion systems with dimensional splitting JF - Journal of Computational Physics N2 - A second-order L-stable exponential time-differencing (ETD) method is developed by combining an ETD scheme with approximating the matrix exponentials by rational functions having real distinct poles (RDP), together with a dimensional splitting integrating factor technique. A variety of non-linear reaction-diffusion equations in two and three dimensions with either Dirichlet, Neumann, or periodic boundary conditions are solved with this scheme and shown to outperform a variety of other second-order implicit-explicit schemes. An additional performance boost is gained through further use of basic parallelization techniques. KW - Exponential time differencing KW - Real distinct pole KW - Dimensional splitting KW - Reaction-diffusion systems KW - Matrix exponential Y1 - 2020 U6 - https://doi.org/10.1016/j.jcp.2020.109490 SN - 0021-9991 N1 - Corresponding author: Andreas Kleefeld VL - 415 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Breuß, Michael A1 - Kleefeld, Andreas T1 - Implicit monotone difference methods for scalar conservation laws with source terms JF - Acta Mathematica Vietnamica N2 - In this article, a concept of implicit methods for scalar conservation laws in one or more spatial dimensions allowing also for source terms of various types is presented. This material is a significant extension of previous work of the first author (Breuß SIAM J. Numer. Anal. 43(3), 970–986 2005). Implicit notions are developed that are centered around a monotonicity criterion. We demonstrate a connection between a numerical scheme and a discrete entropy inequality, which is based on a classical approach by Crandall and Majda. Additionally, three implicit methods are investigated using the developed notions. Next, we conduct a convergence proof which is not based on a classical compactness argument. Finally, the theoretical results are confirmed by various numerical tests. KW - Entropy solution KW - Source term KW - Monotone methods KW - Implicit methods KW - Finite difference methods KW - Conservation laws Y1 - 2020 U6 - https://doi.org/10.1007/s40306-019-00354-1 SN - 2315-4144 N1 - Corresponding author: Andreas Kleefeld VL - 45 SP - 709 EP - 738 PB - Springer Singapore CY - Singapore ER - TY - JOUR A1 - Ketelhut, Maike A1 - Brügge, G. M. A1 - Göll, Fabian A1 - Braunstein, Bjoern A1 - Albracht, Kirsten A1 - Abel, Dirk T1 - Adaptive iterative learning control of an industrial robot during neuromuscular training JF - IFAC PapersOnLine N2 - To prevent the reduction of muscle mass and loss of strength coming along with the human aging process, regular training with e.g. a leg press is suitable. However, the risk of training-induced injuries requires the continuous monitoring and controlling of the forces applied to the musculoskeletal system as well as the velocity along the motion trajectory and the range of motion. In this paper, an adaptive norm-optimal iterative learning control algorithm to minimize the knee joint loadings during the leg extension training with an industrial robot is proposed. The response of the algorithm is tested in simulation for patients with varus, normal and valgus alignment of the knee and compared to the results of a higher-order iterative learning control algorithm, a robust iterative learning control and a recently proposed conventional norm-optimal iterative learning control algorithm. Although significant improvements in performance are made compared to the conventional norm-optimal iterative learning control algorithm with a small learning factor, for the developed approach as well as the robust iterative learning control algorithm small steady state errors occur. KW - Iterative learning control KW - Robotic rehabilitation KW - Adaptive control Y1 - 2020 U6 - https://doi.org/10.1016/j.ifacol.2020.12.741 SN - 2405-8963 VL - 53 IS - 2 SP - 16468 EP - 16475 PB - Elsevier CY - Amsterdam ER - TY - THES A1 - Achtsnicht, Stefan T1 - Multiplex-Magnetdetektion von superparamagnetischen Beads zur Identifizierung von Trinkwasserkontaminationen N2 - Die qualitative und quantitative Detektion von Zielsubstanzen innerhalb einer wässrigen Probe ist für viele Fragestellungen von Interesse, etwa bei der Detektion von Kontaminationen in Trinkwasser in Krisensituationen. Hierbei ist es nicht nur wichtig, dass Pathogene möglichst sensitiv detektiert werden können, sondern auch, dass die Analyse schnell erfolgt, um Betroffenen im Katastrophenfall zügig sicheres Trinkwasser zu Verfügung stellen zu können. Da bei einem solchen Szenario nicht von einer in der Nähe befindlichen funktionierenden Laborinfrastruktur ausgegangen werden kann, ist es wichtig, dass die Messung direkt vor Ort erfolgen kann. Im Rahmen dieser Arbeit wurde untersucht, ob eine derartige Schnellanalytik mithilfe von superparamagnetischen Beads (MBs) und der magnetischen Frequenzmischtechnik möglich ist. Dabei werden die MBs mit Hilfe von primären Antikörpern an die Zielsubstanz gebunden und mit sekundären Antikörpern an die Poren-Oberfläche eines Polyethylen-Filters fixiert (Sandwich-Immunoassay). So kann die Quantifizierung der Zielsubstanz auf eine magnetische Messung der immobilisierten MB-Marker zurückgeführt werden. Die magnetische Frequenzmischtechnik basiert auf der Anregung der Probe mit Magnetfeldern zweier verschiedener Frequenzen. Die durch die nichtlineare Magnetisierungsform der superparamagnetischen MBs entstehenden Mischfrequenzen werden typischerweise mithilfe einer zweistufigen Lock-in-Detektion analysiert (analoge Demodulation), die in einem Magnetreader als Handheldgerät realisiert wurde. Zusätzlich zu dieser Technik wurde das Prinzip der direkten Digitalisierung des gesamten Antwortsignals mit anschließender Fourier-Analyse der erzeugten Mischfrequenzen experimentell umgesetzt, um die Amplituden und Phasen mehrerer Mischfrequenzen simultan zu erfassen. Eine Möglichkeit zur Sensitivitätssteigerung ist die magnetische Aufkonzentration, indem vor der magnetischen Analyse eine Separation der MBs aus einem größeren Probenvolumen mittels magnetischem Feldgradienten durchgeführt wird. Zur Charakterisierung verschiedener kommerzieller MBs hinsichtlich ihrer magnetischen Separierbarkeit wurde ein Aufbau zur Messung ihrer magnetophoretischen Beweglichkeiten realisiert und ihre Geschwindigkeiten im Gradientenfeld mikroskopisch gemessen.Da eine Probe oftmals nicht nur auf eine einzige Zielsubstanz, sondern simultan auf mehrere verschiedene Pathogene hin untersucht werden soll, wurden verschiedene Ansätze entwickelt und getestet, die einen solchen multiparametrischen magnetischen Immunoassay ermöglichen. Einerseits wurde eine räumliche Separation der Bindungsbereiche für verschiedene Zielsubstanzen realisiert, die sequentiell ausgewertet werden können. Andererseits wurde die Unterscheidung von verschiedenen Zielsubstanzen anhand der Charakteristika der an sie gebundenen, verschieden funktionalisierten MB-Typen untersucht. Für eine solche Unterscheidung wurde zum einen die Anregefrequenz der magnetischen Frequenzmischtechnik während einer Messung variiert. Damit konnte gezeigt werden, dass sich verschiedene MB-Sorten anhand der Phase ihrer Frequenzmischsignale voneinander unterscheiden lassen. Weiterhin wurde gezeigt, dass sich der Signalverlauf einer binären Mischung zweier verschiedener MB-Typen als gradueller Übergang der Verläufe der beiden reinen MB-Lösungen ergibt. Eine weitere Analysemethode für einen multiparametrischen Immunoassay besteht darin, ein zusätzliches einstellbares statisches magnetisches Offsetfeld zu verwenden. Hierfür wurden mehrere Aufbauten auf Basis von Permanent- und Elektromagneten simuliert, konstruiert und charakterisiert. Mithilfe von Simulationen konnte gezeigt werden, dass eine auf diesem Verfahren beruhende Unterscheidung für MBs mit unterschiedlichen magnetischen Partikelmomenten möglich ist. Als direkte Anwendung des hier entwickelten Magnetreaders in Zusammenspiel mit der digitalen Demodulation wurde ein magnetischer Assay gegen die B-Untereinheit des Choleratoxins in Trinkwasser mit einem niedrigen Detektionslimit von 0,2 ng/ml demonstriert. N2 - The qualitative and quantitative detection of target substances in an aqueous sample is of interest for many questions, for example in the detection of contaminations in drinking water in crisis situations. It is not only important that pathogens can be detected with highest possible sensitivity, but also that the analysis is carried out quickly so that safe drinking water can be provided in the event of a disaster. During such a scenario one cannot rely on a functioning laboratory infrastructure nearby. Therefore it is important that the measurement can be carried out directly on site. Within the scope of this work, it was investigated whether such a quick analysis can be performed using superparamagnetic beads (MBs) and the magnetic frequency mixing technique. The MBs are bound to the target substance with the aid of primary antibodies and fixed to the pore surfaces of a polyethylene filter with secondary antibodies (sandwich immunoassay). The quantification of the target substance can thus be traced back to a magnetic measurement of the immobilized MB markers. The magnetic frequency mixing technique is based on the excitation of the sample with magnetic fields of two different frequencies. The mixing frequencies generated due to the non-linear shape of the magnetization of the superparamagnetic MBs are typically analyzed using a two-stage Lock-in detection (analog demodulation), which was implemented in a magnetic reader as a handheld device. In addition to this technique, the principle of direct digitization of the entire response signal with subsequent Fourier analysis of the generated mixing frequencies was experimentally implemented in order to simultaneously record the amplitudes and phases of several mixing frequencies. One possibility for increasing the sensitivity is magnetic concentration. In that case, the MBs are separated from a larger sample volume by means of a magnetic field gradient before the magnetic analysis. To characterize various commercial MBs with regard to their magnetic separability, a setup for measuring their magnetophoretic mobility was implemented and their velocities in the gradient field were measured with an optical microscope.Often, a sample has to be examined not only for a single target substance, but for several different pathogens simultaneously. Various approaches have been developed and tested which enable such a multiparametric magnetic immunoassay. On the one hand, a spatial separation of the binding areas for different target substances was realized, which can be evaluated sequentially. On the other hand, a distinction among different target substances based on the magnetic characteristics of their attached different MB types was examined. For this discrimination, the excitation frequency of the magnetic frequency mixing technology was varied during measurement. It is shown that different MB types can be distinguished from one another based on the phase of their frequency mixing signals. The signal curve of a binary mixture of two different MB types is obtained as a gradual transition of the curves of the two pure MB solutions. Another method of analysis for a multiparametric immunoassay is based on an additional adjustable static magnetic offset field. For this purpose, several setups based on permanent magnets and electromagnets were simulated, designed and characterized. The simulations show that a distinction based on this method is possible for MBs with different magnetic particle moments. As a direct application of the developed magnetic reader in conjunction with digital demodulation, a magnetic assay against the B subunit of cholera toxin in drinking water was demonstrated, and a low detection limit of 0.2 ng/ml was achieved. KW - Choleratoxin B KW - Trinkwassersicherheit KW - cholera toxin B KW - drinking water safety KW - magnetic frequency mixing technique Y1 - 2020 U6 - https://doi.org/10.18154/RWTH-2020-12052 N1 - Dissertation, RWTH Aachen University, 2020 ER - TY - JOUR A1 - Gaigall, Daniel T1 - Rothman–Woodroofe symmetry test statistic revisited JF - Computational Statistics & Data Analysis N2 - The Rothman–Woodroofe symmetry test statistic is revisited on the basis of independent but not necessarily identically distributed random variables. The distribution-freeness if the underlying distributions are all symmetric and continuous is obtained. The results are applied for testing symmetry in a meta-analysis random effects model. The consistency of the procedure is discussed in this situation as well. A comparison with an alternative proposal from the literature is conducted via simulations. Real data are analyzed to demonstrate how the new approach works in practice. Y1 - 2020 U6 - https://doi.org/10.1016/j.csda.2019.106837 SN - 0167-9473 VL - 2020 IS - 142 SP - Artikel 106837 PB - Elsevier CY - Amsterdam ER -