Choice of the regularization parameter for the Cauchy problem for the Laplace equation problem for the Laplaceequation

Purpose – In this paper, the Cauchy-type problem for the Laplace equation was solved in the rectangular domain with the use of the Chebyshev polynomials. The purpose of this paper is to present an optimal choice of the regularization parameter for the inverse problem, which allows determining the stable distribution of temperatureononeof theboundariesof therectangle domainwith therequired accuracy. Design/methodology/approach – The Cauchy-type problem is ill-posed numerically, therefore, it has been regularized with the use of the modi ﬁ ed Tikhonov and Tikhonov – Philips regularization. The in ﬂ uence of the regularization parameter choice on the solution was investigated. To choose the regularization parameter, the Morozov principle,the minimumofenergyintegralcriterion and the L-curve method wereapplied. Findings – Numerical examples for the function with singularities outside the domain were solved in this paper. The values of results change signi ﬁ cantly within the calculation domain. Next, results of the sought temperature distributions, obtained withthe use of different methods of choosing the regularizationparameter, werecompared. Methodsofchoosing the regularizationparameter wereevaluated bythenorm N max . Practical implications – Calculation model described in this paper can be applied to determine temperature distribution on the boundary of the heated wall of, for instance, a boiler or a body of the turbine, that is,everywherethe temperature measurement is impossibleto be performed ona part of theboundary. Originality/value – The paper presents a new method for solving the inverse Cauchy problem with the use of the Chebyshev polynomials. The choice of the regularization parameter was analyzed to obtain a solutionwith thelowestpossible sensitivity to input data disturbances.


Nomenclature a
= multinomial coefficient of the function of distribution of temperatureT w ð Þ; c = multinomial coefficient of the function of distribution of temperature T (x, y); E (a) = functional, energy integral; J a = regularizing functional (J a = J þ a 2 I ); m = number of Chebyshev nodes on the y-axis; n = number of Chebyshev nodes on the x-axis; N 1 -1 = degree of the polynomial describing unknown distribution of temperature on the C 1 boundary; N max = norm; q = heat flux density, W/m 2 K; T = temperature, K; T = temperature, function dependent on the Chebyshev node; w = Chebyshev node; W i = Chebyshev polynomial of the first kind of i-th degree; x, y = Cartesian coordinates; [x] n = integer part of the division of number x by n; and x mod n = remainder of the division of number x by n.
Greek symbols a = regularization parameter; d = error; d M = error of measurement data (Morozov principle); g = multinomial coefficient, pertains to the sought temperature distribution on the boundary C 1 ; C = boundary of the domain X, (C = C 1 | C 2 | C 3 | C 4 ); and X = calculation domain.
Subscript A = analytical solution; c = calculated value; m = measured value; ran = randomly disturbed value; and C i = on the boundary C i (for i = 1, 2, 3 and 4).

Introduction
Inverse problems are ill-posed in the Hadamard (1902) sense. It means that a slight disturbance to measurement data results in significant errors of the obtained results (Joachimiak and Ciałkowski, 2017, 2018Nowak, 2017). Therefore, problems of such type need to be regularized. There are many methods used to regularize inverse problems. Among them, there is the Tikhonov regularization (Beck and Woodbury, 2016;Chen et al., 2019;Djerrar et al., 2017;Frąckowiak et al., 2019a;Laneev, 2018;Marin, 2010Marin, , 2016Niu et al., 2014;Sun, 2016;Tikhonov and Arsenin, 1977;Yaparova, 2016), the Tikhonov-Philips regularization (Joachimiak et al., 2019a), the discrete Fourier transform (Frąckowiak and Ciałkowski, 2018;Wr oblewska et al., 2015) and SVD algorithm (Hasanov and Mukanova, 2015). In her article, Cheruvu (2017) applied the wavelet regularization of Laplace's equation in the arbitrarily shaped domain. The solution to the Cauchy problem for the Laplace's equation was also sought with the use of the iterative Tikhonov-type method (Delvare and Cimetière, 2017). Han et al. (2011) in the article presented numerical tests concerning the solution to the Cauchy problem for Laplace's equation with the use of the energy regularization method. Obtained results were compared with the Tikhonov regularization for which the regularization parameter was chosen based on the Morozov principle. In the paper of , the Cauchy problem for the Laplace's equation was solved with the use of the method of fundamental solutions and the energy regularization technique to choose the source points. The Laplace's equation was also solved with the use of iterative algorithms (Frąckowiak et al., 2015a(Frąckowiak et al., , 2015b, of the Trefftz method (Ciałkowski and Frąckowiak, 2002;Ciałkowski and Grysa, 2010;Grysa et al., 2012;Ho_ zejowski, 2016;Lin et al., 2018), of the method of fundamental solution (Kołodziej and Mierzwiczak, 2008;Mierzwiczak et al., 2015;Mierzwiczak and Kołodziej, 2011) and of the collocation method (Joachimiak et al., 2016). In many cases, the regularization of the inverse problem concerns the HFF 30,10 problem of choosing the regularization parameter. The regularization parameter can be chosen based on the Morozov principle (Chen et al., 2019;Han et al., 2011;Joachimiak et al., 2019a;Marin, 2016;Morozov, 1984;Sun, 2016) or using the L-curve method (Jin and Zheng, 2006;Marin and Munteanu, 2010;Marin, 2005). In the study of Marin (2011), the optimal regularization parameter was sought based on the generalized cross-validation criterion. Currently, research work focuses on finding new methods of regularization (Cheng and Feng, 2014;Zhuang and Chen, 2017) and on the modification of already known and used methods (Yang et al., 2015;Zheng and Zhang, 2018). Because of a wide application of inverse problems in engineering problems, such as the cooling of the blades in gas turbines (Frąckowiak et al., 2017;Frąckowiak et al., 2019b;Frąckowiak et al., 2011), analysis of the boiling heat transfer in minichannels (Ho_ zejowska et al., 2009;Maciejewska and Piasecka, 2017), analysis of thermal and thermo-chemical treatment (Joachimiak et al., 2019b) or monitoring of power boilers operation (Taler et al., 2016(Taler et al., , 2017, developing methods for regularization of these problems and investigating the process of choosing the regularization parameter are very significant. In this article, the solution to the Cauchy problem for the Laplace's equation was investigated with the use of the Chebyshev polynomials. To regularize the solution to the inverse problem, the modification of the Tikhonov and of the Tikhonov-Philips regularizations, described in the article (Joachimiak et al., 2019a) was applied. The choice of the regularization parameter was made based on the Morozov principle, the minimum energy integral criterion and the L-curve method.

Calculation model
In many technical applications, it is impossible to measure temperature on the boundary of the heated component of the device or machine, such as the combustion chamber, the inner side of the body of a turbine or a boiler. Then, the distribution of temperature can be determined by finding the solution to the inverse problem. Based on the distribution of temperature on the part of the boundary T C 2 ; T C 3 ; T C 4 ; fig: 1 À Á and, additionally, knowing the heat flux density on the part of the boundary q C 3 ; fig: 1 À Á one can determine the distribution of temperature on the boundary, where it is impossible to measure this temperature T C 1 ; fig: 1 À Á . Such a posed problem is the Cauchy problem, particularly sensitive to errors in measurement and in the calculation. In the stationary thermal field, the heat equation is reduced to the Laplace's equation (for the non-linear case, the Kirchhoff's substitution transforms the heat equation into the Laplace's equation).
with the following boundary conditions ( Figure 1): Cauchy problem for the Laplace equation It was assumed that the solution can be noted as the linear combination of the Chebyshev polynomials (Paszkowski, 1975).
To solve the Cauchy problem, the collocation method was used. It was assumed that there are n points along the x-axis and m points along the y-axis (including points on the boundary). Collocation points being inside the interval (À1, 1) are the Chebyshev nodes (Paszkowski, 1975). Nodes were renumbered, which enables the temperature function to be noted in the following equation (7): where the coefficients a k (k = 1, 2, . . ., mn) are unknown. Sought temperature distribution on the boundary C 1 was assumed as the linear combination of the Chebyshev polynomials (Paszkowski, 1975).
Coefficients a k (k = 1, 2, . . ., mn) are expressed by the values of coefficients g i (i = 1, 2, . . ., N 1 ). Hence, the determination of the temperature distribution is reduced to the determination of coefficients g i . To do so, the functional of the following form was minimized.
where c in subscript denotes the calculated value, while m denotes the measured value. The integral in equation (9) on the boundary C 2 can be noted in the following equation (10): Calculation domain HFF 30,10 Applying numerical integration we have: Having inserted the equation (7) into the equation (11), we obtained: Solving the direct problem, where the temperature on boundaries C 1 , C 2 , C 3 and C 4 was known, was reduced to solving the matrix equation what was described in detail in the paper (Joachimiak et al., 2019a). Based on the solution to the direct problem, constants a k [Equation (12)] are of the following equation (14): where H h;k ¼ X mÀ1 j¼2Ã k;jn W hÀ1 y j ð Þ , whileÃ k;jn (k = 1, 2, . . ., mn; j = 2, 3, . . ., m À 1) are elements of the matrix A À1 . After substituting equation (14) into equation (12) we obtained: We would like the integral J C 2 to have a value equal to zero or as close to zero as possible, hence, we equate the squared expression [Equation (15) We obtain n linear equations of the following equation (17): for i = 1, 2, . . ., n and h = 1, 2, . . ., N 1 . It can be reduced to the matrix equation.
Similarly, for other integrals [Equation (9)] we obtained: After J C 3 , J q;C 3 and J C 4 had been equated to zero, equations of the following forms were obtained: Based on equations (16) and (22)-(24), an oversized system of linear equations was obtained as the matrix equation, which would be solved with the use of the SVD algorithm: what can be noted in the shorter form: Because of a great sensitivity of results to disturbances to measurement data, the Cauchy problem was regularized. The regularizing functional of the following form was assumed: Regularization term can be noted as the sum of integrals.
C 1i . Performing numerical integration using the trapezoidal rule, we obtained: On the boundary C 1 we have: Each of the components [Equation (32)] was equated to zero. The equation of the following form was obtained: Equations (33)-(36) can be reduced to the following system of equations.
where a is the regularization parameter. When the equation (26) and regularization [Equation (37)] are included, the following system of equations is obtained: The solution to the system of equations (38) was sought in the least-squares sense with the use of the SVD algorithm.

HFF 30,10
3. Choice of the regularization parameter To determine unknown regularization parameter a, the Morozov principle, the minimum of energy integral criterion and L-curve method were applied. For the solution obtained with the use of the Morozov principle, the mean (Morozov_A) and the maximal (Morozov_B) errors of the heat flux density d M on the boundary C 3 were evaluated. Interval halving method was used to determine zero of the function F M (a) defined by the following equation (39): Unknown regularization parameter a was also sought based on the minimization of the functional (energy integral) of the following equation (40): The minimum of the energy integral corresponds to satisfying the Laplace's equation (with respective boundary conditions), which is discussed in the paper (Gelfand and Fomin, 1979). Therefore, we choose the parameter a for which min a E a ð Þ occurs, i.e. the derivative E 0 ¼ dE da reverses the sign. Domain X was divided into rectangular domains with the use of equidistant nodes and next into domains being right-angled triangles. The integral value was calculated with the use of the finite element method. Value rT was determined based on the form of the solution equation (6). Values g h were obtained by solving the equation (38).
On the basis of the solution of the equation (38), the L-curve was drawn as the correlation between k[BM]{g } À {BW}k on the x-axis and k[CM]{g } À {CW}k on the y-axis. We sought for the regularization parameter a with which corresponded the point of the L-curve locating on the curvature of this line. To evaluate the choice of the regularization parameter a, the following norm was defined:

Numerical examples
Calculations were made in the domain X for the function.
and  (26)], without disturbance (d ran = 0) and with random disturbance to the heat flux density up to 0.01q (d ran = 0.01) and to 0.02q (d ran = 0.02) are summarized in Table I. Disturbance q is an additive function with the uniform distribution. A slight disturbance to measurement data results in a significant error of the sought temperature on the boundary C 1 . Hence, it is necessary to regularize the inverse problem [Equation (38)] and to choose the regularization parameter a properly.

Example 1
Calculations were made for the function f 1 [Equation (42)]. Heat flux density was disturbed randomly to 0.02q (d ran = 0.02). Regularization parameter a was chosen with the use of the Morozov principle, the minimum of energy integral criterion and the Lcurve method. Figure 2 presents the course of the energy integral and its derivative depending on the parameter a. To solve the Cauchy problem, the authors applied such value of the regularization parameter a for which the energy integral E (a) took the minimal value, which meant that the derivative of the energy integral E 0 ¼ dE da reversed the sign. The value of a was 5.13 Â 10 À4 (Table II). The Cauchy problem was also regularized for Table I. Values of the norm N max for calculations without regularization (a = 0), without disturbance (d ran = 0) and with random disturbance to the heat flux density up to 0.01q (d ran = 0.01) and 0.02q (d ran = 0.02) Error of heat flux density disturbed randomly HFF 30,10 the regularization parameter a amounting to 4.64 Â 10 À1 , which was determined based on the L-curve course (Figure 3). To choose the regularization parameter with the Morozov principle, the values of mean and maximal error d M were evaluated for the heat flux on the boundary C 3 . The respective values were obtained: 0.008 and 0.02 (Table II). Next, zero of the function F M (a) was calculated as per the equation (39). The respective values of the regularization parameter were obtained: 1.98 Â 10 À6 and 4.0799 Â 10 À2 (Table II).
The lowest value of the norm N max amounting to 6.18 Â 10 À2 (Table II) for the function f 1 was obtained for the case of choosing the regularization parameter with the use of the Morozov principle for the maximal error of the heat flux d M (Morozov_B). This criterion brought satisfying results, as did the choice of the regularization parameter made with the use of the minimum energy integral criterion (N max = 9.796 Â 10 À2 ). When the L-curve method was used, the obtained results were considerably worse (N max = 2.22 Â 10 À1 ). For the Morozov principle, for the mean error d M of the heat flux (Morozov_A), the highest value of the norm N max amounting to 50.42 was obtained. Distributions of temperature on the boundary C 1 resulting from the analytical solution (AS) and from the solution of the Cauchy problem are presented in Figure 4.

Example 2
Calculations were made for the function f 2 [Equation (43)]. Heat flux density was disturbed randomly to 0.02q (d ran = 0.02). The best results were obtained for the choice of the regularization parameter made with the use of the minimum energy integral criterion   Table III. Distributions of temperature on the boundary C 1 for the AS and for the solution to the Cauchy problem with regularization are presented in Figure 5. Distribution of temperature on the boundary C 1 obtained with the use of the minimum of energy integral criterion slightly diverges from the AS. To examine thoroughly the criterion for the regularization parameter selection with the use of the minimum of energy integral, calculations were performer also for the following functions: f 6 ¼ e x 2 Ày 2 sin2xy (47) which were chosen based on publications Conde Mones et al., 2017;Fu et al., 2013;Sun, 2017). Values of the regularization parameter and of the norm N max for functions f 3 À f 6 are summarized in Table IV. For the function f 6 and the disturbance to the heat flux density d ran = 0.02 and d ran = 0.05 the minimum of energy integral was not achieved. Distributions of temperature on the boundary C 1 being sought are presented in Figure 6. For functions f 3 À f 5 the disturbance d ran = 0.05 was taken into account, and for the function f 6 it was d ran = 0.01.

Conclusion
This paper presents the solution of the Cauchy problem for Laplace's equation.
Obtained distributions of temperature on the boundary C 1 were analyzed in terms of the dependence on the method for choosing the regularization parameter. The best results were obtained for the choice of the regularization parameter made with the use of the minimum of energy integral criterion and the Morozov principle (d M is the maximal error for the heat flux on the boundary C 3 ). The advantage of the application of the minimum energy integral criterion is a unique determination of the regularization parameter a for which E (a) has minimal value. Regularization made with the use of the minimum energy integral criterion gives satisfying results. However, its disadvantage is the fact that not for all calculation examples the minimum energy integral was determined. For the Morozov principle, the obtained  and of the L-curve method (L-curve) for the function f 2 Cauchy problem for the Laplace equation results (distribution of temperature on the boundary C 1 ) depend on calculation or evaluation of the value of the heat flux d M error, what is a disadvantage of this method. Inexact evaluation of the measurement data error can result in obtaining the distribution of temperature on the boundary C 1 , which is subject to great uncertainty. The choice of the parameter a made with the use of the L-curve method gave the worst results. Smooth L-curve course was obtained, what was related to the problem with the unique determination of the regularization parameter a using this method. 3.46 Â 10 À2 1.0 Â 10 À3 3.98 Â 10 À2 1.0 Â 10 À3 6.406 Â 10 À2 f 5 4.99 Â 10 À4 3.46 Â 10 À2 1.0 Â 10 À3 3.98 Â 10 À2 1.0 Â 10 À2 2.16 Â 10 À1 f 6 1.0 Â 10 À3 2.33 Â 10 À1 No minimum No minimum Figure 6. Distribution of temperature on the boundary C 1 obtained from the AS and form the solution to the Cauchy problem when the regularization parameter a was chosen with the use of the minimum of energy integral criterion (E) for functions f 2 À f 6