The numerical simulation of convection delayed dominated diffusion equation
MATEC Web of Conferences 57, 05007 (2016)
DOI: 10.1051/ matecconf/20165705007
ICAET- 2016
The numerical simulation of convection delayed dominated diffusion
equation
P. Murali Mohan Kumar
1, a
, A.S.V. Ravi Kanth
1
Department of Mathematics, National Institute of Technology Kurukshetra, Haryana – 136119, India
1
Abstract. In this paper, we propose a fitted numerical method for solving convection delayed dominated diffusion
equation. A fitting factor is introduced and the model equation is discretized by cubic spline method. The error
analysis is analyzed for the consider problem. The numerical examples are solved using the present method and
compared the result with the exact solution.
1 Introduction
Consider the following convection delayed dominated
diffusion equation
y ''( x) a( x) y '( x) b( x) y( x ) 0 on [0,1] (1)
subject to the interval conditions
(2)
y( x) ( x) on x 0 , y(1)
where 0
1 is a perturbation parameter and is a
small shifting parameter of order . It is also assumed
is a
that a( x), b( x), ( x) are smooth functions and
constant. These convection diffusion delayed types with
dominated convection term problems play an important
role in the mathematical modelling of various practical
phenomena in the engineering and environmental
sciences, for examples include high Reynold’s number
flow in fluid dynamics, heat transport problems with
large Pecklet number, modelling the problems in
mathematical biology and semi-conductor devices etc.
It is challenging to develop efficient numerical methods
for solving convection diffusion with dominated
convection term due to the existence of boundary layers.
Standard discretization methods for solving such kind of
problems are unstable and fails to give accurate results
when perturbation parameter is small. Therefore, it is
challenging to develop suitable numerical methods to
these problems, whose accuracy does not depend on the
parameter . Lange and Miura[1] initiated the singular
perturbation analysis of boundary value problems for
differential difference equations with small shifts. The
numerical study of second order singularly perturbed
delay differential equations has been given in [2-3] and
references therein.
In this paper, we present an exponentially fitted method
on uniform mesh based on cubic spline method for the
convection delayed dominated diffusion equation. The
layout of the paper is organized as follows: Continuation
of the problem is presented in next section and follows
a
description of the method. In Section 3, the error analysis
of the method is discussed. Section 4 ends with the
Numerical results.
2 Continuous problem
An application of Taylor series in (1) yields
y ''( x) (a( x) b( x)) y '( x) b( x) y( x) 0,
xi 1
x
xi 1
y( x) 0 , y(1)
(3)
a( x) b( x) M 0 throughout the
interval [0,1] , where M is positive constant, then the
problem (3) exhibits boundary layer at x 0 .
From the theory of singular perturbations [4],
We assume that
y x y0 x y0 0 e
a 0 b 0
x
o
(4)
where y0 x is the solution of the reduced problem given
by
a x b x y0 x b x y0 x 0 with y0 1
From Eq.(4) as h 0 , we obtain
lim y ih y0 0 0 y0 0 e
a 0 b 0
ih
h 0
Let
h
, then
lim y ih y0 0 0 y0 0 e
h0
a 0 b 0 i
(5)
Now introducing an exponentially fitting factor to
the Eq.(3), we get
y x a x b x y x b x y x 0
(6)
with
y(0)=0 , y(1)
(7)
Corresponding author:
© The Authors, published by EDP Sciences. This is an open access article distributed under the terms of the Creative Commons Attribution
License 4.0 (http://creativecommons.org/licenses/by/4.0/).
MATEC Web of Conferences 57, 05007 (2016)
DOI: 10.1051/ matecconf/20165705007
ICAET- 2016
The fitting factor is to be determined in such a way
n Eq.(10) and using the following three
approximations for first order derivatives
y yi 1
yi i 1
2h
3 yi 1 4 yi yi 1
yi1
2h
yi 1 4 yi 3 yi 1
yi1
2h
We get the following difference scheme
Ei yi 1 Fi yi Gi yi 1 0,i 1 1 N 1
that the solution of Eq.(6) converges uniformly to the
solution Eq.(3) .
Lemma 1. Let u(x) be a smooth function satisfing the
and
then
.
Prof We can prove the above lemma by method of
contradiction.
Let
be
such
that
and assume that
.
and
and
Now
Clearly
consider
(10)
Where
"
3"
Ei 1 ai 1 bi 1 "2 ai bi 2 ai 1 bi 1 "1hbi 1
2
2
2
Fi
2"1 ai 1 bi 1 2" 2 ai 1 bi 1 2" 2 hbi
3"
"
Gi 1 ai 1 bi 1 " 2 ai bi 2 ai 1 bi 1 "1hbi 1
2
2
Eq.(11) gives a system of N 1 equations with
N 1 unknowns. These N 1 equations together with the
Eq.(7) are sufficient to solve the system by using Thomas
algorithm.
Which is contradiction to our assumption. Hence
.
Lemma 2. Let u(x) is the solution of the boundary value
problem (3), then
Proof. Let
point
be two barrier functions defined by
Then, we have
and
as
.
Using this inquality in the above inequality, we get
2.2 Determination of fitting factor
Taking the limit as h 0 in Eq.(11), we obtain
2
yi 1
lim yi
a 0 b 0 "1 " 2 hlim
0
h 0
!
Therefore by the maximum principle[5], we obtain
which gives the required
estimate.
(11)
yi 1 0
a 0 b 0 "1 " 2 hlim
! 0
Substituting Eq.(5) in Eq.(12) and then simplifying, we
get the variable fitting factor as follows
2.1 Description of the method
The spline function
Let
S ( x, ) S ( x) satisfying in the interval [ xi , xi 1 ] and the
differential equation
x x
S '' ( x) S ( x) S '' ( xi ) S ( xi ) i 1
h !
i
a xi b xi
"1 " 2 a xi b xi coth
2
2
!
( 12)
x xi
S ( xi 1 ) S ( xi 1 )
h !
''
(8)
Where S ( xi ) yi and 0 is termed as cubic spline in
compression. Following Aziz and Khan [6], we obtain the
tridiagonal system
h2 "1M i 1 2"2 M i "1M i 1 yi 1 2 yi yi 1
(9)
Where
"
1
1
"1 2 1
, " 2 2 " cot " 1 ,
sin " !
"
"
is a constant fitting factor for left end boundary layer.
3 Error analysis
Substituting M j a j b j y 'j b j y j , j i, i # 1 in
(10), we obtain
M i y xi , i 1 1 N 1.
Substituting
M j a x j b x j
a 0 b 0
"1 " 2 a 0 b 0 coth
2
2
!
y x j b x j y x j , i
j i,i # 1,
2
MATEC Web of Conferences 57, 05007 (2016)
DOI: 10.1051/ matecconf/20165705007
ICAET- 2016
yi 1 2 yi yi 1
h 2 (2)
ai ('1 )
2!
h2
ai 1 ai hai' ai(2) ('2 )
2!
2
h
bi 1 bi hbi' bi(2) ('3 )
2!
2
h
bi 1 bi hbi' bi(2) (' 4 )
2!
Where xi 1 '1 ,'2 ,'3 ,'4 xi 1
Using these expansions and Eq.(21), we have
( "1h 2bi 1 )ei 1 (2 2"2 h 2bi )ei
ai (...truncated)
