**Applied Mathematics** Vol.4 No.5A(2013), Article ID:31909,8 pages DOI:10.4236/am.2013.45A009

Convergence and Error of Some Numerical Methods for Solving a Convection-Diffusion Problem

Applied Mathematics Department, Babes-Bolyai University, Cluj Napoca, Romania

Email: djantai@yahoo.com, ioana@math.ubbcluj.ro, blaga@math.ubbcluj.ro

Copyright © 2013 Gabriela Nut et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received March 5, 2013; revised April 17, 2013; accepted April 24, 2013

**Keywords:** Time Dependent Convection Diffusion; Multigrid Method; Finite Element and Finite Differences Discretization

ABSTRACT

We use the local Fourier analysis to determine the properties of the multigrid method when used in modeling the skin penetration of a drug. The analyses of these properties can be very in designing an efficient structure of the multigrid method and in comparing the element and finite difference discretization techniques. After the theoretical results obtained, we also present some numerical results for a problem for which the solution is known.

1. Introduction

In this paper we present an eoretical study of the smoothing, convergence and error reduction properties of the multigrid method for a time dependent convection diffusion equation. This is an equation that arises in the mathematical modeling of many physical phenomena, which makes the efficient numerical solution very important.

The equation studied here models the transport of molecules through the layers of the skin, until it reacheas the blood stream. The parameters used for the diffusion coefficients are smaller by several order than those of the convection, thus the equation is a convection dominated one.

The discretization of the differential equation is realized by two different methods: the finite difference method [1] with Euler backward discretization and the Galerkin finite element method [2,3].

The system obtained after the discretization process is solved using the multigrid method. This method was first introduced by Fedorenko [4,5]. The first practical results and efficiency of the method were given by Brandt [6,7].

The theory of multigrid convergence is well established for the Poisson equation [8-10]. In more recent articles the convergence has been studied for the convectiondiffusion equation [11,12].

The novelty in this paper is that we study the smoothing factor, asimptotic convergence factor and the error reduction factor of the multigrid method for a time dependent convection-diffusion equation, on a domain comprising three layers with different physical properties. The analyse is performed using the local Fourier techniques [9,13] which represents a good tool for constructing efficient multigrid methods for a given differential equation.

We also determined the error obtained for a given solution of the model problem, using the multigrid method on different number of grid levels, for both discretization methods mentioned above.

2. Mathematical Model

represents the concentration of the substance transported through the blood stream, is the vector of convection coefficients and is the vector of diffusion coefficients.

The substance which is transported through the skin is applied at the surface on an area with a radius of a few centimeters. The depth to which the active substance is transported by the diffusion and convection process is of the order of nanometers, thus much smaller than the radius of the surface where it is appplied. As a consequence, the problem can be reduced to the unidimensional case. From this point on, the variable will represent the depth where the concentration is to be calculated, and the vectors and are the coefficients in different layers of the skin.

The concentration of the substance applied on the skin is known, and the amount of it is sufficiently large to be constant at any moment of time:

(1)

this being the initial condition of the problem.

On the frontiers between the skin layers the law of flux conservation gives:

(2)

is the number of layers where the diffusion takes place and:

After the discretization process, the system obtained from the Equation (1) has the form:

(3)

.

If the finite elements method [14] is used, the weak formulation of Equation (1) gives:

is a function that has a derivative of order 1 and is square-integrable on. The functions and are approximated using the continuous functions, being the number of interior points of the grid on level, through the relations:. Replacing the functions

and with these approximates and using the standard integration-by-parts formula the equation becomes:

or:

(4)

Computing the integrals from (4) for

the coefficients in the system (3) are:

(5)

For the finite differences method using the explicit backward Euler scheme:

the coefficients for the system (3) will be:

(6)

In the nodes that are on the frontiers between different layers of the skin, the law of flux conservation (2) becomes:

In these points the system (3) has the coefficients:

(7)

3. The Components of the Multigrid Method

For the components of the multigrid method we give in the following the matrices associated to their operators, needed for the local Fourier analysis of the convergence. The essential property used by this method is the fact that the discretiztion of the problem leads to a system that has the eigenvectors equals to the Fourier modes and when the multigrid components have a block structure when computed in the Fourier basis, the analysis of the multigrid method is reduced to the one of diagonal blocks of small size.

3.1. The Matrix of an Operator

If is an operator that can be described by a difference stencil:

meaning that:

then the functions are the eigen functions of:

(8)

and:

are the eigenvalues of.

As, it is sufficient to take.

If, the set of low frequencies, then:

is the set of high frequencies.

Using the above notations, for an arbitrary function

, the operator applied to gives:

and represents the matrix associated with the operator.

3.2. The Operator of the Discretized System

from system (3) applied to a function will give:

where.

Thus the matrix of is:

(9)

3.3. Preand PostSmoother

The Gauss-Seidel red-black method is used before and after the coarse grid correction, and reduces well the high frequency components of the error. The smoothing operator has two components of Jacobi type:

In the relations above:

(10)

As:

the matrix of the smoother will be:

(11)

3.4. Restriction of the Defect

Full-weighting restriction is used as a fine to coarse grid transfer operator :

As

(12)

the restriction operator applied to the function will give:

(13)

with.

Thus the restriction operator has the matrix:

(14)

3.5. Solution on the Coarse Grid

In the two-grid method, the exact solution on the coarse grid is required. After the restriction of the defect, the function on the grid has the form

(12,13). Thus:

wherefrom the matrix of the operator is:

(15)

3.6. Prolongation

The coarse to fine interpolation operator used is the bilinear interpolation:

From this relation it follows that:

and the matrix of the prolongation operator is:

(16)

3.7. Two-Grid Operator

The multigrid method [8,15] is a combination between a relaxation method (that reduces very well the high frequency components of the error, but is slowly convergent because of the low frequency components) and the coarse grid correction (which has complementary properties to the smoother).

The matrices from (9), (11), (14), (15) and (16) are used to create the two-grid operator for the multigrid method:

(17)

where:

(18)

is the matrix of the coarse grid correction operator.

It has been proven [9] that it is sufficient to derive the convergence properties for the two-grid method and the multigrid method will have similar properties. As a consequence, the following factors are defined for the two-grid operator.

Asimptotic convergence factor

(19)

Error reduction factor

(20)

Here: denotes the spectral norm associated with the Euclidian vector norm in, and

(21)

Smoothing factor

(22)

where is the spectral radius of the matrix

.

, introduced in [9], is an "ideal" coarse grid operator that anihilates the low frequency error components and leaves the high frequency components unchanged:

and:

Since,

(23)

4. Local Fourier Analysis Results for the Studied Problem

4.1. Smoothing Factor

If then the matrix (11) of the smoother after steps is:

and has the eigenvalues:

and being given in (10),

For, the eigen value attains its maximum absolute value:

(24)

for. The smoothing factor for the problem (1) is:

(25)

Here, are the coefficients given in (5)-(7).

For a = 10^{−}^{4}, ct = 1, d_{1} = 1 × 10^{−}^{12}, d_{2} = 1 × 10^{−}^{10}, d_{3} = 3 × 10^{−}^{10}, v_{1} = 1 × 10^{−}^{9}, v_{2} = 1 × 10^{−}^{6}, v_{3} = 1 × 10^{−}^{6}, the smoothing factors for the Gauss-Seidel relaxation method are presented in Table 1.

The data from Table 1 show that:

Table 1. The smoothig factor as a function of and.

Ÿ the Gauss-Seidel red-black relaxation method is a very good smoother for this problem as the smoothing factors in the cases presented here are;

Ÿ both the discretization methods lead to good smoothing factors. The finite element method seems slightly more appropriate when the number of grids used in the multigrid method is bigger;

Ÿ the number of relaxation steps before and after the coarse grid correction should not be too big as the smoothing factor increases with.

4.2. Asimptotic Convergence Factor and Error Reduction Factor

For the multigrid method having the components described in (9-16), the matrix of the two-grid operator for the problem (1) is:

(26)

For the corresponding asimptotic convergence factor and error reduction factor have been computed from the matrix (26) and are given in Table 2 for different numbers of preand postsmoothing steps.

The data from Table 2 show that the multigrid method is very rapidly convergent: if at least one smoothing step is performed before and after the coarse grid correction, then the error is reduced by at least a factor per multigrid cycle.

5. Numerical Results

The problem (1) has been solved on a domain containing tree layers with different diffusion and convection coefficients ([16,17]).

The error was computed for the exact solution

The time step in the discretization process has been. Figure 1, Figure 2 and the Table 3 represent the error after eight multigrid cycles, with two smoothing steps before and two after the coarse grid correction.

Table 2. Asimptotic convergence factor and error reduction factor.

(a)(b)

Figure 1. The multigrid error at ad = 100 nm in the skin for v_{1} = 1.0 × 10^{−}^{10}, v_{2} = 1.0 × 10^{−}^{7}, v_{3} = 1.0 × 10^{−}^{7}, d_{1} = 1 × 10^{−}^{12}, d_{2} = 1 × 10^{−}^{10}, d_{3} = 3 × 10^{−}^{10}; c = 10^{4}; a = 0.

(a)(b)

Figure 2. The multigrid error at ad = 100 nm in the skin for v_{1} = 1.0 × 10^{−}^{9}, v_{2} = 1.0 × 10^{−}^{6}, v_{3} = 1.0 × 10^{−}^{6}, d_{1} = 1 × 10^{−}^{12}, d_{2} = 1 × 10^{−}^{10}, d_{3} = 3 × 10^{−}^{10}; c = 10^{3}; a = 0.

Table 3. Multigrid error for FD and FEM.

Table 3 shows the maximum absolute value of the error and the norm of the error vector corresponding to Figures 1(a) and (b), for the finite differences discretization method (FD) and finite element method (FEM).

6. Conclusion

We have presented a convergence and error analysis for the multigrid method applied to a time dependent diffusion-convection problem that is convection dominated. The mathematical model is applied to the study of the concentration of a solute that is transported by the blood, or the penetration of a substance through the skin layers.

The convergence analysis showed that the discretization process is better realized by the finite element method than the finite differences. Also the red-black GaussSeidel is a good smoother for the problem presented here, and needs not to be applied more than two or three times in the construction of the multigrid method.The numerical results in the previous paragraph confirmed the good convergence and error reduction as predicted by the coefficients computed with the local Fourier analysis.

REFERENCES

- K. W. Morton and D. F. Mayers, “Numerical Solution of Partial Differential Equations, An Introduction,” Cambridge University Press, Cambridge, 2005. doi:10.1017/CBO9780511812248
- E. Becker, G. Carey and J. Oden, “Finite Elements. An Introduction,” Prentice-Hall, Englewood Cliffs, 1981.
- H. C. Elman, D. J. Silvester and A. J. Wathen, “Finite Elements and Fast Iterative Solvers,” Oxford University Press, Oxford, 2005.
- R. P. Fedorenko, “A Relaxation Method for Solving Elliptic Difference Equations,” USSR Computational Mathematics and Mathematical Physics, Vol. 1, 1962, pp. 1092- 1096.
- R. P. Fedorenko, “The Speed of Convergence of One Iterative Process,” USSR Computational Mathematics and Mathematical Physics, Vol. 4, 1964, pp. 227-235.
- A. Brandt, “Multilevel Adaptive Solutions to Boundary Value Problems,” Mathematics of Computation, Vol. 31, 1977, pp. 333-390. doi:10.1090/S0025-5718-1977-0431719-X
- A. Brandt, “Multigrid Techniques: 1984 Guide with Applications to Fluid Dynamics,” GMD-Studien Nr. 85, Gesellschaft für Matematik und Datenverarbeitung, St. Augustin, Bonn, 1984.
- W. Hackbush, “Elliptic Differential Equations,” SpringerVerlag, New York, 1992. doi:10.1007/978-3-642-11490-8
- U. Trottenberg, C. Oosterlee and A. Schuller, “Multigrid,” Elsevier Academic Press, London, 2001.
- P. Weseling, “An Introduction to Multigrid Method,” John Wiley & Sons, New York, 1991.
- M. A. Olshanski and A. Reusken, “On a Robust Multigrid Method for Connection-Diffusion Finite Element Problems.” http://www.math.uh.edu/ molshan/ftp/pub/proceed_cd.pdf
- A. Reusken, “Convergence Analysis of a Multigrid Method for Convection-Diffusion Equations,” Numerische Mathematik, Vol. 91, No. 2, 2002, pp. 323-349. doi:10.1007/s002110100312
- R. Wienands and W. Joppich, “Practical Fourier Analysis for Multigrid Methods,” Chapman& Hall/CRC Press, Boca Raton, 2005.
- R. W. Lewis, P. Nithiarasu and K. N. Seetharamu, “Fundamentals of the Finite Element Method for Heat and Fluid Flow,” John Wiley & Sons Ltd, The Atrium, 2004. doi:10.1002/0470014164
- W. L. Briggs, V. E. Henson and S. McCormick, “A Multigrid Tutorial,”2nd Edition, Siam, Philadelphia, 2000. doi:10.1137/1.9780898719505
- D. Neumann, “Modeling Transdermal Absorption,” Biotechnology: Pharmaceutical Aspects, Springer, New York, 2008.
- B. Al-Qallaf, D. Bhusan Das, D. Mori and Z. Cui, “Modelling Transdermal Delivery of High Molecular Weight Drugs from Microneedle Systems,” Philosophical Transactions of the Royal Society A, Vol. 365, 2007, pp. 2951- 2967. doi:10.1098/rsta.2007.0003