# Stability of Exact and Discrete Energy for Non-Fickian Reaction-Diffusion Equations with a Variable Delay

Authors: Dongfang Li; Chao Tong; Jinming Wen
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

This paper is concerned with the stability of non-Fickian reaction-diffusion equations with a variable delay. It is shown that the perturbation of the energy function of the continuous problems decays exponentially, which provides a more accurate and convenient way to express the rate of decay of energy. Then, we prove that the proposed numerical methods are sufficient to preserve energy stability of the continuous problems. We end the paper with some numerical experiments on a biological model to confirm the theoretical results.


---

## Body

## 1. Introduction

Reaction-diffusion equations with delay are widely proposed as models for biological, physical, and control systems [1–5]. Over the past several years, such equations have been extensively investigated and several important properties such as existence and stability of the travelling wavefront are well understood (e.g., [6–9]). It is known that traveling waves describe transition processes from the physical point of view. If the reaction is very fast, the speed of propagation will become rather large (e.g., [10–12]). This pathologic behavior is not presented in the physical phenomena but introduced by some mathematical models (cf. [13]). In order to avoid the unphysical behavior, certain memory effects are taken into account in the mathematical models. In [13–15], the time memory flux is proposed by assuming that a flux observed at some time should have something to do with the gradient of the solution at some past time; that is,
(1)
J
(
x
,
t
+
δ
)
=
-
D
2
∂
u
∂
x
(
x
,
t
)
.
The first order approximation of the above formula gives
(2)
J
(
x
,
t
)
+
δ
∂
J
∂
t
(
x
,
t
)
=
-
D
2
∂
u
∂
x
(
x
,
t
)
.
Its solution satisfies
(3)
J
(
x
,
t
)
=
-
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
∂
x
(
x
,
s
)
d
s
.
In fact, the memory term is presented to avoid the infinite propagation speed in flux definition. These ideas yield the study on the following non-Fickian reaction-diffusion equation with a constant delay:
(4)
∂
u
∂
t
(
x
,
t
)
=
D
1
∂
2
u
∂
x
2
(
x
,
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
∂
x
2
(
x
,
s
)
d
s
+
f
(
t
,
u
(
x
,
t
)
,
u
(
x
,
t
-
τ
)
)
,
where D
2
>
0 and δ
>
0.Recently, many researchers have investigated these kinds of equations. The existence and uniqueness of the solutions were derived by Chang [16]. The well posedness of the integrodifferential models was studied by Araújo et al. in [17]. Later, they further studied the effect of the integrodifferential term in the qualitative behavior of the solution of the diffusion equations in [18]. Zhang and Vandewalle [19] studied the general linear methods for solving the nonlinear integrodifferential equations with memory. Wang et al. investigated the asymptotic stability of exact and discrete solutions of nonlinear integrodifferential equation in [20–22]. Wen et al. considered dissipativity of the functional differential equations [23]. For a more detailed description of this subject as well as its open problems, we refer the readers to the books [1–3, 24, 25], the papers [26–32], and the references therein.In particular, some continuous and discrete energy estimates, which are determined by theL
2-norm of the solution and the estimates for the L
2-norm of the past in time of its gradient, were derived by Branco et al. in [13]. However, the derived energy estimates rely on time. When time goes to infinity, the energy estimates seem to become boundless. Later, inspired by their instructive work, Li et al. further derived some asymptotic stability of problem (
1.1
) in [33]. The stability result implies that the perturbation of the energy estimate tends to zero as t approaches infinity. However, it does not imply anything about how long it takes to decay. Moreover, when referring to the impact of various factors (e.g., environment, temperature, and other potential effects) in some real-world problems, the model should be modified to an equation with a variable delay, which makes the numerical simulation and its analysis more complicated. Hence, some further stability results for the non-Fickian equations with a variable delay are much more significant and challenging from both the physical and mathematical points of view.One aim of the paper is to study the exponential stability of some more general models. Consider the non-Fickian reaction-diffusion equation with a variable delay(5)
∂
u
∂
t
(
x
,
t
)
=
D
1
∂
2
u
∂
x
2
(
x
,
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
∂
x
2
(
x
,
s
)
d
s
+
f
(
t
,
u
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
,
where (
x
,
t
)
∈
(
x
a
,
x
b
)
×
(
0
,
+
∞
), 0
≤
τ
(
t
)
≤
τ
*, the boundary conditions
(6)
u
(
x
a
,
t
)
=
u
(
x
b
,
t
)
=
0
,
t
>
0
,
and initial condition
(7)
u
(
x
,
t
)
=
u
0
(
x
,
t
)
,
x
∈
(
x
a
,
x
b
)
,
t
∈
[
-
τ
*
,
0
]
.
We find that, the perturbation of the energy estimate decays exponentially. This exponential stability result provides a more accurate and convenient way to express the rate of decay. Meanwhile, the stability result can also be applied to the equations investigated in [13–15]. The given results imply that the perturbations of the energy of the equations are controlled by the initial perturbations from the systems.The other aim is to investigate the discrete energy stability of (5). Since finite element method can be easily designed for high order of accuracy in space (see, e.g., [34, 35]), we introduce the method to solve the equation. Exponential stability of the semidiscrete system is derived. Then, the implicit Euler method is further applied to discretize the equation. And the integral is approximated by the right rectangular rule and the delay term is approximated by a linear interpolation. It is proved that the proposed numerical methods are sufficient to preserve stability of the underlying systems. Besides, if we discretize the diffusion term by using the centered difference operator, the energy stability of full discrete non-Fickian equation with a variable delay can also be derived. The given result indicates that the perturbations of the numerical solutions vanish eventually. Finally, a numerical simulation on a biological model is presented to confirm the theoretical and numerical results.The rest of the paper is organized as follows. In Section2, we discuss energy stability of non-Fickian delay reaction-diffusion equations. Section 3 describes in detail the stability of the numerical methods for the equation. Section 4 shows experimental studies for verifying the proposed results. Finally, conclusions and discussions for this paper are summarized in Section 5.


## 2. Exponential Stability

The following lemma can be derived from Lemma 2.1 in [36], where it was used to investigate the dissipativity of delay differential equations. Here it will play a key role in discussing the energy stability of the problem (5).Lemma 1.
Suppose that(8)
Y
′
(
t
)
≤
-
α
(
t
)
Y
(
t
)
+
β
(
t
)
sup
⁡
t
-
τ
≤
ξ
≤
t
Y
(
ξ
)
,
t
≥
t
0
.
Here τ
≥
0 and α
(
t
), β
(
t
) are continuous functions such that α
(
t
)
≥
α
0
>
0, and 0
≤
β
(
t
)
≤
q
α
(
t
) for all t
≥
t
0 with constants α
0
>
0 and 0
≤
q
<
1. Then
(9)
Y
(
t
)
≤
G
exp
⁡
(
-
μ
*
(
t
-
t
0
)
)
,
t
≥
t
0
,
where G
=
sup
⁡
t
0
-
τ
≤
t
≤
t
0
|
Y
(
t
)
| and μ
*
>
0 is defined as
(10)
μ
*
=
inf
⁡
t
≥
t
0
{
μ
(
t
)
:
μ
(
t
)
-
α
(
t
)
+
β
(
t
)
e
μ
(
t
)
τ
=
0
}
.For the stability analysis, we need to consider the perturbed problem, that is, non-Fickian delay reaction-diffusion equation (5) with a different initial condition u
~
(
x
,
t
)
=
u
~
0
(
x
,
t
). Its solution is denoted by u
~
(
x
,
t
) and satisfies the following equation:
(11)
∂
u
~
∂
t
(
x
,
t
)
=
D
1
∂
2
u
~
∂
x
2
(
x
,
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
~
∂
x
2
(
x
,
s
)
d
s
+
f
(
t
,
u
~
(
x
,
t
)
,
u
~
(
x
,
t
-
τ
(
t
)
)
)
,
t
>
0
,
u
~
(
x
,
t
)
=
u
~
0
(
x
,
t
)
,
(
x
,
t
)
∈
(
x
a
,
x
b
)
×
[
-
τ
*
,
0
]
,
u
~
(
x
a
,
t
)
=
u
~
(
x
b
,
t
)
=
0
,
t
>
0
.Besides, let〈
·
,
·
〉 denote the inner product in L
2
(
x
a
,
x
b
) and by ∥
·
∥ we denote the corresponding induced norm. To simplify the notation, we represent v
(
·
,
t
) by v
(
t
) if v is defined in [
x
a
,
x
b
]
×
[
-
τ
,
+
∞
). The energy function is defined by
(12)
E
v
(
t
)
=
∥
v
(
t
)
∥
2
+
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
.
Then, we have the following result.Theorem 2.
Assume that for eacht
∈
[
0
,
+
∞
)
(13)
∂
u
∂
x
(
t
)
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
∂
x
(
s
)
d
s
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
0
∂
x
(
s
)
d
s
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
~
0
∂
x
(
s
)
d
s
∈
L
2
(
x
a
,
x
b
)
.
If
(14)
〈
u
1
-
u
2
,
f
(
t
,
u
1
,
v
1
)
-
f
(
t
,
u
2
,
v
1
)
〉
≤
α
1
(
t
)
∥
u
1
-
u
2
∥
2
,
∥
f
(
t
,
u
,
v
1
)
-
f
(
t
,
u
,
v
2
)
∥
≤
β
1
(
t
)
∥
v
1
-
v
2
∥
,
where α
1
(
t
) and β
1
(
t
)
≥
0 are continuous and there exists 0
≤
q
1
<
1 such that
(15)
q
1
max
⁡
{
-
2
δ
,
2
α
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
+
β
1
(
t
)
}
+
β
1
(
t
)
≤
0
,
f
o
r
t
≥
0
,
then
(16)
E
v
(
t
)
≤
G
1
exp
⁡
(
-
μ
1
*
t
)
,
where G
1
=
sup
⁡
-
τ
≤
t
≤
0
|
E
v
(
t
)
| and μ
1
*
>
0 is defined as
(17)
μ
1
*
=
inf
⁡
t
≥
t
0
{
{
-
2
δ
,
2
α
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
+
β
1
(
t
)
}
μ
(
t
)
:
μ
(
t
)
=
inf
⁡
t
≥
t
0
+
max
⁡
{
-
2
δ
,
2
α
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
+
β
1
(
t
)
}
=
inf
⁡
t
≥
t
0
+
β
1
(
t
)
e
μ
(
t
)
τ
=
0
{
-
2
δ
,
2
α
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
+
β
1
(
t
)
}
}
.Proof.
Letv
(
x
,
t
)
=
u
(
x
,
t
)
-
u
~
(
x
,
t
); we can obtain the following equation:
(18)
∂
v
∂
t
(
t
)
=
D
1
∂
2
v
∂
x
2
(
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
v
∂
x
2
(
s
)
d
s
+
f
(
t
,
u
(
t
)
,
u
(
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
t
)
,
u
~
(
t
-
τ
(
t
)
)
)
.
Multiplying each member of (18) by v
(
x
,
t
) with respect to 〈
·
,
·
〉 and integrating by part, we find
(19)
1
2
d
d
t
∥
v
(
t
)
∥
2
=
-
D
1
∥
∂
v
∂
x
(
t
)
∥
2
-
D
2
δ
〈
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
,
∂
v
∂
x
(
t
)
〉
+
〈
f
(
t
,
u
(
t
)
,
u
(
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
t
)
,
u
~
(
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
.
As in [15],
(20)
〈
∂
v
∂
x
(
t
)
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
〉
=
1
2
d
d
t
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
+
1
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
,
we have
(21)
d
d
t
(
∥
v
(
t
)
∥
2
+
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
)
=
-
2
D
1
∥
∂
v
∂
x
(
t
)
∥
2
-
2
δ
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
+
2
〈
f
(
t
,
u
(
t
)
,
u
(
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
t
)
,
u
~
(
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
≤
-
2
D
1
(
x
b
-
x
a
)
2
∥
v
(
t
)
∥
2
-
2
δ
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
+
2
〈
f
(
t
,
u
(
t
)
,
u
(
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
t
)
,
u
~
(
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
,
where we have used the Poincaré inequality.
Now, using assumption (14), we have
(22)
〈
f
(
t
,
u
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
x
,
t
)
,
u
~
(
x
,
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
=
〈
f
(
t
,
u
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
+
〈
f
(
t
,
u
~
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
-
f
(
t
,
u
~
(
x
,
t
)
,
u
~
(
x
,
t
-
τ
(
t
)
)
)
,
v
(
t
)
〉
≤
α
1
(
t
)
∥
v
(
t
)
∥
2
+
β
1
(
t
)
∥
v
(
t
-
τ
(
t
)
)
∥
·
∥
v
(
t
)
∥
≤
(
α
1
(
t
)
+
β
1
(
t
)
2
)
∥
v
(
t
)
∥
2
+
β
1
(
t
)
2
∥
v
(
t
-
τ
(
t
)
)
∥
2
.
Therefore,
(23)
d
d
t
(
∥
v
(
t
)
∥
2
+
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
)
≤
-
2
δ
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
+
(
2
α
1
(
t
)
+
β
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
)
∥
v
(
t
)
∥
2
+
β
1
(
t
)
∥
v
(
t
-
τ
(
t
)
)
∥
2
≤
max
⁡
⁡
{
-
2
δ
,
2
α
1
(
t
)
+
β
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
}
×
(
∥
v
(
t
)
∥
2
+
D
2
δ
∥
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
v
∂
x
(
s
)
d
s
∥
2
)
+
β
1
(
t
)
∥
v
(
t
-
τ
(
t
)
)
∥
2
.
Then, we derive the following inequality:(24)
d
d
t
E
v
(
t
)
≤
max
⁡
{
-
2
δ
,
2
α
1
(
t
)
+
β
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
}
×
E
v
(
t
)
+
β
1
(
t
)
E
v
(
t
-
τ
(
t
)
)
≤
max
⁡
{
-
2
δ
,
2
α
1
(
t
)
+
β
1
(
t
)
-
2
D
1
(
x
b
-
x
a
)
2
}
×
E
v
(
t
)
+
β
1
(
t
)
sup
⁡
t
-
τ
*
≤
ξ
≤
t
E
v
(
ξ
)
.
Now, applying Lemma1, we finally have the stability result.It is remarkable that formula (16) implies that
(25)
∥
v
(
t
)
∥
2
≤
G
1
exp
⁡
(
-
μ
1
*
t
)
.
Hence, the perturbation of the solution of the problem also decays exponentially.When the right-hand side function of the problem (5) does not possess the delay term, the problem degenerates into a non-Fickian reaction-diffusion equation:
(26)
∂
u
∂
t
(
x
,
t
)
=
D
1
∂
2
u
∂
x
2
(
x
,
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
∂
x
2
(
x
,
s
)
d
s
+
f
(
t
,
u
(
x
,
t
)
)
,
t
>
0
,
which is also discussed in [13]. We have the following conclusion.Corollary 3.
Assume that for eacht
∈
[
0
,
+
∞
)
(27)
∂
u
∂
x
(
t
)
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
∂
x
(
s
)
d
s
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
0
∂
x
(
s
)
d
s
∈
L
2
(
x
a
,
x
b
)
.
If
(28)
〈
f
(
t
,
u
)
-
f
(
t
,
v
)
,
u
-
v
〉
≤
α
2
(
t
)
∥
u
-
v
∥
2
,
where α
2
(
t
) are continuous and 2
α
2
(
t
)
-
2
D
1
/
(
x
b
-
x
a
)
2
<
0, for t
≥
0. Then there exist constants G
2 and μ
2
* such that
(29)
E
v
(
t
)
≤
G
2
exp
⁡
(
-
μ
2
*
t
)
.WhenD
1
=
0, the problem (5) degenerates into an integrodifferential delay equation,
(30)
∂
u
∂
t
(
x
,
t
)
=
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
∂
x
2
(
x
,
s
)
d
s
+
f
(
t
,
u
(
x
,
t
)
,
u
(
x
,
t
-
τ
(
t
)
)
)
,
t
>
0
,
which is investigated in [14, 15]. We have the following conclusion.Corollary 4.
Assume that for eacht
∈
[
0
,
+
∞
)
(31)
∂
u
∂
x
(
t
)
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
∂
x
(
s
)
d
s
,
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
0
∂
x
(
s
)
d
s
∈
L
2
(
x
a
,
x
b
)
.
If
(32)
〈
u
1
-
u
2
,
f
(
t
,
u
1
,
v
1
)
-
f
(
t
,
u
2
,
v
1
)
〉
≤
α
3
(
t
)
∥
u
1
-
u
2
∥
2
,
∥
f
(
t
,
u
,
v
1
)
-
f
(
t
,
u
,
v
2
)
∥
≤
β
3
(
t
)
∥
v
1
-
v
2
∥
,
where α
3
(
t
) and β
3
(
t
)
≥
0 are continuous and there exists 0
≤
q
3
<
1 such that
(33)
q
1
max
⁡
{
-
2
δ
,
2
α
1
(
t
)
+
β
1
(
t
)
}
+
β
1
(
t
)
<
0
,
for
t
≥
0
,
then there exist constants G
3 and μ
3
* such that
(34)
E
v
(
t
)
≤
G
3
exp
⁡
(
-
μ
3
*
t
)
.If we assume that(35)
〈
f
(
t
,
u
,
v
)
,
u
〉
≤
α
1
(
t
)
∥
u
∥
2
+
β
1
(
t
)
∥
v
∥
2
+
γ
1
(
t
)
,
where β
1
(
t
)
≥
0 and γ
1
(
t
)
≥
0 and α
1
(
t
) are continuous functions, there may exist a bounded set ℬ and a time t
*, such that for any given initial function u
0
(
x
,
t
)
,
(
x
,
t
)
∈
[
x
a
,
x
b
]
×
[
-
τ
,
0
], the corresponding energy ∥
E
u
(
t
)
∥ of the problem is contained in ℬ for all t
≥
t
*. The analogous results and their conclusions for non-Fickian equations with a variable delay can be derived without any difficulty. We do not list them here for brevity.


## 3. Stability of the Numerical Approximation

In this section, we will concentrate on the stability of two kinds of numerical approximation. Here and below, for the discretization of system (5), we divide the interval I
=
[
x
a
,
x
b
] with a mesh: x
a
=
x
0
<
x
1
<
⋯
<
x
N
=
x
b with the space stepsize h
=
(
x
b
-
x
a
)
/
N, where N is a positive integer. And k
=
1
/
N
k is a time stepsize, where N
k is a positive integer. We will numerically solve the problem at time n
k, n
=
1,2
,
….


### 3.1. Stability of the Finite Element Approximation

Firstly, the finite element method is introduced to discrete the problems. LetS
h be a finite dimensional subspace of H
0
1
(
I
); then the semidiscrete finite element method is to find U
∈
S
h, such that, for all test functions, χ
i
∈
S
h and 1
≤
i
≤
N, satisfying
(36)
(
∂
u
h
∂
t
,
χ
)
+
(
D
1
∂
u
h
∂
x
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
h
∂
x
(
s
)
d
s
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
(
x
,
t
)
,
u
h
(
x
,
t
-
τ
(
t
)
)
)
,
χ
)
.The exponential stability of semidiscrete approximation is shown as follows.Theorem 5.
Assume that conditions in Theorem2 hold; then the difference E
v
h
(
t
) satisfies the following formulae:
(37)
∥
E
v
h
(
t
)
∥
2
≤
G
1
exp
⁡
(
-
μ
*
(
t
-
t
0
)
)
,
t
≥
t
0
,
where the parameters G
1 and μ
* are defined in Theorem 2 and v
h
=
u
h
-
u
~
h with u
~
h satisfying the semidiscrete finite element approximation to the perturbed problem.Proof.
The proof is similar to that of Theorem2.If we further apply the right rectangular rule to approximate the integral and implicit Euler method to discrete (36), then the full discrete numerical approximation is to find u
h
n
∈
S
h, such that, for all test functions, χ
∈
S
h and 1
≤
i
≤
N, satisfying for n
=
0,1
,
2
,
…
(38)
(
u
h
n
+
1
-
u
h
n
k
,
χ
)
+
(
∑
i
=
1
n
+
1
‍
D
1
∂
∂
x
u
h
n
+
1
+
D
2
k
δ
(
u
h
n
+
1
-
u
h
n
k
,
χ
)
+
×
∑
i
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
i
)
/
δ
∂
∂
x
u
h
i
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
,
χ
)
.
Here u
¯
h
n
+
1 denotes an approximation to u
(
x
,
t
n
+
1
-
τ
(
t
n
+
1
)
), which is obtained by the following linear interpolation procedure at the point t
n
+
1
-
τ
(
t
n
+
1
):
(39)
u
¯
h
n
+
1
=
θ
n
u
h
n
+
1
-
m
n
+
(
1
-
θ
n
)
u
h
n
-
m
n
,
where τ
(
t
)
=
(
m
n
-
θ
n
)
k with integer m
n
≥
1 and θ
n
∈
[
0,1
].Now, we can establish the discrete version of the energy:(40)
E
~
(
u
h
n
)
=
∥
u
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
j
=
1
n
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
u
h
l
∥
2
.Then, we have the following result.Theorem 6.
Assume that conditions (14) and (15) hold. Then there exists a constant 0
≤
q
2
<
1 such that
(41)
E
~
(
v
h
n
+
1
)
≤
q
2
E
~
(
v
h
n
)
,
where v
h
n
=
u
h
n
-
u
~
h
n and u
~
h
n
(
x
i
) is the full discrete approximations to the perturbed problem.Proof.
The differencev
h
n satisfies the following equation:
(42)
(
v
h
n
+
1
-
v
h
n
k
,
χ
)
+
(
∑
i
=
1
n
+
1
‍
D
1
∂
∂
x
v
h
n
+
1
+
D
2
k
δ
(
v
h
n
+
1
-
v
h
n
k
,
χ
)
+
×
∑
i
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
i
)
/
δ
∂
∂
x
v
h
i
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
χ
)
.
Taking χ
=
v
h
n
+
1 in the above formula, we find that
(43)
∥
v
h
n
+
1
∥
2
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
v
h
n
+
1
)
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
¯
h
n
+
1
)
,
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
v
h
n
+
1
)
≤
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
(
α
n
+
1
+
β
n
+
1
2
)
k
∥
v
h
n
+
1
∥
2
+
β
n
+
1
2
k
∥
v
¯
h
n
+
1
∥
2
,
where we have used condition (14).
As in [15], we have
(44)
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
=
1
2
∥
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
-
1
2
e
-
2
(
k
/
δ
)
∥
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
1
2
∥
∂
∂
x
v
h
n
+
1
∥
2
.
Substituting (44) into (43) yields
(45)
∥
v
h
n
+
1
∥
2
+
D
2
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
+
D
2
2
δ
e
-
2
(
k
/
δ
)
×
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
-
D
2
k
2
2
δ
∥
∂
∂
x
v
h
n
+
1
∥
2
+
(
α
n
+
1
+
β
n
+
1
2
)
k
∥
v
h
n
+
1
∥
2
+
β
n
+
1
2
k
∥
v
¯
h
n
+
1
∥
2
.
Applying the following discrete Poincaré inequality(46)
-
∥
∂
∂
x
v
h
n
+
1
∥
2
≤
-
1
(
x
b
-
x
a
)
2
∥
v
h
n
+
1
∥
2
and (
v
h
n
,
v
h
n
+
1
)
≤
(
1
/
2
)
∥
v
h
n
+
1
∥
2
+
(
1
/
2
)
∥
v
h
n
∥
2 to (45), we arrive at
(47)
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
(
x
b
-
x
a
)
2
)
∥
v
h
n
+
1
∥
2
+
D
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
≤
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
∥
v
¯
h
n
+
1
∥
2
=
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
∥
θ
n
v
n
n
+
1
-
m
n
+
(
1
-
θ
n
)
v
h
n
-
m
n
∥
2
=
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
(
θ
n
2
∥
v
n
n
+
1
-
m
n
∥
2
+
(
1
-
θ
n
)
2
∥
v
h
n
-
m
n
∥
2
+
2
θ
n
(
1
-
θ
n
)
∥
v
n
n
+
1
-
m
n
∥
∥
v
h
n
-
m
n
∥
∥
v
n
n
+
1
-
m
n
∥
2
)
≤
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
(
θ
n
∥
v
n
n
+
1
-
m
n
∥
2
+
(
1
-
θ
n
)
∥
v
h
n
-
m
n
∥
2
)
≤
E
~
(
v
h
n
)
+
k
β
n
+
1
max
⁡
{
∥
v
h
n
+
1
-
m
n
∥
2
,
∥
v
h
n
-
m
n
∥
2
}
≤
E
~
(
v
h
n
)
+
k
β
n
+
1
max
⁡
⁡
{
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
≤
(
1
+
k
β
n
+
1
)
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
,
where we have used the inequality 2
α
β
≤
α
2
+
β
2.
It follows from (47) that
(48)
min
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
E
~
(
v
h
n
+
1
)
≤
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
(
x
b
-
x
a
)
2
)
×
∥
v
h
n
+
1
∥
2
+
D
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
≤
(
1
+
k
β
n
+
1
)
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
,
where C
n
+
1
=
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
/
(
x
b
-
x
a
)
2
).
Therefore,(49)
E
~
(
v
h
n
+
1
)
≤
(
1
+
k
β
n
+
1
)
min
⁡
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
×
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
.
In fact, in view of condition (15), there exists a constant ϱ, such that
(50)
1
C
n
+
1
+
k
β
n
+
1
C
n
+
1
≤
ϱ
<
1
,
(51)
e
-
2
(
k
/
δ
)
(
1
+
k
β
n
+
1
)
≤
e
-
k
β
n
+
1
(
1
+
k
β
n
+
1
)
=
(
1
-
k
β
n
+
1
+
𝒪
(
k
2
)
)
×
(
1
+
k
β
n
+
1
)
≤
ϱ
<
1
.
Hence,
(52)
1
+
k
β
n
+
1
min
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
≤
ϱ
<
1
,
n
=
1,2
,
…
Together with (49) and (52), we have the final conclusion.Theorem5 shows that, for any given initial perturbation, the perturbation of the corresponding discrete energy of the problem decreases in time. Meantime, noting that ∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
(
∂
/
∂
x
)
v
h
l
∥
2
>
0 in the energy function, the difference ∥
v
h
n
∥ tends to zero very quickly.


### 3.2. Stability of the Finite Difference Approximation

In this subsection, we will apply the central finite difference method to discrete the diffusion term and discuss the stability of the full discrete problems.Letk be the time stepsize and U
h
n
(
x
i
) denote the numerical approximation to U
(
x
i
,
t
n
) with t
n
=
n
k; then the full discrete finite difference approximation for (5) can be written as
(53)
U
h
n
+
1
(
x
i
)
=
U
h
n
(
x
i
)
+
k
D
1
D
2
,
x
U
h
n
+
1
(
x
i
)
+
k
2
D
2
δ
×
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
2
,
x
U
h
l
(
x
i
)
+
k
f
(
t
n
+
1
,
U
h
n
+
1
(
x
i
)
,
U
¯
h
n
+
1
(
x
i
)
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
n
=
0,1
,
2
,
…
,
where
(54)
D
2
,
x
U
h
(
x
i
)
=
U
h
(
x
i
+
1
)
-
2
U
h
(
x
i
)
+
U
h
(
x
i
-
1
)
h
2
,
iiiiiiiiiiiiiiiiiiiiiiiiiiii
i
=
1
,
…
,
N
-
1
,
and the argument U
¯
h
n
+
1
(
x
) denotes an approximation to u
(
x
,
t
n
+
1
-
τ
(
t
n
+
1
)
), which is obtained by the linear interpolation procedure at the point t
n
+
1
-
τ
(
t
n
+
1
)
(55)
U
¯
h
n
+
1
=
θ
n
U
h
n
+
1
-
m
n
+
(
1
-
θ
n
)
U
h
n
-
m
n
,
where τ
(
t
)
=
(
m
n
-
θ
n
)
k with integer m
n
≥
1 and θ
n
∈
[
0,1
].The discrete version of the energy of the finite difference approximation is defined by(56)
E
~
~
(
U
h
n
+
1
)
=
∥
U
h
n
+
1
∥
L
2
(
I
h
)
2
+
e
-
2
(
k
/
δ
)
D
2
δ
×
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
U
h
l
∥
L
2
(
I
h
,
+
)
2
,
where D
-
x is the backward finite difference operator:
(57)
∥
U
h
∥
L
2
(
I
h
)
2
=
(
U
h
,
U
h
)
h
=
(
U
h
,
U
h
)
h
=
h
∑
i
=
1
N
-
1
‍
U
h
2
(
x
i
)
,
∥
V
∥
L
2
(
I
h
,
+
)
2
=
(
V
,
V
)
h
,
+
=
h
∑
i
=
1
N
‍
V
(
x
i
)
V
(
x
i
)
,
V
=
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
U
h
l
.Theorem 7.
Assume that conditions (14) and (15) hold. Then there exists a constant 0
≤
q
3
<
1 such that
(58)
E
~
~
(
V
h
n
+
1
)
≤
q
3
E
~
~
(
V
h
n
)
,
where V
h
n
=
U
h
n
-
U
~
h
n and U
~
h
n
(
x
i
) is the full discrete finite difference approximations to the perturbed problem.Proof.
The proof is similar to the one done in Theorem6 with the following relation:
(59)
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
j
)
/
δ
D
-
x
v
h
j
,
D
-
x
v
h
n
+
1
)
h
,
+
=
1
2
∥
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
v
h
l
∥
L
2
(
I
h
,
+
)
2
-
1
2
e
-
2
(
k
/
δ
)
∥
∑
l
=
1
n
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
v
h
l
∥
L
2
(
I
h
,
+
)
2
+
1
2
∥
D
-
x
v
h
n
+
1
∥
L
2
(
I
h
,
+
)
2
.


## 3.1. Stability of the Finite Element Approximation

Firstly, the finite element method is introduced to discrete the problems. LetS
h be a finite dimensional subspace of H
0
1
(
I
); then the semidiscrete finite element method is to find U
∈
S
h, such that, for all test functions, χ
i
∈
S
h and 1
≤
i
≤
N, satisfying
(36)
(
∂
u
h
∂
t
,
χ
)
+
(
D
1
∂
u
h
∂
x
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
u
h
∂
x
(
s
)
d
s
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
(
x
,
t
)
,
u
h
(
x
,
t
-
τ
(
t
)
)
)
,
χ
)
.The exponential stability of semidiscrete approximation is shown as follows.Theorem 5.
Assume that conditions in Theorem2 hold; then the difference E
v
h
(
t
) satisfies the following formulae:
(37)
∥
E
v
h
(
t
)
∥
2
≤
G
1
exp
⁡
(
-
μ
*
(
t
-
t
0
)
)
,
t
≥
t
0
,
where the parameters G
1 and μ
* are defined in Theorem 2 and v
h
=
u
h
-
u
~
h with u
~
h satisfying the semidiscrete finite element approximation to the perturbed problem.Proof.
The proof is similar to that of Theorem2.If we further apply the right rectangular rule to approximate the integral and implicit Euler method to discrete (36), then the full discrete numerical approximation is to find u
h
n
∈
S
h, such that, for all test functions, χ
∈
S
h and 1
≤
i
≤
N, satisfying for n
=
0,1
,
2
,
…
(38)
(
u
h
n
+
1
-
u
h
n
k
,
χ
)
+
(
∑
i
=
1
n
+
1
‍
D
1
∂
∂
x
u
h
n
+
1
+
D
2
k
δ
(
u
h
n
+
1
-
u
h
n
k
,
χ
)
+
×
∑
i
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
i
)
/
δ
∂
∂
x
u
h
i
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
,
χ
)
.
Here u
¯
h
n
+
1 denotes an approximation to u
(
x
,
t
n
+
1
-
τ
(
t
n
+
1
)
), which is obtained by the following linear interpolation procedure at the point t
n
+
1
-
τ
(
t
n
+
1
):
(39)
u
¯
h
n
+
1
=
θ
n
u
h
n
+
1
-
m
n
+
(
1
-
θ
n
)
u
h
n
-
m
n
,
where τ
(
t
)
=
(
m
n
-
θ
n
)
k with integer m
n
≥
1 and θ
n
∈
[
0,1
].Now, we can establish the discrete version of the energy:(40)
E
~
(
u
h
n
)
=
∥
u
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
j
=
1
n
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
u
h
l
∥
2
.Then, we have the following result.Theorem 6.
Assume that conditions (14) and (15) hold. Then there exists a constant 0
≤
q
2
<
1 such that
(41)
E
~
(
v
h
n
+
1
)
≤
q
2
E
~
(
v
h
n
)
,
where v
h
n
=
u
h
n
-
u
~
h
n and u
~
h
n
(
x
i
) is the full discrete approximations to the perturbed problem.Proof.
The differencev
h
n satisfies the following equation:
(42)
(
v
h
n
+
1
-
v
h
n
k
,
χ
)
+
(
∑
i
=
1
n
+
1
‍
D
1
∂
∂
x
v
h
n
+
1
+
D
2
k
δ
(
v
h
n
+
1
-
v
h
n
k
,
χ
)
+
×
∑
i
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
i
)
/
δ
∂
∂
x
v
h
i
,
∂
χ
∂
x
)
=
(
f
(
t
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
χ
)
.
Taking χ
=
v
h
n
+
1 in the above formula, we find that
(43)
∥
v
h
n
+
1
∥
2
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
v
h
n
+
1
)
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
¯
h
n
+
1
)
,
v
h
n
+
1
)
+
k
(
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
¯
h
n
+
1
)
-
f
(
t
n
+
1
,
u
~
h
n
+
1
,
u
~
¯
h
n
+
1
)
,
v
h
n
+
1
)
≤
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
-
D
2
k
2
δ
×
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
+
(
α
n
+
1
+
β
n
+
1
2
)
k
∥
v
h
n
+
1
∥
2
+
β
n
+
1
2
k
∥
v
¯
h
n
+
1
∥
2
,
where we have used condition (14).
As in [15], we have
(44)
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
j
)
/
δ
∂
∂
x
v
h
j
,
∂
∂
x
v
h
n
+
1
)
=
1
2
∥
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
-
1
2
e
-
2
(
k
/
δ
)
∥
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
1
2
∥
∂
∂
x
v
h
n
+
1
∥
2
.
Substituting (44) into (43) yields
(45)
∥
v
h
n
+
1
∥
2
+
D
2
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
=
(
v
h
n
,
v
h
n
+
1
)
-
k
D
1
∥
∂
∂
x
v
h
n
+
1
∥
2
+
D
2
2
δ
e
-
2
(
k
/
δ
)
×
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
-
D
2
k
2
2
δ
∥
∂
∂
x
v
h
n
+
1
∥
2
+
(
α
n
+
1
+
β
n
+
1
2
)
k
∥
v
h
n
+
1
∥
2
+
β
n
+
1
2
k
∥
v
¯
h
n
+
1
∥
2
.
Applying the following discrete Poincaré inequality(46)
-
∥
∂
∂
x
v
h
n
+
1
∥
2
≤
-
1
(
x
b
-
x
a
)
2
∥
v
h
n
+
1
∥
2
and (
v
h
n
,
v
h
n
+
1
)
≤
(
1
/
2
)
∥
v
h
n
+
1
∥
2
+
(
1
/
2
)
∥
v
h
n
∥
2 to (45), we arrive at
(47)
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
(
x
b
-
x
a
)
2
)
∥
v
h
n
+
1
∥
2
+
D
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
≤
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
∥
v
¯
h
n
+
1
∥
2
=
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
∥
θ
n
v
n
n
+
1
-
m
n
+
(
1
-
θ
n
)
v
h
n
-
m
n
∥
2
=
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
(
θ
n
2
∥
v
n
n
+
1
-
m
n
∥
2
+
(
1
-
θ
n
)
2
∥
v
h
n
-
m
n
∥
2
+
2
θ
n
(
1
-
θ
n
)
∥
v
n
n
+
1
-
m
n
∥
∥
v
h
n
-
m
n
∥
∥
v
n
n
+
1
-
m
n
∥
2
)
≤
∥
v
h
n
∥
2
+
e
-
2
(
k
/
δ
)
D
2
δ
∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
+
k
β
n
+
1
(
θ
n
∥
v
n
n
+
1
-
m
n
∥
2
+
(
1
-
θ
n
)
∥
v
h
n
-
m
n
∥
2
)
≤
E
~
(
v
h
n
)
+
k
β
n
+
1
max
⁡
{
∥
v
h
n
+
1
-
m
n
∥
2
,
∥
v
h
n
-
m
n
∥
2
}
≤
E
~
(
v
h
n
)
+
k
β
n
+
1
max
⁡
⁡
{
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
≤
(
1
+
k
β
n
+
1
)
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
,
where we have used the inequality 2
α
β
≤
α
2
+
β
2.
It follows from (47) that
(48)
min
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
E
~
(
v
h
n
+
1
)
≤
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
(
x
b
-
x
a
)
2
)
×
∥
v
h
n
+
1
∥
2
+
D
2
δ
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
∂
∂
x
v
h
l
∥
2
≤
(
1
+
k
β
n
+
1
)
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
,
where C
n
+
1
=
(
1
-
2
k
α
n
+
1
-
k
β
n
+
1
+
2
k
D
1
/
(
x
b
-
x
a
)
2
).
Therefore,(49)
E
~
(
v
h
n
+
1
)
≤
(
1
+
k
β
n
+
1
)
min
⁡
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
×
max
⁡
{
E
~
(
v
h
n
)
,
E
~
(
v
h
n
+
1
-
m
n
)
,
E
~
(
v
h
n
-
m
n
)
}
.
In fact, in view of condition (15), there exists a constant ϱ, such that
(50)
1
C
n
+
1
+
k
β
n
+
1
C
n
+
1
≤
ϱ
<
1
,
(51)
e
-
2
(
k
/
δ
)
(
1
+
k
β
n
+
1
)
≤
e
-
k
β
n
+
1
(
1
+
k
β
n
+
1
)
=
(
1
-
k
β
n
+
1
+
𝒪
(
k
2
)
)
×
(
1
+
k
β
n
+
1
)
≤
ϱ
<
1
.
Hence,
(52)
1
+
k
β
n
+
1
min
⁡
{
C
n
+
1
,
e
2
(
k
/
δ
)
}
≤
ϱ
<
1
,
n
=
1,2
,
…
Together with (49) and (52), we have the final conclusion.Theorem5 shows that, for any given initial perturbation, the perturbation of the corresponding discrete energy of the problem decreases in time. Meantime, noting that ∥
k
∑
l
=
1
n
‍
e
-
(
t
n
-
t
l
)
/
δ
(
∂
/
∂
x
)
v
h
l
∥
2
>
0 in the energy function, the difference ∥
v
h
n
∥ tends to zero very quickly.


## 3.2. Stability of the Finite Difference Approximation

In this subsection, we will apply the central finite difference method to discrete the diffusion term and discuss the stability of the full discrete problems.Letk be the time stepsize and U
h
n
(
x
i
) denote the numerical approximation to U
(
x
i
,
t
n
) with t
n
=
n
k; then the full discrete finite difference approximation for (5) can be written as
(53)
U
h
n
+
1
(
x
i
)
=
U
h
n
(
x
i
)
+
k
D
1
D
2
,
x
U
h
n
+
1
(
x
i
)
+
k
2
D
2
δ
×
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
2
,
x
U
h
l
(
x
i
)
+
k
f
(
t
n
+
1
,
U
h
n
+
1
(
x
i
)
,
U
¯
h
n
+
1
(
x
i
)
)
,
iiiiiiiiiiiiiiiiiiiiiiiiiiiiii
n
=
0,1
,
2
,
…
,
where
(54)
D
2
,
x
U
h
(
x
i
)
=
U
h
(
x
i
+
1
)
-
2
U
h
(
x
i
)
+
U
h
(
x
i
-
1
)
h
2
,
iiiiiiiiiiiiiiiiiiiiiiiiiiii
i
=
1
,
…
,
N
-
1
,
and the argument U
¯
h
n
+
1
(
x
) denotes an approximation to u
(
x
,
t
n
+
1
-
τ
(
t
n
+
1
)
), which is obtained by the linear interpolation procedure at the point t
n
+
1
-
τ
(
t
n
+
1
)
(55)
U
¯
h
n
+
1
=
θ
n
U
h
n
+
1
-
m
n
+
(
1
-
θ
n
)
U
h
n
-
m
n
,
where τ
(
t
)
=
(
m
n
-
θ
n
)
k with integer m
n
≥
1 and θ
n
∈
[
0,1
].The discrete version of the energy of the finite difference approximation is defined by(56)
E
~
~
(
U
h
n
+
1
)
=
∥
U
h
n
+
1
∥
L
2
(
I
h
)
2
+
e
-
2
(
k
/
δ
)
D
2
δ
×
∥
k
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
U
h
l
∥
L
2
(
I
h
,
+
)
2
,
where D
-
x is the backward finite difference operator:
(57)
∥
U
h
∥
L
2
(
I
h
)
2
=
(
U
h
,
U
h
)
h
=
(
U
h
,
U
h
)
h
=
h
∑
i
=
1
N
-
1
‍
U
h
2
(
x
i
)
,
∥
V
∥
L
2
(
I
h
,
+
)
2
=
(
V
,
V
)
h
,
+
=
h
∑
i
=
1
N
‍
V
(
x
i
)
V
(
x
i
)
,
V
=
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
U
h
l
.Theorem 7.
Assume that conditions (14) and (15) hold. Then there exists a constant 0
≤
q
3
<
1 such that
(58)
E
~
~
(
V
h
n
+
1
)
≤
q
3
E
~
~
(
V
h
n
)
,
where V
h
n
=
U
h
n
-
U
~
h
n and U
~
h
n
(
x
i
) is the full discrete finite difference approximations to the perturbed problem.Proof.
The proof is similar to the one done in Theorem6 with the following relation:
(59)
(
∑
j
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
j
)
/
δ
D
-
x
v
h
j
,
D
-
x
v
h
n
+
1
)
h
,
+
=
1
2
∥
∑
l
=
1
n
+
1
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
v
h
l
∥
L
2
(
I
h
,
+
)
2
-
1
2
e
-
2
(
k
/
δ
)
∥
∑
l
=
1
n
‍
e
-
(
t
n
+
1
-
t
l
)
/
δ
D
-
x
v
h
l
∥
L
2
(
I
h
,
+
)
2
+
1
2
∥
D
-
x
v
h
n
+
1
∥
L
2
(
I
h
,
+
)
2
.


## 4. Application

In this section, we will provide a numerical experiment to illustrate the given results.Consider the following Non-Fickian Mackey-Glass equations [33]:
(60)
∂
u
∂
t
(
x
,
t
)
=
D
1
∂
2
u
∂
x
2
(
x
,
t
)
+
D
2
δ
∫
0
t
‍
e
-
(
t
-
s
)
/
δ
∂
2
u
∂
x
2
(
x
,
s
)
d
s
-
a
u
(
x
,
t
)
+
b
u
(
x
,
t
-
τ
(
t
)
)
1
+
u
m
(
x
,
t
-
τ
(
t
)
)
,
u
(
x
a
,
t
)
=
u
(
x
b
,
t
)
=
0
,
t
>
0
,
u
(
x
,
t
)
=
φ
(
x
,
t
)
,
(
x
,
t
)
∈
[
0
,
π
]
×
[
-
τ
*
,
0
]
,
where a, b, D
1, D
2, and δ are given parameters, m is an even positive integer, and τ
(
t
) is a nonnegative function.There is no difficulty in verifying that condition (14) holds (α
1
(
t
)
=
-
a
,
β
1
(
t
)
=
b). If there exists 0
≤
q
1
<
1 such that q
1
max
⁡
{
-
2
/
δ
,
-
2
a
-
2
D
1
/
(
x
b
-
x
a
)
2
+
b
}
+
b
<
0, the perturbation of the energy estimate decays exponentially. For the numerical simulation of the model, we set the parameter D
1
=
D
2
=
δ
=
1, m
=
2, τ
(
t
)
=
1
-
0.1
|
sin
(
t
)
|, the space stepsize h
=
π
/
10, and time stepsize k
=
0.1 and apply 1-degree finite element method to solve the problem on [
0
,
π
]
×
[
0,20
] with the following different initial conditions:
(61)
u
0
(
x
,
t
)
=
(
2
+
cos
⁡
(
t
)
)
sin
(
x
)
,
u
~
0
(
x
,
t
)
=
10
x
(
π
-
x
)
.
Some statistics of the numerical results with different parameters a and b are shown in Tables 1 and 2, where the discrete energy E
~
(
v
h
n
) and ∥
v
h
n
∥
2 of the problems decreases in time very quickly. Moreover, The numerical differences |
v
h
n
(
x
i
)
| with the parameters a
=
2 and b
=
1 are shown in Figure 1, where x
i
=
i
h with i
=
1,2
,
…
,
9. The numerical differences (the blue lines) are all bounded by the exponential function y
=
20
exp
⁡
(
-
t
) (the red line). Clearly, they all confirm the results in this paper.Table 1
The discrete energyE
~
(
v
h
n
) at different time.
E
~
(
v
h
n
)
t
=
2
t
=
4
t
=
6
t
=
8
t
=
10
a
=
1, b
=
1
1.37
1.01
E
-
1
9.17
E
-
4
6.90
E
-
5
6.53
E
-
6
a
=
2, b
=
1
5.28
E
-
1
2.21
E
-
2
1.72
E
-
4
6.32
E
-
6
2.08
E
-
7
a
=
3, b
=
1
2.88
E
-
1
7.73
E
-
3
8.20
E
-
5
1.55
E
-
6
3.33
E
-
8
a
=
4, b
=
1
1.87
E
-
1
3.94
E
-
3
5.32
E
-
5
8.98
E
-
7
1.69
E
-
8Table 2
The difference∥
v
h
n
∥
2 at different time.
∥
v
h
n
∥
2
t
=
2
t
=
4
t
=
6
t
=
8
t
=
10
a
=
1, b
=
1
9.52
E
-
1
8.43
E
-
2
4.37
E
-
4
5.80
E
-
5
5.77
E
-
6
a
=
2, b
=
1
2.82
E
-
1
4.56
E
-
2
2.86
E
-
5
3.53
E
-
6
1.38
E
-
7
a
=
3, b
=
1
1.19
E
-
1
4.00
E
-
3
8.21
E
-
6
1.92
E
-
7
7.38
E
-
9
a
=
4, b
=
1
6.16
E
-
2
1.39
E
-
3
4.56
E
-
6
3.40
E
-
9
5.04
E
-
10Figure 1
The difference|
v
h
n
(
x
i
)
| at different time.


## 5. Conclusions

In this paper, we have investigated the stability of continuous and discrete non-Fickian delay reaction-diffusion equations with a variable delay. Firstly, we show that the perturbation of the energy estimate decays exponentially. The present stability results for a more general case improved that of our previous paper. Secondly, the numerical analysis indicates that the finite element method or central finite difference method, combined with implicit Euler method and a linear interpolation procedure for the delay term, can preserve the stability of the energy function of the equations. All the above findings are confirmed by using the numerical methods on the non-Fickian Mackey-Glass equations.



---

*Source: 840573-2014-03-05.xml*