Abstract

The aim of this work is to study a semidiscrete Crank-Nicolson type scheme in order to approximate numerically the Dirichlet-to-Neumann semigroup. We construct an approximating family of operators for the Dirichlet-to-Neumann semigroup, which satisfies the assumptions of Chernoff’s product formula, and consequently the Crank-Nicolson scheme converges to the exact solution. Finally, we write a finite element scheme for the problem, and we illustrate this convergence by means of a FreeFem++ implementation.

1. Introduction

Let be a bounded smooth domain and let be a real-valued matrix function, which is known as the electrical conductivity matrix. The matrix is symmetric and smooth; that is, . Its positive eigenvalues are uniformly bounded; that is, there exists such that, for all , one has , where is the Euclidean norm of .

In this paper, we denote by , with boundary space . It is well known that, for any , the following Dirichlet problem,has a unique solution in , which is called the -harmonic lifting of ; see [1] for more details.

Define the action of the Dirichlet-to-Neumann operator on as the outward normal derivative of on the boundary; that is, where is the outer normal vector to at .

This operator is defined on its domain: The Dirichlet-to-Neumann semigroup is the trace of the following boundary value problem: In [1], Lax has introduced a semigroup in , where is the unit ball and is the identity matrix. Let be the solution of the problem:where is the trace operator on . The Lax semigroup is defined by where is the lifting operator and .

Therefore, if , then, for any , It has been shown that the Lax semigroup is actually the Dirichlet-to-Neumann semigroup.

In [2], it is shown that Lax’s representation cannot be generalized if is not the unit ball. This was a motivation for many authors to use Chernoff’s formula in order to approximate the solution.

In [3], Emamirad and Shariftabar proved that the Euler Explicit scheme converges to the solution of the (BVP).

Later, Cherif et al. proposed, in [4], the following Implicit Euler Scheme: which converges to the solution of the (BVP).

After studying the implicit and explicit schemes, by the same strategy, we are interested in studying the convergence of the Crank-Nicolson type scheme. This paper is divided into 4 sections, as follows.

First, in Section 2, we introduce the Crank-Nicolson type scheme for evolution equations.

In Section 3, we prove the convergence of a Crank-Nicholson type scheme for the (BVP), by constructing an approximating family, and use Chernoff’s product formula in order to prove the convergence in .

Finally, in Section 4, we present a numerical implementation of the Crank-Nicolson type scheme in the computational framework of the finite element method.

2. The Crank-Nicolson Scheme

In numerical analysis, the Crank-Nicolson method is a finite difference method used for numerically solving the heat equation and similar partial differential equations. It is a second-order method in time and it is numerically stable. The method was developed by John Crank and Phyllis Nicolson in the mid-20th century.

Let be a -accretive operator on the Hilbert space defined as .

Consider the evolution problem: We approximate the time derivative by using the finite difference method as follows: Crank-Nicolson introduced the scheme: The Chernoff product formula holds not only in the case of semigroup of contractions but also if the semigroup is bounded, so that the Crank-Nicholson scheme for the (BVP) can be written as

3. Convergence of the Crank-Nicholson Type Scheme

In order to show the convergence of the (CNS), let us recall the following Chernoff theorem.

Theorem 1 (Chernoff’s product formula). Let X be a Banach space, real or complex, and let be a family of contractions on X with .
Assume that there exists which generates a -semigroup of contractions and that for all there exists Then, for each in , one has .

Remark 2. Note that Chernoff’s product formula holds not only in the case of semigroup of contractions but also if the semigroup is bounded (see [5]).

This theorem was used to prove the convergence of both implicit and explicit schemes.

Now, in order to prove the convergence of the Crank-Nicolson scheme, we will show that the operator , defined below, satisfies the 3 conditions of Chernoff’s theorem.

For any , we have

Using in (CNS) the approximation given in (15), we getTaking , we getwhich can be written as where Now, we define the family as the inverse of in the following sense: and we show that satisfies all the assumptions of Chernoff’s product formula.

(i) . For any , we haveso, .

Moreover, so,Using (6) and (7), we get Therefore, .

(ii) . The derivative of with respect to is At point , we have .

Similarly, the derivative of with respect to is At point , we have .

But , so , and then .

Now, the derivative of with respect to is At , .

Therefore,

(iii) Is a Contraction. In fact,but is an upper bound of the normal derivative , so that this term is positive according to Hopf’s lemma, and, consequently, which implies that Now Thus, is of contraction.

Similarly is also of contraction.

Finally, Therefore, is of contraction.

So, we have proved that satisfies all the assumptions of Chernoff’s product formula, and, consequently, the Crank-Nicholson scheme converges to its exact solution.

Theorem 3. The operator satisfies the Chernoff conditions; hence, the Crank-Nicolson method is convergent.

4. Numerical Results

In the last section of this paper, we present the resolution of (CNS) in the particular case where is the unit matrix and is the open unit disk. The point in dealing with this case is that the exact solution is known to be , where is the solution of the Laplace equation , with on the boundary ( is any regular function). Thus, we can compare this exact solution with the approximate solution obtained via our Crank-Nicolson scheme. For this test, we choose . The variational formulation of the (CNS) can be derived by multiplying both sides of the problem by a test function : Using Green’s formula, we get Multiplying by , we get which is of the form where the bilinear form is the linear form.

Using this variational formulation, the problem can be approximate via -finite elements implemented in Freefem++.

The actual algorithm is as follows.

Choose an initial and, for a given , where , we have the following:

for do(i)(ii)compute via the variational problem(iii)

end do.

As an illustration, let us represent the exact (Figure 1) and the approximate solution (Figure 2), the -error between the exact and the approximate solutions (Figure 3), and the decrement of this -error according to the fineness of the meshing (Figure 4).

5. Conclusion

As a conclusion, we have proved the convergence of the (CNS) by showing that the constructed family of operators satisfies the assumptions of Chernoff’s product formula and consequently approximated numerically the Dirichlet-to-Neumann semigroup.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.