Next Article in Journal
The Allocation of Base Stations with Region Clustering and Single-Objective Nonlinear Optimization
Next Article in Special Issue
Failure Transition and Validity of Brazilian Disc Test under Different Loading Configurations: A Numerical Study
Previous Article in Journal
Parametric Frailty Analysis in Presence of Collinearity: An Application to Assessment of Infant Mortality
Previous Article in Special Issue
epSFEM: A Julia-Based Software Package of Parallel Incremental Smoothed Finite Element Method (S-FEM) for Elastic-Plastic Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving Inverse Conductivity Problems in Doubly Connected Domains by the Homogenization Functions of Two Parameters

1
Nanjing Hydraulic Research Institute, Nanjing 210029, China
2
College of Mechanics and Materials, Hohai University, Nanjing 210098, China
3
Center of Excellence for Ocean Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
4
Department of Mathematics, University of Southern Mississippi, Hattiesburg, MS 39406, USA
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(13), 2256; https://doi.org/10.3390/math10132256
Submission received: 7 April 2022 / Revised: 19 June 2022 / Accepted: 24 June 2022 / Published: 27 June 2022
(This article belongs to the Special Issue Computational Methods and Applications for Numerical Analysis)

Abstract

:
In the paper, we make the first attempt to derive a family of two-parameter homogenization functions in the doubly connected domain, which is then applied as the bases of trial solutions for the inverse conductivity problems. The expansion coefficients are obtained by imposing an extra boundary condition on the inner boundary, which results in a linear system for the interpolation of the solution in a weighted Sobolev space. Then, we retrieve the spatial- or temperature-dependent conductivity function by solving a linear system, which is obtained from the collocation method applied to the nonlinear elliptic equation after inserting the solution. Although the required data are quite economical, very accurate solutions of the space-dependent and temperature-dependent conductivity functions, the Robin coefficient function and also the source function are available. It is significant that the nonlinear inverse problems can be solved directly without iterations and solving nonlinear equations. The proposed method can achieve accurate results with high efficiency even for large noise being imposed on the input data.

1. Introduction

In recent decades, a large number of inverse problems of the nonlinear elliptic-type partial differential equation (PDE) have been well investigated, involving the inverse source problem, inverse conductivity problem as well as inverse Robin problem, which arise in several branches of applications in science and engineering. Analytical solutions to inverse problems are difficult to obtain since some information is missing, such as the boundary conditions or sources compared with the forward problems. Therefore, many numerical approaches have been developed to resolve inverse problems [1]. In the linear elliptic type PDEs, for identifying unknown sources, the regularization methods were advocated in [2,3]. Klose [4] solved an inverse source problem based on the radiative transfer equation arising in optical molecular imaging. In Ref. [5], Hon et al. applied Green’s function for the inverse source identification. Then Li et al. [6] proposed the modified regularization method on the Poisson equation for determining an unknown source. Ahmadabadi and co-workers proposed the method of fundamental solutions for the inverse space-dependent heat source problems by using a new transformation [7]. The source function for a seawater intrusion problem in an unconfined aquifer has been studied by Slimani [8]. The inverse source problems were examined by Alahyane et al. [9] using the regularized optimal control method. Some new regularization methods were proposed for inverse source problems governed by fractional PDEs [10,11]. Nguyen [12] investigated the inverse source problems of the fractional diffusion equations based on the Tikhonov regularization method. Recently, Liu [13] have proposed a new procedure of boundary functions, which preserves the energy identity to identify the sources of 2D elliptic-type nonlinear PDEs. However, the methods proposed in [13] required extra boundary conditions of source function on a rectangle. We will extend the work to any 2D nonlinear elliptic equation without using extra boundary data of the source function in the doubly connected domain.
On the other hand, linear and nonlinear inverse conductivity problems have been studied by many authors. Kwon considered the anisotropic inverse conductivity and scattering problems [14]. The inverse problem of time-dependent thermal conductivity was studied by Huntul and Lesnic by recasting the original problems into the nonlinear least-squares minimization [15]. Isakov and Sever provided an integral equation method for inverse conductivity problems using the linearization method [16]. Based on Calderόn’s linearization method, a new direct algorithm was suggested for the anisotropic conductivities [17]. Liu et al. [18] constructed two-parameter homogenization functions for solving the bending problem of a thin plate in a rectangular domain where the boundary conditions can be exactly satisfied. Using the Lie-group iterative method, Liu and Atluri [19] solved the linear Calderόn inverse problem in a rectangular domain, where the unknown conductivity function is effectively recovered. The linear and nonlinear inverse conductivity problems have also been studied by meshless methods, such as the meshless local Petrov-Galerkin method [20], the singular boundary method [21], and the local radial point interpolation method [22], the method of fundamental solutions [23], etc.
In this paper, based on the previous work in [18,19], we focus on the construction of two-parameter 2D homogenization functions in a doubly connected domain, and take linear equations to identify the space-dependent and temperature-dependent conductivity functions, the Robin coefficient function and also the source function in the 2D nonlinear elliptic equations. The derived homogenization functions are used as the bases. The undetermined expansion coefficients are solved by imposing the extra boundary conditions. In this way, the nonlinear inverse problems can be solved directly with high accuracy and efficiency even when twenty percent of noise is added to the known data.
We arrange the rest of this paper as follows. Section 2 describes some nonlinear inverse problems in a doubly connected domain of a 2D nonlinear elliptic equation, which includes the recovery of conductivity functions α ( x , y ) and α ( u ) , the inverse Robin problem and the inverse source problem. In Section 3, we develop the homogenization functions with two parameters. In Section 4, the two-parameter homogenization functions act as the bases for the solution. In Section 5, the space-dependent conductivities of inverse problems are considered. In Section 6, we solve the temperature-dependent conductivity inverse problems. The inverse Robin problem and one example are given in Section 7, and the inverse source problem is solved in Section 8, where two examples are given. Section 9 makes the conclusions.

2. Nonlinear Inverse Problems

For this part, we briefly sketch the problems to be considered that desire the retrieval of unknown functions in the doubly connected domains.

2.1. Space-Dependent Inverse Conductivity Problem

First a space-dependent conductivity function α ( x , y ) is to be recovered from
· [ α ( x , y ) u ( x , y ) ] = Q ( u , u x , u y ) + S ( x , y ) , ( x , y ) Ω R 2 ,
u ( x , y ) = h 1 ( x , y ) , ( x , y ) Γ o ,
u n ( x , y ) = g ( x , y ) , ( x , y ) Γ o ,
where n is an outward unit normal on Γ o . Besides an unknown conductivity function α ( x , y ) and the unknown solution u ( x , y ) , other functions are given.
Ω is a doubly connected domain with boundary Γ = Γ o Γ i , where Γ o Γ i = . While Γ o : = { ( r , θ ) | r = ρ o ( θ ) , 0 θ 2 π } denotes an outer boundary, Γ i : = { ( r , θ ) | r = ρ i ( θ ) , 0 θ 2 π } is an inner boundary. 0 < ρ i ( θ ) < ρ o ( θ ) are, respectively, the radius functions of inner boundary and outer boundary. In order to recover α ( x , y ) , we over-specify
u ( x , y ) = h 2 ( x , y ) , ( x , y ) Γ i ,
where h 2 ( x , y ) is a given function. In Figure 1, we sketch the inverse conductivity problem.
In the polar coordinates ( r , θ ) , Equation (1) is recast to
α u r r + 1 r u r + 1 r 2 u θ θ + α r u r + 1 r 2 α θ u θ = Q ( u , u r , u θ ) + S ( r , θ ) .
Equation (5) is a first-order PDE for the function α ( r , θ ) = α ( x , y ) with respect to r and θ , where u r , u θ , u r r and u θ θ are the coefficient functions. It is a nontrivial task to determine α even with the known u prescribed inside the Ω unless the boundary information of α on Γ is given in advance. Indeed, the inverse conductivity problem, which is considered in this paper, becomes more difficult and troublesome since the information of u is not given inside the solution domain, and only the boundary information is given according to Equations (2)–(4).

2.2. Temperature-Dependent Inverse Conductivity Problem

Secondly, we attempt to retrieve α ( u ) in
u r r + 1 r u r + 1 r 2 u θ θ α ( u ) + α ( u ) u r 2 + 1 r 2 α ( u ) u θ 2 = Q ( u , u r , u θ ) + S ( r , θ ) ,
when S and Q are given.
The temperature-dependent inverse conductivity problem is to determine the unknown conductivity function α ( u ) considering the above governing equation along with the information from Equations (2)–(4). The problem becomes harder for the reason that Equation (6) is nonlinear for u and linear ODE for α ( u ) with respect to u.

2.3. Inverse Robin Problem to Determine γ ( θ )

In the inner boundary, which is an inaccessible part of the boundary Γ , we cannot directly detect the transfer coefficient γ ( θ ) in
u n ( x , y ) + γ ( θ ) u ( x , y ) = h 3 ( x , y ) , ( x , y ) Γ i .
When the information of u ( x , y ) and u n ( x , y ) on Γ i is unknown, h 3 ( x , y ) is given. Taken into the consideration of Equations (1)–(3), the unknown Robin coefficient γ ( θ ) will be recovered, which is known as the inverse Robin problem.

2.4. Inverse Problem for S ( x , y )

When the function u ( x , y ) is known in advance, we can recover S ( x , y ) by
S ( x , y ) = · [ α ( x , y ) u ( x , y ) ] Q ( u , u x , u y ) ,
where α ( x , y ) and Q ( u , u x , u y ) are given functions. We will show that S ( x , y ) is recoverable from Equations (2)–(4) and (8), without solving nonlinear equations.

3. Two-Parameter Basis Functions

First, we demonstrate the basic idea of homogenization function by starting from a 2D boundary value problem (BVP):
L [ u ( x , y ) ] = S ( x , y ) , ( x , y ) ( 0 , a ) × ( 0 , b ) ,
u ( 0 , y ) = h 1 ( y ) , u ( a , y ) = h 2 ( y ) , u ( x , 0 ) = h 3 ( x ) , u ( x , b ) = h 4 ( x ) ,
where L is a second-order linear differential operator. Let
B 0 ( x , y ) = h 1 ( y ) 1 x a + x a h 2 ( y ) ,
and then,
B 0 ( 0 , y ) = h 1 ( y ) , B 0 ( a , y ) = h 2 ( y )
are apparent.
Upon letting
B ( x , y ) = B 0 ( x , y ) + 1 y b [ h 3 ( x ) B 0 ( x , 0 ) ] + y b [ h 4 ( x ) B 0 ( x , b ) ] ,
and according to the following compatibility conditions:
h 3 ( 0 ) = B 0 ( 0 , 0 ) = h 1 ( 0 ) , h 4 ( 0 ) = B 0 ( 0 , b ) = h 1 ( b ) , h 3 ( a ) = B 0 ( a , 0 ) = h 2 ( 0 ) , h 4 ( a ) = B 0 ( a , b ) = h 2 ( b ) ,
it is easy to verify
B ( 0 , y ) = h 1 ( y ) , B ( a , y ) = h 2 ( y ) , B ( x , 0 ) = h 3 ( x ) , B ( x , b ) = h 4 ( x ) .
Therefore, we can produce the 2D homogenization function for the 2D BVP:
B ( x , y ) = 1 x a h 1 ( y ) 1 y b h 3 ( 0 ) y b h 4 ( 0 ) + x a h 2 ( y ) 1 y b h 3 ( a ) y b h 4 ( a ) + 1 y b h 3 ( x ) + y b h 4 ( x ) .
Due to B ( x , y ) , we can transform the original 2D BVP with non-homogeneous boundary conditions to one with homogeneous boundary conditions:
L [ v ( x , y ) ] = S ( x , y ) L [ B ( x , y ) ] , ( x , y ) ( 0 , a ) × ( 0 , b ) ,
v ( 0 , y ) = v ( a , y ) = v ( x , 0 ) = v ( x , b ) = 0 ,
with the help of the variable transformation from u ( x , y ) to v ( x , y ) = u ( x , y ) B ( x , y ) . Obviously, Equations (17) and (18) are more easy to tackle than Equations (9) and (10). As an extension of B ( x , y ) to a two-parameter family, we have
B ( x , y , j , k ) = 1 x a j h 1 ( y ) 1 y b k h 3 ( 0 ) y b h 4 ( 0 ) + x a j h 2 ( y ) 1 y b k h 3 ( a ) y b k h 4 ( a ) + 1 y b k h 3 ( x ) + y b k h 4 ( x ) .
B ( x , y , j , k ) is indeed a family of 2D polynomials, which are complete bases and satisfy the boundary conditions automatically,
B ( 0 , y , j , k ) = h 1 ( y ) , B ( a , y , j , k ) = h 2 ( y ) , B ( x , 0 , j , k ) = h 3 ( x ) , B ( x , b , j , k ) = h 4 ( x ) .
A function is a so-called homogenization function if it satisfies the boundary conditions on the boundary of a domain. Since the solution u ( x , y ) must satisfy the prescribed boundary conditions, it is a member of homogenization functions.
Continuously, the two-parameter homogenization functions are constructed for developing the present method to solve the inverse problems of Equations (1)–(4).
Definition 1.
B 0 ( r , θ ) C 2 ( Ω ) , with Γ o = { ( r , θ ) | r = ρ o ( θ ) , 0 θ 2 π } , is a homogenization function, if the following conditions:
B 0 ( ρ o , θ ) = h 1 ( θ ) , B n 0 ( ρ o , θ ) = g ( θ )
are fulfilled. h 1 ( θ ) and g ( θ ) read as h 1 ( ρ o ( θ ) cos θ , ρ o ( θ ) sin θ ) and g ( ρ o ( θ ) cos θ , ρ o ( θ ) sin θ ) , respectively, and B n 0 signifying the normal derivative of B 0 ( r , θ ) on Γ o is given by
B n 0 ( ρ o , θ ) = η ( θ ) B 0 ( ρ o , θ ) ρ o ρ o ρ o 2 B 0 ( ρ o , θ ) θ ,
where
η ( θ ) = ρ o ( θ ) ρ o 2 ( θ ) + ρ o ( θ ) 2 .
The following homogenization function has been derived [18]:
B 0 ( r , θ ) = h 1 ( θ ) + [ r ρ o ( θ ) ] u ( ρ o , θ ) ρ o ,
B 0 ( ρ o , θ ) = h 1 ( θ ) , B n 0 ( ρ o , θ ) = g ( θ ) .
Theorem 1.
For the given Cauchy data h 1 ( θ ) and g ( θ ) on Γ o , there exist homogenization functions B ( j , k , r , θ ) in Ω, satisfying Equation (21):
B ( j , k , r , θ ) = 2 r ρ o r 2 ρ o 2 j h 1 ( θ ) + r k k ρ o k 1 ρ o k u ( ρ o , θ ) ρ o ,
where j + 1 , k N are parameters.
Proof. 
By Equation (26), B ( j , k , ρ o , θ ) = h 1 ( θ ) satisfies the first equation in Equation (21). Next, we consider the second equation in Equation (21), for which we need to prove
B ( j , k , ρ o , θ ) ρ o = u ( ρ o , θ ) ρ o , B ( j , k , ρ o , θ ) θ = u ( ρ o , θ ) θ .
It is obvious that
2 r ρ o r 2 ρ o 2 j r = ρ o = 1 , r 2 r ρ o r 2 ρ o 2 j r = ρ o = 0 ,
θ 2 r ρ o r 2 ρ o 2 j r = ρ o = j 2 r ρ o r 2 ρ o 2 j 1 2 r ρ o ρ o 2 + 2 r 2 ρ o ρ o 3 r = ρ o = 0 ,
r k k ρ o k 1 ρ o k r = ρ o = 0 , r r k k ρ o k 1 ρ o k r = ρ o = 1 ,
θ r k k ρ o k 1 ρ o k r = ρ o = r k k ( 1 k ) ρ o k ρ o ρ o k r = ρ o = ρ o .
It follows from Equations (26), (28) and (30) that the first part in Equation (27) holds when B ( j , k , r , θ ) is differentiated to r, and we take r = ρ o ( θ ) on Γ o .
The second part in Equation (27) is proven below. It follows from Equation (2) that
h 1 ( θ ) = u ( ρ o , θ ) ρ o ρ o ( θ ) + u ( ρ o , θ ) θ .
From Equations (26) and (28)–(32), it follows that
B ( j , k , ρ o , θ ) θ = h 1 ( θ ) ρ o ( θ ) u ( ρ o , θ ) ρ o = u ( ρ o , θ ) ρ o ρ o ( θ ) + u ( ρ o , θ ) θ ρ o ( θ ) u ( ρ o , θ ) ρ o = u ( ρ o , θ ) θ .
Due to Equation (27),
B n ( j , k , ρ o , θ ) = η ( θ ) B ( j , k , ρ o , θ ) ρ o ρ o ρ o 2 B ( j , k , ρ o , θ ) θ = η ( θ ) u ( ρ o , θ ) ( ρ o , θ ) ρ o ρ o ρ o 2 u ( ρ o , θ ) θ = g ( θ ) = u n ( x , y ) , ( x , y ) Γ o ,
thus we prove
B ( j , k , ρ o , θ ) = h 1 ( θ ) , B n ( j , k , ρ o , θ ) = g ( θ ) ,
which ends the proof of this theorem. □
In Theorem 1, the numbers ( j , k ) are parameters, and then B ( j , k , r , θ ) is a two-parameter function. In addition to Theorem 1, we also have the following result for another two-parameter function E ( j , k , r , θ ) .
Theorem 2.
On Γ o given the Cauchy data h 1 ( θ ) and g ( θ ) , the two-parameter function E ( j , k , r , θ ) C 2 ( Ω ) satisfies Equation (21):
E ( j , k , r , θ ) = B 0 ( r , θ ) + [ r ρ o ( θ ) ] 2 x j k y k = B 0 ( r , θ ) + [ r ρ o ( θ ) ] 2 r j ( cos θ ) j k ( sin θ ) k ,
where j + 1 , k N are parameters and B 0 was defined by Equation (24).
Proof. 
Let
E 0 ( r , θ ) : = [ r ρ o ( θ ) ] 2 .
When r = ρ o ( θ ) , it is obvious that
E 0 ( r , θ ) = [ r ρ o ( θ ) ] 2 = 0 , E n 0 ( r , θ ) = 2 [ r ρ o ( θ ) ] [ r ρ o ( θ ) ] n = 0 .
Inserting r = ρ o ( θ ) into Equation (36) and using Equations (24), (25) and (37), it follows that
E ( j , k , ρ o , θ ) = h 1 ( θ ) .
Taking the normal derivative of Equation (36) on Γ o and using Equations (24), (25) and (37), we have
E n ( j , k , ρ , θ ) = B n 0 ( ρ o , θ ) + E n 0 ( r , θ ) r j ( cos θ ) j k ( sin θ ) k + E 0 ( r , θ ) [ r j ( cos θ ) j k ( sin θ ) k ] n = B n 0 ( ρ o , θ ) = g ( θ ) , w h e n r = ρ o ( θ ) .
This completes the proof. □

4. A Novel Two-Parameter Homogenization Function Method

Since the set E ( j , k , r , θ ) is generated from the Pascal polynomials x j k y k , it is a complete basis for the problem. By the same token, B ( j , k , r , θ ) is a complete basis. All the homogenization functions consist of a weighted Sobolev space denoted as B : = { v ( x , y ) C 2 ( Ω ) | v ( x , y ) = h 1 ( x , y ) , v n ( x , y ) = g ( x , y ) , ( x , y ) Γ o } , which is a weighted space, because for any two functions v 1 ( x , y ) , v 2 ( x , y ) B with a weighted linear combination w 1 v 1 ( x , y ) + w 2 v 2 ( x , y ) B where w 1 + w 2 = 1 . The Sobolev norm
v ( x , y ) 2 : = 0 2 π [ v 2 ( ρ o cos θ , ρ o sin θ ) + v n 2 ( ρ o cos θ , ρ o sin θ ) ] d θ
is defined in the space B . More importantly, the approximate solution u ( x , y ) B .
In terms of the bases B ( j , k , x , y ) , u ( x , y ) can be expanded by
u ( x , y ) j = 0 m 1 k = 1 m a j k B ( j , k , x , y ) ,
where a j k satisfies
j = 0 m 1 k = 1 m a j k = 1 ,
and guarantees conditions (2) and (3) being satisfied by u ( x , y ) . The number of the coefficients a j k is n 1 = m 2 .
As shown in Equation (4), we suppose that there are N data of u ( x , y ) on the inner boundary Γ i available, and then we can solve a linear system, including Equation (40), to determine a j k :
j = 0 m 1 k = 1 m a j k B ( j , k , x q , y q ) = h 2 ( x q , y q ) ,
where θ q = 2 q π / N , x q = ρ i ( θ q ) cos θ q and y q = ρ i ( θ q ) sin θ q .

5. Numerical Procedure to Determine α ( x , y )

5.1. Numerical Algorithm

Next, when u ( x , y ) is obtained from Equation (39), we recover α ( x , y ) by supposing
α ( x , y ) = i = 0 m 0 j = 0 i b i j x i j y j = i = 0 m 0 j = 0 i b i j r i ( cos θ ) i j ( sin θ ) j ,
where b i j are n : = ( m 0 + 1 ) ( m 0 + 2 ) / 2 unknown weighted parameters to be determined by the proposed numerical algorithm. In order to solve this problem, m 1 × m 2 points of ( x , y ) inside the solution domain Ω are collocated by
x p q = r p cos θ q , y p q = r p sin θ q , θ q = 2 q π / m 1 , r p = ρ i ( θ q ) + p [ ρ o ( θ q ) ρ i ( θ q ) ] / ( m 2 + 1 ) , q = 1 , , m 1 , p = 1 , , m 2 .
Since u ( x , y ) can be approximated by Equation (39), we have
u ( x p q , y p q ) = j = 0 m 1 k = 1 m a j k B ( j , k , x p q , y p q ) .
Then, inserting Equation (42) into Equation (5) and collocating at point ( x p q , y p q ) , the following linear system can be obtained:
Δ u ( x p q , y p q ) i = 0 m 0 j = 0 i b i j x p q i j y p q j + u r ( x p q , y p q ) i = 0 m 0 j = 0 i i b i j r p i 1 ( cos θ q ) i j ( sin θ q ) j + 1 r p 2 u θ ( x p q , y p q ) i = 0 m 0 j = 0 i b i j r p i [ j ( cos θ q ) i j + 1 ( sin θ q ) j 1 ( i j ) ( cos θ q ) i j 1 ( sin θ q ) j + 1 ] = Q ( u ( x p q , y p q ) , u r ( x p q , y p q ) , u θ ( x p q , y p q ) ) + S ( r p , θ q ) , q = 1 , , m 1 , p = 1 , , m 2 ,
from which we can determine b i j easily, and, correspondingly, the α ( x , y ) can be determined from Equation (42).
Therefore, the proposed algorithm for recovering α ( x , y ) consists of two linear systems of equations, Equations (41) and (45). We impose the data by a noise:
h ^ 1 ( θ j ) = h 1 ( θ j ) + s R ( j ) , g ^ ( θ j ) = g ( θ j ) + s R ( j ) ,
where R ( j ) are random numbers between [ 1 , 1 ] , which are used to check the stability of the numerical solution.
To evaluate the accuracy, we consider the maximum error (ME) and a relative error defined by
ME ( α ) : = max | α n ( x p q , y p q ) α ( x p q , y p q ) | ,
e ( α ) = q = 1 N 1 p = 1 N 2 [ α n ( x p q , y p q ) α ( x p q , y p q ) ] 2 q = 1 N 1 p = 1 N 2 α 2 ( x p q , y p q ) ,
upon comparing the numerical solution of α n to the exact one α at N 1 × N 2 grid points ( x i , y j ) inside the domain with
x p q = r p cos θ q , y p q = r p sin θ q , θ q = 2 q π / N 1 , r p = ρ i ( θ q ) + p [ ρ o ( θ q ) ρ i ( θ q ) ] / ( N 2 + 1 ) , q = 1 , , N 1 , p = 1 , , N 2 .
We take N 1 = 50 and N 2 = 10 in all computations.
We must emphasize that all the linear systems we consider are over-determined, which means that the number of linear equations is much larger than the number of unknown coefficients. Therefore, we apply the conjugate gradient method (CGM) to solve the corresponding normal linear system, whose solution is unique in the sense of least squares.

5.2. Example 1

For Equation (5) with Q = 0 , we consider
u = r 2 + r cos θ = x 2 + y 2 + x , α = 10 + r 2 + r 4 cos 2 θ = 10 + ( x 2 + y 2 ) ( x 2 + 1 ) .
The outer boundary of the domain Ω is an ellipse:
ρ o ( θ ) = a b b 2 + ( a 2 b 2 ) sin 2 θ ,
where we take a = 1.5 and b = 0.5 , and the inner boundary is ρ i = 0.2 .
With m = 2 , N = 40 , m 0 = 4 , m 1 = 30 , m 2 = 10 and a noise s = 5 % added into the given data, Figure 2a reveals that the maximum absolute error of u is 2.11 × 10 2 , which is much smaller than max ( u ) = 3.29 . The maximum absolute error of coefficient α ( x , y ) denoted as ME ( α ) is 0.18, which is quite a lot smaller than max ( α ) = 15.56 . The value e ( α ) = 3.8 × 10 3 is smaller than 7.65 × 10 2 from previous studies. For this problem the dimension of the normal matrix of the first linear system of (41) and (40) is n 1 × n 1 = 4 × 4 , and the condition number is small with COND = 570.046. The dimension of the normal matrix of the second linear system (45) is n × n = 15 × 15 , and the condition number is COND = 41,996.23. They show that these two linear systems are stable.
The convergence rate is a central issue in numerical methods and algorithms. In Table 1, we consider a different mesh parameter m 1 × m 2 used in the collocation method to influence the convergence rate as reflected in ME( α ) and e ( α ) . It can be seen that more collocated points lead to a more accurate solution of α .
Although for a nonlinear elliptic equation:
α u r r + 1 r u r + 1 r 2 u θ θ + α r u r + 1 r 2 α θ u θ = u 2 ( r , θ ) + S ( r , θ ) ,
where the exact value of S ( r , θ ) can be obtained by inserting Equation (50) into the above equation, as shown in Figure 2b, the ME of u is 2.11 × 10 2 , the ME of α ( x , y ) is 0.18 and e ( α ) = 3.74 × 10 3 . For the nonlinear problem, we have the same condition numbers because the parameters used are the same.

5.3. Example 2

Consider
α u r r + 1 r u r + 1 r 2 u θ θ + α r u r + 1 r 2 α θ u θ = sin u ( r , θ ) + u 2 ( r , θ ) + S ( r , θ ) ,
u = r 2 + r cos θ = x 2 + y 2 + x , α = 20 + r 2 sin ( 2 θ ) = 20 + 2 x y .
For this problem, the outer boundary is given by Equation (51) with a = 4 and b = 3.5 , and
ρ i ( θ ) = 1.5 + cos θ
is the inner boundary. Feeding Equation (54) into Equation (53), S ( r , θ ) can be obtained.
With m = 2 , N = 50 , m 0 = 2 , m 1 = 10 , m 2 = 10 and a noise s = 20 % , as shown in Figure 3, the ME of u is 1.44 × 10 2 , which is much smaller than max ( u ) = 18.79 . The ME of α ( x , y ) is 2.47, which is much smaller than max ( α ) = 32.96 . The value e ( α ) = 3.85 × 10 2 is small. For this problem the dimension of the normal matrix of the first linear system (41) and (40) is n 1 × n 1 = 4 × 4 , and the condition number is small with COND = 1844.258. The dimension of the normal matrix of the second linear system (45) is n × n = 6 × 6 , and the condition number is COND = 1013.03. They show that these two linear systems are stable.

6. Numerical Algorithm to Determine α ( u )

6.1. Numerical Algorithm

From the last section, we have already recovered the coefficients preceding α ( u ) and u r 2 + u θ 2 / r 2 preceding α ( u ) in Equation (6), if Q and S are prescribed in advance. Indeed u ( x , y ) can be derived from Equation (39). In this case, Δ u before α ( u ) , and u r 2 + u θ 2 / r 2 before α ( u ) can be obtained numerically from Equation (39).
Suppose that
α ( u ) = i = 0 m 0 c i u i ,
where c i are under-determined weighted parameters to be determined.
Similar to Equation (43) in the last section, we arrange m 1 × m 2 points of ( x , y ) inside the solution domain. Then, inserting Equation (56) into Equation (6) and collocating ( x p q , y p q ) , we come to
i = 0 m 0 c i u i ( x p q , y p q ) Δ u ( x p q , y p q ) + i = 0 m 0 c i i u i 1 ( x p q , y p q ) [ u r 2 ( x p q , y p q ) + u θ 2 ( x p q , y p q ) / r p 2 ] = Q ( u ( x p q , y p q ) , u r ( x p q , y p q ) , u θ ( x p q , y p q ) ) + S ( x p q , y p q ) , p = 1 , , m 2 , q = 1 , , m 1 .
from which we can obtain c i , and then, α ( u ) is recovered from Equation (56).
It should be noted here that that even for the highly nonlinear inverse problems for coefficient α ( u ) , solving nonlinear equations is not needed.

6.2. Example 3

For a quadratic nonlinear Poisson equation:
Δ u α ( u ) + α ( u ) u r 2 + 1 r 2 α ( u ) u θ 2 + u 2 = S ( r , θ ) ,
u = r 2 = x 2 + y 2 ,
α ( u ) = 10 + u 2 + u is to be recovered. The outer boundary is Equation (51) with a = 3.5 and b = 2.5 , and ρ i = 1 is a unit circle.
With m = 2 , N = 40 , m 0 = 2 , m 1 = 20 , m 2 = 5 and s = 5 % , as shown in Figure 4a, the ME of u is 5.68 × 10 3 , which is much smaller than max ( u ) = 10.7 . The ME of α ( u ) is 0.082, which is much smaller than max ( α ) = 135.43 . The value e ( α ) = 6.03 × 10 4 is quite small. Figure 4b compares the numerical and exact α ( u ) in the range of u 11 . These two curves almost coincide.
For this problem, the dimension of the normal matrix of the first linear system (41) and (40) is n 1 × n 1 = 4 × 4 , and the condition number is very small with COND = 33.959. The dimension of the normal matrix of the second linear system (57) is n × n = 3 × 3 , and the condition number is COND = 177,571.79. They show that these two linear systems are stable.

6.3. Example 4

For a quadratic and cubic nonlinear Poisson equation:
Δ u α ( u ) + α ( u ) u r 2 + 1 r 2 α ( u ) u θ 2 + u 2 + u 3 = S ( r , θ ) ,
u = r 2 = x 2 + y 2 ,
α ( u ) = 10 + u 2 + sin u is to be recovered. The boundaries of Ω are given by
ρ o ( θ ) = 3 cos ( 2 θ ) + 1.5 sin 2 ( 2 θ ) ,
ρ i ( θ ) = exp ( sin θ ) sin 2 ( 2 θ ) + exp ( cos θ ) cos 2 ( 2 θ ) .
With m = 2 , N = 40 , m 0 = 3 , m 1 = 15 , m 2 = 10 and s = 5 % , as shown in Figure 5a, the ME of u is 2.78 × 10 2 , which is much smaller than max ( u ) = 18.6 . The ME of α ( u ) is 2.23, which is much smaller than max ( α ) = 356.45 . The value e ( α ) = 9.35 × 10 3 is quite small. Figure 5b compares the numerical and exact α ( u ) in the range of u 19 . These two curves almost coincide.
For this problem, the dimension of the normal matrix of the first linear system (41) and (40) is n 1 × n 1 = 4 × 4 , and the condition number is very small with COND = 6.3185. The dimension of the normal matrix of the second linear system (57) is n × n = 4 × 4 , and the condition number is COND = 2.15 × 10 8 .
To reduce the condition number for the second linear system (57), we choose m 0 = 2 , m 1 = 10 and m 2 = 10 , such that the dimension of the normal matrix reduces to n × n = 3 × 3 , and the condition number reduces to COND = 297,651.317. Meanwhile, ME ( α ) and e ( α ) are slightly increased to 2.38 and 1.08 × 10 2 , respectively.

7. Numerical Method to Detect γ ( θ )

7.1. Numerical Method

Now, we detect γ ( θ ) by using the data in Equations (2) and (3). Basically, we need to solve Equations (1)–(3) in Ω . For this purpose, we take
u ( x , y ) = j = 0 m k = 0 j a j k E ( j , k , x , y ) ,
where the number of a j k , j , k = 1 , , m is n 1 = ( m + 1 ) ( m + 2 ) / 2 , which are to be determined. Instead of B ( j , k , x , y ) used in Equation (39), we employ E ( j , k , x , y ) from Equation (36) as the bases of u ( x , y ) . The reason is that the order of r j + 2 in E ( j , k , x , y ) is much lower than the order of r 2 j in B ( j , k , x , y ) .
Like Equation (43), we arrange m 1 × m 2 points of ( x , y ) and collocating which comes to:
α ( x p q , y p q ) j = 0 m k = 0 j a j k Δ E ( j , k , x p q , y p q ) + α x ( x p q , y p q ) j = 0 m k = 0 j a j k E x ( j , k , x p q , y p q ) + α y ( x p q , y p q ) j = 0 m k = 0 j a j k E y ( j , k , x p q , y p q ) = Q j = 0 m k = 0 j a j k E ( j , k , x p q , y p q ) , j = 0 m k = 0 j a j k E x ( j , k , x p q , y p q ) , j = 0 m k = 0 j a j k E y ( j , k , x p q , y p q ) + S ( x p q , y p q ) , p = 1 , , m 2 , q = 1 , , m 1 .
from which we can compute a j k , and then u ( x , y ) is obtained from Equation (64), which is inserted in Equation (7) to find γ ( θ ) along the inner boundary Γ i .

7.2. Example 5

We give a solution of a linear diffusion-convection equation:
α Δ u ( x , y ) + α x u x ( x , y ) + α y u y ( x , y ) = S ( x , y ) ,
u = x 2 + y 2 + x ,
where α = 1 + x 2 + y 2 , which is defined in a domain Ω by Equation (51) with a = 2.5 and b = 1.5 and by
ρ i ( θ ) = 1.3 + 0.1 cos θ .
With m = 2 , m 1 = 20 , m 2 = 20 and s = 20 % , as shown in Figure 6a, the ME of u is 4.58 × 10 3 , which is much smaller than max ( u ) = 8.16 . Figure 6b compares the numerical and exact γ ( θ ) , of which these two curves almost coincide with the ME being 1.28 × 10 2 . For this problem, we merely solve the linear system (65) and (40), whose dimension of the normal matrix is n 1 × n 1 = 6 × 6 , and the condition number is small with COND = 446.16.

8. Numerical Method to Recover S ( x , y )

The numerical method to recover S ( x , y ) is very simple, which is obtained by merely inserting the numerical solution of u ( r , θ ) in Equation (39) into the following equation:
S ( r , θ ) = α u r r + 1 r u r + 1 r 2 u θ θ + α r u r + 1 r 2 α θ u θ Q ( u ( r , θ ) , u r ( r , θ ) , u θ ( r , θ ) ) ,
where α ( r , θ ) and Q ( u ( r , θ ) , u r ( r , θ ) , u θ ( r , θ ) ) are given functions.

8.1. Example 6

For Equation (69) with Q = u 2 , we consider
u = r 2 + r cos θ = x 2 + y 2 + x , α = 10 + r 2 + r 4 cos 2 θ = 10 + ( x 2 + y 2 ) ( x 2 + 1 ) ,
with a = 1.5 and b = 0.5 in Equation (51), and we take ρ i = 0.2 .
In this case, we have m = 2 , N = 30 and a noise s = 5 % added into the given data, as shown in Figure 7a, the maximum absolute error of u is 1.35 × 10 2 , which is much smaller than max ( u ) = 3.29 . The maximum absolute error of S ( x , y ) , denoted as ME ( S ) , is 1.22, which is much more accurate than max S = 101.5 . The value e ( S ) = 5.2 × 10 3 is small. For this problem, we merely solve the linear system (41) and (40), whose dimension of the normal matrix is n 1 × n 1 = 4 × 4 , and the condition number is small with COND = 476.623.
In Table 2, we consider different mesh parameters N used in the collocation method for u ( x , y ) to influence the convergence rate as reflected in ME(S) and e ( S ) . It can be seen that more collocated points lead to a more accurate solution of S.

8.2. Example 7

For Equation (69) with Q = u 2 + cos u , we consider
u = r 2 = x 2 + y 2 , α = 10 + r 2 + r 4 cos 2 θ = 10 + ( x 2 + y 2 ) ( x 2 + 1 ) .
The outer and inner boundaries are given by Equations (62) and (63), respectively.
For this example, we have m = 2 , N = 20 and a noise s = 5 % , as shown in Figure 7b, the maximum absolute error of u is 2.03 × 10 2 , which is more accurate than max ( u ) = 18.62 . The ME of S ( x , y ) is 1.04, which is much smaller than max S = 4001.42 . The value e ( S ) = 2.99 × 10 4 is quite small. The dimension of the normal matrix of the linear system (41) and (40) is n 1 × n 1 = 4 × 4 , and the condition number is small with COND = 19.858.
In Table 3, we consider different mesh parameters of N used in the collocation method for u ( x , y ) to influence the convergence rate as reflected in ME(S) and e ( S ) . It can be seen that more collocated points lead to a more accurate solution of S and even for a small N = 3 the accuracy is good.

9. Conclusions

In the paper, we have constructed a category of two-parameter homogenization functions in the 2D doubly connected domain for automatically satisfying the outer Dirichlet and Neumann boundary conditions of the nonlinear elliptic equation. A new numerical method was developed for solving the inverse problems through the technique of two-parameter homogenization functions, which include the recovery of the space-dependent and temperature-dependent conductivity functions and also the source function. We first determine u ( x , y ) in terms of the bases and then a linear system to satisfy the inner boundary condition by the method of collocation is solved. Back-substituting the solution into the nonlinear elliptic equation, we recovered the unknown space-dependent and temperature-dependent conductivity functions by collocating points inside the domain and solving the derived linear equations. The basis B ( j , k , x , y ) has good behavior used in the interpolation for u ( x , y ) in a weighted Sobolev space, such that we can recover u ( x , y ) very well; hence, after the back substitution of u ( x , y ) into the governing equation, the source function was directly recovered with high accuracy. It maintains the same advantages of accuracy and efficiency for solving the inverse conductivity problems and inverse Robin problems, even for large noise.

Author Contributions

Conceptualization, C.-S.L.; Investigation, L.S.; Methodology, J.L.; Writing—review & editing, C.S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program (Nos. 2021YFB2600700, 2021YFC3090100), the National Natural Science Foundation of China (No. 52171272), and Key Special Projects of the Science and Technology Help Economy 2020 (No. 1).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tuan, N.H.; Khoa, V.A.; Tran, T. On an inverse boundary value problem of a nonlinear elliptic equation in three dimensions. J. Math. Anal. Appl. 2015, 426, 1232–1261. [Google Scholar] [CrossRef]
  2. Farcas, A.; Elliott, L.; Ingham, D.B.; Lesnic, D.; Mera, N.S. A dual reciprocity boundary element method for the regularized numerical solution of the inverse source problem associated to the Poisson equation. Inv. Prob. Sci. Eng. 2003, 11, 123–139. [Google Scholar] [CrossRef]
  3. Jin, B.T.; Marin, L. The method of fundamental solutions for inverse source problems associated with the steady-state heat conduction. Int. J. Numer. Meth. Eng. 2007, 69, 1570–1589. [Google Scholar] [CrossRef]
  4. Klose, A.D.; Ntziachristos, V.; Hielscher, A.H. The inverse source problem based on the radiative transfer equation in optical molecular imaging. J. Comput. Phys. 2005, 202, 323–345. [Google Scholar] [CrossRef]
  5. Hon, Y.C.; Li, M.; Melnikov, Y.A. Inverse source identification by Green’s function. Eng. Anal. Bound. Elem. 2010, 34, 352–358. [Google Scholar] [CrossRef]
  6. Li, X.X.; Guo, H.Z.; Wan, S.M.; Yang, F. Inverse source identification by the modified regularization method on Poisson equation. J. Appl. Math. 2012, 2012, 971952. [Google Scholar] [CrossRef]
  7. Ahmadabadi, M.N.; Arab, M.; Ghaini, F. The method of fundamental solutions for the inverse space-dependent heat source problem. Eng. Anal. Bound. Elem. 2009, 33, 1231–1235. [Google Scholar] [CrossRef]
  8. Slimani, S.; Medarhri, I.; Najib, K.; Zine, A. Identification of the source function for a seawater intrusion problem in unconfined aquifer. Numer. Algor. 2020, 84, 1565–1587. [Google Scholar] [CrossRef]
  9. Alahyane, M.; Boutaayamou, I.; Chrifi, A.; Echarroudi, Y.; Ouakrim, Y. Numerical study of inverse source problem for internal degenerate parabolic equation. Int. J. Comput. Meth. 2020, 18, 2050032. [Google Scholar] [CrossRef]
  10. Djennadi, S.; Shawagfeh, N.; Arqub, O.A. A fractional Tikhonov regularization method for an inverse backward and source problems in the time-space fractional diffusion equations. Chaos Solitons Fractals 2021, 150, 111127. [Google Scholar] [CrossRef]
  11. Ma, Y.K.; Prakash, P.; Deiveegan, A. Generalized Tikhonov methods for an inverse source problem of the time-fractional diffusion equation. Chaos Solitons Fractals 2018, 108, 39–48. [Google Scholar] [CrossRef]
  12. Nguyen, H.T.; Le, D.L. Regularized solution of an inverse source problem for a time fractional diffusion equation. Appl. Math. 2016, 40, 8244–8264. [Google Scholar] [CrossRef]
  13. Liu, C.-S. An energetic boundary functional method for solving the inverse source problems of 2D nonlinear elliptic equations. Eng. Anal. Bound. Elem. 2020, 118, 204–215. [Google Scholar] [CrossRef]
  14. Kwon, K. Identification of anisotropic anomalous region in inverse problems. Inverse Probl. 2004, 20, 1117–1136. [Google Scholar] [CrossRef]
  15. Huntul, M.J.; Lesnic, D. An inverse problem of finding the time-dependent thermal conductivity from boundary data. Int. Commun. Heat Mass Transf. 2017, 85, 147–154. [Google Scholar] [CrossRef] [Green Version]
  16. Isakov, V.; Sever, A. Numerical implementation of an integral equation method for the inverse conductivity problem. Inverse Probl. 1996, 12, 939–951. [Google Scholar] [CrossRef]
  17. Murthy, R.; Lin, Y.H.; Shin, K.; Mueller, J.L. A direct reconstruction algorithm for the anisotropic inverse conductivity problem based on Calderόn method in the plane. Inverse Probl. 2020, 12, 125008. [Google Scholar] [CrossRef]
  18. Liu, C.-S.; Qiu, L.; Lin, J. Simulating thin plate bending problems by a family of two-parameter homogenization functions. Appl. Math. Model. 2020, 79, 284–299. [Google Scholar] [CrossRef]
  19. Liu, C.-S.; Atluri, S.N. An iterative and adaptive Lie-group method for solving the Calderόn inverse problem. Comput. Model. Eng. Sci. 2010, 64, 299–326. [Google Scholar]
  20. Sladek, J.; Sladek, V.; Hon, Y.C. Inverse heat conduction problems by meshless local Petrov–Galerkin method. Eng. Anal. Boundary Elem. 2006, 30, 650–661. [Google Scholar] [CrossRef]
  21. Gu, Y.; Chen, W.; Zhang, C.; He, X. A meshless singular boundary method for three-dimensional inverse heat conduction problems in general anisotropic media. Int. J. Heat Mass Transf. 2015, 84, 91–102. [Google Scholar] [CrossRef]
  22. Shivanian, E.; Khodabandehlo, H.R. Application of meshless local radial point interpolation (MLRPI) on a one-dimensional inverse heat conduction problem. Ain Shams Eng. J. 2016, 7, 993–1000. [Google Scholar] [CrossRef] [Green Version]
  23. Sun, Y.; He, S. A meshless method based on the method of fundamental solution for three-dimensional inverse heat conduction problems. Int. J. Heat Mass Transf. 2017, 108, 945–960. [Google Scholar] [CrossRef]
Figure 1. A schematic plot to show a doubly connected domain and for identification.
Figure 1. A schematic plot to show a doubly connected domain and for identification.
Mathematics 10 02256 g001
Figure 2. For example 1, showing the errors in the numerical recovery of u and α under a noise s = 0.05 , (a) Q = 0 and (b) Q = u 2 .
Figure 2. For example 1, showing the errors in the numerical recovery of u and α under a noise s = 0.05 , (a) Q = 0 and (b) Q = u 2 .
Mathematics 10 02256 g002
Figure 3. For example 2, showing the errors in the numerical recovery of u and α under a noise s = 0.2 with Q = u 2 + sin u .
Figure 3. For example 2, showing the errors in the numerical recovery of u and α under a noise s = 0.2 with Q = u 2 + sin u .
Mathematics 10 02256 g003
Figure 4. For example 3, showing (a) the errors in the numerical recovery of u and α and (b) the numerical recovery of α under a noise s = 0.05 .
Figure 4. For example 3, showing (a) the errors in the numerical recovery of u and α and (b) the numerical recovery of α under a noise s = 0.05 .
Mathematics 10 02256 g004
Figure 5. For example 4, showing (a) the errors in the numerical recovery of u and α and (b) the numerical recovery of α under a noise s = 0.05 .
Figure 5. For example 4, showing (a) the errors in the numerical recovery of u and α and (b) the numerical recovery of α under a noise s = 0.05 .
Mathematics 10 02256 g005
Figure 6. For example 5, showing (a) the error in the numerical recovery of u and (b) comparing the numerical recovery of γ under a noise s = 0.2 .
Figure 6. For example 5, showing (a) the error in the numerical recovery of u and (b) comparing the numerical recovery of γ under a noise s = 0.2 .
Mathematics 10 02256 g006
Figure 7. Showing the errors in the numerical recovery of u and S under a noise s = 0.05 , (a) example 6 and (b) example 7.
Figure 7. Showing the errors in the numerical recovery of u and S under a noise s = 0.05 , (a) example 6 and (b) example 7.
Mathematics 10 02256 g007
Table 1. For example 1, the influence of mesh parameter m 1 × m 2 on ME( α ) and e ( α ) .
Table 1. For example 1, the influence of mesh parameter m 1 × m 2 on ME( α ) and e ( α ) .
m 1 × m 2 5 × 5 7 × 5 10 × 10 10 × 5 30 × 10
ME( α )7.73691.84300.40170.21470.18364
e ( α ) 1.46 × 10 1 2.34 × 10 2 1.83 × 10 2 5.15 × 10 2 3.79 × 10 3
Table 2. For example 6, the influence of mesh parameter N on ME(S) and e ( S ) .
Table 2. For example 6, the influence of mesh parameter N on ME(S) and e ( S ) .
N510202530
ME(S)57.3443.5945.9337.061.22
e ( S ) 1.18 × 10 1 9.55 × 10 2 8.70 × 10 2 6.52 × 10 2 5.24 × 10 3
Table 3. For example 7, the influence of mesh parameter N on ME(S) and e ( S ) .
Table 3. For example 7, the influence of mesh parameter N on ME(S) and e ( S ) .
N38141820
ME(S)12.0910.988.988.511.04
e ( S ) 3.70 × 10 3 2.87 × 10 3 2.57 × 10 3 2.17 × 10 3 2.99 × 10 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, J.; Shi, L.; Liu, C.-S.; Chen, C.S. Solving Inverse Conductivity Problems in Doubly Connected Domains by the Homogenization Functions of Two Parameters. Mathematics 2022, 10, 2256. https://doi.org/10.3390/math10132256

AMA Style

Lu J, Shi L, Liu C-S, Chen CS. Solving Inverse Conductivity Problems in Doubly Connected Domains by the Homogenization Functions of Two Parameters. Mathematics. 2022; 10(13):2256. https://doi.org/10.3390/math10132256

Chicago/Turabian Style

Lu, Jun, Lianpeng Shi, Chein-Shan Liu, and C. S. Chen. 2022. "Solving Inverse Conductivity Problems in Doubly Connected Domains by the Homogenization Functions of Two Parameters" Mathematics 10, no. 13: 2256. https://doi.org/10.3390/math10132256

APA Style

Lu, J., Shi, L., Liu, C. -S., & Chen, C. S. (2022). Solving Inverse Conductivity Problems in Doubly Connected Domains by the Homogenization Functions of Two Parameters. Mathematics, 10(13), 2256. https://doi.org/10.3390/math10132256

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop