liu.seSearch for publications in DiVA
Change search

Achieng, Pauline
Publications (4 of 4) Show all publications
Achieng, P. (2023). Reconstruction of solutions of Cauchy problems for elliptic equations in bounded and unbounded domains using iterative regularization methods. (Doctoral dissertation). Linköping: Linköping University Electronic Press
Open this publication in new window or tab >>Reconstruction of solutions of Cauchy problems for elliptic equations in bounded and unbounded domains using iterative regularization methods
2023 (English)Doctoral thesis, comprehensive summary (Other academic)
Abstract [en]

Cauchy problems for elliptic equations arise in applications in science and engineering. These problems often involve finding important information about an elliptical system from indirect or incomplete measurements. Cauchy problems for elliptic equations are known to be disadvantaged in the sense that a small pertubation in the input can result in a large error in the output. Regularization methods are usually required in order to be able to find stable solutions. In this thesis we study the Cauchy problem for elliptic equations in both bounded and unbounded domains using iterative regularization methods. In Paper I and II, we focus on an iterative regularization technique which involves solving a sequence of mixed boundary value well-posed problems for the same elliptic equation. The original version of the alternating iterative technique is based on iterations alternating between Dirichlet-Neumann and Neumann-Dirichlet boundary value problems. This iterative method is known to possibly work for Helmholtz equation. Instead we study a modified version based on alternating between Dirichlet-Robin and Robin-Dirichlet boundary value problems. First, we study the Cauchy problem for general elliptic equations of second order with variable coefficients in a limited domain. Then we extend to the case of unbounded domains for the Cauchy problem for Helmholtz equation. For the Cauchy problem, in the case of general elliptic equations, we show that the iterative method, based on Dirichlet-Robin, is convergent provided that parameters in the Robin condition are chosen appropriately. In the case of an unbounded domain, we derive necessary, and sufficient, conditions for convergence of the Robin-Dirichlet iterations based on an analysis of the spectrum of the Laplacian operator, with boundary conditions of Dirichlet and Robin types.

In the numerical tests, we investigate the precise behaviour of the Dirichlet-Robin iterations, for different values of the wave number in the Helmholtz equation, and the results show that the convergence rate depends on the choice of the Robin parameter in the Robin condition. In the case of unbounded domain, the numerical experiments show that an appropriate truncation of the domain and an appropriate choice of Robin parameter in the Robin condition lead to convergence of the Robin-Dirichlet iterations.

In the presence of noise, additional regularization techniques have to implemented for the alternating iterative procedure to converge. Therefore, in Paper III and IV we focus on iterative regularization methods for solving the Cauchy problem for the Helmholtz equation in a semi-infinite strip, assuming that the data contains measurement noise. In addition, we also reconstruct a radiation condition at infinity from the given Cauchy data. For the reconstruction of the radiation condition, we solve a well-posed problem for the Helmholtz equation in a semi-infinite strip. The remaining solution is obtained by solving an ill-posed problem. In Paper III, we consider the ordinary Helmholtz equation and use seperation of variables to analyze the problem. We show that the radiation condition is described by a non-linear well-posed problem that provides a stable oscillatory solution to the Cauchy problem. Furthermore, we show that the ill–posed problem can be regularized using the Landweber’s iterative method and the discrepancy principle. Numerical tests shows that the approach works well.

Paper IV is an extension of the theory from Paper III to the case of variable coefficients. Theoretical analysis of this Cauchy problem shows that, with suitable bounds on the coefficients, can iterative regularization methods be used to stabilize the ill-posed Cauchy problem.

Series
Linköping Studies in Science and Technology. Dissertations, ISSN 0345-7524 ; 2352
National Category
Computational Mathematics
Identifiers
urn:nbn:se:liu:diva-199096 (URN)10.3384/9789180753715 (DOI)9789180753708 (ISBN)9789180753715 (ISBN)
Public defence
2023-12-08, BL32 (Nobel), B Building, Campus Valla, Linköping, 13:15 (English)
Note

Funding agencies: The International Science Programme (ISP) and the Eastern Africa Universities Mathematics Programme (EAUMP).

Available from: 2023-11-13 Created: 2023-11-13 Last updated: 2023-11-13Bibliographically approved
Achieng, P., Berntsson, F. & Kozlov, V. (2023). Robin-Dirichlet alternating iterative procedure for solving the Cauchy problem for Helmholtz equation in an unbounded domain. Journal of Inverse and Ill-Posed Problems, 31(5)
Open this publication in new window or tab >>Robin-Dirichlet alternating iterative procedure for solving the Cauchy problem for Helmholtz equation in an unbounded domain
2023 (English)In: Journal of Inverse and Ill-Posed Problems, ISSN 0928-0219, E-ISSN 1569-3945, Vol. 31, no 5Article in journal (Refereed) Published
Abstract [en]

We consider the Cauchy problem for the Helmholtz equation with a domain in $\mathbb{R}^{d}, d \geq 2$ with N cylindrical outlets to infinity with bounded inclusions in $\mathbb{R}^{d-1}$. Cauchy data are prescribed on the boundary of the bounded domains and the aim is to find solution on the unbounded part of the boundary. In 1989, Kozlov and Mazya proposed an alternating iterative method for solving Cauchy problems associated with elliptic, selfadjoint and positive-definite operators in bounded domains. Different variants of this method for solving Cauchy problems associated with Helmholtz-type operators exists. We consider the variant proposed by Berntsson, Kozlov, Mpinganzima and Turesson (2018) for bounded domains and derive the necessary conditions for the convergence of the procedure in unbounded domains. For the numerical implementation, a finite difference method is used to solve the problem in a simple rectangular domain in R-2 that represent a truncated infinite strip. The numerical results shows that by appropriate truncation of the domain and with appropriate choice of the Robin parameters mu(0) and mu(1), the Robin-Dirichlet alternating iterative procedure is convergent.

Place, publisher, year, edition, pages
WALTER DE GRUYTER GMBH, 2023
Keywords
Helmholtz equation; Cauchy problem; inverse problem ill-posed problem
National Category
Computational Mathematics
Identifiers
urn:nbn:se:liu:diva-192481 (URN)10.1515/jiip-2020-0133 (DOI)000940871600001 ()
Available from: 2023-03-21 Created: 2023-03-21 Last updated: 2024-03-18Bibliographically approved
Achieng, P., Berntsson, F., Chepkorir, J. & Kozlov, V. (2021). Analysis of Dirichlet–Robin Iterations for Solving the Cauchy Problem for Elliptic Equations. Bulletin of the Iranian Mathematical Society, 47, 1681-1699
Open this publication in new window or tab >>Analysis of Dirichlet–Robin Iterations for Solving the Cauchy Problem for Elliptic Equations
2021 (English)In: Bulletin of the Iranian Mathematical Society, ISSN 1735-8515, Vol. 47, p. 1681-1699Article in journal (Refereed) Published
Abstract [en]

The Cauchy problem for general elliptic equations of second order is considered. In a previous paper (Berntsson et al. in Inverse Probl Sci Eng 26(7):1062–1078, 2018), it was suggested that the alternating iterative algorithm suggested by Kozlov and Maz’ya can be convergent, even for large wavenumbers k2, in the Helmholtz equation, if the Neumann boundary conditions are replaced by Robin conditions. In this paper, we provide a proof that shows that the Dirichlet–Robin alternating algorithm is indeed convergent for general elliptic operators provided that the parameters in the Robin conditions are chosen appropriately. We also give numerical experiments intended to investigate the precise behaviour of the algorithm for different values of k2 in the Helmholtz equation. In particular, we show how the speed of the convergence depends on the choice of Robin parameters.

Springer, 2021
Keywords
Helmholtz equation, Cauchy problem, Inverse problem, Ill-posed problem
National Category
Mathematical Analysis
Identifiers
urn:nbn:se:liu:diva-170834 (URN)10.1007/s41980-020-00466-7 (DOI)000575739300001 ()2-s2.0-85092146699 (Scopus ID)
Available from: 2020-10-26 Created: 2020-10-26 Last updated: 2024-02-22Bibliographically approved
Achieng, P. (2020). Analysis of the Robin-Dirichlet iterative procedure for solving the Cauchy problem for elliptic equations with extension to unbounded domains. (Licentiate dissertation). Linköping: Linköping University Electronic Press
Open this publication in new window or tab >>Analysis of the Robin-Dirichlet iterative procedure for solving the Cauchy problem for elliptic equations with extension to unbounded domains
2020 (English)Licentiate thesis, comprehensive summary (Other academic)
Abstract [en]

In this thesis we study the Cauchy problem for elliptic equations. It arises in many areas of application in science and engineering as a problem of reconstruction of solutions to elliptic equations in a domain from boundary measurements taken on a part of the boundary of this domain. The Cauchy problem for elliptic equations is known to be ill-posed.

We use an iterative regularization method based on alternatively solving a sequence of well-posed mixed boundary value problems for the same elliptic equation. This method, based on iterations between Dirichlet-Neumann and Neumann-Dirichlet mixed boundary value problems was first proposed by Kozlov and Maz’ya [13] for Laplace equation and Lame’ system but not Helmholtz-type equations. As a result different modifications of this original regularization method have been proposed in literature. We consider the Robin-Dirichlet iterative method proposed by Mpinganzima et.al [3] for the Cauchy problem for the Helmholtz equation in bounded domains.

We demonstrate that the Robin-Dirichlet iterative procedure is convergent for second order elliptic equations with variable coefficients provided the parameter in the Robin condition is appropriately chosen. We further investigate the convergence of the Robin-Dirichlet iterative procedure for the Cauchy problem for the Helmholtz equation in a an unbounded domain. We derive and analyse the necessary conditions needed for the convergence of the procedure.

In the numerical experiments, the precise behaviour of the procedure for different values of k2 in the Helmholtz equation is investigated and the results show that the speed of convergence depends on the choice of the Robin parameters, μ0 and μ1. In the unbounded domain case, the numerical experiments demonstrate that the procedure is convergent provided that the domain is truncated appropriately and the Robin parameters, μ0 and μ1 are also chosen appropriately.

Series
Linköping Studies in Science and Technology. Licentiate Thesis, ISSN 0280-7971 ; 1891
National Category
Computational Mathematics
Identifiers
urn:nbn:se:liu:diva-170835 (URN)10.3384/lic.diva-170835 (DOI)9789179297565 (ISBN)
Presentation
2020-12-03, Planck, F-building, entrance 57, Campus Valla, Linköping, 10:00 (English)
Note

Funding Agencies: International Science Programme (ISP) and the Eastern Africa Universities Mathematics Programme (EAUMP).

Available from: 2020-10-26 Created: 2020-10-26 Last updated: 2020-10-28Bibliographically approved

Search in DiVA

Show all publications