Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 878109, 13 pages

doi:10.1155/2012/878109

Research Article

A Tikhonov-Type Regularization Method for Identifying the Unknown Source in the Modified Helmholtz Equation

Jinghuai Gao,¹ Dehua Wang,¹ and Jigen Peng²

¹Institute of Wave and Information, National Engineering Laboratory for Offshore Oil Exploration, School of Electronic and Information Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China

²School of Mathematics and Statistics, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China

Received 30 September 2012; Revised 16 November 2012; Accepted 22 November 2012

Academic Editor: Carla Roque

Copyright © 2012 Jinghuai Gao et al. 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.

Abstract.

An inverse source problem in the modified Helmholtz equation is considered. We give a Tikhonov-type regularization method and set up a theoretical frame to analyze the convergence of such method. A priori and a posteriori choice rules to find the regularization parameter are given. Numerical tests are presented to illustrate the effectiveness and stability of our proposed method.

1. Introduction

The Helmholtz equation often arises in the study of physical problems involving partial differential equations (PDEs) in both space and time. It has a wide range of applications, for example, radar, sonar, geographical exploration, and medical imaging. A kind of important equations similar with the Helmholtz equation in science and engineering is

where the constant is the wave number and is the source term. This equation is called the modified Helmholtz equation. It appears, for example, in the semi-implicit temporal discretization of the heat or the Navier-Stokes equations [1] and in the linearized Poisson-Boltzmann equation.

Inverse source problems have attracted great attention of many researchers over recent years because of their applications to many practical problems such as crack determination [2, 3], heat source determination [4–6], inverse heat conduction [7–11], pollution source identification [12], electromagnetic source identification [13], Stefan design problems [14], sound source reconstruction [15], and identification of current dipolar sources in the so-called inverse electroencephalography/magnetoencephalography (EEG/MEG) problems [16, 17]. Theoretical investigation on the inverse source identification problems can be found in the works of [18–23].

The main difficulty of inverse source identification problems is that they are typically ill posed in the sense of Hadamard [24]. In other words, any small error in the scattered measurement data may induce enormous error to the solution. In general, the unknown source can only be recovered from boundary measurements if some a priori knowledge is assumed. For instances, if one of the products in the separation of variables is known [25, 26], the base area of a cylindrical source is known [25], or a nonseparable type is in the form of a moving front [26], then the boundary condition plus some scattered boundary measurements can uniquely determine the unknown source term. Furthermore, when the unknown source term is relatively smooth, some regularization techniques can be employed, see [5, 27–29] for more details. In addition, due to the complexity and ill posedness of the inverse source identification problems, some of the variational methods [6, 30] are also employed to deal with them.

In this paper, we will consider the following problem (see [29]):

for determining the source term such that the solution of the modified Helmholtz equation satisfies the given supplementary condition , where the constant is the wave number. In practice, the data is usually obtained through measurement and the measured data is denoted by .

To determine the source term , we require the following assumptions:

 (A) and ; (B) there exists a relation between the function and the measured data :

 where denotes the norm in the space and is the noise level. (C) The source term satisfies the a priori bound

 where is a positive constant, and denotes the norm in Sobolev space which is defined by [31] as follows:

We can refer to [31] for the details of the Sobolev space .

Using the separation of variables, we can obtain the explicit solution of the modified Helmholtz equation:

where

is an orthogonal basis in , and

According to (1.6) and the supplementary condition , we have

Based on this relation, we can define an operator on the space as

Thus, our inverse source problem is formulated as follows: give the data and the operator and then determine the unknown source term such that (1.10) holds.

It is straightforward that the operator is invertible and then the exact solution of (1.2) is

with

Note that the factor increases rapidly and tends to infinity as , so a small perturbation in the data may cause a dramatically large error in the solution . Therefore, this inverse source problem is mildly ill posed. It is impossible to obtain the unknown source using classical methods as above.

In [29], the simplified Tikhonov regularization method was given for (1.2). In this method, the regularization parameter is a priori chosen. It is well known that the ill-posed problem is usually sensitive to the regularization parameter and the a priori bound is usually difficult to be obtained precisely in practice. So the a priori choice rule of the regularization parameter is unreliable in practical problems. In this paper, we will present a Tikhonov-type regularization method to deal with (1.2) and show that the regularization parameter can be chosen by an a posteriori rule based on the discrepancy principle in [27].

The rest of this paper is organized as follows. In Section 2, we establish a quasinormal equation, which is crucial for proving the convergence of the Tikhonov-type regularization method. In Section 3, we give the Tikhonov-type regularization method and then prove the convergence of such method. Also, we give a priori and a posteriori choice rules to find the regularization parameter in the regularization method. In Section 4, we demonstrate a numerical example to illustrate the effectiveness of the method. In Section 5, we give some conclusions.

2. Preparation

In this section, we give an auxiliary result which will be used in this paper.

We first define an operator on as follows:

Let us observe that the operator is well defined and is a self-adjoint linear operator.

Next we give a lemma, which is important for discussing the regularization method. For simplicity, we denote the spaces and by and , respectively.

Lemma 2.1. Let be a linear and bounded operator between two Hilbert spaces, and let be defined as in (2.1), . Then for any the following Tikhonov functional

has a unique minimum , and this minimum is the unique solution of the quasinormal equation .

Proof. We divide the proof into three steps.

Step 1. The existence of a minimum of is proved. Let be a minimizing sequence; that is, as . We first need to show that is a Cauchy sequence. According to the definition of , we have

This implies that . Since the left-hand side converges to as tend to infinity. This shows that is a Cauchy sequence and thus convergent. Let , noting that . From the continuity of , we conclude that , that is, . This proves the existence of a minimum of .

Step 2. The equivalence of the quasinormal equation with the minimization problem for is shown. According to the definition and (2.1), we can obtain the following formula:

for all . If satisfies , then , that is, minimizes .

Conversely, if minimizes , then we substitute for any and , and then we can arrive at

Dividing both sides of the above inequality by and taking , we get

for all . This implies that . It follows that solves the quasinormal equation. From this, the equivalence of the quasinormal equation with the minimization problem for is shown exactly.

Step 3. We show that the operator is one-one for every . Let . Multiplication by yields , that is, .

3. A Tikhonov-Type Regularization Method

In this section, we first present a Tikhonov-type regularization method to obtain the approximate solution of (1.2) and then consider an a priori strategy and a posteriori choice rule to find the regularization parameter. Under each choice of the regularization parameter, the corresponding estimate can be obtained.

Since (1.2) is an ill-posed problem, we give its regularized solution which minimizes the Tiknonov functional

where the operator is defined as in (1.10), and is a regularization parameter.

According to Lemma 2.1, this minimum is the unique solution of the quasinormal equation , that is, . Because is a linear self-adjoint operator, that is, , we have the equivalent form of as

Further, the function can be reduced to

Now we are ready to formulate the main results of this paper. Before proceeding, the following lemmas are needed.

Lemma 3.1. For any , , it holds .

Proof. The proof is elementary and is omitted.

Lemma 3.2. For , , and the operator defined in (1.10), one has

Proof. By the Hölder inequality and Lemma 3.1, we have

The proof is completed.

In the following we give the corresponding convergence results for an a priori choice rule and an a posteriori choice rule.

3.1. An A Priori Choice Rule

Choose the regularization parameter as

The next theorem shows that the choice (3.6) is valid under suitable assumptions.

Theorem 3.3. Let be the minimizer of defined by (3.1) and be the exact solution of (1.2), and let assumptions (A), (B), and (C) hold. If is chosen by (3.6), then is convergent to the exact solution as the noise level tends to zero. Furthermore, one has the following estimate:

Proof. Since is the minimizer of defined by (3.1), we can obtain

Furthermore, we get

By Lemma 3.2, we have

The proof is completed.

3.2. An A Posteriori Choice Rule

Choose the regularization parameter as the solution of the equation

where the operator is defined by (1.10) and .

In the following theorem, an a posteriori rule based on the discrepancy principle [27] is considered in the convergence estimate.

Theorem 3.4. Let be the minimizer of defined by (3.1) and be the exact solution of (1.2), and let assumptions (A), (B), and (C) hold. If is chosen as the solution of (3.11), then is convergent to the exact solution as the noise level tends to zero. Furthermore, one has the following estimate:

Proof. Since is the minimizer of defined by (3.1), we can obtain

Consequently, it has

This leads to . It follows from Lemma 3.2 that the assertion of this theorem is true.

4. Numerical Tests

In this section, we present an example to illustrate the effectiveness and stability of our proposed method. The numerical results verify the validity of the theoretical results for the two cases of the a priori and a posteriori parameter choice rules.

Substituting (1.8) into (1.9), and then using trapezoid's rule to discretize (1.9) can result in the following discrete form:

where , , and .

We conduct two tests, and the tests are performed in the following way: first, from (1.6) and (1.11), we can select the source term and then . Consequently, the data function , and

We choose . Next, we add a random distributed perturbation to each data function, giving the vector

The function randn(·) generates arrays of random numbers whose elements are normally distributed with mean and variance . Thus, the total noise level can be measured in the sense of root mean square error according to

And can be obtained according to (1.12). Our error estimates use the relative error, which is given as follows:

where is given by (4.4).

Test 1. In the case of the a priori choice rule, we, respectively, compute with different parameters , , , , and . Tables 1 and 2 show that and have small influence on when they become larger. So, we always take and in this test. Table 3 shows for and10 with the perturbation , and 0.001. Table 4 shows for and8 with the perturbation and. In conclusion, the regularized solution well converges to the exact solution when tends to zero.

Table 1: Relative errors with , , , and for different .

10
50
100
200
400
800
1600

0.0471
0.0141
0.0116
0.0036
0.0039
0.0020
0.0018

Table 2: Relative errors with , , and for different .

1
5
10
20
40
80
100

0.0055
0.0040
0.0081
0.0067
0.0054
0.0046

Table 3: Relative errors with , , and , for different and .

0.1660
0.1345
0.0564
0.0241
0.0097
0.0076
0.0223

0.0283
0.0308
0.0221
0.0040
0.0063
0.0033
0.0021

0.0026
0.0048
0.0030
0.0018
0.0011

Table 4: Relative errors with , , and , for different and .

0.0412
0.0548
0.0304
0.0245

0.0133
0.0174
0.0064
0.0036

0.0032
0.0031

Test 2. In the case of the a posteriori choice rule (3.11), by taking , we also give the corresponding results as described in Test 1. The results can be easily seen from Tables 5, 6, 7, and 8.

Table 5: Relative errors with , , , and for different .

10
50
100
200
400
800
1600

0.0364
0.0239
0.0194
0.0093
0.0074
0.0060
0.0045

Table 6: Relative errors with , , , and for different .

1
5
10
20
40
80
100

1.0000
0.0041
0.0071
0.0086
0.0098
0.0094
0.0090

Table 7: Relative errors with , , and , for different and .

0.1160
0.1849
0.1986
0.1230
0.0179
0.0170
0.0337

0.0292
0.0128
0.0227
0.0026
0.0031
0.0040

0.0017
0.0027
0.0010
0.0025

0.0028

Table 8: Relative errors with , , and , for different and .

0.0250
0.0190
0.0315
0.0516

0.0076
0.0074
0.0030
0.0536

0.0047
0.0013

0.0071

From Tests 1 and 2, we conclude that the proposed regularization method is effective and stable.

5. Conclusion

In this paper, we proposed a Tikhonov-type regularization method to deal with the inverse source identification for the modified Helmholtz equation and set up a theoretical frame to analyze the convergence of such method. For instance, we provided the quasinormal equation to obtain the regularized solution. Moreover, besides the a priori parameter choice rule we studied an a posteriori rule for choosing the regularization parameter. Finally, we presented a numerical example whose results seem to be in excellent agreement with the convergence estimates of the method.

Acknowledgments

The authors would like to thank the referee for his (her) valuable comments and suggestions which improved this work to a great extent. This work was supported by the National Natural Science Foundation of China (NSFC) under Grant 40730424, the National Science and Technology Major Project under Grant 2011ZX05023-005, the National Science and Technology Major Project under Grant 2011ZX05044, the National Basic Research Program of China (973 Program) under Grant 2013CB329402, and the National Natural Science Foundation of China (NSFC) under Grant 11131006.

References

L. Greengard and M. C. Kropinski, “An integral equation approach to the incompressible Navier-Stokes equations in two dimensions,” SIAM Journal on Scientific Computing, vol. 20, no. 1, pp. 318–336, 1998.

C. J. S. Alves, J. B. Abdallah, and M. Jaoua, “Recovery of cracks using a point-source reciprocity gap function,” Inverse Problems in Science and Engineering, vol. 12, no. 5, pp. 519–534, 2004.

S. Andrieux and A. Ben Abda, “Identification of planar cracks by complete overdetermined data: inversion formulae,” Inverse Problems, vol. 12, no. 5, pp. 553–563, 1996.

J. M. Barry, Heat Source Determination in Waste Rock Dumps, World Scientific Publishing Co., River Edge, NJ, USA, 1998.

F. Yang and C.-L. Fu, “The method of simplified Tikhonov regularization for dealing with the inverse time-dependent heat source problem,” Computers & Mathematics with Applications, vol. 60, no. 5, pp. 1228–1236, 2010.

F. Geng and Y. Lin, “Application of the variational iteration method to inverse heat source problems,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2098–2102, 2009.

Z. Yi and D. A. Murio, “Source term identification in 1-D IHCP,” Computers & Mathematics with Applications, vol. 47, no. 12, pp. 1921–1933, 2004.

Z. Yi and D. A. Murio, “Identification of source terms in 2-D IHCP,” Computers & Mathematics with Applications, vol. 47, no. 10-11, pp. 1517–1533, 2004.

W. Cheng, L.-L. Zhao, and C.-L. Fu, “Source term identification for an axisymmetric inverse heat conduction problem,” Computers & Mathematics with Applications, vol. 59, no. 1, pp. 142–148, 2010.

J. I. Frankel, “Residual-minimization least-squares method for inverse heat conduction,” Computers & Mathematics with Applications, vol. 32, no. 4, pp. 117–130, 1996.

Y. C. Hon and T. Wei, “The method of fundamental solution for solving multidimensional inverse heat conduction problems,” Computer Modeling in Engineering and Sciences, vol. 7, no. 2, pp. 119–132, 2005.

M. Andrle and A. E. Badia, “Identification of multiple moving pollution sources in surface waters or atmospheric media with boundary observations,” Inverse Problems, vol. 28, no. 7, Article ID 075009, 2012.

N. Magnoli and G. A. Viano, “The source identification problem in electromagnetic theory,” Journal of Mathematical Physics, vol. 38, no. 5, pp. 2366–2388, 1997.

J. I. Frankel, “Constraining inverse Stefan design problems,” Zeitschrift für Angewandte Mathematik und Physik, vol. 47, no. 3, pp. 456–466, 1996.

A. Schuhmacher, J. Hald, K. B. Rasmussen, and P. C. Hansen, “Sound source reconstruction using inverse boundary element calculations,” Journal of the Acoustical Society of America, vol. 113, no. 1, pp. 114–127, 2003.

M. Hämäläinen, R. Hari, R. J. Ilmoniemi, J. Knuutila, and O. V. Lounasmaa, “Magnetoencephalography—theory, instrumentation, and applications to noninvasive studies of the working human brain,” Reviews of Modern Physics, vol. 65, no. 2, pp. 413–497, 1993.

S. He and V. G. Romanov, “Identification of dipole sources in a bounded domain for Maxwell's equations,” Wave Motion, vol. 28, no. 1, pp. 25–40, 1998.

Yu. E. Anikonov, B. A. Bubnov, and G. N. Erokhin, Inverse and Ill-Posed Sources Problems, Inverse and Ill-Posed Problems Series, VSP, Utrecht, The Netherlands, 1997.

A. E. Badia and T. Nara, “Inverse source problem for Helmholtz's equation from the Cauchy data with a single wave number,” Inverse Problems, vol. 27, Article ID 105001, 2011.

Y. C. Hon, M. Li, and Y. A. Melnikov, “Inverse source identification by Green's function,” Engineering Analysis with Boundary Elements, vol. 34, no. 4, pp. 352–358, 2010.

V. Isakov, Inverse Source Problems, vol. 34 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 1990.

M. Chapouly and M. Mirrahimi, “Distributed source identification for wave equations: an offline observer-based approach,” IEEE Transactions on Automatic Control, vol. 57, pp. 2067–2073, 2012.

S. Acosta, S. Chow, J. Taylor, and V. Villamizar, “On the multi-frequency inverse source problem in heterogeneous media,” Inverse Problems, vol. 28, Article ID 075013, 2012.

R. A. Adams, Sobolev Spaces, Pure and Applied Mathematics, Academic Press, New York, NY, USA, 1975.

A. El Badia and T. Ha Duong, “Some remarks on the problem of source identification from boundary measurements,” Inverse Problems, vol. 14, no. 4, pp. 883–891, 1998.

G. A. Kriegsmann and W. E. Olmstead, “Source identification for the heat equation,” Applied Mathematics Letters, vol. 1, no. 3, pp. 241–245, 1988.

H. W. Engl, M. Hanke, and A. Neubauer, Regularization of Inverse Problems, vol. 375 of Mathematics and its Applications, Kluwer Academic, Dordrecht, The Netherlands, 1996.

A. N. Tiknonov and V. Y. Arsenin, Solutions of Ill-Posed Problems, John Wiley & Sons, New York, NY, USA, 1977.

F. Yang, H. Z. Guo, and X. X. Li, “The simplified Tikhonov regularization method for identifying the unknown source for the modified helmholtz equation,” Mathematical Problems in Engineering, vol. 2011, Article ID 953492, 14 pages, 2011.

Y.-J. Ma, C.-L. Fu, and Y.-X. Zhang, “Identification of an unknown source depending on both time and space variables by a variational method,” Applied Mathematical Modelling, vol. 36, no. 10, pp. 5080–5090, 2012.

A. Kirsch, An Introduction to the Mathematical Theory of Inverse Problems, vol. 120, Springer, New York, NY, USA, 1996.