@article{KleefeldReissel2011, author = {Kleefeld, Andreas and Reißel, Martin}, title = {The Levenberg-Marquardt method applied to a parameter estimation problem arising from electrical resistivity tomography}, series = {Applied Mathematics and Computation}, volume = {217}, journal = {Applied Mathematics and Computation}, number = {9}, publisher = {Elsevier}, address = {Amsterdam}, isbn = {0096-3003}, pages = {4490 -- 4501}, year = {2011}, language = {en} } @article{GrajewskiKleefeld2023, author = {Grajewski, Matthias and Kleefeld, Andreas}, title = {Detecting and approximating decision boundaries in low-dimensional spaces}, series = {Numerical Algorithms}, volume = {93}, journal = {Numerical Algorithms}, number = {4}, publisher = {Springer Science+Business Media}, address = {Dordrecht}, issn = {1572-9265}, pages = {35 Seiten}, year = {2023}, abstract = {A method for detecting and approximating fault lines or surfaces, respectively, or decision curves in two and three dimensions with guaranteed accuracy is presented. Reformulated as a classification problem, our method starts from a set of scattered points along with the corresponding classification algorithm to construct a representation of a decision curve by points with prescribed maximal distance to the true decision curve. Hereby, our algorithm ensures that the representing point set covers the decision curve in its entire extent and features local refinement based on the geometric properties of the decision curve. We demonstrate applications of our method to problems related to the detection of faults, to multi-criteria decision aid and, in combination with Kirsch's factorization method, to solving an inverse acoustic scattering problem. In all applications we considered in this work, our method requires significantly less pointwise classifications than previously employed algorithms.}, language = {en} } @article{AyalaHarrisKleefeldetal.2023, author = {Ayala, Rafael Ceja and Harris, Isaac and Kleefeld, Andreas and Pallikarakis, Nikolaos}, title = {Analysis of the transmission eigenvalue problem with two conductivity parameters}, series = {Applicable Analysis}, journal = {Applicable Analysis}, publisher = {Taylor \& Francis}, issn = {0003-6811}, doi = {10.1080/00036811.2023.2181167}, pages = {37 Seiten}, year = {2023}, abstract = {In this paper, we provide an analytical study of the transmission eigenvalue problem with two conductivity parameters. We will assume that the underlying physical model is given by the scattering of a plane wave for an isotropic scatterer. In previous studies, this eigenvalue problem was analyzed with one conductive boundary parameter whereas we will consider the case of two parameters. We prove the existence and discreteness of the transmission eigenvalues as well as study the dependence on the physical parameters. We are able to prove monotonicity of the first transmission eigenvalue with respect to the parameters and consider the limiting procedure as the second boundary parameter vanishes. Lastly, we provide extensive numerical experiments to validate the theoretical work.}, language = {en} } @article{PieronekKleefeld2024, author = {Pieronek, Lukas and Kleefeld, Andreas}, title = {On trajectories of complex-valued interior transmission eigenvalues}, series = {Inverse problems and imaging : IPI}, volume = {18}, journal = {Inverse problems and imaging : IPI}, number = {2}, publisher = {AIMS}, address = {Springfield, Mo}, issn = {1930-8337 (Print)}, doi = {10.3934/ipi.2023041}, pages = {480 -- 516}, year = {2024}, abstract = {This paper investigates the interior transmission problem for homogeneous media via eigenvalue trajectories parameterized by the magnitude of the refractive index. In the case that the scatterer is the unit disk, we prove that there is a one-to-one correspondence between complex-valued interior transmission eigenvalue trajectories and Dirichlet eigenvalues of the Laplacian which turn out to be exactly the trajectorial limit points as the refractive index tends to infinity. For general simply-connected scatterers in two or three dimensions, a corresponding relation is still open, but further theoretical results and numerical studies indicate a similar connection.}, language = {en} } @article{HarrisKleefeld2018, author = {Harris, Isaac and Kleefeld, Andreas}, title = {The inverse scattering problem for a conductive boundary condition and transmission eigenvalues}, series = {Applicable Analysis}, volume = {99}, journal = {Applicable Analysis}, number = {3}, publisher = {Taylor \& Francis}, address = {London}, issn = {1563-504X}, doi = {10.1080/00036811.2018.1504028}, pages = {508 -- 529}, year = {2018}, abstract = {In this paper, we consider the inverse scattering problem associated with an inhomogeneous media with a conductive boundary. In particular, we are interested in two problems that arise from this inverse problem: the inverse conductivity problem and the corresponding interior transmission eigenvalue problem. The inverse conductivity problem is to recover the conductive boundary parameter from the measured scattering data. We prove that the measured scatted data uniquely determine the conductivity parameter as well as describe a direct algorithm to recover the conductivity. The interior transmission eigenvalue problem is an eigenvalue problem associated with the inverse scattering of such materials. We investigate the convergence of the eigenvalues as the conductivity parameter tends to zero as well as prove existence and discreteness for the case of an absorbing media. Lastly, several numerical and analytical results support the theory and we show that the inside-outside duality method can be used to reconstruct the interior conductive eigenvalues.}, language = {en} } @article{KleefeldPieronek2020, author = {Kleefeld, Andreas and Pieronek, J.}, title = {Elastic transmission eigenvalues and their computation via the method of fundamental solutions}, series = {Applicable Analysis}, volume = {100}, journal = {Applicable Analysis}, number = {16}, publisher = {Taylore \& Francis}, address = {London}, issn = {1563-504X}, doi = {10.1080/00036811.2020.1721473}, pages = {3445 -- 3462}, year = {2020}, abstract = {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.}, language = {en} } @article{BreussKleefeld2020, author = {Breuß, Michael and Kleefeld, Andreas}, title = {Implicit monotone difference methods for scalar conservation laws with source terms}, series = {Acta Mathematica Vietnamica}, volume = {45}, journal = {Acta Mathematica Vietnamica}, publisher = {Springer Singapore}, address = {Singapore}, issn = {2315-4144}, doi = {10.1007/s40306-019-00354-1}, pages = {709 -- 738}, year = {2020}, abstract = {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.}, language = {en} } @article{AsanteAsamaniKleefeldWade2020, author = {Asante-Asamani, E.O. and Kleefeld, Andreas and Wade, B.A.}, title = {A second-order exponential time differencing scheme for non-linear reaction-diffusion systems with dimensional splitting}, series = {Journal of Computational Physics}, volume = {415}, journal = {Journal of Computational Physics}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0021-9991}, doi = {10.1016/j.jcp.2020.109490}, year = {2020}, abstract = {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.}, language = {en} } @incollection{Kleefeld2020, author = {Kleefeld, Andreas}, title = {Numerical calculation of interior transmission eigenvalues with mixed boundary conditions}, series = {Computational and Analytic Methods in Science and Engineering}, booktitle = {Computational and Analytic Methods in Science and Engineering}, editor = {Constanda, Christian}, publisher = {Birkh{\"a}user}, address = {Cham}, isbn = {978-3-030-48185-8 (Hardcover)}, doi = {10.1007/978-3-030-48186-5_9}, pages = {173 -- 195}, year = {2020}, abstract = {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.}, language = {en} } @article{MartinVaqueroKleefeld2020, author = {Mart{\´i}n-Vaquero, J. and Kleefeld, Andreas}, title = {Solving nonlinear parabolic PDEs in several dimensions: Parallelized ESERK codes}, series = {Journal of Computational Physics}, journal = {Journal of Computational Physics}, number = {423}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0021-9991}, doi = {10.1016/j.jcp.2020.109771}, year = {2020}, abstract = {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.}, language = {en} }