SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research Article

A Spectral Regularization Method for a Cauchy Problem of the Modified Helmholtz Equation

Ailin Qian*, Jianfeng Mao and Lianghua Liu

Author Affiliations

School of Mathematics and Statistics, Xianning College, Xianning, Hubei 437100, China

For all author emails, please log on.

Boundary Value Problems 2010, 2010:212056  doi:10.1155/2010/212056


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2010/1/212056


Received:15 December 2009
Accepted:9 May 2010
Published:13 June 2010

© 2010 The Author(s)

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We consider the Cauchy problem for a modified Helmholtz equation, where the Cauchy data is given at and the solution is sought in the interval . A spectral method together with choice of regularization parameter is presented and error estimate is obtained. Combining the method of lines, we formulate regularized solutions which are stably convergent to the exact solutions.

1. Introduction

The Cauchy problem for Helmholtz equation arises from inverse scattering problems. Specific backgrounds can be seen in the existing literature; we can refer to [16] and so forth. A number of numerical methods for stabilizing this problem are developed. Several boundary element methods combined with iterative, conjugate gradient, Tikhonov regularization, and singular value decomposition methods are compared in [6]. Cauchy problem for elliptic equations is well known to be severely ill-posed [7]; that is, the solution does not depend continuously on the boundary data, and small errors in the boundary data can amplify the numerical solution infinitely; hence it is impossible to solve Cauchy problem of Helmholtz equation by using classical numerical methods and it requires special techniques, for example, regularization methods. Although theoretical concepts and computational implementation related to the Cauchy problem of Helmholtz equation have been discussed by many authors [811], there are many open problems deserved to be solved. For example, many authors have considered the following Cauchy problem of the Helmholtz equation [811]:

(11)

However, the boundary condition is very strict. If the boundary condition is replaced by , where is a function of , then their methods cannot be applied easily. Therefore, in this paper, we consider the following Cauchy problem for the Helmholtz equation:

(12)

We want to seek the solution in the interval from the Cauchy data pairs located at . Of course, since and are assumed to be measured, there must be measurement errors, and we would actually have noisy data function , for which the measurement errors and are small. Here and in the following sections, denotes the norm. Thus (1.2) is a noncharacteristic Cauchy problem with appropriate Cauchy data given on the line .

However, the ill-posedness is caused by high frequency. By introducing a "cutoff'' frequency we can obtain a well-posed problem. This method has been used for solving inverse heat conduction problem [12], and sideways heat equation [13], sideways parabolic equation [14]. An error estimate for the proposed method can be found in Section 2. The implementation of the numerical method is explained.

2. Regularization and Error Estimate

First, let be the Fourier transform of :

(21)

Taking Fourier transformation for (1.2), we have a family of problems parameterized by :

(22)

The solution can easily be verified to be

(23)

Following the idea of Fourier method, we consider (2.2) only for by cutting off high frequency and define a regularized solution:

(24)

where is the characteristic function of the interval . The solution can be found by using the inverse Fourier transform. Define the regularized solution with measured data by . The difference between the exact solution and the regularized solution can be divided into two parts:

(25)

where

(26)

We rewrite (2.2) as a system of ordinary differential equations:

(27)

Letting , we can rewrite (2.7) as

(28)

where the matrix is the one in (2.7). The reason for using instead of in the definition of is that with this choice, the matrix is normal and hence diagonalized by a unitary matrix. The eigenvalues of are . Thus we can factorize as

(29)

where is a unitary matrix, it follows that the solution of the system (2.8) can be written as

(210)

Since is unitary, and therefore

(211)

where denotes both the Euclidean norm in the complex vector space and the subordinate matrix norm, remembering in (2.11), we obtain

(212)

This inequality is valid for all and we can integrate over and use the Parseval theorem to obtain estimates for and in the -norm. First we will prove a bound on the difference between any two regularized solutions (2.4). We have errors in the measured and . These two cases are treated separately.

Lemma 2.1.

Assume that one has two regularized solutions and defined by (2.4), with the Cauchy data and . Then

(213)

Proof.

The function satisfies the differential equation; thus solves (2.2), with initial data given by

(214)

and thus satisfies inequality(2.12). We have

(215)

If , then

(216)

holds. By inserting (2.15) into (2.16) we get

(217)

If we insert and integrate over , then we obtain

(218)

Thus (2.13) holds.

Lemma 2.2.

Assume that one has two regularized solutions and defined by (2.4), with the Cauchy data and . Then

(219)

Proof.

The function satisfies the differential equation, with initial data given by

(220)

and thus inequality (2.12) holds. It follows that

(221)

and thus, for ,

(222)

By integrating over the interval , we get

(223)

Since is equal to zero outside the interval , we can extend the integrals. Inserting we get

(224)

This is precisely (2.19).

Next we prove that, for the regularized problem, we have stable dependence of the data. By using the two previous lemmas, we get the following.

Lemma 2.3 (stability).

Assume that is the regularized solution (2.4), with exact data and is the regularized solution with noisy data ; then

(225)

Proof.

Let be a regularized solution defined by (2.4), with data . Then by Lemma 2.1,

(226)

Using Lemma 2.2, we have

(227)

By the triangle inequality,

(228)

This completes the proof.

By Lemma 2.3 the regularized solution depends continuously on the data. Next we derive a bound on the truncation error when we neglect frequencies in (2.3). So far we have not used any information about for we, assume that the Helmholtz equation is valid in a large interval . By imposing a priori bounds on the solution at and at , we obtain an estimate of the difference between the exact solution and a regularized solution with "cutoff'' level . This is a convergence result in the sense that as for the case of exact data . The following estimate holds.

Lemma 2.4 (convergence).

Suppose that is the solution of the problem (2.2), and that the Helmholtz equation is valid for . Then the difference between and a regularized solution can be estimated:

(229)

where the constant is defined by

(230)

Proof.

The solution of (2.2) can be written in the floowing form:

(231)

where and can be determined from the boundary conditions:

(232)

Solving for and we find that the solution can be written:

(233)

where

(234)

We make the observation that

(235)

and that

(236)

Using the expression (2.33) and the triangle inequality, we get

(237)

The first term on the right-hand side satisfies

(238)

Similarly, the second term can be estimated:

(239)

By combining these two expressions, we have

(240)

Thus the proof is complete.

Remark 2.5.

The constant is well defined, but its value has to be estimated. From a numerical computation we conclude that .

Remark 2.6.

When solving (2.2) numerically we need Cauchy data , along the line . The most natural way to obtain this is to use two thermocouples, located at , and , and compute by solving a well-posed problem for the Helmholtz equation in the interval . Hence it is natural to assume knowledge about at a second point.

Let us summarize what we have so far. The constant in Lemma 2.4 is unchanged. The propagated data error is estimated using Lemma 2.3:

(241)

and the truncation error is estimated using Lemma 2.4:

(242)

These two results can be combined into an error estimate for the spectral method. This is demonstrated in two examples.

Example 2.7.

Suppose that , and that we have an estimate of the noise level, . If and if we choose , with this choice, then by expression (2.42),

(243)

where we have assume that and used the bound . By expression (2.41),

(244)

Thus we obtain an error estimate:

(245)

Note that, under these assumptions, can be used as a rule for selecting the regularization parameter.

Example 2.8.

Suppose that the Helmholtz equation is valid in the interval and that we have a priori bounds and . Furthermore, we assume that the measured data satisfies and . Then we have the estimates:

(246)

By balancing these two components and , we can find a suitable value for the regularization parameter.

3. Numerical Implementation

In this section, we use the method of Lines. Let

(31)

and then we can obtain

(32)

where

(33)

The matrix approximates the second-order derivative by spectral cutoff method (see [2]). This ordinary differential equation system can easily be solved by various numerical methods.

Example 3.1.

It is easy to verify that the function is the exact solution of (1.2). Now we need to seek the solution , where from the Cauchy data pairs and .

In our numerical experiment, we give the comparison between the exact solution of problem (2.2) and its approximation by spectral cutoff method for different noise levels and different locations . Please see Figures 1, 2, 3, and 4. In the experiment, the regularization parameter is chosen according to the Remark of Remark 2.5, where we take . From the above results, it is easy to see that the numerical effect of the spectral method works well. Moreover, we can also see that the lower the noise level is, the better the approximate effect is; the closer to the boundary the location is, the better the approximate effect is. These accord with the theory in Section 2.

Acknowledgments

The author wants to express his thanks to the referee for many valuable comments. This work is supported by the Educational Commission of Hubei Province of China (Q20102804, T201009).

References

  1. Kirsch, A: An Introduction to the Mathematical Theory of Inverse Problems, Applied Mathematical Sciences,p. x+282. Springer, New York, NY, USA (1996)

  2. Harari, I, Barbone, PE, Slavutin, M, Shalom, R: Boundary infinite elements for the Helmholtz equation in exterior domains. International Journal for Numerical Methods in Engineering. 41, 1105–1131 (1998). Publisher Full Text OpenURL

  3. Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: An alternating iterative algorithm for the Cauchy problem associated to the Helmholtz equation. Computer Methods in Applied Mechanics and Engineering. 192(5-6), 709–722 (2003). Publisher Full Text OpenURL

  4. Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: Conjugate gradient-boundary element solution to the Cauchy problem for Helmholtz-type equations. Computational Mechanics. 31(3-4), 367–377 (2003)

  5. Regińska, T, Regiński, K: Approximate solution of a Cauchy problem for the Helmholtz equation. Inverse Problems. 22(3), 975–989 (2006). Publisher Full Text OpenURL

  6. Marin, L, Elliott, L, Heggs, PJ, Ingham, DB, Lesnic, D, Wen, X: Comparison of regularization methods for solving the Cauchy problem associated with the Helmholtz equation. International Journal for Numerical Methods in Engineering. 60(11), 1933–1947 (2004). Publisher Full Text OpenURL

  7. Hadamard, J: Lectures on the Cauchy Problem in Linear Partial Differential Equations, Yale University Press, New Haven, Conn, USA (1923)

  8. Xiong, XT, Fu, CL: Two approximate methods of a Cauchy problem for the Helmholtz equation. Computational & Applied Mathematics. 26(2), 285–307 (2007). PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Fu, CL, Feng, XL, Qian, Z: The Fourier regularization for solving the Cauchy problem for the Helmholtz equation. Applied Numerical Mathematics. 59(10), 2625–2640 (2009). Publisher Full Text OpenURL

  10. Qin, HH, Wei, T: Modified regularization method for the Cauchy problem of the Helmholtz equation. Applied Mathematical Modelling. 33(5), 2334–2348 (2009). Publisher Full Text OpenURL

  11. Qian, AL, Xiong, XT, Wu, YJ: On a quasi-reversibility regularization method for a Cauchy problem of the Helmholtz equation. Journal of Computational and Applied Mathematics. 233(8), 1969–1979 (2010). Publisher Full Text OpenURL

  12. Eldén, L, Berntsson, F: Spectral and wavelet methods for solving an inverse heat conduction problem. Inverse Problems in Engineering Mechanics, pp. 3–10. Elsevier, Oxford, UK (1998)

  13. Berntsson, F: A spectral method for solving the sideways heat equation. Inverse Problems. 15(4), 891–906 (1999). Publisher Full Text OpenURL

  14. Xiong, XT, Fu, CL, Cheng, J: Spectral regularization methods for solving a sideways parabolic equation within the framework of regularization theory. Mathematics and Computers in Simulation. 79(5), 1668–1678 (2009). Publisher Full Text OpenURL