On the Coupled of NBEM and FEM for an Anisotropic Quasilinear Problem in Elongated Domains ()
1. Introduction
Based on the Green’s function and Green’s formula, natural boundary element method (NBEM) reduces the boundary value problem of partial differential equation into a hypersingular integral equation on the boundary, and then solves the latter numerically [1,2]. It has advantages over the usual boundary reduction methods: such as the diminution of the number of space dimensions by 1, the conservation of energy functional, the preservation of self-adjointness and coerciveness. But it also has evident limitations, it’s difficult to obtain Green’s functions for solving problem in general domains. Therefore, the coupling of NBEM which is also called artificial boundary condition [3,4] or DtN method [5,6] and finite element method (FEM) [2] is useful and necessary for general cases.
The standard procedure of the coupling method can be described as follows. We introduce an artificial boundary to divide the original domain into two subregions, a bounded inner region and an unbounded one with a special boundary, such as circle, ellipse, and spherical surface, on which the boundary element method and finite element method are used respectively. This technique has been used to solve many linear problems [1,2,4-6] and it has also been successfully generalized to solve nonlinear boundary value problems [7-9] or quasilinear problems [3,10,11]. The problems were discussed in [3,10,11] take circle as artificial boundary, but for the problems with elongated domains, an elliptic boundary that leads to a smaller computational domain is obviously better than the circle one. The purpose of the paper is to study the coupling of NBEM and FEM to solve the anisotropic quasilinear problems with an elliptic artificial boundary.
Let
be a elongated, bounded and simple connected domain in
with sufficiently smooth boundary
.
. We consider the numerical solution of the exterior anisotropic quasilinear problem
(1.1)
With
or
,
,
and
are given functions which will be ranked as below. Following [3,12], suppose that the given function
satisfies
(1.2)
where two positive constants
, and
(1.3)
with a constant
. We also assume that
,
are continuous. In the following, we suppose that the function
has compact support, i.e., there exists a constant
, such that
(1.4)
We also assume that
(1.5)
Now, we introduce an elliptic artificial boundary
![](https://www.scirp.org/html/12-1100100\1d474555-b1e7-4dc5-a900-96233dd4143e.jpg)
divide
into two regions, a bounded domain
and an unbounded domain
with elliptic artificial boundary. Then the problem (1.1) can be rewritten in the coupled form:
(1.6)
(1.7)
(1.8)
where
is the unit exterior normal vector on
. Particularly, when
which is independent of
and
, [13-15] have obtained the natural integral equation. We introduce the so-called Kirichhoff transformation [16]
(1.9)
then we have
(1.10)
and
(1.11)
From equation (1.7) we have that
satisfies the following problem
(1.12)
The rest of the paper is organized as follows. In section 2, we obtain the natural integral equation for elliptic unbounded domain cases. In section 3, we give the equivalent variational problems and the finite element approximations. The reduced problem’s well-posedness, the convergence results and error estimate are also discussed. At last, in section 4, we present some numerical examples to illuminate the efficiency and feasibility of our method.
2. Natural Boundary Reduction
In this section, by virtue of the Poisson integral formula and natural integral equation for the linear problem, we shall obtain the corresponding results for the quasilinear problem in
. For this purpose, we need to discuss some properties between elliptic coordinates
and Cartesian coordinates
first. The relationship between the two coordinates can be expressed as below
(2.1)
where
,
,
. Following from [15], we have
Theorem 2.1 The transformation between elliptic coordinates and Cartesian coordinates (2.1) possesses the following property.
1) The Jacobi determinant of equation (2.1) is
(2.2)
if and only if
;
2)
(2.3)
for
;
3) For the exterior domain ![](https://www.scirp.org/html/12-1100100\a0993046-0230-4b14-bc55-736af30b1f15.jpg)
(2.4)
where
refers to the unit exterior normal vector on
(regarded as the inner boundary of
).
Proof The conclusions 1 and 2 can be obtained by direct computation. And 3 follows from the property
![](https://www.scirp.org/html/12-1100100\e5bc8500-8862-446e-88c7-a5a50d86a417.jpg)
2.1. Natural Integral Equation for α = β = 1
Assume that
is the solution of the problem (1.12), and the value
is given, namely
![](https://www.scirp.org/html/12-1100100\eff56d76-33bf-48c6-a005-655a96b4b6a8.jpg)
Then based on the natural boundary reduction, there are the Poisson integral formulas
(2.5)
or
(2.6)
And the natural integral equation
(2.7)
or
(2.8)
the definition of
can be found in the following. The Poisson integral formulas (2.5) and (2.6) and the natural integral equations (2.7) and (2.8) can also be expressed in the Fourier series forms
(2.9)
(2.10)
where
,
and
![](https://www.scirp.org/html/12-1100100\7bf0e5ad-8bd5-4c39-815b-95a7709c2e0b.jpg)
From (1.10), we obtain
(2.11)
Combining (1.9), (2.10) and (2.11), we get the exact artificial boundary condition of
on
,
(2.12)
where
,
,
. Then by (1.6)-(1.8) and (2.12), the original problem with
confines in
can be defined as follows
(2.13)
Therefore, the solution of problem (2.13) is the solution of the problem (1.1) with
confining in the bounded domain
.
2.2. Natural Integral Equation for β > α > 0
Now we assume that
can be expressed in the form:
, with
. We also assume that
is the solution of the problem
(1.12), and the value
is given, namely
![](https://www.scirp.org/html/12-1100100\afccba8a-7abf-4726-9c9f-c521116e5fad.jpg)
Let
,
, then the boundary
is changed by the elliptic boundary
the unit exterior normal vector on
is
![](https://www.scirp.org/html/12-1100100\33254048-8753-481b-bf1d-b14f8f469b77.jpg)
By the above transformation, the problem (1.12) changes into
(2.14)
This is the right problem we talked in section 2.1. Similar with equation (2.1), we let
![](https://www.scirp.org/html/12-1100100\5a40de20-8288-4b6a-92e9-8aa84443b4c2.jpg)
where
![](https://www.scirp.org/html/12-1100100\4583a35a-97dd-402c-84a8-f91b40cf6938.jpg)
Then just the same as the problem discussed in Section 2.1, we have the natural integral equation on ![](https://www.scirp.org/html/12-1100100\7e29a235-bef8-4bcb-9bb4-cc82a6d0a8fa.jpg)
(2.15)
where
is the unit exterior normal vector on
. From (1.11), we obtain
(2.16)
Combining (1.9), (2.15) and (2.16), we obtain the exact artificial boundary condition of
on
,
(2.17)
Then by (1.6)-(1.8) and (2.17), the original problem with
confines in
can be defined as follows
(2.18)
Therefore, the solution of problem (2.18) is the solution of the problem (1.1) with
confining in the bounded domain
.
3. Variational Problem and Finite Element Approximation
3.1. The Equivalent Variational Problems
Now we consider the problems (2.13) and (2.18). We shall use
denoting the standard Sobolev spaces,
and
referring to the corresponding norms and semi-norms. Especially, we define
,
and
. Let us introduce the space
(3.1)
and the corresponding norms
![](https://www.scirp.org/html/12-1100100\32aa4b35-edf9-4496-ac2f-e72528d7325b.jpg)
The boundary value problems (2.13) and (2.18) are equivalent to the following variational problem
(3.2)
where
(3.3)
(3.4)
where
is gotten from Green’s formula, (2.7) and (2.8) with
and (2.17) with
.
(3.5)
For any real number
, we let
(3.6)
with
and
,
.
Lemma 3.1 There exists a constant
which has different meaning in different place and is related to
and
, such that
![](https://www.scirp.org/html/12-1100100\5e382fc4-4d12-4a25-b19d-57164ebf3d6f.jpg)
![](https://www.scirp.org/html/12-1100100\e587d8ec-9015-4d55-8dd2-6fa6a5385789.jpg)
In practice, we need to truncate the series in (2.12) and (2.17) for some nonnegative integer
, that is
(3.7)
with
(3.8)
when
, and
(3.9)
when
. So we only use the summation of the first
terms in (2.13) and (2.18). We will consider the following approximate problem
(3.10)
(3.11)
Both (3.10) and (3.11) are equivalent to the following variational problem
(3.12)
where
(3.13)
Similar with Lemma 3.1, we have
Lemma 3.2 There exists a constant
which has different meaning in different place, such that
![](https://www.scirp.org/html/12-1100100\6b61d11a-b9b0-4610-889f-941d1ff60458.jpg)
3.2. Finite Element Approximation
Divide the arc
into
parts and take a finite element subdivision in
such that their nodes on
are coincident. That is, we make a regular and quasiuniform triangulation
on
, such that
(3.14)
with
is a (curved) triangle;
the maximum side of the triangles. Let
(3.15)
Then the approximate problem of (3.12) can be written as
(3.16)
Some existence and uniqueness results for this type of problem are given in [12,17,18] under some conditions on the coefficients
, so by the constraint conditions
(1.2) and (1.3) we have Lemma 3.3 Problems (3.2), (3.12) and (3.16) have unique solvability.
3.2.1. Convergence Theorems
In this section, we obtain the convergence result of the problems discussed above. We let
and
be the solution of problems (3.2), (3.12), (3.16) respectively. We also assume that
(3.17)
And we require that
is a family of finitedimensional subspaces of
, which satisfies for any
(3.18)
(3.19)
where
is independent of
.
The continuous piecewise polynomial spaces, such as (3.15), satisfy the condition (3.17). And if we let
, where
is the interpolation operator, then by (3.19), we have
![](https://www.scirp.org/html/12-1100100\bad787bf-5e02-470b-bbd5-9fedc97fb19b.jpg)
And we can also obtain the following result.
Lemma 3.4
.
Proof From the (1.2), (3.12) and Lemma 3.2, we have
![](https://www.scirp.org/html/12-1100100\ffd72ee6-00a2-41c1-bbf8-b3fd75facf70.jpg)
For
, we assume that
![](https://www.scirp.org/html/12-1100100\806ed442-aa09-4ece-9d54-7f405c1a623d.jpg)
![](https://www.scirp.org/html/12-1100100\1c26411f-2e52-4ce7-bdb6-f59944366abb.jpg)
with
and
.
Then we have
![](https://www.scirp.org/html/12-1100100\5b4a0cd2-e61a-4465-9cb0-b2c967f4598d.jpg)
From
, we obtain that
is bounded in
. Therefore, there exists a subsequence
such that
. Then similar with the proof of Lemma 3.4 of [3], we obtain
![](https://www.scirp.org/html/12-1100100\be14b40d-f425-4f30-a41d-c77bead75674.jpg)
By the above lemmas, we get the following convergence result.
Theorem 3.1 Let
, and the assumptions (3.17)-( 3.19) be satisfied, then we have
(3.20)
3.2.2. Error Analysis
In the following, we shall get error estimates for the approximate solution obtained from a FEM-NBEM discrete scheme in the cases
. We assume that the solution
of problem (1.1) satisfies
![](https://www.scirp.org/html/12-1100100\d1d3eb59-de86-4d78-9e02-95d5bb2129c4.jpg)
For simplicity let us define the following notation
![](https://www.scirp.org/html/12-1100100\2d109906-0dc8-4353-a1bd-09c83d043af9.jpg)
![](https://www.scirp.org/html/12-1100100\096b958c-ada6-4f1a-9d4f-1028c3e74d0f.jpg)
![](https://www.scirp.org/html/12-1100100\ef6a5ff9-5b0e-4b8d-b282-1b9e9a35376c.jpg)
Then (3.2), (3.12), (3.16) can be replaced by the corresponding simple forms respectively.
Now we introduce the bilinear form
and
defined by
![](https://www.scirp.org/html/12-1100100\7579a280-f8b7-473a-ad53-fbc6be549cd4.jpg)
![](https://www.scirp.org/html/12-1100100\39171db0-b354-4142-84c1-199202228406.jpg)
Let
be the dual space of
. By (1.2) and continuity of
, we obtain that
is bounded in
. Then there exists an operator
such that
(3.21)
Similar with the proof of [10], we have the lemma as follows
Lemma 3.5 The bilinear form
defined by
satisfies the following inequality
(3.22)
where
is a sufficient large constant and
.
We assume that
(3.23)
Let
be the canonical injection. Since
is compactly embedded in
, we have that the operator
defined by
is also compact. By (3.21) and (3.23) and
satisfies the property of
, we obtain that
is an isomorphism.
By the conditions (3.2), (3.22), (3.23) and Theorem 10.1.2 of [20], one can get that there exists
, such that the following inequality is satisfied
(3.24)
for some constant
independent of
.
We define the Galerkin projection with respect to
, ![](https://www.scirp.org/html/12-1100100\b4a86b70-30a4-4d3f-818e-767ce29eb062.jpg)
![](https://www.scirp.org/html/12-1100100\5d78d24c-326d-492f-a076-f3278542edfa.jpg)
Then the operator
satisfies
(3.25)
We define the set
![](https://www.scirp.org/html/12-1100100\b0aa0c7e-f841-4e27-b498-8d64fd58c7b2.jpg)
Lemma 3.6
is a solution of (3.14) if and only if the following equation
![](https://www.scirp.org/html/12-1100100\fd22b9ef-f456-4e46-bc32-88d36cb9f17a.jpg)
holds, where
![](https://www.scirp.org/html/12-1100100\fc5a0532-18a2-4f8a-b501-b607af829f54.jpg)
with
,
.
Proof. Let
, then by
![](https://www.scirp.org/html/12-1100100\c6b02ebd-0cbb-4665-85f7-c2499f5c4604.jpg)
and
![](https://www.scirp.org/html/12-1100100\4224f604-93d3-4e57-b94d-a4be59a92e5a.jpg)
We can get the desired result.
Let
. Then following [10,11], we have
Lemma 3.7 There exists a positive constant C independent of h, such that
![](https://www.scirp.org/html/12-1100100\5123ccc0-cdc4-4ab5-b92a-6fec7e3c6427.jpg)
We also have the following result.
Lemma 3.8 ![](https://www.scirp.org/html/12-1100100\50eec7ec-3753-42d4-8918-ffecd18f1104.jpg)
Proof For any
, we only need to show that
.
![](https://www.scirp.org/html/12-1100100\00ce74f3-1a7c-41d0-a28f-22a3c9a88ea5.jpg)
![](https://www.scirp.org/html/12-1100100\c7b68bee-c591-42a4-9550-c2afd2ddac1d.jpg)
![](https://www.scirp.org/html/12-1100100\2aadbcdd-df30-4633-a97d-2df6ff2adf28.jpg)
Since
is regular and quasi-uniform, referring to [19], we obtain the following inverse inequality
![](https://www.scirp.org/html/12-1100100\5059d90c-aa71-4f8b-a7f0-614a4842dd5e.jpg)
Combining the above inequalities with the definition of
and (3.26), we obtain
![](https://www.scirp.org/html/12-1100100\e22049ac-552a-4809-ad1d-50cf6eaa60e1.jpg)
By the definition of
, we get the desired result.
Theorem 3.2 Assume
be the solution of (1), with
,
, and we also assume that
and
satisfies (3.23). With sufficiently small
, the finite element equation (3.16) has the approximate solution
such that
![](https://www.scirp.org/html/12-1100100\d2f13d7e-5ee8-4c55-9d75-ae61f2e6dbce.jpg)
where C is a constant independent of h and N.
Proof Firstly, for any
, we have
![](https://www.scirp.org/html/12-1100100\638987f3-1ac1-4d8e-b187-a08a916574d7.jpg)
Then by (3.12), we have
![](https://www.scirp.org/html/12-1100100\d95a61fe-0c35-46cb-8334-f0d502350cd7.jpg)
Let
, we have
![](https://www.scirp.org/html/12-1100100\d1ddc151-6de9-4739-9f27-29e87348ff57.jpg)
From (3.2), (3.22), (3.23) and [20], we obtain
(3.26)
We denote a nonlinear mapping
, which satisfies that for any given
,
is the unique solution of
(3.27)
Therefore, we have
![](https://www.scirp.org/html/12-1100100\ba4ca3c2-3492-4f4e-9870-e1ba3fdad867.jpg)
Combining the above equation with (3.25), we obtain the operator
is continuous, i.e.,
![](https://www.scirp.org/html/12-1100100\9a5e0741-0572-4504-87ae-b1c7e4b56199.jpg)
Next, we assume that
, then by Lemma 3.8, we have that
. By the definition of
, (3.27) can be rewritten as
![](https://www.scirp.org/html/12-1100100\94d60cfe-56b7-47a1-bf66-23fc87b43e84.jpg)
Then, from (3.24), Lemma 3.6 and Lemma 3.7, we have
![](https://www.scirp.org/html/12-1100100\e38a709d-f436-4e19-9a54-8c15e42ee0d0.jpg)
This implies that
. And since
is also continuous, following from Brouwer’s fixed theorem, one can obtain that there exists
, such that
. From
Lemma 3.6, we deduce that
is the solution of (3.16). What’s more, by (3.25) and the fact
, we obtain
(3.28)
Combining (3.26) with (3.28), one can obtain
![](https://www.scirp.org/html/12-1100100\3acaa508-6068-43ec-9114-06219d7b7083.jpg)
This completes the proof.
4. Numerical Examples
In this section, we shall give some examples to confirm our theoretical results. In the following, we choose the finite element space as given in (3.16). For simplicity, we let
![](https://www.scirp.org/html/12-1100100\36863a52-ff7c-459a-a72c-c93b726ad36b.jpg)
Example 4.1 We assume the exterior domain
with elliptical boundary
![](https://www.scirp.org/html/12-1100100\d9ede1ca-7640-47b5-bb56-72549947bf5c.jpg)
Now we consider the problem
(4.1)
when
,
and
.
The exact solution of Example 4.1 is
.
The numerical results are given in Figures 1 and 2 and Table 1.
Example 4.2 Similar with Example 4.1,
and
are replaced by
![](https://www.scirp.org/html/12-1100100\b25e2d87-f392-4df2-9e73-9c185af1d504.jpg)
and
respectively.
Figure 1. Example 4.1 with N = 16, µ1 = 1.7.
Table 1. The errors with N = 16 for Example 4.1.
The exact solution of Example 4.2 is
.
The numerical results are given in Figure 3, Figure 4 and Table 2.
Figure 3. Example 4.2 with N = 6, µ1 = 1.0.
Table 2. The errors with N = 6 for Example 4.2.
Example 4.3 We assume the exterior domain
with elliptical boundary
![](https://www.scirp.org/html/12-1100100\d4ab7209-16bb-4a21-99ce-a9b9a912a512.jpg)
Now we consider the problem
(4.2)
when
,
and
![](https://www.scirp.org/html/12-1100100\c3e6666c-cf4e-40ca-9a59-7ef626c4c6ed.jpg)
The exact solution of Example 4.3 is
.
The numerical results are given in Figures 5 and 6 and Table 3.
Example 4.4 Similar with Example 4.3,
and
are replaced by
Figure 5. Example 4.3 with N = 10, ε = 0.005.
![](https://www.scirp.org/html/12-1100100\72ab6dcb-ee26-4f77-a00b-95785723f169.jpg)
and
respectively. And we take
![](https://www.scirp.org/html/12-1100100\958fdf3f-0e3f-4606-a20d-de23c6549c7c.jpg)
The exact solution of Example 4.3 is
.
The numerical results are given in Figures 7 and 8 and Table 4.
From the numerical results, one obtains that the numerical errors can be affected by the order of artificial boundary condition, the mesh of the domain and the location of the artificial boundary, and it can be reduced by increasing the order of the artificial boundary condition and refining the mesh. What’s more, the convergence rate of anisotropic problems can also be affected by the choice of
as it is shown in Tables 3 and 4. The numerical results are in agreement with the error analysis we obtain and show the efficiency of the coupling method.
Table 3. The errors with N = 10 for Example 4.3.
Figure 7. Example 4.4 with N = 5, ε = 0.05.
Table 4. The errors with N = 10 for Example 4.4.
5. Acknowledgements
This research was partly supported by the National Natural Science Foundation of China (contact/grant number: 10871100,11071109). The computations in this paper have been carried out in the Jiangsu Provincial Key Laboratory for Numerical Simulation of Large Scale Complex Systems. The authors express their thanks to them.
NOTES