\newsiamremark

remarkRemark \newsiamremarkhypothesisHypothesis \newsiamthmclaimClaim \headersConvexification for the 3D Problem of Travel Time Tomography M. V. Klibanov, J. Li, V. G. Romanov, and Z. Yang

Convexification for the 3D Problem of Travel Time Tomography

Michael V. Klibanov Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC, 28223, USA ([email protected])    Jingzhi Li Department of Mathematics & National Center for Applied Mathematics Shenzhen & SUSTech International Center for Mathematics, Southern University of Science and Technology, Shenzhen 518055, P. R. China ([email protected])    Vladimir G. Romanov Sobolev Institute of Mathematics, Novosibirsk 630090, Russian Federation ([email protected])    Zhipeng Yang School of Mathematics and Statistics, Lanzhou University, Lanzhou 730000, P. R. China ([email protected])
Abstract

The travel time tomography problem is a coefficient inverse problem for the eikonal equation. This problem has well known applications in seismic. The eikonal equation is considered here in the circular cylinder, where point sources run along its axis and measurements of travel times are conductes on the whole surface of this cylinder. A new version of the globally convergent convexification numerical method for this problem is developed. Results of numerical studies are presented.

Key Words: eikonal equation, geodesic lines, coefficient inverse problem, travel time tomography in 3d, globally convergent numerical method, Carleman estimate, numerical solution

2020 MSC codes: 35R30.

1 Introduction

The travel time tomography problem (TTTP) has a long history since it has broad applications in seismic studies of Earth, see, e.g. [21, 24]. First, we refer here to the pioneering works of Herglotz in 1905 [7] and Wiechert and Zoeppritz in 1907 [25], who have solved this problem in the 1d case. There are some important uniqueness and stability results for the TTTP. Since we are focused here on the numerical aspect, we cite only a few of those [17, 18, 19, 21]. We also refer to [22, 26] for numerical studies of the TTTP by the methods, which are significantly different from the convexification method of this paper.

Our goal here is to present the theory and computational experiments for a new version of the globally convergent convexification numerical method for the TTTP in the 3d case. The theory of the first version of the convexification for this problem was published in [12], [14, Chapter 11] and numerical studies were published in [15]. However, it was assumed in these references that the domain of interest ΩΩ\Omegaroman_Ω is a vertically elongated rectangular prism, and point sources run along an interval of a straight line, which is located in a horizontal plane “below” Ω.Ω\Omega.roman_Ω . As a result, only horizontally oriented abnormalities were imaged in [15].

Since we do not want to be restricted to only horizontally oriented abnormalities, then we consider in this paper the case when ΩΩ\Omegaroman_Ω is a vertically oriented circular cylinder with point sources running along its axis and measurements of travel times being conducted on the lateral surface, top and bottom boundaries of this cylinder, see Figure 1. This source/detector configuration has applications in seismic imaging of wells. Naturally, our source/detector configuration requires some changes in the theory of [12, 15].

Refer to caption
Figure 1: A schematic diagram of the source/detectors of this paper. The large circle is the outer boundary of the domain ΩΩ\Omegaroman_Ω defined in (3), the small circle is the inner boundary of domain ΩεsubscriptΩ𝜀\Omega_{\varepsilon}roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT defined in (4 ), and the red line with red star is position of the source Lsourcesubscript𝐿sourceL_{\text{source}}italic_L start_POSTSUBSCRIPT source end_POSTSUBSCRIPT defined in (5).

Below we first formulate and prove a new theorem about the monotonicity of the solution of the eikonal equation with respect to the radius of the cylindrical coordinate system. Next, using this theorem, we present a new version of the convexification method, which is adapted to the sources/detectors configuration of Figure 1. Next, we formulate some theorems of the convergence analysis. Since their close analogs were proven in [15], then we do not prove them here. Finally, we describe our numerical studies and present their results.

Definition. We call a numerical method for a Coefficient Inverse Problem (CIP) globally convergent, if there exists a theorem claiming that this method delivers points in a sufficiently small neighborhood of the true solution of that problem without any advanced knowledge of this neighborhood. The size of this neighborhood depends only on the noise level in the input data for this CIP. In other words, a good first guess about the true solution is not required.

The first theoretical results about the convexification method for CIPs were obtained in [9, 10]. More recent works on this method for various CIPs contain both theory and numerical results, see, e.g. [1, 2, 13, 14, 15, 16] and references cited therein. The development of the convexification is caused by the desire to avoid the phenomenon of local minima and ravines of conventional least squares cost functionals for CIPs, see, e.g. [3, 4, 5]. In the convexification, a globally strongly cost functional is constructed. The key element of this functional is the presence of the Carleman Weight Function (CWF) in it. This is the function, which is used as the weight in the Carleman estimate for an associated operator. Usually this is a partial differential operator. In the case of the TTTP, however, this is a Volterra-like integral operator. Finally, it is proven that the gradient descent method of the minimization of this functional converges globally to the true solution of the corresponding CIP as long as the level of noise in the data tends to zero.

All functions considered below are real valued ones. In section 2 we pose our inverse problem. In section 3 we prove the above mentioned important theorem about the monotonicity of the travel time function with respect to the radius of the cylindrical coordinate system. In section 4 we derive an integral differential equation, which does not contain the unknown coefficient. In section 5 we present an important element of our method: the orthonormal basis of [11], [14, section 6.2.3]. In section 6 we derive a boundary value problem and also prove a Carleman estimate for the Volterra integral operator. In section 7 we introduce partial finite differences for our problem. In section 8 we present the convexification functional and formulate theorems of its global convergence analysis. Section 9 is devoted to numerical studies.

2 Statement of the Problem

Below 𝐱=(x,y,z)𝐱𝑥𝑦𝑧\mathbf{x}=\left(x,y,z\right)bold_x = ( italic_x , italic_y , italic_z ) denotes points in 3.superscript3\mathbb{R}^{3}.blackboard_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT . Let c(𝐱)𝑐𝐱c\left(\mathbf{x}\right)italic_c ( bold_x ) be the speed of waves propagation, c(𝐱)=1/n(𝐱)𝑐𝐱1𝑛𝐱c\left(\mathbf{x}\right)=1/n\left(\mathbf{x}\right)italic_c ( bold_x ) = 1 / italic_n ( bold_x ), where n(𝐱)𝑛𝐱n\left(\mathbf{x}\right)italic_n ( bold_x ) is the refractive index. The function n(𝐱)𝑛𝐱n\left(\mathbf{x}\right)italic_n ( bold_x ) generates the Riemannian metric [21, Chapter 3]

(1) dτ=n(𝐱)(dx)2+(dy)2+(dz)2.𝑑𝜏𝑛𝐱superscript𝑑𝑥2superscript𝑑𝑦2superscript𝑑𝑧2d\tau=n\left(\mathbf{x}\right)\sqrt{\left(dx\right)^{2}+\left(dy\right)^{2}+% \left(dz\right)^{2}}.italic_d italic_τ = italic_n ( bold_x ) square-root start_ARG ( italic_d italic_x ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_d italic_y ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_d italic_z ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

Let 𝐱0subscript𝐱0\mathbf{x}_{0}bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT be a source of waves and 𝐱𝐱\mathbf{x}bold_x be an observation point. The time, which the wave needs to propagate from a source 𝐱0subscript𝐱0\mathbf{x}_{0}bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to a point 𝐱,𝐱\mathbf{x,}bold_x , is called the “first arrival time” or “travel time”. In the case of a heterogeneous medium with n(𝐱)const.,𝑛𝐱𝑐𝑜𝑛𝑠𝑡n\left(\mathbf{x}\right)\neq const.,italic_n ( bold_x ) ≠ italic_c italic_o italic_n italic_s italic_t . , the first arriving signal, which arrives at the arrival time to the point 𝐱𝐱\mathbf{x}bold_x, travels along the so-called “geodesic line” Γ(𝐱,𝐱0)Γ𝐱subscript𝐱0\Gamma\left(\mathbf{x},\mathbf{x}_{0}\right)roman_Γ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) generated by metric (1). The first arrival time τ(𝐱,𝐱0)𝜏𝐱subscript𝐱0\tau\left(\mathbf{x},\mathbf{x}_{0}\right)italic_τ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is

τ(𝐱,𝐱0)=Γ(𝐱,𝐱0)n(𝐲(s))𝑑σ,𝜏𝐱subscript𝐱0subscriptΓ𝐱subscript𝐱0𝑛𝐲𝑠differential-d𝜎\tau\left(\mathbf{x},\mathbf{x}_{0}\right)=\int\limits_{\Gamma\left(\mathbf{x}% ,\mathbf{x}_{0}\right)}n\left(\mathbf{y}\left(s\right)\right)d\sigma,italic_τ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∫ start_POSTSUBSCRIPT roman_Γ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_POSTSUBSCRIPT italic_n ( bold_y ( italic_s ) ) italic_d italic_σ ,

where dσ𝑑𝜎d\sigmaitalic_d italic_σ is the element of the euclidean length. The function τ(𝐱,𝐱0)𝜏𝐱subscript𝐱0\tau\left(\mathbf{x},\mathbf{x}_{0}\right)italic_τ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) satisfies the following eikonal equation [21, Chapter 3]:

(2) (𝐱τ)2=n2(𝐱),τ(𝐱,𝐱0)=O(|𝐱𝐱0|) as 𝐱𝐱0.formulae-sequencesuperscriptsubscript𝐱𝜏2superscript𝑛2𝐱𝜏𝐱subscript𝐱0𝑂𝐱subscript𝐱0 as 𝐱subscript𝐱0\begin{split}&\hskip 28.45274pt\left(\nabla_{\mathbf{x}}\tau\right)^{2}=n^{2}% \left(\mathbf{x}\right),\\ &\tau\left(\mathbf{x},\mathbf{x}_{0}\right)=O\left(\left|\mathbf{x}-\mathbf{x}% _{0}\right|\right)\text{ as }\mathbf{x}\rightarrow\mathbf{x}_{0}.\end{split}start_ROW start_CELL end_CELL start_CELL ( ∇ start_POSTSUBSCRIPT bold_x end_POSTSUBSCRIPT italic_τ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( bold_x ) , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL italic_τ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_O ( | bold_x - bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT | ) as bold_x → bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT . end_CELL end_ROW

Consider two numbers R,B>0𝑅𝐵0R,B>0italic_R , italic_B > 0 and let ε(0,R)𝜀0𝑅\varepsilon\in\left(0,R\right)italic_ε ∈ ( 0 , italic_R ) be a small number. We define the domains Ω,Ωε3ΩsubscriptΩ𝜀superscript3\Omega,\Omega_{\varepsilon}\subset\mathbb{R}^{3}roman_Ω , roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ⊂ blackboard_R start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT as:

(3) Ω={𝐱=(x,y,z):x2+y2<R,z(0,B)},Ωconditional-set𝐱𝑥𝑦𝑧formulae-sequencesuperscript𝑥2superscript𝑦2𝑅𝑧0𝐵\displaystyle\Omega=\left\{\mathbf{x}=\left(x,y,z\right):\sqrt{x^{2}+y^{2}}<R,% z\in\left(0,B\right)\right\},roman_Ω = { bold_x = ( italic_x , italic_y , italic_z ) : square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG < italic_R , italic_z ∈ ( 0 , italic_B ) } ,
(4) Ωε={𝐱=(x,y,z):x2+y2<ε,z(0,B)},subscriptΩ𝜀conditional-set𝐱𝑥𝑦𝑧formulae-sequencesuperscript𝑥2superscript𝑦2𝜀𝑧0𝐵\displaystyle\Omega_{\varepsilon}=\left\{\mathbf{x}=\left(x,y,z\right):\sqrt{x% ^{2}+y^{2}}<\varepsilon,z\in\left(0,B\right)\right\},roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT = { bold_x = ( italic_x , italic_y , italic_z ) : square-root start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG < italic_ε , italic_z ∈ ( 0 , italic_B ) } ,
(5) Lsource={𝐱0=(x0,y0,z0):x0=y0=0,z0(0,B)}.subscript𝐿sourceconditional-setsubscript𝐱0subscript𝑥0subscript𝑦0subscript𝑧0formulae-sequencesubscript𝑥0subscript𝑦00subscript𝑧00𝐵\displaystyle\hskip 5.69046ptL_{\text{source}}=\left\{\mathbf{x}_{0}=\left(x_{% 0},y_{0},z_{0}\right):x_{0}=y_{0}=0,z_{0}\in\left(0,B\right)\right\}.italic_L start_POSTSUBSCRIPT source end_POSTSUBSCRIPT = { bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) : italic_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_y start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) } .

Hence, ΩΩ\Omegaroman_Ω and ΩεsubscriptΩ𝜀\Omega_{\varepsilon}roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT are vertically oriented cylinders with the radii R𝑅Ritalic_R and ε𝜀\varepsilonitalic_ε respectively, ΩεΩ,subscriptΩ𝜀Ω\Omega_{\varepsilon}\subset\Omega,roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ⊂ roman_Ω , and Lsourcesubscript𝐿sourceL_{\text{source}}italic_L start_POSTSUBSCRIPT source end_POSTSUBSCRIPT is the vertical axis of both cylinders, along which sources 𝐱0subscript𝐱0\mathbf{x}_{0}bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT run.

Consider the cylindrical coordinates 𝐱=(x,y,z)=(rcosφ,rsinφ,z).𝐱𝑥𝑦𝑧𝑟𝜑𝑟𝜑𝑧\mathbf{x}=\left(x,y,z\right)=\left(r\cos\varphi,r\sin\varphi,z\right).bold_x = ( italic_x , italic_y , italic_z ) = ( italic_r roman_cos italic_φ , italic_r roman_sin italic_φ , italic_z ) . Therefore by (3) and (4) in cylindrical coordinates

(6) Ω={(r,φ,z):r(0,R),φ[0,2π),z(0,B)},Ωconditional-set𝑟𝜑𝑧formulae-sequence𝑟0𝑅formulae-sequence𝜑02𝜋𝑧0𝐵\displaystyle\Omega=\left\{\left(r,\varphi,z\right):r\in\left(0,R\right),% \varphi\in[0,2\pi),z\in\left(0,B\right)\right\},roman_Ω = { ( italic_r , italic_φ , italic_z ) : italic_r ∈ ( 0 , italic_R ) , italic_φ ∈ [ 0 , 2 italic_π ) , italic_z ∈ ( 0 , italic_B ) } ,
(7) Ωε={(r,φ,z):r(0,ε),φ[0,2π),z(0,B)},subscriptΩ𝜀conditional-set𝑟𝜑𝑧formulae-sequence𝑟0𝜀formulae-sequence𝜑02𝜋𝑧0𝐵\displaystyle\Omega_{\varepsilon}=\left\{\left(r,\varphi,z\right):r\in\left(0,% \varepsilon\right),\varphi\in[0,2\pi),z\in\left(0,B\right)\right\},roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT = { ( italic_r , italic_φ , italic_z ) : italic_r ∈ ( 0 , italic_ε ) , italic_φ ∈ [ 0 , 2 italic_π ) , italic_z ∈ ( 0 , italic_B ) } ,
(8) D=ΩΩε={(r,φ,z):r(ε,R),φ[0,2π),z(0,B)}.𝐷ΩsubscriptΩ𝜀conditional-set𝑟𝜑𝑧formulae-sequence𝑟𝜀𝑅formulae-sequence𝜑02𝜋𝑧0𝐵\displaystyle\hskip 14.22636ptD=\Omega\diagdown\Omega_{\varepsilon}=\left\{% \left(r,\varphi,z\right):r\in\left(\varepsilon,R\right),\varphi\in[0,2\pi),z% \in\left(0,B\right)\right\}.italic_D = roman_Ω ╲ roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT = { ( italic_r , italic_φ , italic_z ) : italic_r ∈ ( italic_ε , italic_R ) , italic_φ ∈ [ 0 , 2 italic_π ) , italic_z ∈ ( 0 , italic_B ) } .

We consider below the eikonal equation (2) only in the domain Ω.Ω\Omega.roman_Ω . We assume that the refractive index n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) satisfies the following conditions:

(9) n(𝐱)𝑛𝐱\displaystyle n(\mathbf{x})italic_n ( bold_x ) =1 in Ωε,absent1 in subscriptΩ𝜀\displaystyle=1\text{ in }\Omega_{\varepsilon},= 1 in roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ,
(10) n(𝐱)𝑛𝐱\displaystyle n(\mathbf{x})italic_n ( bold_x ) C1(Ω¯),absentsuperscript𝐶1¯Ω\displaystyle\in C^{1}\left(\overline{\Omega}\right),∈ italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( over¯ start_ARG roman_Ω end_ARG ) ,
(11) nr(r,φ,z)subscript𝑛𝑟𝑟𝜑𝑧\displaystyle n_{r}\left(r,\varphi,z\right)italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z ) 0 in Ω.absent0 in Ω\displaystyle\geq 0\text{ in }\Omega.≥ 0 in roman_Ω .

In cylindrical coordinates eikonal equation (2) has the form

(12) (τr)2+1r2(τφ)2+(τz)2=n2(r,φ,z).superscriptsubscript𝜏𝑟21superscript𝑟2superscriptsubscript𝜏𝜑2superscriptsubscript𝜏𝑧2superscript𝑛2𝑟𝜑𝑧(\tau_{r})^{2}+\frac{1}{r^{2}}(\tau_{\varphi})^{2}+(\tau_{z})^{2}=n^{2}(r,% \varphi,z).( italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ( italic_τ start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z ) .

By (5) denote

(13) τ(𝐱,𝐱0)=τ(𝐱,z0), 𝐱0Lsource, z0(0,B).formulae-sequence𝜏𝐱subscript𝐱0𝜏𝐱subscript𝑧0formulae-sequence subscript𝐱0subscript𝐿source subscript𝑧00𝐵\tau\left(\mathbf{x,x}_{0}\right)=\tau\left(\mathbf{x,}z_{0}\right),\text{ }% \mathbf{x}_{0}\in L_{\text{source}},\text{ }z_{0}\in\left(0,B\right).italic_τ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_τ ( bold_x , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ italic_L start_POSTSUBSCRIPT source end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) .

Travel Time Tomography Problem (TTTP). Let conditions (3)-(13) hold. Find the refractive index n(𝐱)𝑛𝐱n\left(\mathbf{x}\right)italic_n ( bold_x ) for 𝐱D,𝐱𝐷\mathbf{x}\in D,bold_x ∈ italic_D , assuming that the following functions p(φ,z,z0)𝑝𝜑𝑧subscript𝑧0p\left(\varphi,z,z_{0}\right)italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ), p0(r,φ,z0),pB(r,φ,z0)subscript𝑝0𝑟𝜑subscript𝑧0subscript𝑝𝐵𝑟𝜑subscript𝑧0p_{0}\left(r,\varphi,z_{0}\right),p_{B}\left(r,\varphi,z_{0}\right)italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) are known:

(14) p(φ,z,z0)𝑝𝜑𝑧subscript𝑧0\displaystyle p\left(\varphi,z,z_{0}\right)italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =τ(R,φ,z,z0), φ(0,2π),z,z0(0,B),formulae-sequenceabsent𝜏𝑅𝜑𝑧subscript𝑧0formulae-sequence for-all𝜑02𝜋for-all𝑧subscript𝑧00𝐵\displaystyle=\tau\left(R,\varphi,z,z_{0}\right),\text{ }\forall\varphi\in% \left(0,2\pi\right),\forall z,z_{0}\in\left(0,B\right),= italic_τ ( italic_R , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ∀ italic_φ ∈ ( 0 , 2 italic_π ) , ∀ italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) ,
(15) p0(r,φ,z0)subscript𝑝0𝑟𝜑subscript𝑧0\displaystyle p_{0}\left(r,\varphi,z_{0}\right)italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =τ(r,φ,0,z0),r(ε,R), φ(0,2π),z0(0,B),formulae-sequenceabsent𝜏𝑟𝜑0subscript𝑧0formulae-sequencefor-all𝑟𝜀𝑅formulae-sequence for-all𝜑02𝜋for-allsubscript𝑧00𝐵\displaystyle=\tau\left(r,\varphi,0,z_{0}\right),\forall r\in\left(\varepsilon% ,R\right),\text{ }\forall\varphi\in\left(0,2\pi\right),\forall z_{0}\in\left(0% ,B\right),= italic_τ ( italic_r , italic_φ , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ∀ italic_r ∈ ( italic_ε , italic_R ) , ∀ italic_φ ∈ ( 0 , 2 italic_π ) , ∀ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) ,
(16) pB(r,φ,z0)subscript𝑝𝐵𝑟𝜑subscript𝑧0\displaystyle p_{B}\left(r,\varphi,z_{0}\right)italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =τ(r,φ,B,z0),r(ε,R),φ(0,2π),z0(0,B).formulae-sequenceabsent𝜏𝑟𝜑𝐵subscript𝑧0formulae-sequencefor-all𝑟𝜀𝑅formulae-sequencefor-all𝜑02𝜋for-allsubscript𝑧00𝐵\displaystyle=\tau\left(r,\varphi,B,z_{0}\right),\forall r\in\left(\varepsilon% ,R\right),\forall\varphi\in\left(0,2\pi\right),\forall z_{0}\in\left(0,B\right).= italic_τ ( italic_r , italic_φ , italic_B , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ∀ italic_r ∈ ( italic_ε , italic_R ) , ∀ italic_φ ∈ ( 0 , 2 italic_π ) , ∀ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) .

Everywhere below we assume that the following condition is valid:

Assumption of the Regularity of Geodesic Lines. For any pair of points (𝐱,𝐱0)Ω¯×Lsource𝐱subscript𝐱0¯Ωsubscript𝐿source\left(\mathbf{x},\mathbf{x}_{0}\right)\in\overline{\Omega}\times L_{\text{% source}}( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ over¯ start_ARG roman_Ω end_ARG × italic_L start_POSTSUBSCRIPT source end_POSTSUBSCRIPT there exists unique geodesic line Γ(𝐱,𝐱0)Γ𝐱subscript𝐱0\Gamma\left(\mathbf{x},\mathbf{x}_{0}\right)roman_Γ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) connecting these two points.

3 The Monotonicity Property of the Function τ(r,φ,z,z0)𝜏𝑟𝜑𝑧subscript𝑧0\tau\left(r,\varphi,z,z_{0}\right)italic_τ ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT )


Theorem 3.1. Assume that conditions (3)-(11) are in place. Then the following inequality is valid:

(17) τr(r,φ,z,z0)εε2+B2=c, (r,φ,z)D¯.formulae-sequencesubscript𝜏𝑟𝑟𝜑𝑧subscript𝑧0𝜀superscript𝜀2superscript𝐵2𝑐 for-all𝑟𝜑𝑧¯𝐷\tau_{r}\left(r,\varphi,z,z_{0}\right)\geq\frac{\varepsilon}{\sqrt{\varepsilon% ^{2}+B^{2}}}=c,\text{ }\forall\left(r,\varphi,z\right)\in\overline{D}.italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≥ divide start_ARG italic_ε end_ARG start_ARG square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG = italic_c , ∀ ( italic_r , italic_φ , italic_z ) ∈ over¯ start_ARG italic_D end_ARG .

Proof. Denote τr=q1subscript𝜏𝑟subscript𝑞1\tau_{r}=q_{1}italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, τφ=q2subscript𝜏𝜑subscript𝑞2\tau_{\varphi}=q_{2}italic_τ start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, τz=q3subscript𝜏𝑧subscript𝑞3\tau_{z}=q_{3}italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT. Then the eikonal equation (12) becomes: 

(18) q12+1r2q22+q32=n2(r,φ,z).superscriptsubscript𝑞121superscript𝑟2superscriptsubscript𝑞22superscriptsubscript𝑞32superscript𝑛2𝑟𝜑𝑧q_{1}^{2}+\frac{1}{r^{2}}q_{2}^{2}+q_{3}^{2}=n^{2}(r,\varphi,z).italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z ) .

Differentiating (18) with respect to r,φ,z𝑟𝜑𝑧r,\varphi,zitalic_r , italic_φ , italic_z and using the equalities

(q2)r=(q1)φ,(q3)r=(q1)z,(q1)φ=(q2)r,(q3)φ=(q2)z,(q2)z=(q3)φ,formulae-sequencesubscriptsubscript𝑞2𝑟subscriptsubscript𝑞1𝜑formulae-sequencesubscriptsubscript𝑞3𝑟subscriptsubscript𝑞1𝑧formulae-sequencesubscriptsubscript𝑞1𝜑subscriptsubscript𝑞2𝑟formulae-sequencesubscriptsubscript𝑞3𝜑subscriptsubscript𝑞2𝑧subscriptsubscript𝑞2𝑧subscriptsubscript𝑞3𝜑\left.\begin{array}[]{c}(q_{2})_{r}=(q_{1})_{\varphi},\quad(q_{3})_{r}=(q_{1})% _{z},(q_{1})_{\varphi}=(q_{2})_{r},(q_{3})_{\varphi}=(q_{2})_{z},(q_{2})_{z}=(% q_{3})_{\varphi},\end{array}\right.start_ARRAY start_ROW start_CELL ( italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = ( italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , ( italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = ( italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , ( italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = ( italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT , ( italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = ( italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , ( italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = ( italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , end_CELL end_ROW end_ARRAY

we obtain that along the geodesic line Γ(𝐱,𝐱0)Γ𝐱subscript𝐱0\Gamma\left(\mathbf{x},\mathbf{x}_{0}\right)roman_Γ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) with 𝐱0=(0,0,z0)subscript𝐱000subscript𝑧0\mathbf{x}_{0}=(0,0,z_{0})bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( 0 , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT )

(19) drds=q1n2, dφds=q2r2n2, dzds=q3n2,formulae-sequence𝑑𝑟𝑑𝑠subscript𝑞1superscript𝑛2formulae-sequence 𝑑𝜑𝑑𝑠subscript𝑞2superscript𝑟2superscript𝑛2 𝑑𝑧𝑑𝑠subscript𝑞3superscript𝑛2\frac{dr}{ds}=\frac{q_{1}}{n^{2}},\text{ }\frac{d\varphi}{ds}=\frac{q_{2}}{r^{% 2}n^{2}},\text{ }\frac{dz}{ds}=\frac{q_{3}}{n^{2}},divide start_ARG italic_d italic_r end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , divide start_ARG italic_d italic_φ end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , divide start_ARG italic_d italic_z end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,
(20) dq1ds=nrn+q22r3n, dq2ds=nφn, dq3ds=nzn, dτds=1,formulae-sequence𝑑subscript𝑞1𝑑𝑠subscript𝑛𝑟𝑛superscriptsubscript𝑞22superscript𝑟3𝑛formulae-sequence 𝑑subscript𝑞2𝑑𝑠subscript𝑛𝜑𝑛formulae-sequence 𝑑subscript𝑞3𝑑𝑠subscript𝑛𝑧𝑛 𝑑𝜏𝑑𝑠1\frac{dq_{1}}{ds}=\frac{n_{r}}{n}+\frac{q_{2}^{2}}{r^{3}n},\text{ }\frac{dq_{2% }}{ds}=\frac{n_{\varphi}}{n},\text{ }\frac{dq_{3}}{ds}=\frac{n_{z}}{n},\text{ % }\frac{d\tau}{ds}=1,divide start_ARG italic_d italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_ARG start_ARG italic_n end_ARG + divide start_ARG italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_n end_ARG , divide start_ARG italic_d italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_n start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT end_ARG start_ARG italic_n end_ARG , divide start_ARG italic_d italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_n start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG italic_n end_ARG , divide start_ARG italic_d italic_τ end_ARG start_ARG italic_d italic_s end_ARG = 1 ,

where ds𝑑𝑠dsitalic_d italic_s is the element of the Riemannian length of Γ(𝐱,𝐱0)Γ𝐱subscript𝐱0\Gamma\left(\mathbf{x},\mathbf{x}_{0}\right)roman_Γ ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ).

Equations (19), (20) form a coupled system of ordinary differential equations. We need to figure out initial conditions for this system. Consider the geodesic line Γ(r,φ,z,𝐱0)Γ𝑟𝜑𝑧subscript𝐱0\Gamma\left(r,\varphi,z,\mathbf{x}_{0}\right)roman_Γ ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) that connects the point 𝐱0=(0,0,z0)subscript𝐱000subscript𝑧0\mathbf{x}_{0}=(0,0,z_{0})bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( 0 , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) with the point (r,φ,z).𝑟𝜑𝑧(r,\varphi,z).( italic_r , italic_φ , italic_z ) . The tangent vector to this line at the point 𝐱0subscript𝐱0\mathbf{x}_{0}bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is the vector

ν0=(sinθ0cosφ0\nu_{0}=(\sin\theta_{0}\cos\varphi_{0}italic_ν start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = ( roman_sin italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_cos italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, sinθ0sinφ0subscript𝜃0subscript𝜑0\sin\theta_{0}\sin\varphi_{0}roman_sin italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_sin italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, cosθ0)\cos\theta_{0})roman_cos italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). Hence, the initial conditions for this geodesic line are:

(21) r|s=0=0,evaluated-at𝑟𝑠00\displaystyle r|_{s=0}=0,\quaditalic_r | start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = 0 , φ|s=0=φ0,z|s=0=z0,formulae-sequenceevaluated-at𝜑𝑠0subscript𝜑0evaluated-at𝑧𝑠0subscript𝑧0\displaystyle\varphi|_{s=0}=\varphi_{0},\quad z|_{s=0}=z_{0},italic_φ | start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_z | start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ,
(22) (q1)s=0=sinθ0,subscriptsubscript𝑞1𝑠0subscript𝜃0\displaystyle(q_{1})_{s=0}=\sin\theta_{0},\quad( italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = roman_sin italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , (q2)s=0=0,(q3)s=0=cosθ0.formulae-sequencesubscriptsubscript𝑞2𝑠00subscriptsubscript𝑞3𝑠0subscript𝜃0\displaystyle(q_{2})_{s=0}=0,\quad(q_{3})_{s=0}=\cos\theta_{0}.( italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = 0 , ( italic_q start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT = roman_cos italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT .

Thus, solving Cauchy problem (19)-(22), we obtain equations of Γ((r,φ,z),(0,0,z0))Γ𝑟𝜑𝑧00subscript𝑧0\Gamma\left(\left(r,\varphi,z\right),\left(0,0,z_{0}\right)\right)roman_Γ ( ( italic_r , italic_φ , italic_z ) , ( 0 , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) in the parametric form

r=r(s,z0,θ0,φ0),φ=φ(s,z0,θ0,φ0),z=z(s,z0,θ0,φ0).formulae-sequence𝑟𝑟𝑠subscript𝑧0subscript𝜃0subscript𝜑0formulae-sequence𝜑𝜑𝑠subscript𝑧0subscript𝜃0subscript𝜑0𝑧𝑧𝑠subscript𝑧0subscript𝜃0subscript𝜑0r=r(s,z_{0},\theta_{0},\varphi_{0}),\quad\varphi=\varphi(s,z_{0},\theta_{0},% \varphi_{0}),\quad z=z(s,z_{0},\theta_{0},\varphi_{0}).italic_r = italic_r ( italic_s , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_φ = italic_φ ( italic_s , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_z = italic_z ( italic_s , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) .

Due to the regularity of geodesic lines, there exists a one-to-one correspondence between points (0,0,z0)00subscript𝑧0\left(0,0,z_{0}\right)( 0 , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) and (r,φ,z).𝑟𝜑𝑧(r,\varphi,z).( italic_r , italic_φ , italic_z ) . In particular, this means that for any fixed z0subscript𝑧0z_{0}italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT there exists a one-to-one correspondence between (r,φ,z)𝑟𝜑𝑧(r,\varphi,z)( italic_r , italic_φ , italic_z ) and parameters (s,θ0,φ0)𝑠subscript𝜃0subscript𝜑0(s,\theta_{0},\varphi_{0})( italic_s , italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ).  We denote this correspondence as:

s=s(r,φ,z,z0)τ(r,φ,z,z0),θ0=θ0(r,φ,z,z0),φ0=φ0(r,φ,z,z0).formulae-sequence𝑠𝑠𝑟𝜑𝑧subscript𝑧0𝜏𝑟𝜑𝑧subscript𝑧0formulae-sequencesubscript𝜃0subscript𝜃0𝑟𝜑𝑧subscript𝑧0subscript𝜑0subscript𝜑0𝑟𝜑𝑧subscript𝑧0s=s(r,\varphi,z,z_{0})\equiv\tau(r,\varphi,z,z_{0}),\quad\theta_{0}=\theta_{0}% (r,\varphi,z,z_{0}),\quad\varphi_{0}=\varphi_{0}(r,\varphi,z,z_{0}).italic_s = italic_s ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≡ italic_τ ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_θ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) .

Moreover, we have along the geodesic line Γ((r,φ,z),(0,0,z0))::Γ𝑟𝜑𝑧00subscript𝑧0absent\Gamma\left(\left(r,\varphi,z\right),\left(0,0,z_{0}\right)\right):roman_Γ ( ( italic_r , italic_φ , italic_z ) , ( 0 , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) :

(23) dq1ds=dτrds=nrn+q22r3n20.𝑑subscript𝑞1𝑑𝑠𝑑subscript𝜏𝑟𝑑𝑠subscript𝑛𝑟𝑛superscriptsubscript𝑞22superscript𝑟3superscript𝑛20\frac{dq_{1}}{ds}=\frac{d\tau_{r}}{ds}=\frac{n_{r}}{n}+\frac{q_{2}^{2}}{r^{3}n% ^{2}}\geq 0.divide start_ARG italic_d italic_q start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_d italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_ARG start_ARG italic_d italic_s end_ARG = divide start_ARG italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_ARG start_ARG italic_n end_ARG + divide start_ARG italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ≥ 0 .

For any point (r,φ,z)Ω¯Ωε𝑟𝜑𝑧¯ΩsubscriptΩ𝜀\left(r,\varphi,z\right)\in\overline{\Omega}\diagdown\Omega_{\varepsilon}( italic_r , italic_φ , italic_z ) ∈ over¯ start_ARG roman_Ω end_ARG ╲ roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT the geodesic line Γ(r,φ,z,𝐱0)Γ𝑟𝜑𝑧subscript𝐱0\Gamma\left(r,\varphi,z,\mathbf{x}_{0}\right)roman_Γ ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is:

Γ(r,φ,z,𝐱0)=Γε(r,φ,z,𝐱0)Γ1(r,φ,z,𝐱0),Γ𝑟𝜑𝑧subscript𝐱0subscriptΓ𝜀𝑟𝜑𝑧subscript𝐱0subscriptΓ1𝑟𝜑𝑧subscript𝐱0\Gamma\left(r,\varphi,z,\mathbf{x}_{0}\right)=\Gamma_{\varepsilon}\left(r,% \varphi,z,\mathbf{x}_{0}\right)\cup\Gamma_{1}\left(r,\varphi,z,\mathbf{x}_{0}% \right),roman_Γ ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = roman_Γ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∪ roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ,

where

(24) Γε(r,φ,z,𝐱0)=Γ(r,φ,z,𝐱0)Ωε and Γ1(r,φ,z,𝐱0)D¯.subscriptΓ𝜀𝑟𝜑𝑧subscript𝐱0Γ𝑟𝜑𝑧subscript𝐱0subscriptΩ𝜀 and subscriptΓ1𝑟𝜑𝑧subscript𝐱0¯𝐷\Gamma_{\varepsilon}\left(r,\varphi,z,\mathbf{x}_{0}\right)=\Gamma\left(r,% \varphi,z,\mathbf{x}_{0}\right)\cap\Omega_{\varepsilon}\text{ and }\Gamma_{1}\left(r,\varphi,z,\mathbf{x}_{0}\right)\subset\overline{D}.roman_Γ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = roman_Γ ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∩ roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT and roman_Γ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ⊂ over¯ start_ARG italic_D end_ARG .

By (9) Γε(r,φ,z,𝐱0)subscriptΓ𝜀𝑟𝜑𝑧subscript𝐱0\Gamma_{\varepsilon}\left(r,\varphi,z,\mathbf{x}_{0}\right)roman_Γ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is a part of a straight line. Let (ε,φ,z)𝜀superscript𝜑superscript𝑧\left(\varepsilon,\varphi^{\ast},z^{\ast}\right)( italic_ε , italic_φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ) be the intersection of Γε(𝐱,𝐱0)subscriptΓ𝜀𝐱subscript𝐱0\Gamma_{\varepsilon}(\mathbf{x,x}_{0})roman_Γ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT ( bold_x , bold_x start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) with the lateral boundary {r=ε,z(0,B)}formulae-sequence𝑟𝜀𝑧0𝐵\left\{r=\varepsilon,z\in\left(0,B\right)\right\}{ italic_r = italic_ε , italic_z ∈ ( 0 , italic_B ) } of the cylinder Ωε.subscriptΩ𝜀\Omega_{\varepsilon}.roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT . Then τ(r,φ,z,z0)=r2+(zz0)2𝜏𝑟superscript𝜑superscript𝑧subscript𝑧0superscript𝑟2superscriptsuperscript𝑧subscript𝑧02\tau\left(r,\varphi^{\ast},z^{\ast},z_{0}\right)=\sqrt{r^{2}+\left(z^{\ast}-z_% {0}\right)^{2}}italic_τ ( italic_r , italic_φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT - italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG for r(0,ε).𝑟0𝜀r\in\left(0,\varepsilon\right).italic_r ∈ ( 0 , italic_ε ) . Hence,

(25) τr(ε,φ,z,z0)=εε2+(zz0)2εε2+B2.subscript𝜏𝑟𝜀superscript𝜑superscript𝑧subscript𝑧0𝜀superscript𝜀2superscriptsuperscript𝑧subscript𝑧02𝜀superscript𝜀2superscript𝐵2\tau_{r}\left(\varepsilon,\varphi^{\ast},z^{\ast},z_{0}\right)=\frac{% \varepsilon}{\sqrt{\varepsilon^{2}+\left(z^{\ast}-z_{0}\right)^{2}}}\geq\frac{% \varepsilon}{\sqrt{\varepsilon^{2}+B^{2}}}.italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_ε , italic_φ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = divide start_ARG italic_ε end_ARG start_ARG square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT - italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ≥ divide start_ARG italic_ε end_ARG start_ARG square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG .

Hence, by (11), (20) and (25) for (r,φ,z)Ω¯Ωε𝑟𝜑𝑧¯ΩsubscriptΩ𝜀\left(r,\varphi,z\right)\in\overline{\Omega}\diagdown\Omega_{\varepsilon}( italic_r , italic_φ , italic_z ) ∈ over¯ start_ARG roman_Ω end_ARG ╲ roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT

τr(r,φ,z,z0)=ε2+(zz0)2ssubscript𝜏𝑟𝑟𝜑𝑧subscript𝑧0superscriptsubscriptsuperscript𝜀2superscriptsuperscript𝑧subscript𝑧02𝑠\displaystyle\tau_{r}\left(r,\varphi,z,z_{0}\right)=\int\limits_{\sqrt{% \varepsilon^{2}+\left(z^{\ast}-z_{0}\right)^{2}}}^{s}italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∫ start_POSTSUBSCRIPT square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT - italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_s end_POSTSUPERSCRIPT (nrn+q22r3n2)ds+εε2+(zz0)2subscript𝑛𝑟𝑛superscriptsubscript𝑞22superscript𝑟3superscript𝑛2𝑑𝑠𝜀superscript𝜀2superscriptsuperscript𝑧subscript𝑧02absent\displaystyle\left(\frac{n_{r}}{n}+\frac{q_{2}^{2}}{r^{3}n^{2}}\right)ds+\frac% {\varepsilon}{\sqrt{\varepsilon^{2}+\left(z^{\ast}-z_{0}\right)^{2}}}\geq( divide start_ARG italic_n start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT end_ARG start_ARG italic_n end_ARG + divide start_ARG italic_q start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ) italic_d italic_s + divide start_ARG italic_ε end_ARG start_ARG square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_z start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT - italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG ≥
εε2+B2=c. formulae-sequenceabsent𝜀superscript𝜀2superscript𝐵2𝑐 \displaystyle\geq\frac{\varepsilon}{\sqrt{\varepsilon^{2}+B^{2}}}=c.\text{ \ }\square≥ divide start_ARG italic_ε end_ARG start_ARG square-root start_ARG italic_ε start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_B start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG = italic_c . □

4 Integral Differential Equation

Denote

(26) u(r,φ,z,z0)=τr2(r,φ,z,z0),(r,φ,z)D,z0(0,B).formulae-sequence𝑢𝑟𝜑𝑧subscript𝑧0superscriptsubscript𝜏𝑟2𝑟𝜑𝑧subscript𝑧0formulae-sequence𝑟𝜑𝑧𝐷subscript𝑧00𝐵\left.\begin{array}[]{c}u\left(r,\varphi,z,z_{0}\right)=\tau_{r}^{2}\left(r,% \varphi,z,z_{0}\right),\ \left(r,\varphi,z\right)\in D,z_{0}\in\left(0,B\right% ).\end{array}\right.start_ARRAY start_ROW start_CELL italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ( italic_r , italic_φ , italic_z ) ∈ italic_D , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY

Since by Theorem 3.1 τr(r,φ,z,z0)c>0,subscript𝜏𝑟𝑟𝜑𝑧subscript𝑧0𝑐0\tau_{r}\left(r,\varphi,z,z_{0}\right)\geq c>0,italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≥ italic_c > 0 , then  by (26)

(27) τr(r,φ,z,z0)=u(r,φ,z,z0).subscript𝜏𝑟𝑟𝜑𝑧subscript𝑧0𝑢𝑟𝜑𝑧subscript𝑧0\tau_{r}\left(r,\varphi,z,z_{0}\right)=\sqrt{u\left(r,\varphi,z,z_{0}\right)}.italic_τ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = square-root start_ARG italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG .

Hence, (14) and (27) imply

(28) τ(r,φ,z,z0)=rRu(t,φ,z,z0)𝑑t+p(φ,z,z0),τ(r,φ,0,z0)=p0(r,φ,z0),τ(r,φ,B,z0)=pB(r,φ,z0),(r,φ,z,z0)Q=(ε,R)×(0,2π)×(0,B)×(0,B).𝜏𝑟𝜑𝑧subscript𝑧0superscriptsubscript𝑟𝑅𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡𝑝𝜑𝑧subscript𝑧0formulae-sequence𝜏𝑟𝜑0subscript𝑧0subscript𝑝0𝑟𝜑subscript𝑧0𝜏𝑟𝜑𝐵subscript𝑧0subscript𝑝𝐵𝑟𝜑subscript𝑧0𝑟𝜑𝑧subscript𝑧0𝑄𝜀𝑅02𝜋0𝐵0𝐵\left.\begin{array}[]{c}\tau\left(r,\varphi,z,z_{0}\right)=-\int\limits_{r}^{R% }\sqrt{u\left(t,\varphi,z,z_{0}\right)}dt+p\left(\varphi,z,z_{0}\right),\\ \tau\left(r,\varphi,0,z_{0}\right)=p_{0}\left(r,\varphi,z_{0}\right),\ \tau% \left(r,\varphi,B,z_{0}\right)=p_{B}\left(r,\varphi,z_{0}\right),\\ \left(r,\varphi,z,z_{0}\right)\in Q=\left(\varepsilon,R\right)\times\left(0,2% \pi\right)\times\left(0,B\right)\times\left(0,B\right).\end{array}\right.start_ARRAY start_ROW start_CELL italic_τ ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = - ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG italic_d italic_t + italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL italic_τ ( italic_r , italic_φ , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_τ ( italic_r , italic_φ , italic_B , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ italic_Q = ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) × ( 0 , italic_B ) × ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY

By the first line of (28)

(29) τφ(r,φ,z,z0)subscript𝜏𝜑𝑟𝜑𝑧subscript𝑧0\displaystyle\tau_{\varphi}\left(r,\varphi,z,z_{0}\right)italic_τ start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =12rRuφ(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pφ(φ,z,z0),absent12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝜑𝜑𝑧subscript𝑧0\displaystyle=-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{\varphi}\left(t,\varphi,% z,z_{0}\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_{\varphi}\left(% \varphi,z,z_{0}\right),= - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ,
(30) τz(r,φ,z,z0)subscript𝜏𝑧𝑟𝜑𝑧subscript𝑧0\displaystyle\tau_{z}\left(r,\varphi,z,z_{0}\right)italic_τ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =12rRuz(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pz(φ,z,z0).absent12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝑧𝜑𝑧subscript𝑧0\displaystyle=-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{z}\left(t,\varphi,z,z_{0% }\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_{z}\left(\varphi,z,z_{0}% \right).= - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) .

Substituting (28)-(30) in the eikonal equation (12), we obtain the following equation in the domain Q𝑄Qitalic_Q defined in the last line of (28):

u(r,φ,z,z0)+1r2[12rRuφ(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pφ(φ,z,z0)]2+𝑢𝑟𝜑𝑧subscript𝑧0limit-from1superscript𝑟2superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝜑𝜑𝑧subscript𝑧02u\left(r,\varphi,z,z_{0}\right)+\frac{1}{r^{2}}\left[-\frac{1}{2}\int\limits_{% r}^{R}\frac{u_{\varphi}\left(t,\varphi,z,z_{0}\right)}{\sqrt{u\left(t,\varphi,% z,z_{0}\right)}}dt+p_{\varphi}\left(\varphi,z,z_{0}\right)\right]^{2}+italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) + divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT +
(31) +[12rRuz(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pz(φ,z,z0)]2=n2(r,φ,z).superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝑧𝜑𝑧subscript𝑧02superscript𝑛2𝑟𝜑𝑧+\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{z}\left(t,\varphi,z,z_{0}\right% )}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_{z}\left(\varphi,z,z_{0}\right)% \right]^{2}=n^{2}(r,\varphi,z).+ [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z ) .

Differentiating (31) with respect to z0subscript𝑧0z_{0}italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and using z0[n2(r,φ,z)]0,subscriptsubscript𝑧0delimited-[]superscript𝑛2𝑟𝜑𝑧0\partial_{z_{0}}\left[n^{2}(r,\varphi,z)\right]\equiv 0,∂ start_POSTSUBSCRIPT italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT [ italic_n start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z ) ] ≡ 0 , we obtain

uz0(r,φ,z,z0)+1r2z0[12rRuφ(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pφ(φ,z,z0)]2+subscript𝑢subscript𝑧0𝑟𝜑𝑧subscript𝑧0limit-from1superscript𝑟2subscript𝑧0superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝜑𝜑𝑧subscript𝑧02u_{z_{0}}\left(r,\varphi,z,z_{0}\right)+\frac{1}{r^{2}}\frac{\partial}{% \partial z_{0}}\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{\varphi}\left(t,% \varphi,z,z_{0}\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_{\varphi}% \left(\varphi,z,z_{0}\right)\right]^{2}+italic_u start_POSTSUBSCRIPT italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) + divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG divide start_ARG ∂ end_ARG start_ARG ∂ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT +
(32) +z0[12rRuz(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pz(φ,z,z0)]2=0.subscript𝑧0superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝑧𝜑𝑧subscript𝑧020+\frac{\partial}{\partial z_{0}}\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{% z}\left(t,\varphi,z,z_{0}\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_% {z}\left(\varphi,z,z_{0}\right)\right]^{2}=0.+ divide start_ARG ∂ end_ARG start_ARG ∂ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 0 .

In addition, (15), (16) and (26) imply:

(33) u(r,φ,0,z0)=(rp0(r,φ,z0))2, u(r,φ,B,z0)=(rpB(r,φ,z0))2.formulae-sequence𝑢𝑟𝜑0subscript𝑧0superscriptsubscript𝑟subscript𝑝0𝑟𝜑subscript𝑧02 𝑢𝑟𝜑𝐵subscript𝑧0superscriptsubscript𝑟subscript𝑝𝐵𝑟𝜑subscript𝑧02u\left(r,\varphi,0,z_{0}\right)=\left(\partial_{r}p_{0}\left(r,\varphi,z_{0}% \right)\right)^{2},\text{ }u\left(r,\varphi,B,z_{0}\right)=\left(\partial_{r}p% _{B}\left(r,\varphi,z_{0}\right)\right)^{2}.italic_u ( italic_r , italic_φ , 0 , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ( ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_u ( italic_r , italic_φ , italic_B , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ( ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .

In (31)-(33)

(34) (r,φ,z,z0)Q¯.𝑟𝜑𝑧subscript𝑧0¯𝑄\left(r,\varphi,z,z_{0}\right)\in\overline{Q}.( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ over¯ start_ARG italic_Q end_ARG .

Thus, we came up with the following boundary value problem (BVP).

Boundary Value Problem 1 (BVP1). Find the function u(r,φ,z,z0)C1(Q¯)𝑢𝑟𝜑𝑧subscript𝑧0superscript𝐶1¯𝑄u\left(r,\varphi,z,z_{0}\right)\in C^{1}\left(\overline{Q}\right)italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( over¯ start_ARG italic_Q end_ARG )  satisfying conditions (32), (33) in the domain (34).

Suppose that we have solved this problem. Then we substitute its solution in the left hand side of equation (31) and find the target function n(r,φ,z).𝑛𝑟𝜑𝑧n(r,\varphi,z).italic_n ( italic_r , italic_φ , italic_z ) . Clearly BVP1 is a very complicated one. Therefore, we construct below an approximate mathematical model for its solution, and then confirm the validity of this model computationally.

5 A Special Orthonormal Basis in L2(0,B)subscript𝐿20𝐵L_{2}\left(0,B\right)italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 , italic_B )

The basis described in this section was first introduced in [11]. Then it was used for the convexification method in some other works, see, e.g. [12, 14, 15]. Consider the set of functions {z0nez0}n=0L2(0,B).superscriptsubscriptsuperscriptsubscript𝑧0𝑛superscript𝑒subscript𝑧0𝑛0subscript𝐿20𝐵\{z_{0}^{n}e^{z_{0}}\}_{n=0}^{\infty}\subset L_{2}(0,B).{ italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT } start_POSTSUBSCRIPT italic_n = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ⊂ italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 , italic_B ) . Obviously, this set is complete in L2(0,B),subscript𝐿20𝐵L_{2}(0,B),italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 , italic_B ) , and these functions are linearly independent. Apply the Gram-Schmidt orthonormalization procedure to this set. Then we obtain the orthonormal basis {Ψn(z0)}n=0superscriptsubscriptsubscriptΨ𝑛subscript𝑧0𝑛0\{\Psi_{n}\left(z_{0}\right)\}_{n=0}^{\infty}{ roman_Ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) } start_POSTSUBSCRIPT italic_n = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT in L2(0,B)subscript𝐿20𝐵L_{2}(0,B)italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( 0 , italic_B ). The function Ψn(z0)subscriptΨ𝑛subscript𝑧0\Psi_{n}(z_{0})roman_Ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) has the form Ψn(z0)=Pn(z0)ez0subscriptΨ𝑛subscript𝑧0subscript𝑃𝑛subscript𝑧0superscript𝑒subscript𝑧0\Psi_{n}(z_{0})=P_{n}(z_{0})e^{z_{0}}roman_Ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_P start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT, n0,for-all𝑛0\forall n\geq 0,∀ italic_n ≥ 0 , where Pn(z0)subscript𝑃𝑛subscript𝑧0P_{n}(z_{0})italic_P start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is a polynomial of the degree n𝑛nitalic_n. Denote

amn=0BΦm(z0),Φn(z0)dz0.subscript𝑎𝑚𝑛superscriptsubscript0𝐵subscriptΦ𝑚subscript𝑧0superscriptsubscriptΦ𝑛subscript𝑧0𝑑subscript𝑧0a_{mn}=\int\limits_{0}^{B}\Phi_{m}\left(z_{0}\right),\Phi_{n}^{\prime}\left(z_% {0}\right)dz_{0}.italic_a start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT roman_Φ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , roman_Φ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) italic_d italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT .

It was proven in [11], [14, Theorem 6.2.1] that

(35) amn={1 if m=n,0 if m>n.subscript𝑎𝑚𝑛cases1 if 𝑚𝑛0 if 𝑚𝑛a_{mn}=\left\{\begin{array}[]{c}1\text{ if }m=n,\\ 0\text{ if }m>n.\end{array}\right.italic_a start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT = { start_ARRAY start_ROW start_CELL 1 if italic_m = italic_n , end_CELL end_ROW start_ROW start_CELL 0 if italic_m > italic_n . end_CELL end_ROW end_ARRAY

Let N1𝑁1N\geq 1italic_N ≥ 1 be an integer. Consider the N×N𝑁𝑁N\times Nitalic_N × italic_N matrix AN=(amn)m,n=0N.subscript𝐴𝑁superscriptsubscriptsubscript𝑎𝑚𝑛𝑚𝑛0𝑁A_{N}=\left(a_{mn}\right)_{m,n=0}^{N}.italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = ( italic_a start_POSTSUBSCRIPT italic_m italic_n end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT italic_m , italic_n = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT . By (35) detAN=1.subscript𝐴𝑁1\det A_{N}=1.roman_det italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = 1 . Hence, the matrix ANsubscript𝐴𝑁A_{N}italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT is invertible. We observe that neither classical orthonormal polynomials nor the basis of trigonometric functions do not provide a corresponding invertible matrix ANsubscript𝐴𝑁A_{N}italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT. This is because the first function in these two cases is an identical constant, implying that the first column of that analog of ANsubscript𝐴𝑁A_{N}italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT is formed only by zeros.

6 Boundary Value Problem 2

The first step of our approximate mathematical model mentioned in section 4 is the assumption that the function u𝑢uitalic_u can be represented as truncated Fourier-like series,

(36) u(r,φ,z,z0)=s=0N1us(r,φ,z)Φs(z0), (r,φ,z,z0)Q.formulae-sequence𝑢𝑟𝜑𝑧subscript𝑧0superscriptsubscript𝑠0𝑁1subscript𝑢𝑠𝑟𝜑𝑧subscriptΦ𝑠subscript𝑧0 𝑟𝜑𝑧subscript𝑧0𝑄u\left(r,\varphi,z,z_{0}\right)=\sum\limits_{s=0}^{N-1}u_{s}\left(r,\varphi,z% \right)\Phi_{s}\left(z_{0}\right),\text{ }\left(r,\varphi,z,z_{0}\right)\in Q.italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z ) roman_Φ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ italic_Q .

We do not prove convergence of our method as N.𝑁N\rightarrow\infty.italic_N → ∞ . We point out that the absence of such proofs is a rather common place in the theory of Inverse Problems, basically due to their ill-posedness, see, e.g. [6, 8].

Consider the Nlimit-from𝑁N-italic_N -D vector function V(r,φ,z)𝑉𝑟𝜑𝑧V\left(r,\varphi,z\right)italic_V ( italic_r , italic_φ , italic_z ) of coefficients of the series (36)

(37) V(r,φ,z)=(u0,,uN1)T(r,φ,z), (r,φ,z)(ε,R)×(0,2π)×(0,B).𝑉𝑟𝜑𝑧superscriptsubscript𝑢0subscript𝑢𝑁1𝑇𝑟𝜑𝑧 𝑟𝜑𝑧𝜀𝑅02𝜋0𝐵\left.\begin{array}[]{c}V\left(r,\varphi,z\right)=\left(u_{0},...,u_{N-1}% \right)^{T}\left(r,\varphi,z\right),\text{ }\\ \left(r,\varphi,z\right)\in\left(\varepsilon,R\right)\times\left(0,2\pi\right)% \times\left(0,B\right).\end{array}\right.start_ARRAY start_ROW start_CELL italic_V ( italic_r , italic_φ , italic_z ) = ( italic_u start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , italic_u start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z ) , end_CELL end_ROW start_ROW start_CELL ( italic_r , italic_φ , italic_z ) ∈ ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) × ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY

We want to find the vector function V(r,φ,z)𝑉𝑟𝜑𝑧V\left(r,\varphi,z\right)italic_V ( italic_r , italic_φ , italic_z ) in (37). Naturally, we assume that for functions in (14)-(16)

(38) p(φ,z,z0)=s=0N1ps(φ,z)Φs(z0),(φ,z,z0)(0,2π)×(0,B)×(0,B).𝑝𝜑𝑧subscript𝑧0superscriptsubscript𝑠0𝑁1subscript𝑝𝑠𝜑𝑧subscriptΦ𝑠subscript𝑧0𝜑𝑧subscript𝑧002𝜋0𝐵0𝐵\left.\begin{array}[]{c}p\left(\varphi,z,z_{0}\right)=\sum\limits_{s=0}^{N-1}p% _{s}\left(\varphi,z\right)\Phi_{s}\left(z_{0}\right),\\ \left(\varphi,z,z_{0}\right)\in\left(0,2\pi\right)\times\left(0,B\right)\times% \left(0,B\right).\end{array}\right.start_ARRAY start_ROW start_CELL italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_φ , italic_z ) roman_Φ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ ( 0 , 2 italic_π ) × ( 0 , italic_B ) × ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY
(39) (rp0)2(r,φ,z0)=s=0N1p~0,s(r,φ)Φs(z0),(rpB)2(r,φ,z0)=s=0N1p~B,s(r,φ)Φs(z0),(r,φ,z0)(ε,R)×(0,2π)×(0,B).superscriptsubscript𝑟subscript𝑝02𝑟𝜑subscript𝑧0superscriptsubscript𝑠0𝑁1subscript~𝑝0𝑠𝑟𝜑subscriptΦ𝑠subscript𝑧0superscriptsubscript𝑟subscript𝑝𝐵2𝑟𝜑subscript𝑧0superscriptsubscript𝑠0𝑁1subscript~𝑝𝐵𝑠𝑟𝜑subscriptΦ𝑠subscript𝑧0𝑟𝜑subscript𝑧0𝜀𝑅02𝜋0𝐵\left.\begin{array}[]{c}\left(\partial_{r}p_{0}\right)^{2}\left(r,\varphi,z_{0% }\right)=\sum\limits_{s=0}^{N-1}\widetilde{p}_{0,s}\left(r,\varphi\right)\Phi_% {s}\left(z_{0}\right),\\ \left(\partial_{r}p_{B}\right)^{2}\left(r,\varphi,z_{0}\right)=\sum\limits_{s=% 0}^{N-1}\widetilde{p}_{B,s}\left(r,\varphi\right)\Phi_{s}\left(z_{0}\right),\\ \left(r,\varphi,z_{0}\right)\in\left(\varepsilon,R\right)\times\left(0,2\pi% \right)\times\left(0,B\right).\end{array}\right.start_ARRAY start_ROW start_CELL ( ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 0 , italic_s end_POSTSUBSCRIPT ( italic_r , italic_φ ) roman_Φ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL ( ∂ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_B , italic_s end_POSTSUBSCRIPT ( italic_r , italic_φ ) roman_Φ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) × ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY

Denote

(40) G(φ,z)=(p0,,pN1)T(φ,z),(φ,z)(0,2π)×(0,B), G0(r,φ)=(p~0,0,p~0,N1)T(r,φ), (r,φ)(ε,R)×(0,2π)GB(r,φ)=(p~B,0,p~B,N1)T(r,φ), (r,φ)(ε,R)×(0,2π).formulae-sequence𝐺𝜑𝑧superscriptsubscript𝑝0subscript𝑝𝑁1𝑇𝜑𝑧𝜑𝑧02𝜋0𝐵formulae-sequence subscript𝐺0𝑟𝜑superscriptsubscript~𝑝00subscript~𝑝0𝑁1𝑇𝑟𝜑 𝑟𝜑𝜀𝑅02𝜋formulae-sequencesubscript𝐺𝐵𝑟𝜑superscriptsubscript~𝑝𝐵0subscript~𝑝𝐵𝑁1𝑇𝑟𝜑 𝑟𝜑𝜀𝑅02𝜋\left.\begin{array}[]{c}G\left(\varphi,z\right)=\left(p_{0},...,p_{N-1}\right)% ^{T}\left(\varphi,z\right),\quad\left(\varphi,z\right)\in\left(0,2\pi\right)% \times\left(0,B\right),\\ \text{ }G_{0}\left(r,\varphi\right)=\left(\widetilde{p}_{0,0},...\widetilde{p}% _{0,N-1}\right)^{T}\left(r,\varphi\right),\text{ }\left(r,\varphi\right)\in% \left(\varepsilon,R\right)\times\left(0,2\pi\right)\\ G_{B}\left(r,\varphi\right)=\left(\widetilde{p}_{B,0},...\widetilde{p}_{B,N-1}% \right)^{T}\left(r,\varphi\right),\text{ }\left(r,\varphi\right)\in\left(% \varepsilon,R\right)\times\left(0,2\pi\right).\end{array}\right.start_ARRAY start_ROW start_CELL italic_G ( italic_φ , italic_z ) = ( italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , italic_p start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_φ , italic_z ) , ( italic_φ , italic_z ) ∈ ( 0 , 2 italic_π ) × ( 0 , italic_B ) , end_CELL end_ROW start_ROW start_CELL italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ ) = ( over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 0 , 0 end_POSTSUBSCRIPT , … over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT 0 , italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_r , italic_φ ) , ( italic_r , italic_φ ) ∈ ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) end_CELL end_ROW start_ROW start_CELL italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ ) = ( over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_B , 0 end_POSTSUBSCRIPT , … over~ start_ARG italic_p end_ARG start_POSTSUBSCRIPT italic_B , italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_r , italic_φ ) , ( italic_r , italic_φ ) ∈ ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) . end_CELL end_ROW end_ARRAY

Substituting (36)-(40) in (32), multiplying the obtained equality sequentially by the functions Ψn(z0),n=0,,N1formulae-sequencesubscriptΨ𝑛subscript𝑧0𝑛0𝑁1\Psi_{n}\left(z_{0}\right),n=0,...,N-1roman_Ψ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_n = 0 , … , italic_N - 1 and integrating with respect to z0(0,B),subscript𝑧00𝐵z_{0}\in\left(0,B\right),italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , italic_B ) , we obtain

(41) ANV+S(Vφ,Vz,G,r,φ,z)=0,V(r,φ,0)= G0(r,φ),V(r,φ,B)= GB(r,φ),V(r,0,z)=V(r,2π,z),(r,φ,z)(ε,R)×(0,2π)×(0,B),subscript𝐴𝑁𝑉𝑆subscript𝑉𝜑subscript𝑉𝑧𝐺𝑟𝜑𝑧0𝑉𝑟𝜑0 subscript𝐺0𝑟𝜑𝑉𝑟𝜑𝐵 subscript𝐺𝐵𝑟𝜑𝑉𝑟0𝑧𝑉𝑟2𝜋𝑧𝑟𝜑𝑧𝜀𝑅02𝜋0𝐵\left.\begin{array}[]{c}A_{N}V+S\left(V_{\varphi},V_{z},G,r,\varphi,z\right)=0% ,\\ V\left(r,\varphi,0\right)=\text{ }G_{0}\left(r,\varphi\right),\\ V\left(r,\varphi,B\right)=\text{ }G_{B}\left(r,\varphi\right),\\ V\left(r,0,z\right)=V\left(r,2\pi,z\right),\\ \left(r,\varphi,z\right)\in\left(\varepsilon,R\right)\times\left(0,2\pi\right)% \times\left(0,B\right),\end{array}\right.start_ARRAY start_ROW start_CELL italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_V + italic_S ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_G , italic_r , italic_φ , italic_z ) = 0 , end_CELL end_ROW start_ROW start_CELL italic_V ( italic_r , italic_φ , 0 ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ ) , end_CELL end_ROW start_ROW start_CELL italic_V ( italic_r , italic_φ , italic_B ) = italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ ) , end_CELL end_ROW start_ROW start_CELL italic_V ( italic_r , 0 , italic_z ) = italic_V ( italic_r , 2 italic_π , italic_z ) , end_CELL end_ROW start_ROW start_CELL ( italic_r , italic_φ , italic_z ) ∈ ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) × ( 0 , italic_B ) , end_CELL end_ROW end_ARRAY

where SN𝑆superscript𝑁S\in\mathbb{R}^{N}italic_S ∈ blackboard_R start_POSTSUPERSCRIPT italic_N end_POSTSUPERSCRIPT is a C1limit-fromsuperscript𝐶1C^{1}-italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT -vector function of its 3N+33𝑁33N+33 italic_N + 3 arguments. The fourth line in (41) means the 2πlimit-from2𝜋2\pi-2 italic_π -periodicity condition with respect to the angle φ.𝜑\varphi.italic_φ . Thus, (41) is the boundary value problem for the nonlinear system of integral differential equations with respect to the vector function V(r,φ,z).𝑉𝑟𝜑𝑧V\left(r,\varphi,z\right).italic_V ( italic_r , italic_φ , italic_z ) . Here, the slimit-from𝑠s-italic_s -th component of the Nlimit-from𝑁N-italic_N -D vector function S=(S0,,SN1)T𝑆superscriptsubscript𝑆0subscript𝑆𝑁1𝑇S=\left(S_{0},...,S_{N-1}\right)^{T}italic_S = ( italic_S start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , italic_S start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT has the form

Ss(Vφ,Vz,r)=Ssφ(Vφ,Vz,r)+Ssz(Vφ,Vz,r).subscript𝑆𝑠subscript𝑉𝜑subscript𝑉𝑧𝑟superscriptsubscript𝑆𝑠𝜑subscript𝑉𝜑subscript𝑉𝑧𝑟superscriptsubscript𝑆𝑠𝑧subscript𝑉𝜑subscript𝑉𝑧𝑟S_{s}\left(V_{\varphi},V_{z},r\right)=S_{s}^{\varphi}\left(V_{\varphi},V_{z},r% \right)+S_{s}^{z}\left(V_{\varphi},V_{z},r\right).italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) = italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_φ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) + italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) .

Here Snφ(Vφ,Vz,Gφ,Gz,r)superscriptsubscript𝑆𝑛𝜑subscript𝑉𝜑subscript𝑉𝑧subscript𝐺𝜑subscript𝐺𝑧𝑟S_{n}^{\varphi}\left(V_{\varphi},V_{z},G_{\varphi},G_{z},r\right)italic_S start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_φ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) has the following form after the integration by parts with respect to z0subscript𝑧0z_{0}italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT:

Ssφ(Vφ,Vz,r)=superscriptsubscript𝑆𝑠𝜑subscript𝑉𝜑subscript𝑉𝑧𝑟absent\displaystyle S_{s}^{\varphi}\left(V_{\varphi},V_{z},r\right)=italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_φ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) = Ψs(B)1r2[12rRuφ(t,φ,z,B)u(t,φ,z,B)𝑑t+pφ(φ,z,B)]2limit-fromsubscriptΨ𝑠𝐵1superscript𝑟2superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧𝐵𝑢𝑡𝜑𝑧𝐵differential-d𝑡subscript𝑝𝜑𝜑𝑧𝐵2\displaystyle\Psi_{s}\left(B\right)\frac{1}{r^{2}}\left[-\frac{1}{2}\int% \limits_{r}^{R}\frac{u_{\varphi}\left(t,\varphi,z,B\right)}{\sqrt{u\left(t,% \varphi,z,B\right)}}dt+p_{\varphi}\left(\varphi,z,B\right)\right]^{2}-roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_B ) divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_B ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_B ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_B ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT -
(42) Ψs(0)1r2subscriptΨ𝑠01superscript𝑟2\displaystyle-\Psi_{s}\left(0\right)\frac{1}{r^{2}}- roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( 0 ) divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG [12rRuφ(t,φ,z,0)u(t,φ,z,0)𝑑t+pφ(φ,z,0)]2limit-fromsuperscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧0𝑢𝑡𝜑𝑧0differential-d𝑡subscript𝑝𝜑𝜑𝑧02\displaystyle\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{\varphi}\left(t,% \varphi,z,0\right)}{\sqrt{u\left(t,\varphi,z,0\right)}}dt+p_{\varphi}\left(% \varphi,z,0\right)\right]^{2}-[ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , 0 ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , 0 ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , 0 ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT -
1r20BΨs1superscript𝑟2superscriptsubscript0𝐵superscriptsubscriptΨ𝑠\displaystyle-\frac{1}{r^{2}}\int\limits_{0}^{B}\Psi_{s}^{\prime}- divide start_ARG 1 end_ARG start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT (z0)[12rRuφ(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pφ(φ,z,z0)]2dz0.subscript𝑧0superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝜑𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝜑𝜑𝑧subscript𝑧02𝑑subscript𝑧0\displaystyle\left(z_{0}\right)\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{% \varphi}\left(t,\varphi,z,z_{0}\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}% }dt+p_{\varphi}\left(\varphi,z,z_{0}\right)\right]^{2}dz_{0}.( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT .

And Ssz(Vφ,Vz,Gφ,Gz,r)superscriptsubscript𝑆𝑠𝑧subscript𝑉𝜑subscript𝑉𝑧subscript𝐺𝜑subscript𝐺𝑧𝑟S_{s}^{z}\left(V_{\varphi},V_{z},G_{\varphi},G_{z},r\right)italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_G start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) has the form after the integration by parts with respect to z0subscript𝑧0z_{0}italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT:

(43) Ssz(Vφ,Vz,r)=superscriptsubscript𝑆𝑠𝑧subscript𝑉𝜑subscript𝑉𝑧𝑟absent\displaystyle S_{s}^{z}\left(V_{\varphi},V_{z},r\right)=italic_S start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_z end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_r ) = Ψs(B)[12rRuz(t,φ,z,B)u(t,φ,z,B)𝑑t+pz(φ,z,B)]2limit-fromsubscriptΨ𝑠𝐵superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧𝐵𝑢𝑡𝜑𝑧𝐵differential-d𝑡subscript𝑝𝑧𝜑𝑧𝐵2\displaystyle\Psi_{s}\left(B\right)\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{% u_{z}\left(t,\varphi,z,B\right)}{\sqrt{u\left(t,\varphi,z,B\right)}}dt+p_{z}% \left(\varphi,z,B\right)\right]^{2}-roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_B ) [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_B ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_B ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_B ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT -
Ψs(0)subscriptΨ𝑠0\displaystyle-\Psi_{s}\left(0\right)- roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( 0 ) [12rRuz(t,φ,z,0)u(t,φ,z,0)𝑑t+pz(φ,z,0)]2limit-fromsuperscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧0𝑢𝑡𝜑𝑧0differential-d𝑡subscript𝑝𝑧𝜑𝑧02\displaystyle\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{z}\left(t,\varphi,z% ,0\right)}{\sqrt{u\left(t,\varphi,z,0\right)}}dt+p_{z}\left(\varphi,z,0\right)% \right]^{2}-[ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , 0 ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , 0 ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , 0 ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT -
0BΨssuperscriptsubscript0𝐵superscriptsubscriptΨ𝑠\displaystyle-\int\limits_{0}^{B}\Psi_{s}^{\prime}- ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT roman_Ψ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT (z0)[12rRuz(t,φ,z,z0)u(t,φ,z,z0)𝑑t+pz(φ,z,z0)]2dz0.subscript𝑧0superscriptdelimited-[]12superscriptsubscript𝑟𝑅subscript𝑢𝑧𝑡𝜑𝑧subscript𝑧0𝑢𝑡𝜑𝑧subscript𝑧0differential-d𝑡subscript𝑝𝑧𝜑𝑧subscript𝑧02𝑑subscript𝑧0\displaystyle\left(z_{0}\right)\left[-\frac{1}{2}\int\limits_{r}^{R}\frac{u_{z% }\left(t,\varphi,z,z_{0}\right)}{\sqrt{u\left(t,\varphi,z,z_{0}\right)}}dt+p_{% z}\left(\varphi,z,z_{0}\right)\right]^{2}dz_{0}.( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) [ - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT divide start_ARG italic_u start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG start_ARG square-root start_ARG italic_u ( italic_t , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG end_ARG italic_d italic_t + italic_p start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT .

In (42), (43) the function u𝑢uitalic_u has the form (36), and the vector function V𝑉Vitalic_V has the form (37). Thus, we have obtained Boundary Value Problem 2 (BVP2).

Boundary Value Problem 2 (BVP2). Let D𝐷Ditalic_D be the domain defined in (8). Find the Nlimit-from𝑁N-italic_N -D vector function VC1(D¯)𝑉superscript𝐶1¯𝐷V\in C^{1}\left(\overline{D}\right)italic_V ∈ italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( over¯ start_ARG italic_D end_ARG ) satisfying conditions (41), where the function u𝑢uitalic_u is connected with V𝑉Vitalic_V via (36), (37).

We now formulate a Carleman estimate for the Volterra integral operator occurring in (41)-(43). First, we introduce the Carleman Weight Function φλ(r)subscript𝜑𝜆𝑟\varphi_{\lambda}\left(r\right)italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) for this Volterra operator, where λ>0𝜆0\lambda>0italic_λ > 0 is a parameter. This function is:

(44) φλ(r)=e2λr.subscript𝜑𝜆𝑟superscript𝑒2𝜆𝑟\varphi_{\lambda}\left(r\right)=e^{2\lambda r}.italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) = italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT .

Theorem 6.1 (Carleman estimate for the Volterra integral operator). The following Carleman estimate is valid:

(45) εR(rRf(s)𝑑s)2e2λr𝑑r1λ2εRf2(r)e2λr𝑑r, λ>0, fL2(ε,R).formulae-sequencesuperscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2superscript𝑒2𝜆𝑟differential-d𝑟1superscript𝜆2superscriptsubscript𝜀𝑅superscript𝑓2𝑟superscript𝑒2𝜆𝑟differential-d𝑟formulae-sequence for-all𝜆0 for-all𝑓subscript𝐿2𝜀𝑅\int\limits_{\varepsilon}^{R}\left(\int\limits_{r}^{R}f\left(s\right)ds\right)% ^{2}e^{2\lambda r}dr\leq\frac{1}{\lambda^{2}}\int\limits_{\varepsilon}^{R}f^{2% }\left(r\right)e^{2\lambda r}dr,\text{ }\forall\lambda>0,\text{ }\forall f\in L% _{2}\left(\varepsilon,R\right).∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT italic_d italic_r ≤ divide start_ARG 1 end_ARG start_ARG italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r ) italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT italic_d italic_r , ∀ italic_λ > 0 , ∀ italic_f ∈ italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_ε , italic_R ) .

Proof. Integrating by parts and using Cauchy-Schwarz inequality, we obtain

εR(rRf(s)𝑑s)2e2λr𝑑r=12λ(εRf(s)𝑑s)2e2λεdr+superscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2superscript𝑒2𝜆𝑟differential-d𝑟limit-from12𝜆superscriptsuperscriptsubscript𝜀𝑅𝑓𝑠differential-d𝑠2superscript𝑒2𝜆𝜀𝑑𝑟\displaystyle\int\limits_{\varepsilon}^{R}\left(\int\limits_{r}^{R}f\left(s% \right)ds\right)^{2}e^{2\lambda r}dr=-\frac{1}{2\lambda}\left(\int\limits_{% \varepsilon}^{R}f\left(s\right)ds\right)^{2}e^{2\lambda\varepsilon}dr+∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT italic_d italic_r = - divide start_ARG 1 end_ARG start_ARG 2 italic_λ end_ARG ( ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_ε end_POSTSUPERSCRIPT italic_d italic_r +
+1λεR[eλrrRf(s)𝑑s][f(r)eλr]𝑑r1𝜆superscriptsubscript𝜀𝑅delimited-[]superscript𝑒𝜆𝑟superscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠delimited-[]𝑓𝑟superscript𝑒𝜆𝑟differential-d𝑟absent\displaystyle\hskip 28.45274pt+\frac{1}{\lambda}\int\limits_{\varepsilon}^{R}% \left[e^{\lambda r}\int\limits_{r}^{R}f\left(s\right)ds\right]\left[f\left(r% \right)e^{\lambda r}\right]dr\leq+ divide start_ARG 1 end_ARG start_ARG italic_λ end_ARG ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT [ italic_e start_POSTSUPERSCRIPT italic_λ italic_r end_POSTSUPERSCRIPT ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ] [ italic_f ( italic_r ) italic_e start_POSTSUPERSCRIPT italic_λ italic_r end_POSTSUPERSCRIPT ] italic_d italic_r ≤
1λ[εR(rRf(s)𝑑s)2e2λr𝑑r]1/2[εRf2(r)e2λr𝑑r]1/2.absent1𝜆superscriptdelimited-[]superscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2superscript𝑒2𝜆𝑟differential-d𝑟12superscriptdelimited-[]superscriptsubscript𝜀𝑅superscript𝑓2𝑟superscript𝑒2𝜆𝑟differential-d𝑟12\displaystyle\leq\frac{1}{\lambda}\left[\int\limits_{\varepsilon}^{R}\left(% \int\limits_{r}^{R}f\left(s\right)ds\right)^{2}e^{2\lambda r}dr\right]^{1/2}% \left[\int\limits_{\varepsilon}^{R}f^{2}\left(r\right)e^{2\lambda r}dr\right]^% {1/2}.≤ divide start_ARG 1 end_ARG start_ARG italic_λ end_ARG [ ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT italic_d italic_r ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT [ ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r ) italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT italic_d italic_r ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT .

Thus, we have proven that

(46) εR(rRf(s)𝑑s)2φλ(r)𝑑r(1/λ)[εR(rRf(s)𝑑s)2φλ(r)𝑑r]1/2[εRf2(r)φλ(r)𝑑r]1/2.superscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2subscript𝜑𝜆𝑟differential-d𝑟absentabsent1𝜆superscriptdelimited-[]superscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2subscript𝜑𝜆𝑟differential-d𝑟12superscriptdelimited-[]superscriptsubscript𝜀𝑅superscript𝑓2𝑟subscript𝜑𝜆𝑟differential-d𝑟12\left.\begin{array}[]{c}\int\limits_{\varepsilon}^{R}\left(\int\limits_{r}^{R}% f\left(s\right)ds\right)^{2}\varphi_{\lambda}\left(r\right)dr\leq\\ \leq\left(1/\lambda\right)\left[\int\limits_{\varepsilon}^{R}\left(\int\limits% _{r}^{R}f\left(s\right)ds\right)^{2}\varphi_{\lambda}\left(r\right)dr\right]^{% 1/2}\left[\int\limits_{\varepsilon}^{R}f^{2}\left(r\right)\varphi_{\lambda}% \left(r\right)dr\right]^{1/2}.\end{array}\right.start_ARRAY start_ROW start_CELL ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) italic_d italic_r ≤ end_CELL end_ROW start_ROW start_CELL ≤ ( 1 / italic_λ ) [ ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) italic_d italic_r ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT [ ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_r ) italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) italic_d italic_r ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT . end_CELL end_ROW end_ARRAY

Dividing both sides of (46) by

[εR(rRf(s)𝑑s)2φλ(r)𝑑r]1/2superscriptdelimited-[]superscriptsubscript𝜀𝑅superscriptsuperscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠2subscript𝜑𝜆𝑟differential-d𝑟12\left[\int\limits_{\varepsilon}^{R}\left(\int\limits_{r}^{R}f\left(s\right)ds% \right)^{2}\varphi_{\lambda}\left(r\right)dr\right]^{1/2}[ ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( ∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_φ start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_r ) italic_d italic_r ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT

and then squaring both sides of the resulting inequality, we obtain (45). \square

A similar estimate was proven in [12]. However, estimate (45) is stronger than the one of [12] due to the presence of the multiplier 1/λ21superscript𝜆21/\lambda^{2}1 / italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT here instead of 1/λ1𝜆1/\lambda1 / italic_λ in [12] and 1/λ2<<1/λmuch-less-than1superscript𝜆21𝜆1/\lambda^{2}<<1/\lambda1 / italic_λ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < < 1 / italic_λ for λ>>1.much-greater-than𝜆1\lambda>>1.italic_λ > > 1 . Note that we use λ>>1much-greater-than𝜆1\lambda>>1italic_λ > > 1 below.

7 Partial Finite Differences

To solve boundary value problem (41), we use partial finite differences with respect to the variables z𝑧zitalic_z and φ.𝜑\varphi.italic_φ . without imposing an additional assumption. This is the second and the last step of our approximate mathematical model.  More precisely, we assume that problem (41) in finite differences with respect to z𝑧zitalic_z and φ,𝜑\varphi,italic_φ , where the grid step sizes hzsubscript𝑧h_{z}italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and hφsubscript𝜑h_{\varphi}italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT with respect to z𝑧zitalic_z and φ𝜑\varphiitalic_φ satisfy the following inequality:

(47) hz,hφh0>0,subscript𝑧subscript𝜑subscript00h_{z},h_{\varphi}\geq h_{0}>0,italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ≥ italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT > 0 ,

where the number h0subscript0h_{0}italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is fixed. Consider finite difference grids with respect to z𝑧zitalic_z and φ,𝜑\varphi,italic_φ ,

(48) 0=z0<z1<<zk1<zk=B,zizi1=hz,formulae-sequence0superscript𝑧0subscript𝑧1subscript𝑧𝑘1subscript𝑧𝑘𝐵subscript𝑧𝑖subscript𝑧𝑖1subscript𝑧\displaystyle 0=z^{0}<z_{1}<...<z_{k-1}<z_{k}=B,z_{i}-z_{i-1}=h_{z},0 = italic_z start_POSTSUPERSCRIPT 0 end_POSTSUPERSCRIPT < italic_z start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < … < italic_z start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT < italic_z start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT = italic_B , italic_z start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_z start_POSTSUBSCRIPT italic_i - 1 end_POSTSUBSCRIPT = italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT ,
(49) 00\displaystyle 0 =φ0<φ1<<φn1<φn=2π,φiφi1=hφ.formulae-sequenceabsentsubscript𝜑0subscript𝜑1subscript𝜑𝑛1subscript𝜑𝑛2𝜋subscript𝜑𝑖subscript𝜑𝑖1subscript𝜑\displaystyle=\varphi_{0}<\varphi_{1}<...<\varphi_{n-1}<\varphi_{n}=2\pi,% \varphi_{i}-\varphi_{i-1}=h_{\varphi}.= italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT < italic_φ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < … < italic_φ start_POSTSUBSCRIPT italic_n - 1 end_POSTSUBSCRIPT < italic_φ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT = 2 italic_π , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT - italic_φ start_POSTSUBSCRIPT italic_i - 1 end_POSTSUBSCRIPT = italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT .

Denote h=(hz,hφ).subscript𝑧subscript𝜑h=\left(h_{z},h_{\varphi}\right).italic_h = ( italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT ) . The vector function V(r,φ,z)𝑉𝑟𝜑𝑧V\left(r,\varphi,z\right)italic_V ( italic_r , italic_φ , italic_z ) is turned now in the semi-discrete function, which is defined on the grid (48), (49) for all r[ε,R],𝑟𝜀𝑅r\in\left[\varepsilon,R\right],italic_r ∈ [ italic_ε , italic_R ] , i.e. r𝑟ritalic_r remains a continuous variable. We denote this function Vh(r)superscript𝑉𝑟V^{h}\left(r\right)italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r ). The semi-discrete form of (37) is:

(50) Vh(r)={(u0,,uN1)T(r,φi,zj)}(i,j)=(0,0)(n,k), (r,φi,zj)Q~=(ε,R)×(0,2π)×(0,B).superscript𝑉𝑟superscriptsubscriptsuperscriptsubscript𝑢0subscript𝑢𝑁1𝑇𝑟subscript𝜑𝑖subscript𝑧𝑗𝑖𝑗00𝑛𝑘 𝑟subscript𝜑𝑖subscript𝑧𝑗~𝑄𝜀𝑅02𝜋0𝐵\left.\begin{array}[]{c}V^{h}\left(r\right)=\left\{\left(u_{0},...,u_{N-1}% \right)^{T}\left(r,\varphi_{i},z_{j}\right)\right\}_{\left(i,j\right)=\left(0,% 0\right)}^{\left(n,k\right)},\text{ }\\ \left(r,\varphi_{i},z_{j}\right)\in\widetilde{Q}=\left(\varepsilon,R\right)% \times\left(0,2\pi\right)\times\left(0,B\right).\end{array}\right.start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r ) = { ( italic_u start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , italic_u start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) } start_POSTSUBSCRIPT ( italic_i , italic_j ) = ( 0 , 0 ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n , italic_k ) end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) ∈ over~ start_ARG italic_Q end_ARG = ( italic_ε , italic_R ) × ( 0 , 2 italic_π ) × ( 0 , italic_B ) . end_CELL end_ROW end_ARRAY

And the semi-discrete form of (36) is:

(51) uh(r,φi,zj,z0)=s=0N1us(r,φi,zj)Φs(z0), (r,φi,zj,z0)Q.formulae-sequencesuperscript𝑢𝑟subscript𝜑𝑖subscript𝑧𝑗subscript𝑧0superscriptsubscript𝑠0𝑁1subscript𝑢𝑠𝑟subscript𝜑𝑖subscript𝑧𝑗subscriptΦ𝑠subscript𝑧0 𝑟subscript𝜑𝑖subscript𝑧𝑗subscript𝑧0𝑄u^{h}\left(r,\varphi_{i},z_{j},z_{0}\right)=\sum\limits_{s=0}^{N-1}u_{s}\left(% r,\varphi_{i},z_{j}\right)\Phi_{s}\left(z_{0}\right),\text{ }\left(r,\varphi_{% i},z_{j},z_{0}\right)\in Q.italic_u start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_s = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT italic_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) roman_Φ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∈ italic_Q .

7.1 Some specifics of partial finite differences

We use the second order approximation accuracy for derivatives zVh,φVh,subscript𝑧superscript𝑉subscript𝜑superscript𝑉\partial_{z}V^{h},\partial_{\varphi}V^{h},∂ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , ∂ start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ,

(52) zVh(r,φi,zj)subscript𝑧superscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑗\displaystyle\partial_{z}V^{h}\left(r,\varphi_{i},z_{j}\right)∂ start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) =Vh(r,φi,zj+1)Vh(r,φi,zj1)2hz, i[0,n],j[1,k1],formulae-sequenceabsentsuperscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑗1superscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑗12subscript𝑧formulae-sequence 𝑖0𝑛𝑗1𝑘1\displaystyle=\frac{V^{h}\left(r,\varphi_{i},z_{j+1}\right)-V^{h}\left(r,% \varphi_{i},z_{j-1}\right)}{2h_{z}},\text{ }i\in\left[0,n\right],j\in\left[1,k% -1\right],= divide start_ARG italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j + 1 end_POSTSUBSCRIPT ) - italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j - 1 end_POSTSUBSCRIPT ) end_ARG start_ARG 2 italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG , italic_i ∈ [ 0 , italic_n ] , italic_j ∈ [ 1 , italic_k - 1 ] ,
(53) φVh(r,φi,zj)subscript𝜑superscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑗\displaystyle\partial_{\varphi}V^{h}\left(r,\varphi_{i},z_{j}\right)∂ start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) =Vh(r,φi+1,zj)Vh(r,φi1,zj)2hφ, i[1,n1],j[0,k].formulae-sequenceabsentsuperscript𝑉𝑟subscript𝜑𝑖1subscript𝑧𝑗superscript𝑉𝑟subscript𝜑𝑖1subscript𝑧𝑗2subscript𝜑formulae-sequence 𝑖1𝑛1𝑗0𝑘\displaystyle=\frac{V^{h}\left(r,\varphi_{i+1},z_{j}\right)-V^{h}\left(r,% \varphi_{i-1},z_{j}\right)}{2h_{\varphi}},\text{ }i\in\left[1,n-1\right],j\in% \left[0,k\right].= divide start_ARG italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i + 1 end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) - italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i - 1 end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) end_ARG start_ARG 2 italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT end_ARG , italic_i ∈ [ 1 , italic_n - 1 ] , italic_j ∈ [ 0 , italic_k ] .

A separate question is about ensuring boundary conditions in the second and third lines of (41) as well as the periodicity condition in the fourth line of (41). To satisfy boundary conditions in the second and third lines of (41), we use in (52):

(54) Vzh(r,φi,z1)superscriptsubscript𝑉𝑧𝑟subscript𝜑𝑖subscript𝑧1\displaystyle V_{z}^{h}\left(r,\varphi_{i},z_{1}\right)italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) =Vh(r,φi,z2) G0(r,φi)2hz, i[0,n],formulae-sequenceabsentsuperscript𝑉𝑟subscript𝜑𝑖subscript𝑧2 subscript𝐺0𝑟subscript𝜑𝑖2subscript𝑧 𝑖0𝑛\displaystyle=\frac{V^{h}\left(r,\varphi_{i},z_{2}\right)-\text{ }G_{0}\left(r% ,\varphi_{i}\right)}{2h_{z}},\text{ }i\in\left[0,n\right],= divide start_ARG italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) - italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) end_ARG start_ARG 2 italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG , italic_i ∈ [ 0 , italic_n ] ,
(55) Vzh(r,φi,zk1)superscriptsubscript𝑉𝑧𝑟subscript𝜑𝑖subscript𝑧𝑘1\displaystyle V_{z}^{h}\left(r,\varphi_{i},z_{k-1}\right)italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_k - 1 end_POSTSUBSCRIPT ) = GB(r,φi) Vh(r,φi,zk2)2hz, i[0,n].formulae-sequenceabsent subscript𝐺𝐵𝑟subscript𝜑𝑖 superscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑘22subscript𝑧 𝑖0𝑛\displaystyle=\frac{\text{ }G_{B}\left(r,\varphi_{i}\right)-\text{ }V^{h}\left% (r,\varphi_{i},z_{k-2}\right)}{2h_{z}},\text{ }i\in\left[0,n\right].= divide start_ARG italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) - italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_k - 2 end_POSTSUBSCRIPT ) end_ARG start_ARG 2 italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG , italic_i ∈ [ 0 , italic_n ] .

In addition, to ensure the periodicity condition in the fourth line of (41), we set

(56) V(r,φ0,zj)=V(r,φn,zj), j[0,k]formulae-sequence𝑉𝑟subscript𝜑0subscript𝑧𝑗𝑉𝑟subscript𝜑𝑛subscript𝑧𝑗 𝑗0𝑘V\left(r,\varphi_{0},z_{j}\right)=V\left(r,\varphi_{n},z_{j}\right),\text{ }j% \in\left[0,k\right]italic_V ( italic_r , italic_φ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) = italic_V ( italic_r , italic_φ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) , italic_j ∈ [ 0 , italic_k ]

with corresponding changes in formulas (52)-(55) for i=0,n.𝑖0𝑛i=0,n.italic_i = 0 , italic_n .

Let H𝐻Hitalic_H be a Hilbert space. Consider the direct product of N𝑁Nitalic_N such spaces HN=H×H×H,subscript𝐻𝑁𝐻𝐻𝐻H_{N}=H\times H\times...H,italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT = italic_H × italic_H × … italic_H , N𝑁Nitalic_N times. If UHsubscriptnorm𝑈𝐻\left\|U\right\|_{H}∥ italic_U ∥ start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT is the norm of UH𝑈𝐻U\in Hitalic_U ∈ italic_H in H𝐻Hitalic_H, then the norm of the vector function U~=(U0,,UN1)THN~𝑈superscriptsubscript𝑈0subscript𝑈𝑁1𝑇subscript𝐻𝑁\widetilde{U}=\left(U_{0},...,U_{N-1}\right)^{T}\in H_{N}over~ start_ARG italic_U end_ARG = ( italic_U start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , … , italic_U start_POSTSUBSCRIPT italic_N - 1 end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT is set below as:

U~HN=[i=0N1UiH2]1/2.subscriptnorm~𝑈subscript𝐻𝑁superscriptdelimited-[]superscriptsubscript𝑖0𝑁1superscriptsubscriptnormsubscript𝑈𝑖𝐻212\left\|\widetilde{U}\right\|_{H_{N}}=\left[\sum\limits_{i=0}^{N-1}\left\|U_{i}% \right\|_{H}^{2}\right]^{1/2}.∥ over~ start_ARG italic_U end_ARG ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT end_POSTSUBSCRIPT = [ ∑ start_POSTSUBSCRIPT italic_i = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_N - 1 end_POSTSUPERSCRIPT ∥ italic_U start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∥ start_POSTSUBSCRIPT italic_H end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT 1 / 2 end_POSTSUPERSCRIPT .

We now introduce some Hilbert spaces of semi-discrete functions, similarly with [15]. Denote 𝐱h={(r,φi,zj)(i,j)=(0,0)(n,k),r(ε,R)}.superscript𝐱superscriptsubscript𝑟subscript𝜑𝑖subscript𝑧𝑗𝑖𝑗00𝑛𝑘𝑟𝜀𝑅\mathbf{x}^{h}=\left\{\left(r,\varphi_{i},z_{j}\right)_{\left(i,j\right)=\left% (0,0\right)}^{\left(n,k\right)},r\in\left(\varepsilon,R\right)\right\}.bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT = { ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT ( italic_i , italic_j ) = ( 0 , 0 ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n , italic_k ) end_POSTSUPERSCRIPT , italic_r ∈ ( italic_ε , italic_R ) } . Along with the above notation Vh(rV^{h}(ritalic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r, φisubscript𝜑𝑖\varphi_{i}italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, zj),z_{j}),italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) , we also use the equivalent one, which is Vh(𝐱h).superscript𝑉superscript𝐱V^{h}\left(\mathbf{x}^{h}\right).italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) . For the domain D𝐷Ditalic_D defined in (8) we denote

Dh={(r,φi,zj)(i,j)=(0,0)(n,k), r(ε,R)}.superscript𝐷superscriptsubscript𝑟subscript𝜑𝑖subscript𝑧𝑗𝑖𝑗00𝑛𝑘 𝑟𝜀𝑅D^{h}=\left\{\left(r,\varphi_{i},z_{j}\right)_{\left(i,j\right)=\left(0,0% \right)}^{\left(n,k\right)},\text{ }r\in\left(\varepsilon,R\right)\right\}.italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT = { ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) start_POSTSUBSCRIPT ( italic_i , italic_j ) = ( 0 , 0 ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n , italic_k ) end_POSTSUPERSCRIPT , italic_r ∈ ( italic_ε , italic_R ) } .

Next, using (48) and (49), we define those Hilbert spaces as:

(57) L2,Nh(Dh)={Vh(𝐱h):Vh(𝐱h)L2,Nh(Dh)2==(i,j)=(0,0)(n,k){εR(Vh(r,φi,zj))2𝑑r}<},HN1,h(Dh)={Vh(𝐱h):VhHN1,h(Dh)2==VzhL2Nh(Dh)2+VφhL2Nh(Dh)2+VhL2Nh(Dh)2<},superscriptsubscript𝐿2𝑁superscript𝐷:superscript𝑉superscript𝐱superscriptsubscriptnormsuperscript𝑉superscript𝐱superscriptsubscript𝐿2𝑁superscript𝐷2absentabsentsuperscriptsubscript𝑖𝑗00𝑛𝑘superscriptsubscript𝜀𝑅superscriptsuperscript𝑉𝑟subscript𝜑𝑖subscript𝑧𝑗2differential-d𝑟superscriptsubscript𝐻𝑁1superscript𝐷:superscript𝑉superscript𝐱superscriptsubscriptnormsuperscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷2absentabsentsuperscriptsubscriptnormsuperscriptsubscript𝑉𝑧superscriptsubscript𝐿2𝑁superscript𝐷2superscriptsubscriptnormsuperscriptsubscript𝑉𝜑superscriptsubscript𝐿2𝑁superscript𝐷2superscriptsubscriptnormsuperscript𝑉superscriptsubscript𝐿2𝑁superscript𝐷2\left.\begin{array}[]{c}L_{2,N}^{h}\left(D^{h}\right)=\left\{\begin{array}[]{c% }V^{h}(\mathbf{x}^{h}):\left\|V^{h}(\mathbf{x}^{h})\right\|_{L_{2,N}^{h}\left(% D^{h}\right)}^{2}=\\ =\sum\limits_{\left(i,j\right)=\left(0,0\right)}^{\left(n,k\right)}\left\{\int% \limits_{\varepsilon}^{R}\left(V^{h}\left(r,\varphi_{i},z_{j}\right)\right)^{2% }dr\right\}<\infty\end{array}\right\},\\ H_{N}^{1,h}\left(D^{h}\right)=\left\{\begin{array}[]{c}V^{h}(\mathbf{x}^{h}):% \left\|V^{h}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}^{2}=\\ =\left\|V_{z}^{h}\right\|_{L_{2N}^{h}\left(D^{h}\right)}^{2}+\left\|V_{\varphi% }^{h}\right\|_{L_{2N}^{h}\left(D^{h}\right)}^{2}+\left\|V^{h}\right\|_{L_{2N}^% {h}\left(D^{h}\right)}^{2}<\infty\end{array}\right\},\end{array}\right.start_ARRAY start_ROW start_CELL italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) : ∥ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = end_CELL end_ROW start_ROW start_CELL = ∑ start_POSTSUBSCRIPT ( italic_i , italic_j ) = ( 0 , 0 ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_n , italic_k ) end_POSTSUPERSCRIPT { ∫ start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_r } < ∞ end_CELL end_ROW end_ARRAY } , end_CELL end_ROW start_ROW start_CELL italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) : ∥ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = end_CELL end_ROW start_ROW start_CELL = ∥ italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ∥ italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ∥ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < ∞ end_CELL end_ROW end_ARRAY } , end_CELL end_ROW end_ARRAY
(58) H0,N1,h(Dh)={Vh(𝐱h)HN1,h(Dh):Vh(r,φi,0)=Vh(r,φi,B)=0,i[0,n]}.superscriptsubscript𝐻0𝑁1superscript𝐷:superscript𝑉superscript𝐱superscriptsubscript𝐻𝑁1superscript𝐷absentformulae-sequencesuperscript𝑉𝑟subscript𝜑𝑖0superscript𝑉𝑟subscript𝜑𝑖𝐵0𝑖0𝑛H_{0,N}^{1,h}\left(D^{h}\right)=\left\{\begin{array}[]{c}V^{h}(\mathbf{x}^{h})% \in H_{N}^{1,h}\left(D^{h}\right):\\ V^{h}\left(r,\varphi_{i},0\right)=V^{h}\left(r,\varphi_{i},B\right)=0,i\in% \left[0,n\right]\end{array}\right\}.italic_H start_POSTSUBSCRIPT 0 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) : end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , 0 ) = italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_B ) = 0 , italic_i ∈ [ 0 , italic_n ] end_CELL end_ROW end_ARRAY } .

In (57) derivatives Vzh,Vφhsuperscriptsubscript𝑉𝑧superscriptsubscript𝑉𝜑V_{z}^{h},V_{\varphi}^{h}italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT are understood as in (52)-(55) being supplied by the periodicity condition (56).

7.2 Completion of the approximate mathematical model

It was stated in the end of section 4 that we construct an approximate mathematical model for our travel time tomography problem, and then we verify this model computationally. The first step of this model is the truncated Fourier-like series (36). We remind that truncated series are used quite often in developments of numerical methods for various inverse problems. On the other hand, convergencies of such methods as the number of terms tends to infinity are usually not proven, see, e.g. [6, 8]. This is caused by the ill-posed nature of inverse problems. The second and the final step of that approximate mathematical model is the assumption of partial finite differences with condition (47). This step is reflected in (48)-(56).

8 Convexification

8.1 The cost functional of the convexification

Let M>0𝑀0M>0italic_M > 0 be an arbitrary number. Let c>0𝑐0c>0italic_c > 0 be the number in (17). Define the set B(M)𝐵𝑀B\left(M\right)italic_B ( italic_M ) of vector functions Vh(r)superscript𝑉𝑟V^{h}\left(r\right)italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r ) of (50) as

(59) B(M)={VhHN1,h(Dh):uh(r,φi,zj,z0)c,z0[0,B]VhHN1,h(Dh)<M,Vh(r,φi,0)= G0(r,φi),Vh(r,φi,B)= GB(r,φi),i[0,n],j[0,k],(54)-(56) hold,}𝐵𝑀:superscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷formulae-sequencesuperscript𝑢𝑟subscript𝜑𝑖subscript𝑧𝑗subscript𝑧0𝑐subscript𝑧00𝐵subscriptnormsuperscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷𝑀formulae-sequencesuperscript𝑉𝑟subscript𝜑𝑖0 subscript𝐺0𝑟subscript𝜑𝑖superscript𝑉𝑟subscript𝜑𝑖𝐵 subscript𝐺𝐵𝑟subscript𝜑𝑖formulae-sequence𝑖0𝑛𝑗0𝑘(54)-(56) hold,\left.\begin{array}[]{c}B\left(M\right)=\left\{\begin{array}[]{c}V^{h}\in H_{N% }^{1,h}\left(D^{h}\right):u^{h}\left(r,\varphi_{i},z_{j},z_{0}\right)\geq c,z_% {0}\in\left[0,B\right]\\ \left\|V^{h}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}<M,\\ V^{h}\left(r,\varphi_{i},0\right)=\text{ }G_{0}\left(r,\varphi_{i}\right),V^{h% }\left(r,\varphi_{i},B\right)=\text{ }G_{B}\left(r,\varphi_{i}\right),\\ i\in\left[0,n\right],j\in\left[0,k\right],\\ \text{(\ref{7.6})-(\ref{7.8}) hold,}\end{array}\right\}\end{array}\right.start_ARRAY start_ROW start_CELL italic_B ( italic_M ) = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) : italic_u start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≥ italic_c , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ [ 0 , italic_B ] end_CELL end_ROW start_ROW start_CELL ∥ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_M , end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , 0 ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_B ) = italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL italic_i ∈ [ 0 , italic_n ] , italic_j ∈ [ 0 , italic_k ] , end_CELL end_ROW start_ROW start_CELL ( )-( ) hold, end_CELL end_ROW end_ARRAY } end_CELL end_ROW end_ARRAY

where the space HN1,h(Dh)superscriptsubscript𝐻𝑁1superscript𝐷H_{N}^{1,h}\left(D^{h}\right)italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) is defined in (57).

Thus, we have replaced BVP2 with BVP3:

Boundary Value Problem 3 (BVP3). Find the vector function VhB(M)¯.superscript𝑉¯𝐵𝑀V^{h}\in\overline{B\left(M\right)}.italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG . 

To solve BVP3, we consider the following minimization problem:

Minimization Problem. Minimize the following functional on the set B(M)¯::¯𝐵𝑀absent\overline{B\left(M\right)}:over¯ start_ARG italic_B ( italic_M ) end_ARG :

(60) Jλ(Vh)=(ANVh+S(Vφh,Vzh,r))eλrL2,Nh(Dh)2.subscript𝐽𝜆superscript𝑉superscriptsubscriptnormsubscript𝐴𝑁superscript𝑉𝑆superscriptsubscript𝑉𝜑superscriptsubscript𝑉𝑧𝑟superscript𝑒𝜆𝑟superscriptsubscript𝐿2𝑁superscript𝐷2J_{\lambda}\left(V^{h}\right)=\left\|\left(A_{N}V^{h}+S\left(V_{\varphi}^{h},V% _{z}^{h},r\right)\right)e^{\lambda r}\right\|_{L_{2,N}^{h}\left(D^{h}\right)}^% {2}.italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) = ∥ ( italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT + italic_S ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_r ) ) italic_e start_POSTSUPERSCRIPT italic_λ italic_r end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .

Recall that e2λrsuperscript𝑒2𝜆𝑟e^{2\lambda r}italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT is the Carleman Weight Function in (44), (45) for the Volterra integral operator

rRf(s)𝑑s.superscriptsubscript𝑟𝑅𝑓𝑠differential-d𝑠\int\nolimits_{r}^{R}f\left(s\right)ds.∫ start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_R end_POSTSUPERSCRIPT italic_f ( italic_s ) italic_d italic_s .

Note that we do not use the Tikhonov penalization term [23] in (60), which is rare in the theory of Inverse Problems.

8.2 The global convergence

We now formulate theorems of the global convergence analysis for the Minimization Problem. Both their formulations and proofs are similar with those of [15]. Therefore, we omit proofs. We now briefly explain the role of the Carleman Weight Function e2λrsuperscript𝑒2𝜆𝑟e^{2\lambda r}italic_e start_POSTSUPERSCRIPT 2 italic_λ italic_r end_POSTSUPERSCRIPT in the proof of our central result, which is Theorem 8.1. Given the set B(M)𝐵𝑀B\left(M\right)italic_B ( italic_M ) in (59), this theorem claims the strong convexity of the functional Jλ(Vh)subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) on B(M)¯¯𝐵𝑀\overline{B\left(M\right)}over¯ start_ARG italic_B ( italic_M ) end_ARG for sufficiently large values of the parameter λ>>1.much-greater-than𝜆1\lambda>>1.italic_λ > > 1 . First, consider for a moment only the quadratic functional J~λ(Vh)=(ANVh)eλrL2,Nh(Dh)2.subscript~𝐽𝜆superscript𝑉superscriptsubscriptnormsubscript𝐴𝑁superscript𝑉superscript𝑒𝜆𝑟superscriptsubscript𝐿2𝑁superscript𝐷2\widetilde{J}_{\lambda}\left(V^{h}\right)=\left\|\left(A_{N}V^{h}\right)e^{% \lambda r}\right\|_{L_{2,N}^{h}\left(D^{h}\right)}^{2}.over~ start_ARG italic_J end_ARG start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) = ∥ ( italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_λ italic_r end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . Since the matrix ANsubscript𝐴𝑁A_{N}italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT is invertible, then this functional is strongly convex on the entire space L2,Nh(Dh)superscriptsubscript𝐿2𝑁superscript𝐷L_{2,N}^{h}\left(D^{h}\right)italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) for any value of λ.𝜆\lambda.italic_λ . However, the presence of the nonlinear term S(Vφh,Vzh,r)𝑆superscriptsubscript𝑉𝜑superscriptsubscript𝑉𝑧𝑟S\left(V_{\varphi}^{h},V_{z}^{h},r\right)italic_S ( italic_V start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_r ) in (60) might destroy the strong convexity property. To dominate this term, we use in the proof of Theorem 8.1 the Carleman estimate of Theorem 6.1, which, roughly speaking, states that the nonlinear term is dominated by the strongly convex quadratic term (ANVh)eλrL2,Nh(Dh)2.superscriptsubscriptnormsubscript𝐴𝑁superscript𝑉superscript𝑒𝜆𝑟superscriptsubscript𝐿2𝑁superscript𝐷2\left\|\left(A_{N}V^{h}\right)e^{\lambda r}\right\|_{L_{2,N}^{h}\left(D^{h}% \right)}^{2}.∥ ( italic_A start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) italic_e start_POSTSUPERSCRIPT italic_λ italic_r end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT .

Another question is on what exactly “sufficiently large λ"𝜆"\lambda"italic_λ " is. Philosophically this is very similar with any asymptotic theory. Indeed, such a theory usually claims that as soon as a certain parameter A𝐴Aitalic_A is sufficiently large, a certain formula B𝐵Bitalic_B is valid with a good accuracy. In a computational practice, however, one always works with specific ranges of many parameters. Therefore, only results of numerical experiments can establish which exactly values of A𝐴Aitalic_A are sufficient to ensure a good accuracy of B𝐵Bitalic_B. In particular, we demonstrate below that λ=3𝜆3\lambda=3italic_λ = 3 is an optimal value of the parameter λ𝜆\lambdaitalic_λ for our computations. We also note that in all previous publications on the convexification of this research group optimal values of λ𝜆\lambdaitalic_λ were λ[1,5]𝜆15\lambda\in\left[1,5\right]italic_λ ∈ [ 1 , 5 ], see, e.g. [13]-[16].

Theorem 8.1 (strong convexity). Assume that conditions (6)-(11) hold, and let c>0𝑐0c>0italic_c > 0 be the number in (17) and (59). Then:

1. The functional Jλ(Vh)subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) has the Fréchet derivative Jλ(Vh)H0,N1,h(Dh)superscriptsubscript𝐽𝜆superscript𝑉superscriptsubscript𝐻0𝑁1superscript𝐷J_{\lambda}^{\prime}\left(V^{h}\right)\in H_{0,N}^{1,h}\left(D^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∈ italic_H start_POSTSUBSCRIPT 0 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) at every point VhB(M)¯superscript𝑉¯𝐵𝑀V^{h}\in\overline{B\left(M\right)}italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG and for all λ>0.𝜆0\lambda>0.italic_λ > 0 .  Furthermore, the Fréchet derivative Jλ(Vh)superscriptsubscript𝐽𝜆superscript𝑉J_{\lambda}^{\prime}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) satisfies Lipschitz continuity condition on B(M)¯,¯𝐵𝑀\overline{B\left(M\right)},over¯ start_ARG italic_B ( italic_M ) end_ARG , i.e. there exists a number

C¯=C¯(h0,c,M,N,Dh,λ)>0 ¯𝐶¯𝐶subscript0𝑐𝑀𝑁superscript𝐷𝜆0 \overline{C}=\overline{C}\left(h_{0},c,M,N,D^{h},\lambda\right)>0\emph{\ }over¯ start_ARG italic_C end_ARG = over¯ start_ARG italic_C end_ARG ( italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_c , italic_M , italic_N , italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_λ ) > 0

depending only on listed parameters such that the following estimate holds:

Jλ(V2h)Jλ(V1h)L2,Nh(Dh)C¯V2hV1hHN1,h(Ωh), V1h,V2h B(M)¯.formulae-sequencesubscriptnormsuperscriptsubscript𝐽𝜆superscriptsubscript𝑉2superscriptsubscript𝐽𝜆superscriptsubscript𝑉1superscriptsubscript𝐿2𝑁superscript𝐷¯𝐶subscriptnormsuperscriptsubscript𝑉2superscriptsubscript𝑉1superscriptsubscript𝐻𝑁1superscriptΩ for-allsuperscriptsubscript𝑉1superscriptsubscript𝑉2 ¯𝐵𝑀\left\|J_{\lambda}^{\prime}\left(V_{2}^{h}\right)-J_{\lambda}^{\prime}\left(V_% {1}^{h}\right)\right\|_{L_{2,N}^{h}\left(D^{h}\right)}\leq\overline{C}\left\|V% _{2}^{h}-V_{1}^{h}\right\|_{H_{N}^{1,h}\left(\Omega^{h}\right)},\text{ }% \forall V_{1}^{h},V_{2}^{h}\in\emph{\ }\overline{B\left(M\right)}.∥ italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) - italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT ≤ over¯ start_ARG italic_C end_ARG ∥ italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( roman_Ω start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT , ∀ italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG .

2. There exist a sufficiently large number λ01subscript𝜆01\lambda_{0}\geq 1italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ≥ 1  and a number C>0,𝐶0C>0,italic_C > 0 ,

(61) λ0=λ0(h0,c,M,N,Dh)1, C=C(h0,c,M,N,Dh)>0,formulae-sequencesubscript𝜆0subscript𝜆0subscript0𝑐𝑀𝑁superscript𝐷1 𝐶𝐶subscript0𝑐𝑀𝑁superscript𝐷0\lambda_{0}=\lambda_{0}\left(h_{0},c,M,N,D^{h}\right)\geq 1,\text{ }C=C\left(h% _{0},c,M,N,D^{h}\right)>0,italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_c , italic_M , italic_N , italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ≥ 1 , italic_C = italic_C ( italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_c , italic_M , italic_N , italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) > 0 ,

both numbers depending only on listed parameters, such that for every λλ0𝜆subscript𝜆0\lambda\geq\lambda_{0}italic_λ ≥ italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT the functional Jλ(Vh)subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) is strongly convex on the set B(M)¯.¯𝐵𝑀\overline{B\left(M\right)}.over¯ start_ARG italic_B ( italic_M ) end_ARG . The strong convexity means that the following estimate holds::::

Jλ(V2h)Jλ(V1h)Jλ(V1h)(V2hV1h)CV2hV1hH2,N1,h(Ωh)2,V1h,V2h B(M)¯.formulae-sequencesubscript𝐽𝜆superscriptsubscript𝑉2subscript𝐽𝜆superscriptsubscript𝑉1superscriptsubscript𝐽𝜆superscriptsubscript𝑉1superscriptsubscript𝑉2superscriptsubscript𝑉1𝐶superscriptsubscriptdelimited-∥∥superscriptsubscript𝑉2superscriptsubscript𝑉1superscriptsubscript𝐻2𝑁1superscriptΩ2for-allsuperscriptsubscript𝑉1superscriptsubscript𝑉2 ¯𝐵𝑀\begin{split}J_{\lambda}\left(V_{2}^{h}\right)-J_{\lambda}\left(V_{1}^{h}% \right)-J_{\lambda}^{\prime}&\left(V_{1}^{h}\right)\left(V_{2}^{h}-V_{1}^{h}% \right)\geq C\left\|V_{2}^{h}-V_{1}^{h}\right\|_{H_{2,N}^{1,h}\left(\Omega^{h}% \right)}^{2},\\ &\forall V_{1}^{h},V_{2}^{h}\in\emph{\ }\overline{B\left(M\right)}.\end{split}start_ROW start_CELL italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) - italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) - italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL start_CELL ( italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ( italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ≥ italic_C ∥ italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( roman_Ω start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL ∀ italic_V start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , italic_V start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG . end_CELL end_ROW

3. Furthermore, for every λλ0𝜆subscript𝜆0\lambda\geq\lambda_{0}italic_λ ≥ italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT there exists unique minimizer Vmin,λhB(M)¯superscriptsubscript𝑉𝜆¯𝐵𝑀V_{\min,\lambda}^{h}\in\ \overline{B\left(M\right)}italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG of the functional Jλ(Vh)subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) on the set B(M)¯¯𝐵𝑀\ \overline{B\left(M\right)}over¯ start_ARG italic_B ( italic_M ) end_ARG and the following inequality holds:

Jλ(Vmin,λh)(VhVmin,λh)0, VhB(M)¯.formulae-sequencesuperscriptsubscript𝐽𝜆superscriptsubscript𝑉𝜆superscript𝑉superscriptsubscript𝑉𝜆0 for-allsuperscript𝑉¯𝐵𝑀J_{\lambda}^{\prime}\left(V_{\min,\lambda}^{h}\right)\left(V^{h}-V_{\min,% \lambda}^{h}\right)\geq 0,\text{ }\forall V^{h}\in\overline{B\left(M\right)}.italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ≥ 0 , ∀ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG .

Below C𝐶Citalic_C denotes different positive numbers depending only on parameters listed in (61). In the regularization theory [23], the minimizer Vmin,λhsuperscriptsubscript𝑉𝜆V_{\min,\lambda}^{h}italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT of functional (60) is called “regularized solution”. It is important to estimate the accuracy of the regularized solution depending on the level of the noise in the data. To do this, we recall first that, following the regularization theory, we need to assume the existence of the “ideal” solution of BVP3, i.e. solution with the noiseless data. The ideal solution is also called “exact” solution. We denote this solution VhHN1,h(Dh).superscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷V^{h\ast}\in H_{N}^{1,h}\left(D^{h}\right).italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) . We denote the noiseless data in the fourth line of (59) as G0(r,φi),GB(r,φi).superscriptsubscript𝐺0𝑟subscript𝜑𝑖superscriptsubscript𝐺𝐵𝑟subscript𝜑𝑖G_{0}^{\ast}\left(r,\varphi_{i}\right),G_{B}^{\ast}\left(r,\varphi_{i}\right).italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) . We assume that the exact solution VhB(M),superscript𝑉superscript𝐵𝑀V^{h\ast}\in B^{\ast}\left(M\right),italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∈ italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M ) , where B(M)superscript𝐵𝑀B^{\ast}\left(M\right)italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M ) is the following analog of the set B(M)𝐵𝑀B\left(M\right)italic_B ( italic_M ) in (59)

(62) B(M)=={VhHN1,h(Dh):uh(r,φi,zj,z0)c,z0[0,B]VhHN1,h(Dh)<M,Vh(r,φi,0)= G0(r,φi),Vh(r,φi,B)= GB(r,φi),i[0,n],j[0,k],(54)-(56) hold.}superscript𝐵𝑀absentabsent:superscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷formulae-sequencesuperscript𝑢𝑟subscript𝜑𝑖subscript𝑧𝑗subscript𝑧0𝑐subscript𝑧00𝐵subscriptnormsuperscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷𝑀formulae-sequencesuperscript𝑉𝑟subscript𝜑𝑖0 superscriptsubscript𝐺0𝑟subscript𝜑𝑖superscript𝑉𝑟subscript𝜑𝑖𝐵 superscriptsubscript𝐺𝐵𝑟subscript𝜑𝑖formulae-sequence𝑖0𝑛𝑗0𝑘(54)-(56) hold.\left.\begin{array}[]{c}B^{\ast}\left(M\right)=\\ =\left\{\begin{array}[]{c}V^{h}\in H_{N}^{1,h}\left(D^{h}\right):u^{h}\left(r,% \varphi_{i},z_{j},z_{0}\right)\geq c,z_{0}\in\left[0,B\right]\\ \left\|V^{h}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}<M,\\ V^{h}\left(r,\varphi_{i},0\right)=\text{ }G_{0}^{\ast}\left(r,\varphi_{i}% \right),V^{h}\left(r,\varphi_{i},B\right)=\text{ }G_{B}^{\ast}\left(r,\varphi_% {i}\right),\\ i\in\left[0,n\right],j\in\left[0,k\right],\\ \text{(\ref{7.6})-(\ref{7.8}) hold.}\end{array}\right\}\end{array}\right.start_ARRAY start_ROW start_CELL italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M ) = end_CELL end_ROW start_ROW start_CELL = { start_ARRAY start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) : italic_u start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ≥ italic_c , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ [ 0 , italic_B ] end_CELL end_ROW start_ROW start_CELL ∥ italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_M , end_CELL end_ROW start_ROW start_CELL italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , 0 ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_B ) = italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL italic_i ∈ [ 0 , italic_n ] , italic_j ∈ [ 0 , italic_k ] , end_CELL end_ROW start_ROW start_CELL ( )-( ) hold. end_CELL end_ROW end_ARRAY } end_CELL end_ROW end_ARRAY

Let δ(0,1)𝛿01\delta\in\left(0,1\right)italic_δ ∈ ( 0 , 1 ) be the level of the noise in the boundary data G0,GBsuperscriptsubscript𝐺0superscriptsubscript𝐺𝐵G_{0}^{\ast},G_{B}^{\ast}italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT. To simplify the presentation, we assume that the noise is not introduced in the function p(φ,z,z0)𝑝𝜑𝑧subscript𝑧0p\left(\varphi,z,z_{0}\right)italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in (14), i.e. we assume that p(φ,z,z0)=p(φ,z,z0),𝑝𝜑𝑧subscript𝑧0superscript𝑝𝜑𝑧subscript𝑧0p\left(\varphi,z,z_{0}\right)=p^{\ast}\left(\varphi,z,z_{0}\right),italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) = italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , although this case can also be included. Suppose that there exist such extensions Ph(𝐱h)HN1,h(Dh)superscript𝑃superscript𝐱superscriptsubscript𝐻𝑁1superscript𝐷P^{h}\left(\mathbf{x}^{h}\right)\in H_{N}^{1,h}\left(D^{h}\right)italic_P start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) and Ph(𝐱h)HN1,h(Dh)superscript𝑃superscript𝐱superscriptsubscript𝐻𝑁1superscript𝐷P^{h\ast}\left(\mathbf{x}^{h}\right)\in H_{N}^{1,h}\left(D^{h}\right)italic_P start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∈ italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) of the pairs of boundary data (G0(r,φi),GB(r,φi))i=0nsuperscriptsubscriptsubscript𝐺0𝑟subscript𝜑𝑖subscript𝐺𝐵𝑟subscript𝜑𝑖𝑖0𝑛\left(G_{0}\left(r,\varphi_{i}\right),G_{B}\left(r,\varphi_{i}\right)\right)_{% i=0}^{n}( italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ) start_POSTSUBSCRIPT italic_i = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT and (G0(r,φi),GB(r,φi))i=0nsuperscriptsubscriptsuperscriptsubscript𝐺0𝑟subscript𝜑𝑖superscriptsubscript𝐺𝐵𝑟subscript𝜑𝑖𝑖0𝑛\left(G_{0}^{\ast}\left(r,\varphi_{i}\right),G_{B}^{\ast}\left(r,\varphi_{i}% \right)\right)_{i=0}^{n}( italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) ) start_POSTSUBSCRIPT italic_i = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT that for i[0,n]𝑖0𝑛i\in\left[0,n\right]italic_i ∈ [ 0 , italic_n ] and r(ε,R)𝑟𝜀𝑅r\in\left(\varepsilon,R\right)italic_r ∈ ( italic_ε , italic_R )

(63) Ph(r,φi,0)=G0(r,φi), Ph(r,φi,B)=GB(r,φi),Ph(r,φi,0)=G0(r,φi), Ph(r,φi,B)=GB(r,φi),PhHN1,h(Dh)<M, PhHN1,h(Dh)<M,PhPhHN1,h(Dh)<δ. formulae-sequencesuperscript𝑃𝑟subscript𝜑𝑖0subscript𝐺0𝑟subscript𝜑𝑖 superscript𝑃𝑟subscript𝜑𝑖𝐵subscript𝐺𝐵𝑟subscript𝜑𝑖formulae-sequencesuperscript𝑃𝑟subscript𝜑𝑖0superscriptsubscript𝐺0𝑟subscript𝜑𝑖 superscript𝑃𝑟subscript𝜑𝑖𝐵superscriptsubscript𝐺𝐵𝑟subscript𝜑𝑖evaluated-atsubscriptnormsuperscript𝑃superscriptsubscript𝐻𝑁1superscript𝐷bra𝑀 superscript𝑃superscriptsubscript𝐻𝑁1superscript𝐷𝑀subscriptnormsuperscript𝑃superscript𝑃superscriptsubscript𝐻𝑁1superscript𝐷𝛿 \left.\begin{array}[]{c}P^{h}\left(r,\varphi_{i},0\right)=G_{0}\left(r,\varphi% _{i}\right),\text{ }P^{h}\left(r,\varphi_{i},B\right)=G_{B}\left(r,\varphi_{i}% \right),\\ P^{h\ast}\left(r,\varphi_{i},0\right)=G_{0}^{\ast}\left(r,\varphi_{i}\right),% \text{ }P^{h\ast}\left(r,\varphi_{i},B\right)=G_{B}^{\ast}\left(r,\varphi_{i}% \right),\\ \left\|P^{h}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}<M,\text{ }\left\|P^{h\ast% }\right\|_{H_{N}^{1,h}\left(D^{h}\right)}<M,\\ \left\|P^{h}-P^{h\ast}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}<\delta.\end{% array}\right.\text{ }start_ARRAY start_ROW start_CELL italic_P start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , 0 ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_P start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_B ) = italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL italic_P start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , 0 ) = italic_G start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , italic_P start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_B ) = italic_G start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_r , italic_φ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) , end_CELL end_ROW start_ROW start_CELL ∥ italic_P start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_M , ∥ italic_P start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_M , end_CELL end_ROW start_ROW start_CELL ∥ italic_P start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_P start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_δ . end_CELL end_ROW end_ARRAY

Theorem 8.2 (an estimate of the accuracy of the regularized solution). Assume that conditions of Theorem 8.1 hold. Let the exact solution VhB(M),superscript𝑉superscript𝐵𝑀V^{h\ast}\in B^{\ast}\left(M\right),italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∈ italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M ) ,where the set B(M)superscript𝐵𝑀B^{\ast}\left(M\right)italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M ) is defined in (62). Assume that conditions (63) hold. Furthermore, assume that VhHN1,h(Ωh)<Mα,subscriptnormsuperscript𝑉superscriptsubscript𝐻𝑁1superscriptΩ𝑀𝛼\left\|V^{h\ast}\right\|_{H_{N}^{1,h}\left(\Omega^{h}\right)}<M-\alpha,∥ italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( roman_Ω start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT < italic_M - italic_α , where the number α(0,M)𝛼0𝑀\alpha\in\left(0,M\right)italic_α ∈ ( 0 , italic_M ) is so small that

(64) α<Cδ.𝛼𝐶𝛿\alpha<C\delta.italic_α < italic_C italic_δ .

Let λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT be the number λ0subscript𝜆0\lambda_{0}italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT of Theorem 8.1 in the case when the number M𝑀Mitalic_M in (61) is replaced with 2M,2𝑀2M,2 italic_M ,

(65) λ1=λ0(h0,c,2M,N,Dh)λ0(h0,c,M,N,Dh).subscript𝜆1subscript𝜆0subscript0𝑐2𝑀𝑁superscript𝐷subscript𝜆0subscript0𝑐𝑀𝑁superscript𝐷\lambda_{1}=\lambda_{0}\left(h_{0},c,2M,N,D^{h}\right)\geq\lambda_{0}\left(h_{% 0},c,M,N,D^{h}\right).italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_c , 2 italic_M , italic_N , italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ≥ italic_λ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_h start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_c , italic_M , italic_N , italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) .

For any λλ1𝜆subscript𝜆1\lambda\geq\lambda_{1}italic_λ ≥ italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT let Vmin,λhB(M)¯superscriptsubscript𝑉𝜆¯𝐵𝑀V_{\min,\lambda}^{h}\in\overline{B\left(M\right)}italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ over¯ start_ARG italic_B ( italic_M ) end_ARG be the minimizer on of the functional Jλ(Vh)subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right)italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) on the set B(M)¯,¯𝐵𝑀\overline{B\left(M\right)},over¯ start_ARG italic_B ( italic_M ) end_ARG , which is found in Theorem 8.1. Then the following accuracy estimate holds:

Vmin,λhVhHN1,h(Dh)Cδ.subscriptnormsuperscriptsubscript𝑉𝜆superscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷𝐶𝛿\left\|V_{\min,\lambda}^{h}-V^{h\ast}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}% \leq C\delta.∥ italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT ≤ italic_C italic_δ .

For λλ1𝜆subscript𝜆1\lambda\geq\lambda_{1}italic_λ ≥ italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT we now construct the gradient descent method of the minimization of the functional Jλ(Vh).subscript𝐽𝜆superscript𝑉J_{\lambda}\left(V^{h}\right).italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) . Let γ(0,1)𝛾01\gamma\in\left(0,1\right)italic_γ ∈ ( 0 , 1 ) be a number and let V0hB(M/3)superscriptsubscript𝑉0𝐵𝑀3V_{0}^{h}\in B\left(M/3\right)italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∈ italic_B ( italic_M / 3 ) be an arbitrary point. The gradient descent method is constructed as the following sequence:

(66) Vnh=Vn1hγJλ(Vn1h), n=1,2,formulae-sequencesuperscriptsubscript𝑉𝑛superscriptsubscript𝑉𝑛1𝛾superscriptsubscript𝐽𝜆superscriptsubscript𝑉𝑛1 𝑛12V_{n}^{h}=V_{n-1}^{h}-\gamma J_{\lambda}^{\prime}\left(V_{n-1}^{h}\right),% \text{ }n=1,2,...italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT = italic_V start_POSTSUBSCRIPT italic_n - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_γ italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_V start_POSTSUBSCRIPT italic_n - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) , italic_n = 1 , 2 , …

Since by Theorem 8.1 Jλ1(Wn1h)H0,N1,h(Ω¯h),superscriptsubscript𝐽subscript𝜆1superscriptsubscript𝑊𝑛1superscriptsubscript𝐻0𝑁1superscript¯ΩJ_{\lambda_{1}}^{\prime}\left(W_{n-1}^{h}\right)\in H_{0,N}^{1,h}\left(% \overline{\Omega}^{h}\right),italic_J start_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_W start_POSTSUBSCRIPT italic_n - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ∈ italic_H start_POSTSUBSCRIPT 0 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( over¯ start_ARG roman_Ω end_ARG start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) , nfor-all𝑛\forall n∀ italic_n, then all vector functions Vnhsuperscriptsubscript𝑉𝑛V_{n}^{h}italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT satisfy the same boundary condition as the ones in the third line of (59), see (58).

Theorem 8.3. Let Cδ(0,M/3)𝐶𝛿0𝑀3C\delta\in\left(0,M/3\right)italic_C italic_δ ∈ ( 0 , italic_M / 3 ) in (64) and let the number β(Cδ,M/3).𝛽𝐶𝛿𝑀3\beta\in\left(C\delta,M/3\right).italic_β ∈ ( italic_C italic_δ , italic_M / 3 ) . Suppose that the exact solution VhB(M/3β)superscript𝑉superscript𝐵𝑀3𝛽V^{h\ast}\in B^{\ast}\left(M/3-\beta\right)italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∈ italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( italic_M / 3 - italic_β ). Let λ=λ1𝜆subscript𝜆1\lambda=\lambda_{1}italic_λ = italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT where λ1subscript𝜆1\lambda_{1}italic_λ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is defined in (65). Then there exists a sufficiently small number γ0(0,1)subscript𝛾001\gamma_{0}\in\left(0,1\right)italic_γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∈ ( 0 , 1 ) such that for any γ(0,γ0)𝛾0subscript𝛾0\gamma\in\left(0,\gamma_{0}\right)italic_γ ∈ ( 0 , italic_γ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) the sequence ( 66) {Vnh}n=0B(M).superscriptsubscriptsuperscriptsubscript𝑉𝑛𝑛0𝐵𝑀\left\{V_{n}^{h}\right\}_{n=0}^{\infty}\subset B\left(M\right).{ italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT } start_POSTSUBSCRIPT italic_n = 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT ⊂ italic_B ( italic_M ) . In addition, there exists a number ρ=ρ(γ)(0,1)𝜌𝜌𝛾01\rho=\rho\left(\gamma\right)\in\left(0,1\right)italic_ρ = italic_ρ ( italic_γ ) ∈ ( 0 , 1 ) such that the following convergence estimate holds

VnhVhHN1,h(Dh)Cδ+ρnV0hVmin,λhHN1,h(Dh).subscriptnormsuperscriptsubscript𝑉𝑛superscript𝑉superscriptsubscript𝐻𝑁1superscript𝐷𝐶𝛿superscript𝜌𝑛subscriptnormsuperscriptsubscript𝑉0superscriptsubscript𝑉𝜆superscriptsubscript𝐻𝑁1superscript𝐷\left\|V_{n}^{h}-V^{h\ast}\right\|_{H_{N}^{1,h}\left(D^{h}\right)}\leq C\delta% +\rho^{n}\left\|V_{0}^{h}-V_{\min,\lambda}^{h}\right\|_{H_{N}^{1,h}\left(D^{h}% \right)}.∥ italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT ≤ italic_C italic_δ + italic_ρ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∥ italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT .

Furthermore, let the unction [nh(𝐱h)]2superscriptdelimited-[]superscript𝑛superscript𝐱2\left[n^{h\ast}\left(\mathbf{x}^{h}\right)\right]^{2}[ italic_n start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT be the exact semi-discrete target function and let [nnh(𝐱h)]2superscriptdelimited-[]superscriptsubscript𝑛𝑛superscript𝐱2\left[n_{n}^{h}\left(\mathbf{x}^{h}\right)\right]^{2}[ italic_n start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT be the semi-discrete target function, which is found via the substitution of Vnhsuperscriptsubscript𝑉𝑛V_{n}^{h}italic_V start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT , sequentially, first in the semi-discrete analog of the first line of (37), then in (36) and finally in the semi-discrete analog of the left hand side of (31). Then the following convergence rate holds:

[nnh(𝐱h)]2[nh(𝐱h)]2L2,Nh(Dh)Cδ+ρnV0hVmin,λhHN1,h(Dh).subscriptnormsuperscriptdelimited-[]superscriptsubscript𝑛𝑛superscript𝐱2superscriptdelimited-[]superscript𝑛superscript𝐱2superscriptsubscript𝐿2𝑁superscript𝐷𝐶𝛿superscript𝜌𝑛subscriptnormsuperscriptsubscript𝑉0superscriptsubscript𝑉𝜆superscriptsubscript𝐻𝑁1superscript𝐷\left\|\left[n_{n}^{h}\left(\mathbf{x}^{h}\right)\right]^{2}-\left[n^{h\ast}% \left(\mathbf{x}^{h}\right)\right]^{2}\right\|_{L_{2,N}^{h}\left(D^{h}\right)}% \leq C\delta+\rho^{n}\left\|V_{0}^{h}-V_{\min,\lambda}^{h}\right\|_{H_{N}^{1,h% }\left(D^{h}\right)}.∥ [ italic_n start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - [ italic_n start_POSTSUPERSCRIPT italic_h ∗ end_POSTSUPERSCRIPT ( bold_x start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 , italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT ≤ italic_C italic_δ + italic_ρ start_POSTSUPERSCRIPT italic_n end_POSTSUPERSCRIPT ∥ italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT - italic_V start_POSTSUBSCRIPT roman_min , italic_λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ∥ start_POSTSUBSCRIPT italic_H start_POSTSUBSCRIPT italic_N end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 1 , italic_h end_POSTSUPERSCRIPT ( italic_D start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) end_POSTSUBSCRIPT .

Since a smallness assumption is not imposed on the number M𝑀Mitalic_M and since the starting point V0hsuperscriptsubscript𝑉0V_{0}^{h}italic_V start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT of the sequence (66) is an arbitrary point of B(M/3),𝐵𝑀3B\left(M/3\right),italic_B ( italic_M / 3 ) , then Definition of section 1 implies that Theorem 8.3 ensures the global convergence of the gradient descent method (66).

9 Numerical Studies

In this section we describe our numerical studies. We specify parameters of the domains ΩΩ\Omegaroman_Ω and ΩεsubscriptΩ𝜀\Omega_{\varepsilon}roman_Ω start_POSTSUBSCRIPT italic_ε end_POSTSUBSCRIPT in (3) and (4) as:

R=1,ε=0.01,B=1.formulae-sequence𝑅1formulae-sequence𝜀0.01𝐵1R=1,\ \varepsilon=0.01,\ B=1.italic_R = 1 , italic_ε = 0.01 , italic_B = 1 .

The refractive index n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) in eikonal equation (2) is taken as

(67) n(𝐱)={ca=const.>1,inside of the tested inclusion,1,outside of the tested inclusion.n(\mathbf{x})=\left\{\begin{array}[]{cc}c_{a}=const.>1,&\text{inside of the % tested inclusion,}\\ 1,&\text{outside of the tested inclusion.}\end{array}\right.italic_n ( bold_x ) = { start_ARRAY start_ROW start_CELL italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = italic_c italic_o italic_n italic_s italic_t . > 1 , end_CELL start_CELL inside of the tested inclusion, end_CELL end_ROW start_ROW start_CELL 1 , end_CELL start_CELL outside of the tested inclusion. end_CELL end_ROW end_ARRAY

To ensure that n(𝐱)C1(Ω¯)𝑛𝐱superscript𝐶1¯Ωn(\mathbf{x})\in C^{1}\left(\overline{\Omega}\right)italic_n ( bold_x ) ∈ italic_C start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT ( over¯ start_ARG roman_Ω end_ARG ) as in (10), we smooth out n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) near the boundaries of our tested inclusions. Then we set:

(68) correct inclusion/background contrast=ca1,correct inclusion/background contrastsubscript𝑐𝑎1\displaystyle\hskip 28.45274pt\mbox{correct inclusion/background contrast}=% \frac{c_{a}}{1},correct inclusion/background contrast = divide start_ARG italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT end_ARG start_ARG 1 end_ARG ,
(69) computed inclusion/background contrast=maxinclusion(ncomp(𝐱))1.computed inclusion/background contrastsubscriptinclusionsubscript𝑛comp𝐱1\displaystyle\text{computed inclusion/background contrast}=\frac{\max_{\text{% inclusion}}\left(n_{\text{comp}}(\mathbf{x})\right)}{1}.computed inclusion/background contrast = divide start_ARG roman_max start_POSTSUBSCRIPT inclusion end_POSTSUBSCRIPT ( italic_n start_POSTSUBSCRIPT comp end_POSTSUBSCRIPT ( bold_x ) ) end_ARG start_ARG 1 end_ARG .

In the numerical tests below, we take ca=1.5,3,5subscript𝑐𝑎1.535c_{a}=1.5,3,5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 , 3 , 5 which means 1.5:1, 3:1 and 5:1 inclusion background contrasts respectively, see (68). To demonstrate that our numerical method can work with inclusions, which have sophisticated non-convex shapes with voids in them, we take shapes inclusions like letters ‘B𝐵Bitalic_B’ and ‘O𝑂Oitalic_O’, This is similar with our previous works [13]-[16] on the convexification.

After choosing the function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ), we use the fast marching toolbox ”Toolbox Fast Marching” [20] in MATLAB to solve eikonal equation (2). Then we convert data from Cartesian coordinates to polar coordinates to gain the observation data in (14)-(16). Then we solve the Minimization Problem (60) to gain the computed solution ncomp(r,φ,z)subscript𝑛comp𝑟𝜑𝑧n_{\text{comp}}(r,\varphi,z)italic_n start_POSTSUBSCRIPT comp end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z ). Finally, we convert the reconstructed solution ncomp(r,φ,z)subscript𝑛comp𝑟𝜑𝑧n_{\text{comp}}(r,\varphi,z)italic_n start_POSTSUBSCRIPT comp end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z ) from polar coordinates to Cartesian coordinates to exhibit results.

To solve eikonal equation (2) for generating the observation data in (14)-(16), we choose hz=hφ=1/40subscript𝑧subscript𝜑140h_{z}=h_{\varphi}=1/40italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = 1 / 40 in (48) and (49), as well as we choose hr=1/40subscript𝑟140h_{r}=1/40italic_h start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = 1 / 40 to generate discrete points along the rlimit-from𝑟r-italic_r -direction. Because r(ε,R)𝑟𝜀𝑅r\in\left(\varepsilon,R\right)italic_r ∈ ( italic_ε , italic_R ) rather than r(0,R)𝑟0𝑅r\in\left(0,R\right)italic_r ∈ ( 0 , italic_R ) in (8), then the first interval along rlimit-from𝑟r-italic_r -direction is hrεsubscript𝑟𝜀h_{r}-\varepsilonitalic_h start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT - italic_ε. To solve the Minimization Problem (60), we choose hz=hφ=hr=1/20subscript𝑧subscript𝜑subscript𝑟120h_{z}=h_{\varphi}=h_{r}=1/20italic_h start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_h start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT = italic_h start_POSTSUBSCRIPT italic_r end_POSTSUBSCRIPT = 1 / 20.

To guarantee that the solution of the Minimization Problem (60) satisfies the boundary conditions in (54)-(56), we adopt the Matlab’s built-in optimization toolbox fmincon. The iterations of fmincon stop when we get

(70) |Jλ(Vh)|<102.subscript𝐽𝜆superscript𝑉superscript102|\nabla J_{\lambda}\left(V^{h}\right)|<10^{-2}.| ∇ italic_J start_POSTSUBSCRIPT italic_λ end_POSTSUBSCRIPT ( italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT ) | < 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT .

The starting point of iterations of fmincon is chosen as Vh=0superscript𝑉0V^{h}=0italic_V start_POSTSUPERSCRIPT italic_h end_POSTSUPERSCRIPT = 0. Although the starting point does not satisfy the boundary conditions in (54)-(56), fmincon makes sure that the boundary conditions (54)-(56) are satisfied on all other iterations of fmincon.

We consider the random noise in observation data in (14)-(16) as follows:

(71) pξ(φ,z,z0)superscript𝑝𝜉𝜑𝑧subscript𝑧0\displaystyle p^{\xi}\left(\varphi,z,z_{0}\right)italic_p start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =p(φ,z,z0)(1+δξp(φ,z,z0)),absent𝑝𝜑𝑧subscript𝑧01𝛿subscript𝜉𝑝𝜑𝑧subscript𝑧0\displaystyle=p\left(\varphi,z,z_{0}\right)\left(1+\delta\xi_{p}\left(\varphi,% z,z_{0}\right)\right),= italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ( 1 + italic_δ italic_ξ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) ,
(72) p0ξ(r,φ,z0)superscriptsubscript𝑝0𝜉𝑟𝜑subscript𝑧0\displaystyle p_{0}^{\xi}\left(r,\varphi,z_{0}\right)italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =p0(r,φ,z0)(1+δξ0(r,φ,z0)),absentsubscript𝑝0𝑟𝜑subscript𝑧01𝛿subscript𝜉0𝑟𝜑subscript𝑧0\displaystyle=p_{0}\left(r,\varphi,z_{0}\right)\left(1+\delta\xi_{0}\left(r,% \varphi,z_{0}\right)\right),= italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ( 1 + italic_δ italic_ξ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) ,
(73) pBξ(r,φ,z0)superscriptsubscript𝑝𝐵𝜉𝑟𝜑subscript𝑧0\displaystyle p_{B}^{\xi}\left(r,\varphi,z_{0}\right)italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) =pB(r,φ,z0)(1+δξB(r,φ,z0)).absentsubscript𝑝𝐵𝑟𝜑subscript𝑧01𝛿subscript𝜉𝐵𝑟𝜑subscript𝑧0\displaystyle=p_{B}\left(r,\varphi,z_{0}\right)\left(1+\delta\xi_{B}\left(r,% \varphi,z_{0}\right)\right).= italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ( 1 + italic_δ italic_ξ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ) .

In (71)-(73), ξp(φ,z,z0)subscript𝜉𝑝𝜑𝑧subscript𝑧0\xi_{p}\left(\varphi,z,z_{0}\right)italic_ξ start_POSTSUBSCRIPT italic_p end_POSTSUBSCRIPT ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is the Gaussian random variable depending on variables φ,z,z0𝜑𝑧subscript𝑧0\varphi,z,z_{0}italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and ξ0(r,φ,z0),ξB(r,φ,z0)subscript𝜉0𝑟𝜑subscript𝑧0subscript𝜉𝐵𝑟𝜑subscript𝑧0\xi_{0}\left(r,\varphi,z_{0}\right),\xi_{B}\left(r,\varphi,z_{0}\right)italic_ξ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) , italic_ξ start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) are the uniformly distributed random variables in the interval [1,1]11[-1,1][ - 1 , 1 ] depending on variables r,φ,z0𝑟𝜑subscript𝑧0r,\varphi,z_{0}italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Also, δ=0.01𝛿0.01\delta=0.01italic_δ = 0.01 and δ=0.03𝛿0.03\delta=0.03italic_δ = 0.03, which correspond to the 1%percent11\%1 % and 3%percent33\%3 % random noise levels respectively. Recall that, to simplify the presentation, we did not consider the noise in the function p(φ,z,z0)𝑝𝜑𝑧subscript𝑧0p\left(\varphi,z,z_{0}\right)italic_p ( italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in our theoretical derivations, see (63). Nevertheless, we still introduce random noise in this function in our numerical studies, see (71).

To calculate the φlimit-from𝜑\varphi-italic_φ -derivative and the zlimit-from𝑧z-italic_z -derivative of the noisy function pξ(φp^{\xi}(\varphiitalic_p start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_φ, z𝑧zitalic_z, z0)z_{0})italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in (29)-(30), as well as the rlimit-from𝑟r-italic_r -derivative of noisy functions p0ξ(r,φ,z0)superscriptsubscript𝑝0𝜉𝑟𝜑subscript𝑧0p_{0}^{\xi}(r,\varphi,z_{0})italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ), pBξ(r,φ,z0)superscriptsubscript𝑝𝐵𝜉𝑟𝜑subscript𝑧0p_{B}^{\xi}(r,\varphi,z_{0})italic_p start_POSTSUBSCRIPT italic_B end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_ξ end_POSTSUPERSCRIPT ( italic_r , italic_φ , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in (33), we firstly use the natural cubic splines to approximate the noisy data in (71)-(73). Then we use the derivatives of those splines to approximate the derivatives of corresponding noisy data.

The parameters λ𝜆\lambdaitalic_λ and N𝑁Nitalic_N are the two key parameters in our numerical method. We firstly choose the optimal λ𝜆\lambdaitalic_λ, when N𝑁Nitalic_N is large enough with N=8𝑁8N=8italic_N = 8. Then, keeping that optimal λ,𝜆\lambda,italic_λ , we find the optimal value of N𝑁Nitalic_N. This is done in Test 1. An important point to make here that, once chosen in Test 1, that optimal pair (λ,N)𝜆𝑁\left(\lambda,N\right)( italic_λ , italic_N ) is kept for all other tests.

Test 1. We test the case when the inclusion in (67) has the shape of the horizontally oriented letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. The goal of this test is to find the optimal values of the parameters λ𝜆\lambdaitalic_λ and N𝑁Nitalic_N. To make sure that the chosen parameter N𝑁Nitalic_N does not impact our results with different values of λ𝜆\lambdaitalic_λ, we set N𝑁Nitalic_N to be large enough, i.e. N=8𝑁8N=8italic_N = 8. Computational results are displayed in Figure 2. We observe that the images have a low quality for λ=0,1,2𝜆012\lambda=0,1,2italic_λ = 0 , 1 , 2. Then the quality is improved with λ=3,4𝜆34\lambda=3,4italic_λ = 3 , 4, and the reconstruction quality significantly deteriorates at λ=10𝜆10\lambda=10italic_λ = 10. Hence, we select λ=3𝜆3\lambda=3italic_λ = 3 as the optimal one.

Refer to caption
Figure 2: Test 1: The reconstructions of n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) of different λ𝜆\lambdaitalic_λ with N=8𝑁8N=8italic_N = 8, when the shape of the inclusion in (67) is horizontal letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5. The value of λ𝜆\lambdaitalic_λ can be seen on the top side of each square. The images have a low quality for λ=0,1,2𝜆012\lambda=0,1,2italic_λ = 0 , 1 , 2. The quality is improved with λ=3,4,5𝜆345\lambda=3,4,5italic_λ = 3 , 4 , 5. Then it starts to deteriorate at λ=6𝜆6\lambda=6italic_λ = 6 and becomes unsatisfactory at λ=10𝜆10\lambda=10italic_λ = 10. Thus, we select λ=3𝜆3\lambda=3italic_λ = 3 as the optimal value of the parameter λ𝜆\lambdaitalic_λ.

Next, we fix the optimal value of λ=3𝜆3\lambda=3italic_λ = 3 and test the influence of the parameter N𝑁Nitalic_N. The results with N=2,4,8𝑁248N=2,4,8italic_N = 2 , 4 , 8 are shown in Figure 3. We can find that the reconstruction for N=2𝑁2N=2italic_N = 2 is of a low quality, and the reconstructions for N=4,8𝑁48N=4,8italic_N = 4 , 8 are almost same. Hence, to reduce the computational cost we select N=4𝑁4N=4italic_N = 4 as the optimal value of N𝑁Nitalic_N, instead of N=8𝑁8N=8italic_N = 8. Furthermore, we also consider the truncated expansion of the function u(r,φ,z,z0)𝑢𝑟𝜑𝑧subscript𝑧0u\left(r,\varphi,z,z_{0}\right)italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) in (26) with N=4𝑁4N=4italic_N = 4, which is denoted as u4(r,φ,z,z0)subscript𝑢4𝑟𝜑𝑧subscript𝑧0u_{4}\left(r,\varphi,z,z_{0}\right)italic_u start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ). We have obtained that

(74) u4(r,φ,z,z0)L2(Ω×(0,B))u(r,φ,z,z0)L2(Ω×(0,B))=99.86%.subscriptnormsubscript𝑢4𝑟𝜑𝑧subscript𝑧0subscript𝐿2Ω0𝐵subscriptnorm𝑢𝑟𝜑𝑧subscript𝑧0subscript𝐿2Ω0𝐵percent99.86\frac{\left\|u_{4}\left(r,\varphi,z,z_{0}\right)\right\|_{L_{2}\left(\Omega% \times\left(0,B\right)\right)}}{\left\|u\left(r,\varphi,z,z_{0}\right)\right\|% _{L_{2}\left(\Omega\times\left(0,B\right)\right)}}=99.86\%.divide start_ARG ∥ italic_u start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( roman_Ω × ( 0 , italic_B ) ) end_POSTSUBSCRIPT end_ARG start_ARG ∥ italic_u ( italic_r , italic_φ , italic_z , italic_z start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) ∥ start_POSTSUBSCRIPT italic_L start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( roman_Ω × ( 0 , italic_B ) ) end_POSTSUBSCRIPT end_ARG = 99.86 % .

It is clear from (74) that N=4𝑁4N=4italic_N = 4 is a quite informative case.

Conclusion: We choose

(75) λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4

as the optimal values of these parameters, and we use these values in Tests 2-5.

Refer to caption
Figure 3: Test 1: Reconstructions of n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) of N=2,4,8𝑁248N=2,4,8italic_N = 2 , 4 , 8 at λ=3𝜆3\lambda=3italic_λ = 3, when the shape of the inclusion in (67) is horizontal letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5. The reconstruction for N=2𝑁2N=2italic_N = 2 has a low quality, and the reconstructions for N=4,8𝑁48N=4,8italic_N = 4 , 8 are almost same. Thus, keeping also in mind (74) and (75) and also to reduce the computational cost, we choose in Tests 2-5 N=4𝑁4N=4italic_N = 4 and λ=3𝜆3\lambda=3italic_λ = 3.

Test 2. We test the case when the inclusion in (67) has the shape of the vertically oriented letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Results are presented on Figure 4. An accurate reconstruction can be observed.

Refer to caption
Figure 4: Test 2: The exact (left) and reconstructed (right) function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ), when the shape of the inclusion in (67) is vertically oriented letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Here λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4 as in (75). The reconstruction is accurate.

Test 3. We test the cases when the inclusion in (67) has the shape of the vertically oriented letter ‘B𝐵Bitalic_B’ with ca=3subscript𝑐𝑎3c_{a}=3italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 3 and ca=5subscript𝑐𝑎5c_{a}=5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 5 in it. Results are displayed on Figures 5 and 6 respectively. The inclusion/background contrasts in (68) are respectively 3:1:313:13 : 1 and 5:1:515:15 : 1. We see that shape of the inclusion is imaged accurately in both cases. In addition, the computed inclusion/background contrasts in (69) are accurate.

Refer to caption
Figure 5: Test 3: The exact (left) and reconstructed (right) function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ), when the shape of the inclusion in (67) is vertically oriented letter ‘B𝐵Bitalic_B’ with ca=3subscript𝑐𝑎3c_{a}=3italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 3 in it. Here λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4 as in (75). The inclusion/background contrast in (68) is 3:1:313:13 : 1. The computed inclusion/background contrast in (69) is accurate.
Refer to caption
Figure 6: Test 3: The exact (left) and reconstructed (right) function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ), when the shape of the inclusion in (67) is vertically oriented letter ‘B𝐵Bitalic_B’ with ca=5subscript𝑐𝑎5c_{a}=5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 5 in it. Here λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4 as in (75). The inclusion/background contrast in ( 68) is 5:1:515:15 : 1. The computed inclusion/background contrast in ( 69) is accurate.

Test 4. We test the case when the inclusion in (67) has the shape of the letter ‘O𝑂Oitalic_O’ elongated along the ylimit-from𝑦y-italic_y -axis and with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Results are presented on Figure 7. We again observe an accurate reconstruction of both the shape of the inclusion and the inclusion/background contrast.

Refer to caption
Figure 7: Test 4: The exact (left) and reconstructed (right) function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ), when the shape of the inclusion in (67) is the letter ‘O𝑂Oitalic_O’ elongated along the y-axis and with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Here λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4 as in (75). The reconstruction is accurate.

Test 5. We consider the case when the random noise is present in the data in (71)-(73) with δ=0.01𝛿0.01\delta=0.01italic_δ = 0.01 and δ=0.03𝛿0.03\delta=0.03italic_δ = 0.03, i.e. with 1% and 3% noise level respectively. We test the reconstruction for the case when the inclusion in (67) has the shape of the vertically oriented letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Results are displayed on Figure 8. The reconstructions of the shape of the inclusion as well as computed inclusion/background contrasts in (69) are accurate.

Refer to caption
Figure 8: Test 5: Reconstructed function n(𝐱)𝑛𝐱n(\mathbf{x})italic_n ( bold_x ) with σ=0.01𝜎0.01\sigma=0.01italic_σ = 0.01 (top) and σ=0.03𝜎0.03\sigma=0.03italic_σ = 0.03 (bottom), i.e. with 1% and 3% noise level, when the shape of the inclusion in (67) is vertically oriented letter ‘B𝐵Bitalic_B’ with ca=1.5subscript𝑐𝑎1.5c_{a}=1.5italic_c start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 1.5 in it. Here λ=3,N=4formulae-sequence𝜆3𝑁4\lambda=3,N=4italic_λ = 3 , italic_N = 4 as in (75). The reconstructions and computed inclusion/background contrasts in (69) are accurate.

Remarks 9.1:

  1. 1.

    Recall that by Theorem 3.1 condition (11) is a sufficient condition for our method to work. Recall also that in the data generation process we smooth out tested inclusions in small neighborhoods of their boundaries. Given these, a careful analysis of correct images of Figures 2-8 indicates that condition (11) is satisfied at least in a major part of the domain ΩΩ\Omegaroman_Ω in each of the above Tests 1-5.

  2. 2.

    We conclude, therefore, that our method works numerically under conditions which are broader than (11). In other words, (11) is not a necessary condition for our method to work. Establishing necessary conditions is outside of the scope of the current publication.

Acknowledgments. The work of Li was partially supported by Guangdong Basic and Applied Basic Research Foundation 2023B1515250005. The work of Romanov was partially supported by the Mathematical Center in Akademgorodok under Agreement 075152022281 with the Ministry of Science and Higher Education of the Russian Federation. The work of Yang was partially supported by Supercomputing Center of Lanzhou University.

References

  • [1] L. Baudouin, M. de Buhan, and S. Ervedoza, Convergent algorithm based on Carleman estimates for the recovery of a potential in the wave equation, SIAM J. Numer. Anal., 55 (2017), pp. 1578–1613.
  • [2] L. Baudouin, M. de Buhan, S. Ervedoza, and A. Osses, Carleman-based reconstruction algorithm for the waves, SIAM J. Numer. Anal., 59 (2021), pp. 998–1039.
  • [3] L. Beilina and E. Lindström, An adaptive finite element/finite difference domain decomposition method for applications in microwave imaging, Electronics, 11 (2022), p. 1359.
  • [4] A. V. Goncharsky, S. Y. Romanov, and S. Y. Seryozhnikov, On mathematical problems of two-coefficient inverse problems of ultrasonic tomography, Inverse Probl., 40 (2024), p. 045026.
  • [5] M. J. Grote and U. Nahum, Adaptive eigenspace for multi-parameter inverse scattering problems, Computers & Mathematics with Applications, 77 (2019), pp. 3264–3280.
  • [6] J. P. Guillement and R. G. Novikov, Inversion of weighted Radon transforms via finite Fourier series weight approximation, Inverse Probl. Sci. En., 22 (2013), pp. 787–802.
  • [7] G. Herglotz, Űber die Elastizitaet der Erde bei Beruecksichtigung ihrer variablen Dichte, Zeitschr. fur Math. Phys., 52 (1905), pp. 275–299.
  • [8] S. I. Kabanikhin, K. K. Sabelfeld, N. S. Novikov, and M. A. Shishlenin, Numerical solution of the multidimensional Gelfand-Levitan equation, J. Inverse Ill-Posed Probl., 23 (2015), pp. 439–450.
  • [9] M. V. Klibanov and O. V. Ioussoupova, Uniform strict convexity of a cost functional for three-dimensional inverse scattering problem, SIAM J. Math. Anal, 26 (1995), pp. 147–179.
  • [10] M. V. Klibanov, Global convexity in a three-dimensional inverse acoustic problem, SIAM J. Math. Anal., 28 (1997), pp. 1371–1388.
  • [11] M. V. Klibanov, Convexification of restricted Dirichlet to Neumann map, J. Inverse Ill-Posed Probl., 25 (2017), pp. 669–685.
  • [12] M. V. Klibanov, Travel time tomography with formally determined incomplete data in 3D, Inverse Problems and Imaging, 13 (2019), pp. 1367–1393.
  • [13] M. V. Klibanov, J. Li, and W. Zhang, Convexification for an inverse parabolic problem, Inverse Probl., 36 (2020), p. 085008.
  • [14] M. V. Klibanov and J. Li, Inverse Problems and Carleman Estimates: Global Uniqueness, Global Convergence and Experimental Data, De Gruyter, Berlin, 2021.
  • [15] M. V. Klibanov, J. Li, and W. Zhang, Numerical solution of the 3-D travel time tomography problem, J. Computational Physics, 476 (2023), p. 111910.
  • [16] M. V. Klibanov, J. Li, and Z. Yang, Convexification numerical method for the retrospective problem of mean field games, Applied Mathematics and Optimization, (2024), 90:6.
  • [17] R. G. Mukhometov, The reconstruction problem of a two-dimensional Riemannian metric and integral geometry, Soviet Math. Dokl., 18 (1977), pp. 32–35.
  • [18] R. G. Mukhometov and V. G. Romanov, On the problem of determining an isotropic Riemannian metric in the n𝑛nitalic_n-dimensional space, Dokl. Acad. Sci. USSR, 19 (1978), pp. 1330–1333.
  • [19] L. Pestov and G. Uhlmann, Two dimensional simple Riemannian manifolds are boundary distance rigid, Annals of Mathematics, 161 (2005), pp. 1093–1110.
  • [20] G. Peyre, Toolbox fast marching, MATLAB Central File Exchange, (2023).
  • [21] V. G. Romanov, Inverse Problems of Mathematical Physics, VNU Press, Utrecht, The Netherlands, 1986.
  • [22] U. Schrőder and T. Schuster, An iterative method to reconstruct the refractive index of a medium from time-off-light measurements, Inverse Problems, 32 (2016), p. 085009.
  • [23] A. N. Tikhonov, A. V. Goncharsky, V. V. Stepanov, and A. G. Yagola,  Numerical methods for the solution of Ill-posed problems, Kluwer, London, 1995.
  • [24] L. Volgyesi and M. Moser, The inner structure of the Earth, Periodica Polytechnica Chemical Engineering, 26 (1982), pp. 155–204.
  • [25] E. Wiechert and K. Zoeppritz, Uber Erdbebenwellen, Nachr. Koenigl. Geselschaft Wiss. Gottingen, (1907), pp. 415–549.
  • [26] H. Zhao and Y. Zhong, A hybrid adaptive phase space method for reflection traveltime tomography, SIAM J. Imaging Sciences, 12 (2019), pp. 28–53.