liu.seSearch for publications in DiVA
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
Reconstruction of the Radiation Condition and Solution for the Helmholtz Equation in a Semi-infinite Strip from Cauchy Data on an Interior Segment
Linköping University, Department of Mathematics, Analysis and Mathematics Education. Linköping University, Faculty of Science & Engineering. Univ Nairobi, Kenya.
Linköping University, Department of Mathematics, Applied Mathematics. Linköping University, Faculty of Science & Engineering.ORCID iD: 0000-0002-2681-8965
Linköping University, Department of Mathematics, Analysis and Mathematics Education. Linköping University, Faculty of Science & Engineering.
2023 (English)In: Computational Methods in Applied Mathematics, ISSN 1609-4840, E-ISSN 1609-9389Article in journal (Refereed) Epub ahead of print
Abstract [en]

We consider an inverse problem for the Helmholtz equation of reconstructing a solution from measurements taken on a segment inside a semi-infinite strip. Homogeneous Neumann conditions are prescribed on both side boundaries of the strip and an unknown Dirichlet condition on the remaining part of the boundary. Additional complexity is that the radiation condition at infinity is unknown. Our aim is to find the unknown function in the Dirichlet boundary condition and the radiation condition. Such problems appear in acoustics to determine acoustical sources and surface vibrations from acoustic field measurements. The problem is split into two sub-problems, a well-posed and an ill-posed problem. We analyse the theoretical properties of both problems; in particular, we show that the radiation condition is described by a stable non-linear problem. The second problem is ill-posed, and we use the Landweber iteration method together with the discrepancy principle to regularize it. Numerical tests show that the approach works well.

Place, publisher, year, edition, pages
WALTER DE GRUYTER GMBH , 2023.
Keywords [en]
Helmholtz Equation; Inverse Problem; Cauchy Problem; Ill-Posed Problem; Well-Posed Problem; Landweber Method
National Category
Fluid Mechanics and Acoustics
Identifiers
URN: urn:nbn:se:liu:diva-196637DOI: 10.1515/cmam-2022-0244ISI: 001035412500001OAI: oai:DiVA.org:liu-196637DiVA, id: diva2:1788778
Available from: 2023-08-17 Created: 2023-08-17 Last updated: 2023-11-13
In thesis
1. Reconstruction of solutions of Cauchy problems for elliptic equations in bounded and unbounded domains using iterative regularization methods
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.

Place, publisher, year, edition, pages
Linköping: Linköping University Electronic Press, 2023. p. 17
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)
Opponent
Supervisors
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

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full text

Search in DiVA

By author/editor
Achieng, PaulineBerntsson, FredrikKozlov, Vladimir
By organisation
Analysis and Mathematics EducationFaculty of Science & EngineeringApplied Mathematics
In the same journal
Computational Methods in Applied Mathematics
Fluid Mechanics and Acoustics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 70 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • oxford
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf