Next Article in Journal
On Strengthened Extragradient Methods Non-Convex Combination with Adaptive Step Sizes Rule for Equilibrium Problems
Previous Article in Journal
Uncertainty of Preferences in the Assessment of Supply Chain Management Systems Using the PROMETHEE Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Meshless Method for Solving 3D Inverse Conductivity Issues of Highly Nonlinear Elliptic Equations

Department of Mechanical Engineering, National United University, Miaoli 360302, Taiwan
Submission received: 18 April 2022 / Revised: 16 May 2022 / Accepted: 18 May 2022 / Published: 19 May 2022
(This article belongs to the Section Mathematics)

Abstract

:
In this research, the 3D inverse conductivity issues of highly nonlinear elliptic partial differential equations (PDEs) are investigated numerically. Even some researchers have utilized several schemes to overcome these multi-dimensional forward issues of those PDEs; nevertheless, an effective numerical algorithm to solve these 3D inverse conductivity issues of highly nonlinear elliptic PDEs is still not available. We apply two sets of single-parameter homogenization functions as the foundations for the answer and conductivity function to cope with the 3D inverse conductivity issue of highly nonlinear PDEs. The unknown conductivity function can be retrieved by working out another linear system produced from the governing equation by collocation skill, while the answer is acquired by dealing with a linear system to gratify over-specified Neumann boundary condition on a fractional border. As this new computational approach is based on a concrete theoretical foundation, it can result in a deeper understanding of 3D inverse conductivity issues with symmetry and asymmetry geometries. The homogenization functions method is rather stable, effective, and accurate in revealing the conductivity function when the over-specified Neumann data with a large level of noise of 28%.

1. Introduction

For inverse conductivity issues of nonlinear elliptic partial differential equations (PDEs), Calderón [1] discussed the inverse boundary value problem (BVP), which originated from the problem of electrical prospection. Furthermore, the Calderón problem is also called the inverse conductivity problem. He was the first to fabricate the phrase of an inverse BVP for studying the electrical impedance tomography (EIT) technique by mathematical modeling because the conductivity occurs as a variable coefficient of diffusion in an elliptic PDE. It is an inverse issue of the Neumann to Dirichlet mapping by realizing the voltage and current on the boundary. Computerized tomography is a standard instrument in the non-destructive testing (NDT) of materials and medical diagnostics right now. In addition, the well-known schemes of magnetic resonance imaging and X-ray, the last four decades have testified to an increasing interest in the new imaging skill of EIT because of its low cost and simple implementation. A comprehensive review of EIT was shown in Borcea [2]. In the medical and industrial fields, we often encounter those difficult ill-posed inverse problems. There have been many papers in the engineering and mathematical literature resolving such various themes as the electrical impedance tomography, the approaches of numerical reestablishment, the stability and identifiability of solutions for conductivity, mathematical modeling of electrodes, schemes of numerical establishment, and the design of measurement devices [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21].
In general, we can utilize the NDT to deal with the corrosion of metals. Identifying the impedance coefficient may be an effective method to reveal the region of corrosion, and one may possibly evaluate the level of corrosion in accordance with the exterior estimations and then employ the solution of the inverse issue by the EIT. The Robin coefficient quantifies the level of corrosion of some inaccessible part of the boundary in a simpler corrosion model [22]. The model issue was to retrieve the Robin coefficient through the present flux and potential surveys on fractionally accessible data in a finite annulus.
Nevertheless, merely a few papers are interested in inverse conductivity issues of nonlinear elliptic PDEs. Liu and Atluri [23] have addressed the Calderón inverse issue to retrieve the unknown conductivity function by utilizing the iterative scheme of Lie-group type together with a finite-strip method. The advantages of the current algorithm were that no a priori information about the functional form of conductivity function was essential, and no added measurement of data inside the area was needed. Later, Liu and Liu [24] established a homogenization function to gratify Cauchy boundary conditions on a smooth 2D boundary of a simply connected area. This issue for choosing the inner boundary of the Poisson equation in an arbitrary doubly-connected plane area was addressed, which retrieved an unknown inner boundary of a rigid inclusion under the over-specified Cauchy data on the accessible outer boundary. The accuracy and robustness of the current homogenization boundary function method (HBFM) were assessed through several numerical experiments by comparing the exact inner boundary to the retrieved one under a large noisy effect. Then, Liu and Wang [25] utilized the HBFM to deal with a nonlinear inverse Cauchy issue of nonlinear elliptic type PDE in an arbitrary doubly-connected plane area. The unknown Dirichlet data on an inner boundary were retrieved by over-specifying the Cauchy data on an outer boundary. A homogenization function was derived from annihilating the Cauchy data on the outer boundary, and then a homogenization skill generated a transformed equation in terms of a transformed variable whose outer Cauchy boundary data were homogeneous. Besides, they examined several numerical instances and obtained good results even with a large, noisy disturbance.
For the mathematical analysis of the Calderón problem, Caro et al. [26] proved the corresponding global evaluations in all dimensions higher than three. The evaluations were based on the construction of solutions to the Schrödinger equation by complicated geometrical optics developed in the anisotropic setting by Dos Santos Ferreira et al. [27] to deal with the Calderón issue in certain admissible geometries. Later, Piiroinen and Simon [28] acquired a probabilistic formulation of Calderón’s inverse conductivity issue. This formulation came in three equivalent versions and each of them may yield both a novel perspective as well as a novel set of (probabilistic) tools when it came to studying questions related to the unique determinability of conductivities from boundary data. After that, Alessandrini et al. [29] contemplated the electrostatic inverse BVP, also known as EIT, for the situation where the conductivity was a piecewise linear function on the area, and they demonstrated that a Lipschitz stability evaluation for the conductivity in terms of the local Dirichlet-to-Neumann (DN) map held true. Apart from that, Muñoz [30] displayed the uniqueness of the conductivity for the quasilinear Calderón’s inverse issue. Under some structural suppositions on the direct issue, a real-valued conductivity allowing a small analytic continuation to the complex plane induce a unique DN map. This approach of proof considered some complex-valued, linear test functions based on the point of the boundary of the area and a linearization of the DN map placed at these special sets of solutions. However, they did not demonstrate some examples to prove their algorithm.
For the applications of EIT, de Moura et al. [31] used EIT in multiphase flow for identifying flow patterns and void-fraction calculation. An image made of the conductivity of the domain was formed by combining an inverse issue and measurements. For that issue, they presented a recursive Gauss–Newton optimization sequential importance resampling (GNOSIR) filter aimed to reconstruct vertical slug flow. Their results demonstrated that the GNOSIR filter improved the evaluation when there was a sudden and fast evolution of the inclusion’s contour between states; however, they only coped with 2D problems. Later, Duran et al. [32] coped with the EIT image reconstruction by using convolutional neural networks (CNN), which was that the resulting image of the convolution process was smaller than the original input image. This issue was usually addressed by introducing padding, which was the addition of layers in the borders of the original input image. Their work presented the use of a doubly periodic padding, which was relevant for toroidal image issues such as the electric potential distribution measured by utilizing the EIT. The CNN was trained by employing a database generated by numerical calculations. Their resulting image reconstructions were presented for different noisy potential inputs; however, they merely dealt with 2D problems. After that, Wagner et al. [33] explored the evaluation of spatial strain distributions by elastoresistive thin-film sensors utilizing 2D Electrical Impedance Tomography. Hence, spatially-limited strains and the subsequent conductivity changes in a thin-film surface sensor were modeled and experimentally tested. Then, the inverse reconstruction of both simulated and experimental estimations by the EIT was compared and analyzed.
This study is classified as follows. Section 2 demonstrates the highly nonlinear issue statement and establishes the single-parameter homogenization functions in Section 3. In Section 4, we interpret the homogenization functions method to apply the inverse conductivity issues. Three numerical experiments of the inverse conductivity issue of highly nonlinear elliptic equations are displayed in Section 5. Finally, we illustrate the conclusions in Section 6.

2. Problem Statement and the Homogenization Functions Method

The purpose of this study is to construct an effective homogenization functions scheme that possesses the merit of easy numerical implementation and owns good flexibility when applied to the 3D inverse conductivity issue of highly nonlinear elliptic PDE for engineering and medical images in a cuboid Ω : = { ( x , y , z ) ( 0 , d ) × ( 0 , e ) × ( 0 , f ) } :
σ ( x , y , z ) [ v x x ( x , y , z ) + v y y ( x , y , z ) + v z z ( x , y , z ) ] + v x ( x , y , z ) σ x ( x , y , z ) + v y ( x , y , z ) σ y ( x , y , z ) + v z ( x , y , z ) σ z ( x , y , z ) + F [ v ( x , y , z ) ] = P ( x , y , z ) ,   ( x , y , z ) Ω ,  
v ( 0 , y , z ) = g 1 ( y , z ) ,       v ( d , y , z ) = g 2 ( y , z ) , v ( x , 0 , z ) = g 3 ( x , z ) ,       v ( x , e , z ) = g 4 ( x , z ) , v ( x , y , 0 ) = g 5 ( x , y ) ,       v ( x , y , f ) = g 6 ( x , y ) ,
where v(x, y, z) is electrical potential, F is a first high-order nonlinear operator, P(x, y, z) is a presented source function and σ(x, y, z) is an unknown electrical conductivity function to be retrieved under the over-specified Neumann data on the bottom:
v z ( x , y , 0 ) = g 7 ( x , y ) ,
Considering this issue, we assume that the Dirichlet data of σ(x, y, z) on boundaries can be estimated:
σ ( 0 , y , z ) = h 1 ( y , z ) ,       σ ( d , y , z ) = h 2 ( y , z ) , σ ( x , 0 , z ) = h 3 ( x , z ) ,       σ ( x , e , z ) = h 4 ( x , z ) , σ ( x , y , 0 ) = h 5 ( x , y ) ,       σ ( x , y , f ) = h 6 ( x , y ) ,
The issue in Equations (1)–(4) is fundamentally highly nonlinear since σ(x, y, z) and v(x, y, z) are both unknown functions. The current inverse coefficient issue is to retrieve the electrical conductivity coefficient function σ(x, y, z) inside the cuboid with the aid of the estimations of the potential on the entire boundary and the present flux on a partially accessible boundary in a finite cuboid.
Illustrating the concept of homogenization function [34], we contemplate a 3D boundary value problem (BVP):
S [ v ( x , y , z ) ] = P ( x , y , z ) ,   ( x , y , z ) Ω ,
v ( 0 , y , z ) = g 1 ( y , z ) ,       v ( d , y , z ) = g 2 ( y , z ) , v ( x , 0 , z ) = g 3 ( x , z ) ,       v ( x , e , z ) = g 4 ( x , z ) , v ( x , y , 0 ) = g 5 ( x , y ) ,       v ( x , y , f ) = g 6 ( x , y ) ,
where S is a second high-order nonlinear operator, P(x, y, z) is a presented source function, and gi, i = 1, …, 6 are presented boundary data.
Through
E 1 ( x , y , z ) = ( 1 x d ) [ g 1 ( y , z ) ( 1 y e ) g 3 ( 0 , z ) y e g 4 ( 0 , z ) ] + x d [ g 2 ( y , z ) ( 1 y e ) g 3 ( d , z ) y e g 4 ( d , z ) ] + ( 1 y e ) g 3 ( x , z ) + y e g 4 ( x , z ) .  
Knowing that
E 1 ( 0 , y , z ) = g 1 ( y , z ) ,       E 1 ( d , y , z ) = g 2 ( y , z ) .
Let
E ( x , y , z ) = E 1 ( x , y , z ) + ( 1 z f ) [ g 5 ( x , y ) E 1 ( x , y , 0 ) ] + z f [ g 6 ( x , y ) E 1 ( x , y , f ) ] .
We can employ the compatibility data of boundary conditions as follows:
g 5 ( 0 , y ) = E 1 ( 0 , y , 0 ) = g 1 ( y , 0 ) ,       g 6 ( 0 , y ) = E 1 ( 0 , y , f ) = g 1 ( y , f ) , g 5 ( d , y ) = E 1 ( d , y , 0 ) = g 2 ( y , 0 ) ,       g 6 ( d , y ) = E 1 ( d , y , f ) = g 2 ( y , f ) ,
it can be verified easily
E 1 ( 0 , y , z ) = g 1 ( y , z ) ,       E 1 ( d , y , z ) = g 2 ( y , z ) , E 1 ( x , 0 , z ) = g 3 ( x , z ) ,       E 1 ( x , e , z ) = g 4 ( x , z ) , E 1 ( x , y , 0 ) = g 5 ( x , y ) ,       E 1 ( x , y , f ) = g 6 ( x , y ) .
Thus, we can acquire a 3D homogenization function:
E ( x , y , z ) = ( 1 x d ) [ g 1 ( y , z ) ( 1 y e ) g 3 ( 0 , z ) y e g 4 ( 0 , z ) ] + x d [ g 2 ( y , z ) ( 1 y e ) g 3 ( d , z ) y e g 4 ( d , z ) ] + ( 1 y e ) g 3 ( x , z ) + y e g 4 ( x , z ) + ( 1 z f ) ( g 5 ( x , y ) ( 1 x d ) [ g 1 ( y , 0 ) ( 1 y e ) g 3 ( 0 , 0 ) y e g 4 ( 0 , 0 ) ] x d [ g 2 ( y , 0 ) ( 1 y e ) g 3 ( d , 0 ) y e g 4 ( d , 0 ) ] ( 1 y e ) g 3 ( x , 0 ) y e g 4 ( x , 0 ) ) + z f ( g 6 ( x , y ) ( 1 x d ) [ g 1 ( y , f ) ( 1 y e ) g 3 ( 0 , f ) y e g 4 ( 0 , f ) ] x d [ g 2 ( y , f ) ( 1 y e ) g 3 ( d , f ) y e g 4 ( d , f ) ] ( 1 y e ) g 3 ( x , f ) y e g 4 ( x , f ) ) .
Because of the attribute of E(x, y, z) and with the aid of the variable transformation from v(x, y, z) to w(x, y, z) = v(x, y, z)-E(x, y, z), we can transform the elemental 3D BVP with non-homogeneous boundary data to a one with homogeneous boundary data for w(x, y, z):
S [ w ( x , y , z ) ] = P ( x , y , z ) S [ E ( x , y , z ) ] ,   ( x , y , z ) Ω ,
w ( 0 , y , z ) = w ( d , y , z ) = w ( x , 0 , z ) = w ( x , e , z ) = w ( x , y , 0 ) = w ( x , y , f ) = 0 .
Hence, E(x, y, z) is a significant part of the homogenization function to demolish the non-zero boundary data.

3. A Single-Parameter Homogenization Function

We address the inverse conductivity issue of highly nonlinear elliptic PDE by utilizing the homogenization functions. For this aim, we bring in the normalized coordinates:
x ¯ = x d , y ¯ = y e , z ¯ = z f .
We can generalize Equation (12) by bringing in the mth order shape functions:
ε m ( x ¯ ) = x ¯ m , ε m ( y ¯ ) = y ¯ m , ε m ( z ¯ ) = z ¯ m ,
where
ε m ( 0 ) = 0 , ε m ( 1 ) = 1 , m N .
To supersede x/d, y/e, and z/f in Equation (12) by ε m ( x ¯ ) ,   ε m ( y ¯ )   and   ε m ( z ¯ ) , we can acquire the mth order homogenization function as follows:
E v ( m , x , y , z ) = [ 1 ε m ( x ¯ ) ] { g 1 ( y , z ) [ 1 ε m ( y ¯ ) ] g 3 ( 0 , z ) ε m ( y ¯ ) g 4 ( 0 , z ) } + ε m ( x ¯ ) { g 2 ( y , z ) [ 1 ε m ( y ¯ ) ] g 3 ( d , z ) ε m ( y ¯ ) g 4 ( d , z ) } + [ 1 ε m ( y ¯ ) ] g 3 ( x , z ) + ε m ( y ¯ ) g 4 ( x , z ) + [ 1 ε m ( z ¯ ) ] { g 5 ( x , y ) [ 1 ε m ( x ¯ ) ] { g 1 ( y , 0 ) [ 1 ε m ( y ¯ ) ] g 3 ( 0 , 0 ) ε m ( y ¯ ) g 4 ( 0 , 0 ) } ε m ( x ¯ ) { g 2 ( y , 0 ) [ 1 ε m ( y ¯ ) ] g 3 ( d , 0 ) ε m ( y ¯ ) g 4 ( d , 0 ) } [ 1 ε m ( y ¯ ) ] g 3 ( x , 0 ) ε m ( y ¯ ) g 4 ( x , 0 ) } + ε m ( z ¯ ) { g 6 ( x , y ) [ 1 ε m ( x ¯ ) ] { g 1 ( y , f ) [ 1 ε m ( y ¯ ) ] g 3 ( 0 , f ) ε m ( y ¯ ) g 4 ( 0 , f ) } ε m ( x ¯ ) { g 2 ( y , f ) [ 1 ε m ( y ¯ ) ] g 3 ( d , f ) ε m ( y ¯ ) g 4 ( d , f ) } [ 1 ε m ( y ¯ ) ] g 3 ( x , f ) ε m ( y ¯ ) g 4 ( x , f ) } .
By utilizing Equation (17), we can verify
E v ( m , 0 , y , z ) = g 1 ( y , z ) ,       E v ( m , d , y , z ) = g 2 ( y , z ) , E v ( m , x , 0 , z ) = g 3 ( x , z ) ,       E v ( m , x , e , z ) = g 4 ( x , z ) , E v ( m , x , y , 0 ) = g 5 ( x , y ) ,       E v ( m , x , y , f ) = g 6 ( x , y ) ,
Thus, we can choose Ev(k, x, y, z) as the foundations of v(x, y, z) and assume
v ( x , y , z ) = m = 1 p 1 a m E v ( m , x , y , z ) ,
where
m = 1 p 1 a m = 1 ,
we can acquire a simple solution of v(x, y, z) and hence the solution of σ(x, y, z) of the nonlinear inverse conductivity issue in the cuboid can be accomplished in the next part. Equation (21) is used to promise that v(x, y, z) in Equation (20) can gratify the boundary conditions in Equation (2).
The data g7(x, y) = vz(x, y, 0) is over-specified in Equation (3), and by taking the z-differential of Equation (20) and introducing z = 0, we can acquire
m = 1 p 1 a m E z v ( m , x , y , 0 ) = g 7 ( x , y ) ,
beginning with which we can produce linear equations by the collocation of points xk and yk, k = 1,…, q1:
m = 1 p 1 a m E z v ( m , x k , y k , 0 ) = g 7 ( x k , y k ) ,       k = 1 , , q 1 .
Dealing with the 2q1 + 1 linear Equations (21) and (23), we can decide the p1 coefficients am, m = 1,…, p1, and then Equation (17) is used to work out v(x, y, z) in the entire area.

4. Numerical Algorithm for Inverse Conductivity Issues

We can produce another mth order 3D homogenization function for σ(x, y, z) by
E σ ( m , x , y , z ) = [ 1 ε m ( x ¯ ) ] { h 1 ( y , z ) [ 1 ε m ( y ¯ ) ] h 3 ( 0 , z ) ε m ( y ¯ ) h 4 ( 0 , z ) } + ε m ( x ¯ ) { h 2 ( y , z ) [ 1 ε m ( y ¯ ) ] h 3 ( d , z ) ε m ( y ¯ ) h 4 ( d , z ) } + [ 1 ε m ( y ¯ ) ] h 3 ( x , z ) + ε m ( y ¯ ) h 4 ( x , z ) + [ 1 ε k ( z ¯ ) ] { h 5 ( x , y ) [ 1 ε k ( x ¯ ) ] { h 1 ( y , 0 ) [ 1 ε m ( y ¯ ) ] h 3 ( 0 , 0 ) ε m ( y ¯ ) h 4 ( 0 , 0 ) } ε m ( x ¯ ) { h 2 ( y , 0 ) [ 1 ε m ( y ¯ ) ] h 3 ( d , 0 ) ε m ( y ¯ ) h 4 ( d , 0 ) } [ 1 ε m ( y ¯ ) ] h 3 ( x , 0 ) ε m ( y ¯ ) h 4 ( x , 0 ) } + ε m ( z ¯ ) { h 6 ( x , y ) [ 1 ε m ( x ¯ ) ] { h 1 ( y , f ) [ 1 ε m ( y ¯ ) ] h 3 ( 0 , f ) ε m ( y ¯ ) h 4 ( 0 , f ) } ε m ( x ¯ ) { h 2 ( y , f ) [ 1 ε m ( y ¯ ) ] h 3 ( d , f ) ε m ( y ¯ ) h 4 ( d , f ) } [ 1 ε m ( y ¯ ) ] h 3 ( x , f ) ε m ( y ¯ ) h 4 ( x , f ) } .
which automatically gratifies the specified data in Equation (4).
By utilizing
σ ( x , y , z ) = m = 1 p 2 b m E σ ( m , x , y , z ) ,
m = 1 p 2 b m = 1 ,
back substituting v(x, y, z) in Equation (20) into Equation (1) and collocating q 2 3 points in the area Ω, we can produce q 2 3 linear equations:
m = 1 p 2 { [ v x ( x i , y j , z r ) E x σ ( m , x i , y j , z r ) + v y ( x i , y j , z r ) E y σ ( m , x i , y j , z r ) + v z ( x i , y j , z r ) E z σ ( m , x i , y j , z r ) + E σ ( m , x i , y j , z r ) [ v x x ( x i , y j , z r ) + v y y ( x i , y j , z r ) + v z z ( x i , y j , z r ) ] } b m = P ( x i , y j , z r ) F [ v ( x i , y j , z r ) ] , i , j , r = 1 , , q 2 .
Coping with q 2 3 + 1 linear Equations (26) and (27), we can decide the p2 coefficients bm, m = 1,…, p2, Equation (25) is used to retrieve σ(x, y, z) in the total area. This is a new algorithm of the 3D inverse conductivity issue of highly nonlinear elliptic PDE by employing the homogenization functions method (HFM). While Ev(m, x, y, z) plays the role of the foundations to broaden the solution σ(x, y, z), we do not require to resolve nonlinear equations.

5. Numerical Experiments

We will utilize the homogenization functions method to solve inverse conductivity issues of highly nonlinear elliptic PDEs through three examples from the second-order nonlinear issue to the highly nonlinear problem. We are interested in the stability of our scheme while the measured data are contaminated by random noise for varied problems. We can estimate the stability by increasing the varied random noise in the measured data:
g ^ 7 ( x i , y j ) = g 7 ( x i , y j ) + s R ( i , j ) ,
where R(i, j) are random numbers between [−1, 1], and s is the intensity of noise.
Besides, we deliberate the relative root-mean-square-error defined by
e m ( σ ) = i = 1 M j = 1 M r = 1 M [ σ m ( x i , y j , z r ) σ ( x i , y j , z r ) ] 2 i = 1 M j = 1 M r = 1 M σ 2 ( x i , y j , z r )
to estimate the accuracy of the retrieved conductivity function, in which we compare the numerically retrieved conductivity σm and the exact conductivity σ at M × M × M grid points (x, y, z), i, j, r = 1,…, M. We also employ the maximum absolute error (MAEm) to estimate the accuracy of the numerical result, which is defined by
MAE m : = max i , j , r = 1 , , M | σ m ( x i , y j , z r ) σ ( x i , y j , z r ) | .
All the computational algorithms were implemented to the Fortran code on the Microsoft Developer Studio platform in OS Windows 10 (64 bit) with i3-4160 3.60 GHz CPU and 16 GB memory.

5.1. Example 1

We deliberate the following 3D second-order nonlinear issue with a large domain
σ ( x , y , z ) [ v x x ( x , y , z ) + v y y ( x , y , z ) + v z z ( x , y , z ) ] + v x ( x , y , z ) σ x ( x , y , z ) + v y ( x , y , z ) σ y ( x , y , z ) + v z ( x , y , z ) σ z ( x , y , z ) + v 2 ( x , y , z ) = P ( x , y , z ) ,   ( x , y , z ) Ω ,  
with an exact solution
v ( x , y , z ) = ( x 5 ) 2 exp ( 2 y ) z 2 ,       σ ( x , y , z ) = ( x 5 ) 2 + y 2 + z 2 ,       ( x , y , z ) ( 0 , 12 ) × ( 0 , 3 ) × ( 0 , 9 ) .
We can obtain P(x, y, z) by inserting Equation (32) into Equation (31). Under the parameters: d = 12, e = 3, f = 9, p1 = 2, q1 = 3, p2 = 3, q1 = 2, and a noise with s = 0.4. We utilize the HFM to obtain the quite accurate solutions, as displayed in Figure 1, with the maximum absolute value of v being 3416.15 and the maximum absolute error (MAEm) of v being 4.69 × 10 2 , and the MAEm of σ is 9.53 × 10 5 . Note that the maximum result of σ is 139.0, and the result em(σ) = 6.86 × 10−7 is very small. For this instance, the CPU time is 3.18 s. From this instance, we find that these results are very accurate, time-saving, and stable even under the noise effect and large domain. Besides, we claim that the proposed algorithm can address the electrical conductivity coefficient function σ(x, y, z) accurately inside the cuboid with the aid of the estimations of the potential on the entire boundary and the present flux on a partially accessible boundary in a bounded cuboid.

5.2. Example 2

We ponder another 3D moderate nonlinear issue with the moderate domain as follows:
σ ( x , y , z ) [ v x x ( x , y , z ) + v y y ( x , y , z ) + v z z ( x , y , z ) ] + v x ( x , y , z ) σ x ( x , y , z ) + v y ( x , y , z ) σ y ( x , y , z ) + v z ( x , y , z ) σ z ( x , y , z ) v 3 ( x , y , z ) + 6 v x ( x , y , z ) v y ( x , y , z ) v z ( x , y , z ) = P ( x , y , z ) ,   ( x , y , z ) Ω ,
with an exact solution
v ( x , y , z ) = sin x + cos y + exp ( z ) + x 3 + y 2 + z 3 , σ ( x , y , z ) = 2 exp ( x + y + z ) ,   ( x , y , z ) ( 0 , 1 ) × ( 0 , 2 ) × ( 0 , 2 ) .
We can acquire P(x, y, z) by inserting Equation (34) into Equation (33). Under the parameters: d = 1, e = f = 2, p1 = 2, q1 = 5, p2 = 2, q2 = 4, and a noise with s = 0.4. We employ the HFM to acquire the accurate solutions, as demonstrated in Figure 2, with the maximum absolute value of v being 19.8, the maximum absolute error (MAEm) of v being 0.218, and the MAEm of σ is 0.237. Note that the maximum result of σ is 296.83, and the result em(σ) = 1.49 × 10−3 is rather small. For this case, the CPU time is only 2.83 s. From this example, we reveal that these numerical results are accurate, effective, and stable even under the noise effect and moderate area.

5.3. Example 3

Next, we contemplate the following 3D highly nonlinear problem
σ ( x , y , z ) [ v x x ( x , y , z ) + v y y ( x , y , z ) + v z z ( x , y , z ) ] + v x ( x , y , z ) σ x ( x , y , z ) + v y ( x , y , z ) σ y ( x , y , z ) + v z ( x , y , z ) σ z ( x , y , z ) + v 7 ( x , y , z ) + v x 2 ( x , y , z ) + v y 2 ( x , y , z ) + v z 2 ( x , y , z ) = P ( x , y , z ) ,   ( x , y , z ) Ω ,
with an exact solution
v ( x , y , z ) = cos x cos y sin z + x 2 y 2 cosh z , σ ( x , y , z ) = sinh z sin ( 2 x ) cosh y + cosh z cos ( 2 x ) sinh y ,   ( x , y , z ) ( 0 , 1 ) × ( 0 , 1 ) × ( 0 , 1 ) .
We can obtain P(x, y, z) by inserting Equation (36) into Equation (35). Under the parameters: d = e = f = 1, p1 = 2, q1 = 9, p2 = 3, q2 = 8, and a noise with s = 0.5, which is a large level of noise 28%. We use the HFM to procure the accurate solutions, as illustrated in Figure 3, with the maximum absolute value of v being 1.79, the maximum absolute error (MAEm) of v being 1.77 × 10 2 , and the MAEm of σ is 0.181 . Note that the maximum result of σ is 23.08, and the result em(σ) = 1.06 × 10−2 is quite small. For this numerical experiment, the CPU time is 4.52 s. From this numerical experiment, we exhibit that these numerical results are effective, accurate, and stable even under the large noise influence.

6. Conclusions

We addressed the 3D inverse conductivity issues of highly nonlinear elliptic equations by the homogenization functions method (HFM) in this article. The significant concept was constructed the single-parameter homogenization functions to produce practical foundations to broaden the solution and the conductivity function in the entire area. Although the 3D inverse conductivity issue is highly nonlinear, we still can effortlessly distinguish the unknown conductivity function by dealing with two sets of linear equations. Furthermore, to the author’s best knowledge, there is no study in the literature that the numerical algorithms for those three issues can supply more accurate results than the present results even under a large level of noise of 28%. The current method can be broadened to deal with the multi-dimensional inverse nonlinear steady-state PDEs, as demonstrated in Figure 4, and will be figured out in the near future.

Funding

This work was funded by National United University grant number [111-NUUPRJ-04] and the APC was funded by National United University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Calderón, A.P. On an Inverse Boundary Value Problem. In Seminar on Numerical Analysis and Its Application to Continuum Physics; Socidade Brasileira de Mathematica: Rio de Janeiro, Brazil, 1980; pp. 65–73. [Google Scholar]
  2. Borcea, L. Electrical impedance tomography. Inverse Probl. 2002, 18, R99–R136. [Google Scholar] [CrossRef]
  3. Adler, A.; Guardo, R. A neural network image reconstruction technique for electrical impedance tomography. IEEE Trans. Med. Imag. 1994, 13, 594–600. [Google Scholar] [CrossRef] [PubMed]
  4. Borcea, L.; Gray, G.A.; Zhang, Y. Variationally constrained numerical solution of electrical impedance tomography. Inverse Probl. 2003, 19, 1159–1184. [Google Scholar] [CrossRef] [Green Version]
  5. Brown, R.M.; Uhlmann, G.A. Uniqueness in the inverse conductivity problem for non-smooth conductivities in two dimensions. Commun. Partial. Differ. Equ. 1997, 22, 1009–1027. [Google Scholar] [CrossRef]
  6. Francini, E. Recoverning a complex coefficient in a planar domain from the Dirichlet-to-Neumann map. Inverse Probl. 2000, 16, 107–119. [Google Scholar] [CrossRef]
  7. Knowles, I. A variational algorithm for electrical impedance tomography. Inverse Probl. 1998, 14, 1513–1526. [Google Scholar] [CrossRef] [Green Version]
  8. Kohn, R.V.; McKenney, A. Numerical implementation of a variational method for electrical impedance tomography. Inverse Probl. 1990, 6, 389–414. [Google Scholar] [CrossRef]
  9. Levy, S.; Adam, D.; Bresler, Y. Electromagnetic impedance tomography (EMIT): A new method for impedance imaging. IEEE Trans. Med. Imaging 2002, 21, 676–687. [Google Scholar] [CrossRef]
  10. Meeson, S.; Killingback, A.L.T.; Blott, B.H. The dependence of EIT images on the assumed initial conductivity distribution: A study of pelvic imaging. Phys. Med. Biol. 1995, 40, 643–657. [Google Scholar] [CrossRef]
  11. Murai, T.; Kagawa, Y. Electrical impedance computed tomography based on a finite element model. IEEE Trans. Biomed. Eng. 1985, 32, 177–184. [Google Scholar] [CrossRef]
  12. Nachman, A.I. Reconstructions from boundary measurements. Ann. Math. 1988, 128, 531–576. [Google Scholar] [CrossRef]
  13. Nachman, A.I. Global uniqueness for a two-dimensional inverse boundary problem. Ann. Math. 1996, 143, 71–96. [Google Scholar] [CrossRef]
  14. Siltanen, S.; Mueller, J.; Isaacson, D. An implementation of the reconstruction algorithm of A Nachman for the 2D inverse conductivity problem. Inverse Probl. 2000, 16, 681–699. [Google Scholar] [CrossRef]
  15. Stasiak, M.; Sikora, J.; Filipowicz, S.F.; Nita, K. Principal component analysis and artificial neural network approach to electrical impedance tomography problems approximated by multi-region boundary element method. Eng. Anal. Bound. Elem. 2007, 31, 713–720. [Google Scholar] [CrossRef]
  16. Sylvester, J. A convergent layer stripping algorithm for radially symmetric impedance tomography problem. Commun. Partial. Differ. Equ. 1992, 17, 1955–1994. [Google Scholar]
  17. Wexler, A.; Fry, B.; Neuman, M.R. Impedance-computed tomography algorithm and system. Appl. Opt. 1985, 24, 3985–3992. [Google Scholar] [CrossRef]
  18. Yorkey, T.J.; Webster, J.G.; Tompkins, W.J. Comparing reconstruction algorithms for electrical impedance tomography. IEEE Trans. Biomed. Eng. 1987, 34, 843–852. [Google Scholar] [CrossRef]
  19. Zlochiver, S.; Rosenfeld, M.; Abboud, S. Induced-current electrical impedance tomography: A 2-D theoretical simulation. IEEE Trans. Med. Imaging 2003, 22, 1550–1560. [Google Scholar] [CrossRef]
  20. Zadehkoochak, M.; Hames, T.K.; Blott, B.H.; George, R.F. A transputer implemented algorithm for electrical impedance tomography. Clin. Phys. Physiol. Meas. 1990, 11, 223–230. [Google Scholar] [CrossRef]
  21. Calderón, A.P. On an inverse boundary value problem. J. Comput. Appl. Math. 2006, 25, 133–138. [Google Scholar] [CrossRef] [Green Version]
  22. Liu, C.-S. Solving the inverse problems of Laplace equation to determine the Robin coefficient/cracks’ position inside a disk. Comput. Model. Eng. Sci. 2009, 40, 1–28. [Google Scholar]
  23. 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]
  24. Liu, C.-S.; Liu, D. A homogenization boundary function method for determining inaccessible boundary of a rigid inclusion for the Poisson equation. Eng. Anal. Bound. Elem. 2018, 86, 56–63. [Google Scholar] [CrossRef]
  25. Liu, C.-S.; Wang, F. A meshless method for solving the nonlinear inverse Cauchy problem of elliptic type equation in a doubly-connected domain. Comput. Math. Appl. 2018, 76, 1837–1852. [Google Scholar] [CrossRef]
  26. Caro, P.; Ferreira, D.D.S.; Ruiz, A. Stability estimates for the Calderón problem with partial data. J. Diff. Eq. 2016, 260, 2457–2489. [Google Scholar] [CrossRef]
  27. Ferreira, D.D.S.; Kenig, C.E.; Salo, M.; Uhlmann, G. Limiting Carleman weights and anisotropic inverse problems. Invent. Math. 2009, 178, 119–171. [Google Scholar] [CrossRef] [Green Version]
  28. Piiroinen, P.; Simon, M. Probabilistic interpretation of the Calderón problem. Inv. Prob. Imag. 2017, 11, 553–575. [Google Scholar] [CrossRef]
  29. Alessandrini, G.; de Hoop, M.V.; Gaburro, R.; Sincich, E. Lipschitz stability for the electrostatic inverse boundary value problem with piecewise linear conductivities. J. Math. Pures Appl. 2017, 107, 638–664. [Google Scholar] [CrossRef] [Green Version]
  30. Muñoz, C.; Uhlmann, G. The Calderón problem for quasilinear elliptic equations. In Annales de l’Institut Henri Poincaré C, Analyse non linéaire; Elsevier Masson: Paris, France, 2020; Volume 37, pp. 1143–1166. [Google Scholar]
  31. De Moura, B.F.; Martins, M.F.; Palma, F.H.S.; da Silva, W.B.; Cabello, J.A.; Ramos, R. Nonstationary bubble shape determination in electrical impedance tomography combining Gauss–Newton optimization with particle filter. Measurement 2021, 186, 110216. [Google Scholar] [CrossRef]
  32. Duran, G.C.; Sato, A.K.; Ueda, E.K.; Takimoto, R.Y.; Martins, T.C.; Tsuzuki, M.S.G. Electrical impedance tomography image reconstruction using convolutional neural network with periodic padding. IFAC-PapersOnLine 2021, 54, 418–423. [Google Scholar] [CrossRef]
  33. Wagner, J.; Kralovec, C.; Schagerl, M. Evaluation of spatial strain distribution by elastoresistive thin-film sensors using 2D Electrical Impedance Tomography. Mater. Today Proc. 2022; in press. [Google Scholar] [CrossRef]
  34. Liu, C.-S.; Chang, C.-W. Solving the inverse conductivity problems of nonlinear elliptic equations by the superposition of homogenization functions method. Appl. Math. Lett. 2019, 94, 272–278. [Google Scholar] [CrossRef]
Figure 1. For example, 1 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with noise effect.
Figure 1. For example, 1 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with noise effect.
Symmetry 14 01044 g001aSymmetry 14 01044 g001b
Figure 2. For example 2 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with noise effect.
Figure 2. For example 2 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with noise effect.
Symmetry 14 01044 g002aSymmetry 14 01044 g002b
Figure 3. For example 3 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with large noise effect.
Figure 3. For example 3 of the 3D inverse conductivity issue of the highly nonlinear elliptic equation, (a) the maximum errors of v and the retrieved of σ, (b) exact solution, and (c) numerical solution with large noise effect.
Symmetry 14 01044 g003
Figure 4. Framework of multi-dimensional inverse nonlinear steady-state PDEs by utilizing the proposed HFM.
Figure 4. Framework of multi-dimensional inverse nonlinear steady-state PDEs by utilizing the proposed HFM.
Symmetry 14 01044 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chang, C.-W. A New Meshless Method for Solving 3D Inverse Conductivity Issues of Highly Nonlinear Elliptic Equations. Symmetry 2022, 14, 1044. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051044

AMA Style

Chang C-W. A New Meshless Method for Solving 3D Inverse Conductivity Issues of Highly Nonlinear Elliptic Equations. Symmetry. 2022; 14(5):1044. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051044

Chicago/Turabian Style

Chang, Chih-Wen. 2022. "A New Meshless Method for Solving 3D Inverse Conductivity Issues of Highly Nonlinear Elliptic Equations" Symmetry 14, no. 5: 1044. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051044

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