Optimal error estimates of a linearized second-order BDF scheme for a nonlocal parabolic problem

Daoussa Haggar, M.S.; Mbehou, M.

doi:10.1108/AJMS-05-2022-0126

Purpose

This paper focuses on the unconditionally optimal error estimates of a linearized second-order scheme for a nonlocal nonlinear parabolic problem. The first step of the scheme is based on Crank–Nicholson method while the second step is the second-order BDF method.

Design/methodology/approach

A rigorous error analysis is done, and optimal L² error estimates are derived using the error splitting technique. Some numerical simulations are presented to confirm the study’s theoretical analysis.

Findings

Optimal L² error estimates and energy norm.

Originality/value

The goal of this research article is to present and establish the unconditionally optimal error estimates of a linearized second-order BDF finite element scheme for the reaction-diffusion problem. An optimal error estimate for the proposed methods is derived by using the temporal-spatial error splitting techniques, which split the error between the exact solution and the numerical solution into two parts, that is, the temporal error and the spatial error. Since the spatial error is not dependent on the time step, the boundedness of the numerical solution in L∞-norm follows an inverse inequality immediately without any restriction on the grid mesh.

1. Introduction

In this paper, we consider the following parabolic problem with nonlocal nonlinearity:

\{\begin{cases} u_{t} - a (l (u)) Δ u + α | u |^{p - 2} u = f (u) i n Ω \times (0, T], \\ u (x, t) = 0 on \partial Ω \times (0, T], \\ u (x, 0) = u_{0} (x) in Ω, \end{cases}

(1.1)

where $Ω \subset R^{d}$ ⁠, d ≥ 1 is again a domain with a smooth boundary ∂Ω, a and f are functions to be defined in the next section and l denote a continuous linear form on L²(Ω) given by

l (u (t)) = \int_{Ω} g (x) u (t, x) d x,

where g is a function on L²(Ω).

The study of nonlocal parabolic problems has received considerable attention in recent years ([1–3] and the references therein). This kind of problems arises in various situations, for instance, u could describe the density of a population (for instance, bacteria) subject to spreading. The diffusion coefficient a is then supposed to depend on the entire population in the domain rather than on the local density, that is, moves are guided by considering the global state of the medium. The problem is nonlocal in the sense that the diffusion coefficient is determined by a global quantity. Besides its mathematical motivation because of the presence of the nonlocal term a(l(u)), such problems come from physical situations related to migration of a population of bacteria in a container in which the velocity of migration v = a∇u depends on the global population in a subdomain Ω′ ⊂ Ω given by a(l(u)).

Simsen and Ferreira [4] have discussed not only the existence and uniqueness of solutions for this problem but also continuity with respect to initial values, the exponential stability of weak solutions and important results on the existence of a global attractor. The numerical methods for the nonlocal problems have been investigated by many authors as like in Refs [5, 6] and the references therein. However, they are restricted to nonlocal reaction terms or nonlocal boundary conditions. Chaudhary et al. [7] studied the convergence analysis of the Crank–Nicolson finite element method for the nonlocal problem involving the Dirichlet energy. Mbehou et al. [8] studied (1.1) using the Crank–Nicolson Galerkin finite element method. The main focus on this paper was to present the exponential decay and vanishing of the solutions in finite time. They also derived the optimal convergence order in L²-norm using P_r with r ≥ 1 finite elements. Yin and Xu [9] applied the finite-volume method to obtain approximate solutions for a nonlocal problem on reactive flows in porous media and derived the optimal convergence order in the L² norm. Almeida et al. [10] presented convergence analysis for a fully discretized approximation to a nonlocal problem involving a parabolic equation with moving boundaries, with the finite element method applied for the space variables and the Crank–Nicolson method for the time. Recently, Yang et al. [11] derived the unconditional optimal error estimate of Galerkin FEMs for the time-dependent Klein–Gordon–Schrodinger equations using the error splitting technique. Also in Ref. [12], Yang and Jiang applied the linearized second-order backward differentiation formulae (BDF) Galerkin Finite element methods (FEMs) for the Landau-Lifshitz equations to derive the unconditional optimal error estimates.

Our goal in this research article is to give and establish the unconditionally optimal error estimates of a linearized second-order BDF finite element scheme for the reaction-diffusion problem (1.1). Using P_r (r ≥ 1) finite element to approximate the solution of (1.1), the optimal error estimates O(Δt² + h^r+1) in L² norm are derived using the error splitting technique.

This paper is organized as follows. In Section 2, we recall few known results and present few regularities, which are used in the proof of the optimal error estimates. To prove the optimal error estimates by the error splitting technique, the temporal errors and the spatial errors are shown in Sections 3 and 4, respectively. Finally numerical results are presented in Section 5 to demonstrate our theoretical analysis.

2. Preliminaries and main results

Let $Ω \subset R^{d}$ (d ≥ 1) be a bounded domain with a smooth boundary ∂Ω = Γ. The standard notations (see for instance Refs [13, 14]) will be used throughout this work. The Lebesgue space is denoted L^p(Ω), 1 ≤ p ≤ ∞, with norms $‖ \cdot ‖_{L^{p}}$ but the L²(Ω)-norm will be denoted by ‖ ⋅‖. For any nonnegative integer m and real number p ≥ 1, the classical Sobolev spaces:

W^{m, p} (Ω) = {v \in L^{p} (Ω); D^{α} v \in L^{p} (Ω) for all | α | \leq m},

equipped with the semi-norm

| v |_{m, p} = {\{\sum_{| α | = m} \int_{Ω} | D^{α} v |^{p} d x\}}^{1 / p},

and the norm

‖ v ‖_{m, p} = {\{\sum_{0 \leq | α | \leq m} \int_{Ω} | D^{α} v |^{p} d x\}}^{1 / p},

(2.1)

with the usual extension when p = ∞. When p = 2, W^m,p(Ω) is the Hilbert space H^m(Ω) with the scalar product:

{((v, w))}_{m} = \sum_{| α | \leq m} (D^{α} v, D^{α} w)

The norm of H^m(Ω) will be denoted by ‖ ⋅‖_m. It should be mentioned that D^α stands for the derivative in the sense of distribution, while α = (α₁, …, α_d) denotes a multi-index of length |α| = α₁ + ⋯ + α_d. We also employ the standard notation of Bochner spaces, such as L^q(0, T, X) with norm

\begin{array}{l} ‖ w ‖_{L^{q} (X)} = {(\int_{0}^{T} ‖ w (t) ‖_{X}^{q} d t)}^{1 / q}, 1 \leq q < \infty, \\ ‖ w ‖_{L^{\infty} (X)} = e s s \sup_{0 \leq t \leq T} ‖ w (t) ‖_{X}, \end{array}

where X is an Hilbert space and ‖ ⋅‖_X the norm of X. For all these notions on Sobolev spaces and Bochner spaces, we refer to Refs [13, 15].

Throughout this paper, the following known inequalities will be frequently used [13].

‖ v ‖_{L^{r}} \leq C ‖ v ‖_{1} (2 \leq r \leq 6) \forall v \in H^{1} (Ω)

(2.2)

‖ v ‖_{\infty} \leq C ‖ v ‖^{1 / 2} ‖ v ‖_{2}^{1 / 2} \forall v \in H^{2} (Ω) .

(2.3)

Let us now suppose that α is a nonnegative constant and p > 1. Simsen and Ferreira [4] proved the existence and uniqueness of global solution under the following hypotheses.

H1.

u₀ ∈ L²(Ω).

H2.

$f : R \to R$ is Lipschitz–continuous function, that is, there exists γ > 0 such that |f(s) − f(t)| ≤ γ|s − t|, for all $s, t \in R$ and f(0) = 0.

H3.

$a : R \to R$ is bounded with 0 < m ≤ a(s) ≤ M, for all $s \in R$ with $λ_{1} > \frac{γ}{m}$ ⁠, where λ₁ is the first eigenvalue of $(- Δ, H_{0}^{1} (Ω))$ ⁠.

H4.

$a : R \to R$ is Lipschitz–continuous with $| a (s_{1}) - a (s_{2}) | \leq A | s_{1} - s_{2} |, \forall s_{1}, s_{2} \in R$ ⁠.

H5.

$l : L^{2} (Ω) \to R$ is a continuous linear form, i.e. there exists g ∈ L²(Ω) such that l(u) = l_g(u) = ∫_Ω g(x)u(x)dx, for all u ∈ L²(Ω).

Theorem 2.1

(Existence and uniqueness of solution, [4]). Assume that p ≥ 2 and if the hypotheses (H1)–(H5) hold, then problem (1.1) possesses a unique solution, that is, there exists a unique function u such that

u \in L^{2} (0, T, H_{0}^{1} (Ω) \cap L^{p} (Ω)) \cap C ([0, T]; L^{2} (Ω)),

(2.4)

u_{t} \in L^{2} (0, T, H^{- 1} (Ω)),

(2.5)

u (x, 0) = u_{0} (x) i n Ω .

(2.6)

(u_{t}, v) + a (l (u)) (\nabla u, \nabla v) + α (| u |^{p - 2} u, v) = (f (u), v), \forall v \in H_{0}^{1} (Ω),

(2.7)

wheremust be understood as an equality in

D^{'} (0, T)

⁠.

Given the hypotheses (H1)–(H5), we will also adopt another hypothesis, namely

H6.

for all r ≥ 1,

‖ u_{0} ‖_{H^{r + 1}} + ‖ u ‖_{L^{\infty} (H^{r + 1} (Ω))} + ‖ u_{t} ‖_{L^{2} (H^{r + 1} (Ω))} + ‖ u_{t t} ‖_{L^{2} (H^{1} (Ω))} + ‖ u_{ttt} ‖_{L^{2} (L^{2} (Ω))} \leq C .

(2.8)

The following lemmas will be useful.

Lemma 2.1

(cf. Ref. [16]). For all p ∈ (1, ∞) and τ ≥ 0, there exists a generic constant C = C(p, d) such that for all $ξ, η \in R^{d}$ with d ≥ 1 we have

‖ ξ |^{p - 2} ξ - | η |^{p - 2} η | \leq C | ξ - η |^{1 - τ} {(| ξ | + | η |)}^{p - 2 + τ} .

(2.9)

(| ξ |^{p - 2} ξ - | η |^{p - 2} η) \cdot (ξ - η) \geq C | ξ - η |^{2 + τ} {(| ξ | + | η |)}^{p - 2 - τ} .

(2.10)

Lemma 2.2

(cf. Ref. [3]). Let a and b be two nonnegaitve numbers. Then for all s ∈ (1, ∞),

| a^{s} - b^{s} | \leq | a - b | {(a + b)}^{s - 1} .

(2.11)

Lemma 2.3

(cf. Ref. [17]). Let a_k, b_k, c_k and γ_k, for integers k ≥ 0, be the positive numbers such that

a_{n} + τ \sum_{k = 0}^{n} b_{k} \leq τ \sum_{k = 0}^{n} γ_{k} a_{k} + τ \sum_{k = 0}^{n} c_{k} + B, for n \geq 0 .

(2.12)

Suppose thatτγ_k < 1, for allk, and set

σ_{k} = \frac{1}{(1 - τ γ_{k})}

⁠. Then

a_{n} + τ \sum_{k = 0}^{n} b_{k} \leq \exp (τ \sum_{k = 0}^{n} γ_{k} σ_{k}) (τ \sum_{k = 0}^{n} c_{k} + B), for n \geq 0 .

(2.13)

Remark. If the first sum on the right hand side of (2.12) extends only up to n − 1, then estimate (2.13) holds for all k > 0 with σ_k = 1.

Lemma 2.4

(H^k-estimate of elliptic equations [18]). Suppose that v is a solution of the boundary value problem

\begin{array}{l} - Δ v = f, in Ω, \\ v = 0, on \partial Ω, \end{array}

where

Ω \subset R^{d}

,d = 2, 3, is a smooth and bounded domain. Then,

‖ v ‖_{H^{k}} \leq C ‖ f ‖_{H^{k - 2}}, k = 2, 3 .

(2.14)

Let $T_{h} = {K}$ be a uniform triangular or tetrahedral partition of Ω into triangles or tetrahedrons. Thus, let $h = \max_{K \in T_{h}} {h_{K}}$ denote the mesh size, where h_K = diam(K) = max{‖x − y‖, x, y ∈ K}, and V_h be the finite dimensional subspace of $H_{0}^{1} (Ω)$ ⁠, which consists of continuous piecewise polynomials of degree r ≥ 1 on $T_{h}$ ⁠.

Let {t_n| t_n = nΔt; 0 ≤ n ≤ N} be a uniform partition of [0, T] with time step Δt = T/N. We write uⁿ = u(x, t_n), Uⁿ ≈ u(x, t_n) and for any sequence of functions ${w^{n}}_{n = 0}^{N}$ define

\begin{array}{l} D_{1} w^{n} & = \frac{w^{n} - w^{n - 1}}{Δ t}, \\ {\bar{w}}^{n} & = \frac{1}{2} (w^{n} + w^{n - 1}), n = 1,2 \dots, N, \\ D_{2} w^{n} & = \frac{3}{2} D_{1} w^{n} - \frac{1}{2} D_{1} w^{n - 1} and {\hat{w}}^{n} = 2 w^{n - 1} - w^{n - 2}, n = 2, \dots, N . \end{array}

The following telescope formula is for n ≥ 2

(D_{2} w^{n}, w^{n}) = \frac{1}{4 Δ t} (‖ w^{n} ‖^{2} - ‖ w^{n - 1} ‖^{2} + ‖ {\hat{w}}^{n + 1} ‖^{2} - ‖ {\hat{w}}^{n} ‖^{2} + ‖ w^{n} - 2 w^{n - 1} + w^{n - 2} ‖^{2}) .

(2.15)

Under the above notations, we propose the following linearized second-order BDF Galerkin finite element scheme associated to (1.1), which is to find $U_{h}^{n} \in V_{h}$ such that

Step 1: For $U_{h}^{0} = Π_{h} u_{0} \in V_{h}$ ⁠, find $U_{h}^{1} \in V_{h}$ such that for all v_h ∈ V_h

(D_{1} U_{h}^{1}, v_{h}) + a (l ({\hat{U}}_{h}^{1})) (\nabla {\bar{U}}_{h}^{1}, \nabla v_{h}) + α (| {\hat{U}}_{h}^{1} |^{p - 2} {\bar{U}}_{h}^{1}, v_{h}) = (f ({\hat{U}}_{h}^{1}), v_{h}),

(2.16)

where ${\hat{U}}_{h}^{1}$ is given by

(\frac{{\hat{U}}_{h}^{1} - U_{h}^{0}}{Δ t / 2}, v_{h}) + a (l (U_{h}^{0}) (\nabla {\hat{U}}_{h}^{1}, \nabla v_{h}) + α (| U_{h}^{0} |^{p - 2} {\hat{U}}_{h}^{1}, v_{h}) = (f (U_{h}^{0}), v_{h}) .

(2.17)

Step 2: For 2 ≤ n ≤ N, find $U_{h}^{n} \in V_{h}$ such that for all v_h ∈ V_h

(D_{2} U_{h}^{n}, v_{h}) + a (l ({\hat{U}}_{h}^{n})) (\nabla U_{h}^{n}, \nabla v_{h}) + α (| {\hat{U}}_{h}^{n} |^{p - 2} U_{h}^{n}, v_{h}) = (f ({\hat{U}}_{h}^{n}), v_{h}) .

(2.18)

Π_h is an interpolation operator from $H_{0}^{1} (Ω)$ to V_h.

Theorem 2.2

Assume that the hypotheses (H1)– (H5) hold. Then the fully discrete system defined in (2.16)–(2.18) has a unique solution $U_{h}^{n}$ which satisfies

‖ U_{h}^{n} ‖^{2} + ‖ {\hat{U}}_{h}^{n + 1} ‖^{2} + C Δ t \sum_{k = 1}^{n} ‖ \nabla U_{h}^{k} ‖^{2} \leq C ‖ u_{0} ‖^{2} .

(2.19)

Proof. 1 For the existence, taking $v_{h} = U_{h}^{1}$ ⁠, $v_{h} = {\hat{U}}_{h}^{1}$ and $v_{h} = U_{h}^{n}$ in (2.16)–(2.18), respectively, the existence and uniqueness of $U_{h}^{1}$ ⁠, ${\hat{U}}_{h}^{1}$ and $U_{h}^{n}$ are from the Lax–Milgram theorem and the hypothesis (H3).

Let $v_{h} = {\bar{U}}_{h}^{n}$ in (2.16), we have

\frac{1}{2 Δ t} (‖ U_{h}^{1} ‖^{2} - ‖ U_{h}^{0} ‖^{2}) + a (l ({\hat{U}}_{h}^{1})) ‖ \nabla {\bar{U}}_{h}^{1} ‖^{2} + α (| {\hat{U}}_{h}^{1} |^{p - 2} {\bar{U}}_{h}^{1}, {\bar{U}}_{h}^{1}) = (f ({\hat{U}}_{h}^{1}), {\bar{U}}_{h}^{1}) \leq ‖ f ({\hat{U}}_{h}^{1}) ‖ ‖ {\bar{U}}_{h}^{1} ‖ .

Drop the third term of the left hand side, use the lower bound of a(⋅) and (H2),

\frac{1}{2 Δ t} (‖ U_{h}^{1} ‖^{2} - ‖ U_{h}^{0} ‖^{2}) + M ‖ \nabla {\bar{U}}_{h}^{1} ‖^{2} \leq L ‖ {\hat{U}}_{h}^{1} ‖ ‖ {\bar{U}}_{h}^{1} ‖ \leq C ‖ {\hat{U}}_{h}^{1} ‖^{2} + \frac{M}{2} ‖ \nabla {\bar{U}}_{h}^{1} ‖^{2} .

and

‖ U_{h}^{1} ‖^{2} + C Δ t ‖ \nabla {\bar{U}}_{h}^{1} ‖^{2} \leq C Δ t ‖ {\hat{U}}_{h}^{1} ‖^{2} + ‖ u_{0} ‖^{2} .

(2.20)

Now, let $v_{h} = {\hat{U}}_{h}^{1}$ in (2.17), the same arguments used above give us

‖ {\hat{U}}_{h}^{1} ‖^{2} + C Δ t ‖ \nabla {\hat{U}}_{h}^{1} ‖^{2} \leq C ‖ u_{0} ‖^{2} .

Taking $v_{h} = U_{h}^{n}$ in (2.18), using the lower bound of a(⋅), (H2) and dropping the third term of the left hand side lead to

(D_{2} U_{h}^{n}, U_{h}^{n}) + M ‖ \nabla U_{h}^{n} ‖^{2} \leq L ‖ {\hat{U}}_{h}^{n} ‖ ‖ U_{h}^{n} ‖ .

From the telescope (2.15), we obtain

\begin{array}{l} \frac{1}{4 Δ t} (‖ U_{h}^{n} ‖^{2} - ‖ U_{h}^{n - 1} ‖^{2} + ‖ {\hat{U}}_{h}^{n + 1} ‖^{2} - ‖ {\hat{U}}_{h}^{n} ‖^{2}) + M ‖ \nabla U_{h}^{n} ‖^{2} \leq C L ‖ {\hat{U}}_{h}^{n} ‖ ‖ \nabla U_{h}^{n} ‖ \\ \leq C ‖ {\hat{U}}_{h}^{n} ‖^{2} + \frac{M}{2} ‖ \nabla U_{h}^{n} ‖^{2} . \end{array}

That is

‖ U_{h}^{n} ‖^{2} - ‖ U_{h}^{n - 1} ‖^{2} + ‖ {\hat{U}}_{h}^{n + 1} ‖^{2} - ‖ {\hat{U}}_{h}^{n} ‖^{2} + \frac{M}{2} Δ t ‖ \nabla U_{h}^{n} ‖^{2} \leq C ‖ {\hat{U}}_{h}^{n} ‖^{2} .

(2.21)

The relation (2.19) is obtained by summing up the above relation (2.21) and using the discrete Gronwall lemma 2.3.

The main result of this work is presented in the following theorem.

Theorem 2.3

Suppose that system (1.1) has a unique solution u satisfying (H6). Then the fully discrete system defined in (2.16)–(2.18) has a unique solution $U_{h}^{n}$ , and

\max_{0 \leq n \leq N} (‖ u^{n} - U_{h}^{n} ‖^{2} + ‖ {\hat{u}}^{n} - {\hat{U}}_{h}^{n} ‖^{2} + Δ t \sum_{k = 0}^{n} ‖ \nabla (u^{n} - U_{h}^{n}) ‖^{2}) \leq C (Δ t^{4} + h^{2 r + 2}),

(2.22)

whereCis a positive constant independent of Δtandh.

The proof of this theorem will be done in the following sections.

3. Error estimates for the semi-discrete problem

Let us introduce the corresponding time discrete system associated with (1.1)

Step 1: for U⁰ = u₀, find U¹ by

\{\begin{cases} D_{1} U^{1} - a (l ({\hat{U}}^{1})) Δ {\bar{U}}^{1} + α | {\hat{U}}^{1} |^{p - 2} {\bar{U}}^{1} = f ({\hat{U}}^{1}) \\ U^{1} = 0 on \partial Ω, \end{cases}

(3.1)

where ${\hat{U}}^{1}$ is the solution to

\{\begin{cases} \frac{{\hat{U}}^{1} - U^{0}}{Δ t / 2} - a (l (U^{0})) Δ {\hat{U}}^{1} + α | U^{0} |^{p - 2} {\hat{U}}^{1} = f (U^{0}) \\ {\hat{U}}^{1} = 0 on \partial Ω . \end{cases}

(3.2)

Step 2: for 2 ≤ n ≤ N, find Uⁿ by

\{\begin{cases} D_{2} U^{n} - a (l ({\hat{U}}^{n})) Δ U^{n} + α | {\hat{U}}^{n} |^{p - 2} U^{n} = f ({\hat{U}}^{n}) \\ U^{n} = 0 on \partial Ω . \end{cases}

(3.3)

The weak formulations of (3.1)–(3.3) are defined as follows: find $U^{n} \in H_{0}^{1} (Ω)$ such that for all $v \in H_{0}^{1} (Ω)$

(D_{1} U^{1}, v) + a (l ({\hat{U}}^{1})) (\nabla {\bar{U}}^{1}, \nabla v) + α (| {\hat{U}}^{1} |^{p - 2} {\bar{U}}^{1}, v) = (f ({\hat{U}}^{1}), v),

(3.4)

and for 2 ≤ n ≤ N

(D_{2} U^{n}, v) + a (l ({\hat{U}}^{n})) (\nabla U^{n}, \nabla v) + α (| {\hat{U}}^{n} |^{p - 2} U^{n}, v) = (f ({\hat{U}}^{n}), v),

(3.5)

with ${\hat{U}}^{1} \in H_{0}^{1} (Ω)$ such that

(\frac{{\hat{U}}^{1} - U^{0}}{Δ t / 2}, v) + a (l (U^{0})) (\nabla {\hat{U}}^{1}, \nabla v) + α (| U^{0} |^{p - 2} {\hat{U}}^{1}, v) = (f (U^{0}), v) .

(3.6)

The existence and uniqueness of the solution to problems (3.4)–(3.6) can be easily proved by using Lax–Milgram theorem.

Let u be the exact solution of (1.1). Then, u satisfies the following equations:

D_{1} u^{1} - a (l ({\hat{u}}^{1})) Δ {\bar{u}}^{1} + α | {\hat{u}}^{1} |^{p - 2} {\bar{u}}^{1} = f ({\hat{u}}^{1}) + R^{1}

(3.7)

D_{2} u^{n} - a (l ({\hat{u}}^{n})) Δ u^{n} + α | {\hat{u}}^{n} |^{p - 2} u^{n} = f ({\hat{u}}^{n}) + R^{n}, n = 2, \dots, N,

(3.8)

where ${\hat{u}}^{1}$ satisfies

\frac{{\hat{u}}^{1} - u_{0}}{Δ t / 2} - a (l (u_{0})) Δ {\hat{u}}^{1} + α | u_{0} |^{p - 2} {\hat{u}}^{1} = f (u_{0}) + {\hat{R}}^{0} .

(3.9)

${\hat{R}}^{0}$ ⁠, R¹ and Rⁿ are, respectively, the truncation errors given by

{\hat{R}}^{0} = (\frac{{\hat{u}}^{1} - u_{0}}{Δ t / 2} - u_{t}^{1 / 2}) - Δ {\hat{u}}^{1} (a (l (u_{0})) - a (l (u^{1 / 2}))) + α | u_{0} |^{p - 2} (u^{1 / 2} - u_{0}) + α (| u_{0} |^{p - 2} u_{0} - | u^{1 / 2} |^{p - 2} u^{1 / 2}) + (f (u_{0}) - f (u^{1 / 2})), R^{1} = (D_{1} 1 u^{1} - u_{t}^{1 / 2}) - Δ {\bar{u}}^{1} (a (l ({\hat{u}}^{1})) - a (l (u^{1 / 2}))) - a (l (u^{1 / 2})) Δ ({\bar{u}}^{1} - u^{1 / 2}) + α | {\hat{u}}^{1} |^{p - 2} ({\bar{u}}^{1} - {\hat{u}}^{1}) + α (| {\hat{u}}^{1} |^{p - 2} {\hat{u}}^{1} - | u^{1 / 2} |^{p - 2} u^{1 / 2}) + (f ({\hat{u}}^{1}) - f (u^{1 / 2})) R^{n} = (D_{2} u^{n} - u_{t}^{n}) - Δ u^{n} (a (l ({\hat{u}}^{n})) - a (l (u^{n}))) + α | {\hat{u}}^{n} |^{p - 2} (u^{n} - {\hat{u}}^{n}) + α (| {\hat{u}}^{n} |^{p - 2} {\hat{u}}^{n} - | u^{n} |^{p - 2} u^{n}) + (f ({\hat{u}}^{n}) - f (u^{n})) .

By Taylor formula and relation (2.9) with τ = 1, it is easy to see that

{(\sum_{n = 1}^{N} Δ t ‖ R^{n} ‖_{1}^{2})}^{1 / 2} + Δ t ‖ {\hat{R}}^{0} ‖_{1} \leq C Δ t^{2} .

(3.10)

Let us denote

{\bar{e}}^{1} = {\bar{u}}^{1} - {\bar{U}}^{1}; {\hat{e}}^{n} = {\hat{u}}^{n} - {\hat{U}}^{n}; e^{n} = u^{n} - U^{n} for 1 \leq n \leq N .

We have the following assumption.

Lemma 3.1

Assume that the exact solution u of (1.1) satisfies the regularities (2.8). Then there exists a positive constant C independent of Δt such that

‖ {\hat{e}}^{1} ‖ + Δ t^{1 / 2} ‖ \nabla {\hat{e}}^{1} ‖ \leq C Δ t^{2} .

(3.11)

Proof. Subtracting (3.6) from (3.9) leads to

2 {\hat{e}}^{1} - Δ t a (l (u_{0})) Δ {\hat{e}}^{1} + α Δ t | u_{0} |^{p - 2} {\hat{e}}^{1} = Δ t {\hat{R}}^{0} .

Testing the above equation by ${\hat{e}}^{1}$ yield

2 ‖ {\hat{e}}^{1} ‖^{2} + Δ t a (l (u_{0})) ‖ \nabla {\hat{e}}^{1} ‖^{2} + α Δ t \int_{Ω} | u_{0} |^{p - 2} | {\hat{e}}^{1} |^{2} d x = Δ t ({\hat{R}}^{0}, {\hat{e}}^{1})

Using the left bound of a(⋅) to the left hand side and Young’s inequality to the right hand side, we obtain

2 ‖ {\hat{e}}^{1} ‖^{2} + m Δ t ‖ \nabla {\hat{e}}^{1} ‖^{2} + α Δ t \int_{Ω} | u_{0} |^{p - 2} | {\hat{e}}^{1} |^{2} d x \leq \frac{1}{2} Δ t^{2} ‖ {\hat{R}}^{0} ‖^{2} + \frac{1}{2} ‖ {\hat{e}}^{1} ‖^{2} .

The proof ended by dropping the third term of the left hand side and applying (3.10) to the right hand side.

Based upon (3.11), we have

Proposition 3.1

Suppose that the solution u of (1.1) satisfies the regularities (2.8). Then there exists a generic constant C that does not dependent on Δt such that

‖ e^{1} ‖^{2} + Δ t ‖ \nabla e^{1} ‖^{2} \leq C Δ t^{5} .

(3.12)

Proof. Subtracting (2.16) from (3.7) and observing that e⁰ = 0 leads to

\begin{array}{l} e^{1} - Δ t a (l (u_{0})) Δ {\bar{e}}^{1} + α Δ t | {\hat{u}}^{1} |^{p - 2} {\bar{e}}^{1} = - α Δ t {\bar{u}}^{1} (| {\hat{u}}^{1} |^{p - 2} - | {\hat{U}}^{1} |^{p - 2}) \\ + Δ t Δ {\bar{u}}^{1} (a (l ({\hat{u}}^{1})) - a (l ({\hat{U}}^{1}))) + Δ t (f ({\hat{u}}^{1}) - f ({\hat{U}}^{1})) + Δ t R^{1} . \end{array}

Testing the above equation by e¹ and using the fact that ${\bar{e}}^{1} = \frac{1}{2} e^{1}$ ⁠, we have

\begin{align} ‖ e^{1} ‖^{2} + \frac{1}{2} m Δ t ‖ \nabla e^{1} ‖^{2} + \frac{1}{2} α Δ t \int_{Ω} | {\hat{u}}^{1} |^{p - 2} | e^{1} |^{2} d x \leq α Δ t ({\bar{u}}^{1} (| {\hat{u}}^{1} |^{p - 2} - | {\hat{U}}^{1} |^{p - 2}), e^{1}) \\ + Δ t (a (l ({\hat{u}}^{1})) - a (l ({\hat{U}}^{1}))) (\nabla {\bar{u}}^{1}, \nabla e^{1}) + Δ t (f ({\hat{u}}^{1}) - f ({\hat{U}}^{1}), e^{1}) + Δ t (R^{1}, e^{1}) \\ = \sum_{i = 1}^{4} I_{i} . \end{align}

(3.13)

We have

\begin{array}{l} I_{1} & = α Δ t ({\bar{u}}^{1} (| {\hat{u}}^{1} |^{p - 2} - | {\hat{U}}^{1} |^{p - 2}), e^{1}) \leq C_{1} (‖ {\bar{u}}^{1} ‖_{L^{\infty}}, ‖ {\hat{u}}^{1} ‖_{L^{\infty}}, ‖ {\hat{U}}^{1} ‖_{L^{\infty}}, p) Δ t ‖ {\hat{e}}^{1} ‖ ‖ e^{1} ‖ \\ \leq C_{1} Δ t^{2} ‖ {\hat{e}}^{1} ‖^{2} + \frac{1}{4} ‖ e^{1} ‖^{2} \leq C Δ t^{6} + \frac{1}{4} ‖ e^{1} ‖^{2} (using (3.11)) . \end{array}

\begin{array}{l} I_{2} & = Δ t (a (l ({\hat{u}}^{1})) - a (l ({\hat{U}}^{1}))) (\nabla {\bar{u}}^{1}, \nabla e^{1}) \leq Δ t A ‖ {\hat{e}}^{1} ‖ ‖ \nabla {\bar{u}}^{1} ‖ ‖ \nabla e^{1} ‖ (using (H 4)) \\ \leq C_{2} (A, ‖ \nabla {\bar{u}}^{1} ‖) Δ t ‖ {\hat{e}}^{1} ‖^{2} + \frac{m}{4} Δ t ‖ \nabla e^{1} ‖^{2} \leq C Δ t^{5} + \frac{m}{4} Δ t ‖ \nabla e^{1} ‖^{2} . \end{array}

\begin{array}{l} I_{3} & = Δ t (f ({\hat{u}}^{1}) - f ({\hat{U}}^{1}), e^{1}) \leq γ Δ t ‖ {\hat{e}}^{1} ‖ ‖ e^{1} ‖ (using (H 2)) \\ \leq C_{3} Δ t^{2} ‖ {\hat{e}}^{1} ‖^{2} + \frac{1}{4} ‖ e^{1} ‖^{2} \leq C Δ t^{6} + \frac{1}{4} ‖ e^{1} ‖^{2} . \end{array}

I_{4} = Δ t (R^{1}, e^{1}) \leq C_{4} ‖ R^{1} ‖^{2} + \frac{1}{4} ‖ e^{1} ‖^{2} \leq C Δ t^{5} + \frac{1}{4} ‖ e^{1} ‖^{2} (using (3.11)) .

Taking these estimates into (3.13), we obtain the desire result.

The main result in this section is as follows.

Theorem 3.1

Suppose that the solution u of (1.1) satisfies the regularities (2.8). Then there exists a generic constant C that does not dependent on Δt such that

\max_{1 \leq n \leq N} (‖ e^{n} ‖^{2} + ‖ {\hat{e}}^{n} ‖^{2} + Δ t \sum_{k = 1}^{n} ‖ \nabla e^{k} ‖^{2}) \leq C Δ t^{4},

(3.14)

\max_{0 \leq n \leq N} ‖ U^{n} ‖_{\infty} \leq C,

(3.15)

whereCis a positive constant independent ofnand Δt.

Proof. The proof of this theorem will be done using the mathematical induction.

In view of (3.11) and (3.12), the inequality (3.14) holds for n = 0, 1. Since U⁰ = u₀, the inquality (3.15) holds for n = 0. Now, let us assume that (3.14) and (3.15) hold for n ≤ m with m ≤ N − 1. Then we need to prove the inequality for n = m + 1. By the definition of ${\hat{U}}^{n}$ and the induction assumption, $‖ {\hat{U}}^{n} ‖_{\infty} \leq C$ ⁠.

Subtracting (2.18) from (3.8) results in the following equation:

\begin{align} D_{2} e^{n} - a (l ({\hat{U}}^{n})) Δ e^{n} + α | {\hat{U}}^{n} |^{p - 2} e^{n} = (a (l ({\hat{u}}^{n})) - a (l ({\hat{U}}^{n}))) Δ u^{n} \\ - α u^{n} (| {\hat{u}}^{n} |^{p - 2} - | {\hat{U}}^{n} |^{p - 2}) + (f ({\hat{u}}^{n}) - f ({\hat{U}}^{n})) + R^{n} . \end{align}

(3.16)

Multiply (3.17) by 4Δteⁿ and integrate it over Ω. The use of the telescope formula to the resulting equation leads to

\begin{array}{l} ‖ e^{n} ‖^{2} & - ‖ e^{n - 1} ‖^{2} + ‖ {\hat{e}}^{n + 1} ‖^{2} - ‖ {\hat{e}}^{n} ‖^{2} + ‖ e^{n} - 2 e^{n - 1} + e^{n - 2} ‖^{2} + 4 Δ t a (l ({\hat{U}}^{n})) ‖ \nabla e^{n} ‖^{2} \\ + 4 Δ t α \int_{Ω} | {\hat{U}}^{n} |^{p - 2} | e^{n} |^{2} d x = - 4 Δ t (a (l ({\hat{u}}^{n})) - a (l ({\hat{U}}^{n}))) (\nabla u^{n}, \nabla e^{n}) \\ - 4 Δ t α (u^{n} (| {\hat{u}}^{n} |^{p - 2} - | {\hat{U}}^{n} |^{p - 2}), e^{n}) + 4 Δ t (f ({\hat{u}}^{n}) - f ({\hat{U}}^{n}), e^{n}) + 4 Δ t (R^{n}, e^{n}) . \end{array}

Use the lower bound of a(⋅) and drop certain positive terms on the left hand side of the above equation leads to

\begin{align} ‖ e^{n} ‖^{2} & - ‖ e^{n - 1} ‖^{2} + ‖ {\hat{e}}^{n + 1} ‖^{2} - ‖ {\hat{e}}^{n} ‖^{2} + 4 m Δ t ‖ \nabla e^{n} ‖^{2} \leq 4 Δ t (a (l ({\hat{u}}^{n})) - a (l ({\hat{U}}^{n}))) (\nabla u^{n}, \nabla e^{n}) \\ + 4 Δ t α (u^{n} (| {\hat{u}}^{n} |^{p - 2} - | {\hat{U}}^{n} |^{p - 2}), e^{n}) + 4 Δ t (f ({\hat{u}}^{n}) - f ({\hat{U}}^{n}), e^{n}) + 4 Δ t (R^{n}, e^{n}) \\ = \sum_{k = 1}^{4} J_{k} . \end{align}

(3.17)

We have

\begin{array}{l} J_{1} = 4 Δ t (a (l ({\hat{u}}^{n})) - a (l ({\hat{U}}^{n}))) (\nabla u^{n}, \nabla e^{n}) & \leq 4 Δ t A ‖ {\hat{e}}^{n} ‖ ‖ \nabla u^{n} ‖ ‖ \nabla e^{n} ‖ (using (H 4)) \\ \leq C Δ t ‖ {\hat{e}}^{n} ‖^{2} + m Δ t ‖ \nabla e^{n} ‖^{2} \end{array}

\begin{array}{l} J_{2} = 4 Δ t α (u^{n} (| {\hat{u}}^{n} |^{p - 2} - | {\hat{U}}^{n} |^{p - 2}), e^{n}) & \leq C (‖ u^{n} ‖_{L^{\infty}}, ‖ {\hat{u}}^{n} ‖_{L^{\infty}}, ‖ {\hat{U}}^{n} ‖_{L^{\infty}}, p) Δ t ‖ {\hat{e}}^{n} ‖ ‖ e^{n} ‖ \\ \leq C Δ t (‖ {\hat{e}}^{n} ‖^{2} + ‖ e^{n} ‖^{2}) \end{array}

\begin{array}{l} J_{3} & = 4 Δ t (f ({\hat{u}}^{n}) - f ({\hat{U}}^{n}), e^{n}) \leq 4 γ Δ t ‖ {\hat{e}}^{n} ‖ ‖ e^{n} ‖ (using (H 2)) \\ \leq C Δ t (‖ {\hat{e}}^{n} ‖^{2} + ‖ e^{n} ‖^{2}) . \end{array}

J_{4} = 4 Δ t (R^{n}, e^{n}) \leq C Δ t (‖ R^{n} ‖^{2} + ‖ e^{n} ‖^{2}) .

Therefore,

‖ e^{n} ‖^{2} - ‖ e^{n - 1} ‖^{2} + ‖ {\hat{e}}^{n + 1} ‖^{2} - ‖ {\hat{e}}^{n} ‖^{2} + Δ t ‖ \nabla e^{n} ‖^{2} \leq C Δ t (‖ e^{n} ‖^{2} + ‖ {\hat{e}}^{n} ‖^{2} + ‖ R^{n} ‖^{2}) .

Summing up the above inequality and using the discrete Gronwall inequality, we get

‖ e^{n} ‖^{2} + ‖ {\hat{e}}^{n} ‖^{2} + Δ t \sum_{k = 1}^{n} ‖ \nabla e^{k} ‖^{2} \leq C Δ t^{4} .

From ‖eⁿ‖ ≤ CΔt², we have

\begin{array}{l} ‖ U^{n} ‖ \leq ‖ u^{n} ‖ + ‖ e^{n} ‖ \leq C, \\ ‖ D_{k} U^{n} ‖ \leq ‖ D_{k} u^{n} ‖ + ‖ D_{k} e^{n} ‖ \leq C, with k = 1, 2 . \end{array}

Applying Lemma 2.4 for the linear elliptic problems (3.1)–(3.3) with the induction assumptions gives the H² estimate

\begin{array}{l} ‖ U^{n} ‖_{2} & \leq C ‖ D_{k} U^{n} ‖ + C ‖ \nabla U^{n} ‖ + C ‖ | {\hat{U}}^{n} |^{p - 2} U^{n} ‖ + L ‖ {\hat{U}}^{n} ‖ \\ \leq C + C ‖ {\hat{U}}^{n} ‖_{\infty}^{p - 2} ‖ U^{n} ‖ \leq C . \end{array}

Using (2.3), we have

‖ U^{n} ‖_{\infty} \leq C ‖ U^{n} ‖^{1 / 2} ‖ U^{n} ‖_{2}^{1 / 2} \leq C

which concludes the proof.

4. Error estimates for the fully discrete problem

In this section, we will prove the optimal spatial error estimates. Let Π_h be an interpolation operator and $R_{h} : H_{0}^{1} (Ω) \to V_{h}$ be a Ritz projection operator defined by

\int_{Ω} \nabla (u - R_{h} u) \cdot \nabla w d x = 0, \forall w \in H_{0}^{1} (Ω) .

(4.1)

Then we have the following lemma.

Lemma 4.1

(cf. Ref. [19]). If $u \in H^{r + 1} (Ω) \cap H_{0}^{1} (Ω)$ , then

‖ u - Π_{h} u ‖ + h ‖ \nabla (u - Π_{h} u) ‖ \leq C h^{r + 1} ‖ u ‖_{H^{r + 1}},

(4.2)

‖ u - R_{h} u ‖ + h ‖ \nabla (u - R_{h} u) ‖ \leq C h^{r + 1} ‖ u ‖_{H^{r + 1}}

(4.3)

‖ v_{h} ‖_{\infty} \leq C h^{- d / 2} ‖ v_{h} ‖, \forall v_{h} \in V_{h} .

(4.4)

whereCis a positive constant that does not depend onhandr.

Let us denote

\begin{array}{l} E_{h}^{0} = u^{0} - U_{h}^{0} \\ {\bar{e}}_{h}^{1} = R_{h} {\bar{U}}^{1} - {\bar{U}}_{h}^{1}; {\hat{e}}_{h}^{n} = R_{h} {\hat{U}}^{n} - {\hat{U}}_{h}^{n}; e_{h}^{n} = R_{h} U^{n} - U_{h}^{n} \\ {\bar{E}}^{1} = R_{h} {\bar{U}}^{1} - {\bar{U}}^{1}; {\hat{E}}^{n} = R_{h} {\hat{U}}^{n} - {\hat{U}}^{n}; E^{n} = R_{h} U^{n} - U^{n} for 1 \leq n \leq N . \end{array}

From lemma 4.1, we have

‖ E_{h}^{0} ‖ + h ‖ \nabla E_{h}^{0} ‖ \leq C h^{r + 1}

(4.5)

‖ {\bar{E}}^{1} ‖ + ‖ {\hat{E}}^{n} ‖ + ‖ E^{n} ‖ + h (‖ \nabla {\bar{E}}^{1} ‖ + ‖ \nabla {\hat{E}}^{n} ‖ + ‖ \nabla E^{n} ‖) \leq C h^{r + 1} .

(4.6)

Lemma 4.2

Assume that the exact solution u of (1.1) satisfies the regularities (2.8). Then there exists a positive constant C independent of Δt and h such that

‖ {\hat{e}}_{h}^{1} ‖^{2} + Δ t ‖ \nabla {\hat{e}}_{h}^{1} ‖^{2} \leq C h^{2 r + 2}

(4.7)

‖ {\hat{U}}_{h}^{1} ‖ \leq C .

(4.8)

Proof. From equations (2.17) and (3.6), ${\hat{e}}_{h}^{1}$ satisfies the following equation:

\begin{array}{l} 2 ({\hat{e}}_{h}^{1}, v_{h}) + Δ t a (l ((U_{h}^{0})) (\nabla {\hat{e}}_{h}^{1}, \nabla v_{h}) + α Δ t (| U_{h}^{0} |^{p - 2} {\hat{e}}_{h}^{1}, v_{h}) = 2 (E_{h}^{0}, v_{h}) - 2 ({\hat{E}}^{1}, v_{h}) \\ - Δ t (a (l ((U^{0})) - a (l ((U_{h}^{0}))) (\nabla {\hat{U}}^{1}, \nabla v_{h}) - α Δ t (| U_{h}^{0} |^{p - 2} {\hat{E}}^{1}, v_{h}) \\ - α Δ t ((| U^{0} |^{p - 2} - | U_{h}^{0} |^{p - 2}) {\hat{U}}^{1}, v_{h}) + Δ t (f (U^{0}) - f (U_{h}^{0}), v_{h}) . \end{array}

Setting $v_{h} = {\hat{e}}_{h}^{1}$ in the above equations leads to

\begin{align} 2 ‖ {\hat{e}}_{h}^{1} ‖^{2} & + m Δ t ‖ \nabla {\hat{e}}_{h}^{1} ‖^{2} + α Δ t \int_{Ω} | U_{h}^{0} |^{p - 2} | {\hat{e}}_{h}^{1} |^{2} d x \leq 2 (E_{h}^{0} + {\hat{E}}^{1}, {\hat{e}}_{h}^{1}) \\ + Δ t (a (l ((U^{0})) - a (l ((U_{h}^{0}))) (\nabla {\hat{U}}^{1}, \nabla {\hat{e}}_{h}^{1}) + α Δ t (| U_{h}^{0} |^{p - 2} {\hat{E}}^{1}, {\hat{e}}_{h}^{1}) \\ + α Δ t ((| U^{0} |^{p - 2} - | U_{h}^{0} |^{p - 2}) {\hat{U}}^{1}, {\hat{e}}_{h}^{1}) + Δ t (f (U^{0}) - f (U_{h}^{0}), {\hat{e}}_{h}^{1}) \\ = \sum_{i = 1}^{5} L_{i} . \end{align}

(4.9)

From (4.5) and (4.6), we have

\begin{array}{l} L_{1} & = 2 (E_{h}^{0} + {\hat{E}}^{1}, {\hat{e}}_{h}^{1}) \leq C h^{2 r + 2} + \frac{1}{2} ‖ {\hat{e}}_{h}^{1} ‖^{2}, \\ L_{2} & = Δ t (a (l ((U^{0})) - a (l ((U_{h}^{0}))) (\nabla {\hat{U}}^{1}, \nabla {\hat{e}}_{h}^{1}) \leq Δ t A ‖ E_{h}^{0} ‖ ‖ \nabla {\hat{U}}^{1} ‖ ‖ \nabla {\hat{e}}_{h}^{1} ‖ \\ \leq C Δ t h^{2 r + 2} + \frac{m}{2} Δ t ‖ \nabla {\hat{e}}_{h}^{1} ‖^{2}, \\ L_{3} & = α Δ t (| U_{h}^{0} |^{p - 2} {\hat{E}}^{1}, {\hat{e}}_{h}^{1}) \leq α Δ t C (‖ U_{h}^{0} ‖_{\infty}) ‖ {\hat{E}}^{1} ‖ ‖ {\hat{e}}_{h}^{1} ‖ \leq C Δ t^{2} h^{2 r + 2} + \frac{1}{4} ‖ {\hat{e}}_{h}^{1} ‖^{2}, \\ L_{4} & = α Δ t ((| U^{0} |^{p - 2} - | U_{h}^{0} |^{p - 2}) {\hat{U}}^{1}, {\hat{e}}_{h}^{1}) \leq C (‖ U^{0} ‖_{\infty}, ‖ U_{h}^{0} ‖_{\infty}, ‖ {\hat{U}}^{1} ‖_{\infty}) Δ t ‖ E_{h}^{0} ‖ ‖ {\hat{e}}_{h}^{1} ‖ \\ \leq C Δ t^{2} h^{2 r + 2} + \frac{1}{4} ‖ {\hat{e}}_{h}^{1} ‖^{2}, \\ L_{5} & = Δ t (f (U^{0}) - f (U_{h}^{0}), {\hat{e}}_{h}^{1}) \leq γ Δ t ‖ E_{h}^{0} ‖ ‖ {\hat{e}}_{h}^{1} ‖ \leq C Δ t^{2} h^{2 r + 2} + \frac{1}{4} ‖ {\hat{e}}_{h}^{1} ‖^{2} . \end{array}

Combining these estimates into (4.9), we get (4.7).

From the inverse inequality, $‖ {\hat{e}}_{h}^{1} ‖_{\infty} \leq C h^{- d / 2} ‖ {\hat{e}}_{h}^{1} ‖ \leq C h^{r / 2}$ and

‖ {\hat{U}}_{h}^{1} ‖_{\infty} \leq ‖ R_{h} U^{1} ‖_{\infty} + ‖ {\hat{e}}_{h}^{1} ‖_{\infty} \leq C .

The main result in this section is as follows.

Theorem 4.1

Suppose that the exact solution u of (1.1) satisfies the regularities (2.8). Then there exists a positive constant C independent of Δt and h such that

\max_{1 \leq n \leq N} (‖ e_{h}^{n} ‖^{2} + ‖ {\hat{e}}_{h}^{n} ‖^{2} + Δ t \sum_{k = 1}^{n} ‖ \nabla e_{h}^{k} ‖^{2}) \leq C h^{2 r + 2},

(4.10)

\max_{0 \leq n \leq N} ‖ U_{h}^{n} ‖ \leq C .

(4.11)

Proof. The proof of this result will be done by mathematical induction. Since $U_{h}^{0} = Π_{h} U^{0}$ ⁠, (4.11) holds for n = 0. To compute the error estimate (4.10) for n = 1, subtract (3.4) from (2.16) and take $v_{h} = {\bar{e}}_{h}^{1} \equiv \frac{1}{2} (e_{h}^{1} + E_{h}^{0})$ ⁠,

\begin{align} ‖ e_{h}^{1} ‖^{2} & + 2 m Δ t ‖ \nabla {\bar{e}}_{h}^{1} ‖^{2} + 2 α Δ t \int_{Ω} | {\hat{U}}_{h}^{1} |^{p - 2} | {\bar{e}}_{h}^{1} |^{2} d x \leq (E^{1} + E_{h}^{0}, {\bar{e}}_{h}^{1}) \\ + 2 Δ t (a (l (({\hat{U}}^{1})) - a (l (({\hat{U}}_{h}^{1}))) (\nabla {\bar{U}}^{1}, \nabla {\bar{e}}_{h}^{1}) + 2 Δ t α (| {\hat{U}}_{h}^{1} |^{p - 2} {\bar{E}}^{1}, {\bar{e}}_{h}^{1}) \\ + 2 Δ t α ((| {\hat{U}}^{1} |^{p - 2} - | {\hat{U}}_{h}^{1} |^{p - 2}) {\bar{U}}^{1}, {\bar{e}}_{h}^{1}) + 2 Δ t (f ({\hat{U}}^{1}) - f ({\hat{U}}_{h}^{1}), {\bar{e}}_{h}^{1}) \\ = \sum_{i = 1}^{5} T_{i} . \end{align}

(4.12)

From (4.5) and (4.6), we have

\begin{array}{l} T_{1} & = (E^{1} + E_{h}^{0}, {\bar{e}}_{h}^{1}) \leq C h^{2 r + 2} + \frac{1}{8} ‖ e_{h}^{1} ‖^{2}, \\ T_{2} & = 2 Δ t (a (l (({\hat{U}}^{1})) - a (l (({\hat{U}}_{h}^{1}))) (\nabla {\bar{U}}^{1}, \nabla {\bar{e}}_{h}^{1}) \leq Δ t A ‖ {\hat{U}}^{1} - {\hat{U}}_{h}^{1} ‖ ‖ \nabla {\bar{U}}^{1} ‖ ‖ \nabla {\bar{e}}_{h}^{1} ‖ \\ \leq C Δ t h^{2 r + 2} + C Δ t ‖ {\hat{e}}_{h}^{1} ‖^{2} + \frac{m}{2} Δ t ‖ \nabla {\bar{e}}_{h}^{1} ‖^{2}, \\ T_{3} & = 2 α Δ t (| {\hat{U}}_{h}^{1} |^{p - 2} {\bar{E}}^{1}, {\bar{e}}_{h}^{1}) \leq α Δ t C (‖ {\hat{U}}_{h}^{1} ‖_{\infty}) ‖ {\bar{E}}^{1} ‖ ‖ {\bar{e}}_{h}^{1} ‖ \leq C Δ t^{2} h^{2 r + 2} + \frac{1}{8} ‖ e_{h}^{1} ‖^{2}, \\ T_{4} & = 2 α Δ t ((| {\hat{U}}^{1} |^{p - 2} - | {\hat{U}}_{h}^{1} |^{p - 2}) {\bar{U}}^{1}, {\bar{e}}_{h}^{1}) \leq C (‖ {\hat{U}}^{1} ‖_{\infty}, ‖ {\hat{U}}_{h}^{1} ‖_{\infty}, ‖ {\bar{U}}^{1} ‖_{\infty}) Δ t ‖ {\hat{U}}^{1} - {\hat{U}}_{h}^{1} ‖ ‖ {\bar{e}}_{h}^{1} ‖ \\ \leq C Δ t^{2} h^{2 r + 2} + C Δ t^{2} ‖ {\hat{e}}_{h}^{1} ‖^{2} + \frac{1}{8} ‖ e_{h}^{1} ‖^{2}, \\ T_{5} & = 2 Δ t (f ({\hat{U}}^{1}) - f ({\hat{U}}_{h}^{1}), {\bar{e}}_{h}^{1}) \leq γ Δ t ‖ {\hat{U}}^{1} - {\hat{U}}_{h}^{1} ‖ ‖ {\bar{e}}_{h}^{1} ‖ \leq C Δ t^{2} h^{2 r + 2} + C Δ t^{2} ‖ {\hat{e}}_{h}^{1} ‖^{2} + \frac{1}{8} ‖ e_{h}^{1} ‖^{2} . \end{array}

Taking these estimates into (4.12) and using Lemma 4.2, we conclude that

‖ e_{h}^{1} ‖^{2} + Δ t ‖ \nabla {\bar{e}}_{h}^{1} ‖^{2} \leq C h^{2 r + 2}

(4.13)

which proves (4.10) for n = 1.

Now, we assume that (4.10) and (4.11) hold for n = m − 1, 2 ≤ m ≤ N, then we need to show it also holds for n = m. By the definition of ${\hat{U}}_{h}^{n}$ and the induction assumption, $‖ {\hat{U}}_{h}^{n} ‖_{\infty} \leq C$ ⁠.

Subtracting (2.18) from (3.5), we obtain

\begin{array}{l} (D_{2} e_{h}^{n}, v_{h}) + a (l (({\hat{e}}_{h}^{n})) (\nabla e_{h}^{n}, \nabla v_{h}) + α (| {\hat{e}}_{h}^{n} |^{p - 2} e_{h}^{n}, v_{h}) = - (D_{2} E^{n}, v_{h}) \\ - (a (l (({\hat{U}}^{n})) - a (l (({\hat{U}}_{h}^{n}))) (\nabla U^{n}, \nabla v_{h}) - α (| {\hat{U}}_{h}^{n} |^{p - 2} E^{n}, v_{h}) \\ - α ((| {\hat{U}}^{n} |^{p - 2} - | {\hat{U}}_{h}^{n} |^{p - 2}) U^{n}, v_{h}) + (f ({\hat{U}}^{n}) - f ({\hat{U}}_{h}^{n}), v_{h}) . \end{array}

If one takes $v_{h} = 4 Δ t e_{h}^{n}$ and uses the telescope formula, one obtains

\begin{array}{l} ‖ e_{h}^{n} ‖^{2} & - ‖ e_{h}^{n - 1} ‖^{2} + ‖ {\hat{e}}_{h}^{n + 1} ‖^{2} - ‖ {\hat{e}}_{h}^{n} ‖^{2} + ‖ e_{h}^{n} - 2 e_{h}^{n - 1} + e_{h}^{n - 2} ‖^{2} + 4 Δ t a (l ({\hat{U}}_{h}^{n})) ‖ \nabla e_{h}^{n} ‖^{2} \\ + 4 Δ t α \int_{Ω} | {\hat{U}}_{h}^{n} |^{p - 2} | e_{h}^{n} |^{2} d x = - 4 Δ t (D_{2} E^{n}, e_{h}^{n}) - 4 Δ t (a (l ({\hat{U}}^{n})) - a (l ({\hat{U}}_{h}^{n}))) (\nabla U^{n}, \nabla e_{h}^{n}) \\ - 4 Δ t α (| {\hat{U}}_{h}^{n} |^{p - 2} E^{n}, e_{h}^{n}) - 4 Δ t α ((| {\hat{U}}^{n} |^{p - 2} - | {\hat{U}}_{h}^{n} |^{p - 2}) U^{n}, e_{h}^{n}) + 4 Δ t (f ({\hat{U}}^{n}) - f ({\hat{U}}_{h}^{n}), e_{h}^{n}) . \end{array}

That is

\begin{align} ‖ e_{h}^{n} ‖^{2} & - ‖ e_{h}^{n - 1} ‖^{2} + ‖ {\hat{e}}_{h}^{n + 1} ‖^{2} - ‖ {\hat{e}}_{h}^{n} ‖^{2} + 4 m Δ t ‖ \nabla e_{h}^{n} ‖^{2} \leq 4 Δ t (D_{2} E^{n}, e_{h}^{n}) \\ + 4 Δ t (a (l ({\hat{U}}^{n})) - a (l ({\hat{U}}_{h}^{n}))) (\nabla U^{n}, \nabla e_{h}^{n}) + 4 Δ t α (| {\hat{U}}_{h}^{n} |^{p - 2} E^{n}, e_{h}^{n}) \\ + 4 Δ t α ((| {\hat{U}}^{n} |^{p - 2} - | {\hat{U}}_{h}^{n} |^{p - 2}) U^{n}, e_{h}^{n}) + 4 Δ t (f ({\hat{U}}^{n}) - f ({\hat{U}}_{h}^{n}), e_{h}^{n}) \\ = \sum_{i = 1}^{5} K_{i} . \end{align}

(4.14)

The quantities K_i, i = 1, …, 5 can be bounded by the similar way T_i, i = 1, …, 5:

\begin{array}{l} K_{1} & = 4 Δ t (D_{2} E^{n}, e_{h}^{n}) \leq C Δ t (‖ D_{2} E^{n} ‖^{2} + ‖ e_{h}^{n} ‖^{2}) \leq C Δ t h^{2 r + 2} + C Δ t ‖ e_{h}^{n} ‖^{2}, \\ K_{2} & = 4 Δ t (a (l (({\hat{U}}^{n})) - a (l (({\hat{U}}_{h}^{n}))) (\nabla U^{n}, \nabla e_{h}^{n}) \leq Δ t A ‖ {\hat{E}}^{n} + {\hat{e}}_{h}^{n} ‖ ‖ \nabla U^{n} ‖ ‖ \nabla e_{h}^{n} ‖ \\ \leq C Δ t h^{2 r + 2} + C Δ t ‖ {\hat{e}}_{h}^{n} ‖^{2} + m Δ t ‖ \nabla e_{h}^{n} ‖^{2}, \\ K_{3} & = 4 α Δ t (| {\hat{U}}_{h}^{n} |^{p - 2} E^{n}, e_{h}^{n}) \leq 4 α Δ t C (‖ {\hat{U}}_{h}^{n} ‖_{\infty}) ‖ E^{n} ‖ ‖ e_{h}^{n} ‖ \leq C Δ t h^{2 r + 2} + C Δ t ‖ e_{h}^{n} ‖^{2}, \\ K_{4} & = 4 α Δ t ((| {\hat{U}}^{n} |^{p - 2} - | {\hat{U}}_{h}^{n} |^{p - 2}) U^{n}, e_{h}^{n}) \leq C (‖ {\hat{U}}^{n} ‖_{\infty}, ‖ {\hat{U}}_{h}^{n} ‖_{\infty}, ‖ U^{n} ‖_{\infty}) Δ t ‖ {\hat{E}}^{n} + {\hat{e}}_{h}^{n} ‖ ‖ e_{h}^{n} ‖ \\ \leq C Δ t h^{2 r + 2} + C Δ t (‖ {\hat{e}}_{h}^{n} ‖^{2} + ‖ e_{h}^{n} ‖^{2}), \\ K_{5} & = 4 Δ t (f ({\hat{U}}^{n}) - f ({\hat{U}}_{h}^{n}), e_{h}^{n}) \leq 4 γ Δ t ‖ {\hat{E}}^{n} + {\hat{e}}_{h}^{n} ‖ ‖ e_{h}^{n} ‖ \leq C Δ t h^{2 r + 2} + C Δ t (‖ {\hat{e}}_{h}^{n} ‖^{2} + ‖ e_{h}^{n} ‖^{2}) . \end{array}

Taking these bounds into (4.14), we obtain

‖ e_{h}^{n} ‖^{2} - ‖ e_{h}^{n - 1} ‖^{2} + ‖ {\hat{e}}_{h}^{n + 1} ‖^{2} - ‖ {\hat{e}}_{h}^{n} ‖^{2} + 4 m Δ t ‖ \nabla e_{h}^{n} ‖^{2} \leq C Δ t h^{2 r + 2} + C Δ t (‖ {\hat{e}}_{h}^{n} ‖^{2} + ‖ e_{h}^{n} ‖^{2}) .

Sum up the above inequality and use the discrete Gronwall Lemma 2.3 leads to

‖ e_{h}^{n} ‖^{2} + ‖ {\hat{e}}_{h}^{n} ‖^{2} + Δ t \sum_{k = 1}^{n} ‖ \nabla e_{h}^{k} ‖^{2} \leq \exp (C T) h^{2 r + 2} .

From the inverse inequality, $‖ e_{h}^{n} ‖_{\infty} \leq C h^{- d / 2} ‖ e_{h}^{n} ‖ \leq C h^{r / 2}$ and

‖ U_{h}^{n} ‖_{\infty} \leq ‖ R_{h} U^{n} ‖_{\infty} + ‖ e_{h}^{n} ‖_{\infty} \leq C .

5. Numerical results

In this section, we present several numerical simulations to illustrate our theoretical analysis. Since the resulting matrix of the linear system (2.16)–(2.18) is sparse, symmetric and positive definite, an incomplete Cholesky factorization is performed and the result is used as preconditioner in the preconditioned conjugate method iterative solver (see for instance Refs [20, 21]).

To analyze the convergence rate, we consider the following problem.

\{\begin{cases} u_{t} - a (l (u)) Δ u + α | u |^{p - 2} u = f (u) + g in Ω \times (0, T] \\ u (x, t) = 0 on \partial Ω \times (0, T] \\ u (x, 0) = u_{0} (x) in Ω, \end{cases}

(5.1)

with Ω = (0,1)², the coefficient α = 1, p = 3.5

f (s) = s (10 - s), a (s) = 3 + \cos (s), l (u) = \int_{Ω} u d x .

g is chosen correspondingly to the exact solution

u (x, y, t) = 2 (1 + t^{2} e x p (- t)) x y (1 - x) (1 - y) .

We simulated the above problem on uniform meshes with a linear finite element approximation (r = 1) and T = 0.1.

For the convergence with respect to the mesh size h, we choose Δt = h² and we solve problem (2.16)–(2.18) with different values of h (h = 1/5; 1/10; 1/15; 1/20; 1/25); from our theoretical analysis, the L²-norm errors are in order O(h² + Δt²) = O(h² + h⁴) ∼ O(h²). H¹-norm errors are in order O(h + Δt²) = O(h + h⁴) ∼ O(h). In Figure 1, we plot the log of errors against log(h). One can see that for L²-norm, the slope is almost 2, and for H¹ − norm, the slope is almost 1, which are in good agreement with our theoretical analysis.

Figure 1

View large Download slide

Convergence rate with respect to the mesh size h in L² and H¹ norm

For the convergence with respect to the time step Δt, h is fixed (h = 0.01), and we solve problem (2.16)–(2.18) with different time steps Δt = 0.1; 0.05; 0.025; 0.0125 (Δt = 0.1 × 2^1−l, l = 1, …, 4), and the L²-norm errors are in order O(h² + Δt²) ∼ O(Δt²). Figure 2 shows the plots of log L²-error norm against log(Δt). Again, one can see that the slope is almost 2. These results are consistent with our theoretical analysis.

Figure 2

View large Download slide

Convergence rate with respect to the time step Δt in L² norm

6. Conclusion

We have presented and analyzed a linearized second-order BDF Galerkin finite element method for the nonlocal parabolic problems. We have proved the L² and energy error estimates using sufficient conditions on the exact solution. We also presented some numerical experiments on Matlab’s environment, and our numerical results confirm the theoretical analysis. The results in this paper lay the foundation for developing finite element based numerical methods for more general and complicated nonlocal problems both stationary and evolutionary.

References

1

Chipot

M

,

Savitska

T

.

Nonlocal p-Laplace equations depending on the L^p norm of the gradient

.

Adv Differ Equ

.

2014

;

19

(

11/12

):

997

-

1020

.

Google Scholar

2.

Robalo

RJ

,

Almeida

RM

,

do Carmo Coimbra

M

,

Ferreira

J

.

A reaction–diffusion model for a class of nonlinear parabolic equations with moving boundaries: existence, uniqueness, exponential decay and simulation

.

Appl Math Model

.

2014

;

38

(

23

):

5609

-

22

.

Google Scholar

Crossref

3.

Mbehou

M

.

The Euler-Galerkin finite element method for nonlocal diffusion problems with a P-Laplace-type operator

.

Appl Anal

.

2019

;

98

(

11

):

2031

-

47

.

Google Scholar

Crossref

4.

Simsen

J

,

Ferreira

J

.

A global attractor for a nonlocal parabolic problem

.

Nonlinear Stud

.

2014

;

21

(

3

):

405

-

16

.

Google Scholar

5.

Mbehou

M

,

Chendjou

G

.

Numerical methods for a nonlocal parabolic problem with nonlinearity of Kirchhoff type

.

Numer Anal Appl

.

2019

;

12

(

3

):

251

-

62

.

Google Scholar

Crossref

6.

Sharma

N

,

Khebchareon

M

,

Sharma

K

,

Pani

AK

.

Finite element g Alerkin approximations to a class of nonlinear and nonlocal parabolic problems

.

Numer Methods Partial Differ Equ

.

2016

;

32

(

4

):

1232

-

64

.

Google Scholar

Crossref

7.

Chaudhary

S

,

Srivastava

V

,

Srinivas Kumar

V

.

Finite element scheme with Crank–Nicolson method for parabolic nonlocal problems involving the Dirichlet energy

.

Int J Comput Methods

.

2016

:

1

-

24

.

Google Scholar

8.

Mbehou

M.

,

Maritz

R.

,

Tchepmo

P

.

Numerical analysis for a nonlocal parabolic problem

.

East Asian J Appl Math

.

2016

;

6

(

4

):

434

-

47

.

Google Scholar

Crossref

9.

Yin

Z.

,

Xu

Q

.

A fully discrete symmetric finite volume element approximation of nonlocal reactive flows in porous media

.

Math Probl Eng

.

2013

;

2013

:

1

-

7

.

Google Scholar

10.

Almeida

RM

,

Duque

JC

,

Ferreira

J

,

Robalo

RJ

.

Finite element schemes for a class of nonlocal parabolic systems with moving boundaries

.

Appl Numer Math

.

2018

;

127

:

226

-

48

.

Google Scholar

Crossref

11.

Yang

Y-B

,

Jiang

Y-L

,

Yu

B-H

.

Unconditional optimal error estimates of linearized, decoupled and conservative Galerkin FEMs for the Klein–Gordon–Schrödinger equation

.

J Sci Comput

.

2021

;

87

(

3

):

1

-

32

.

Google Scholar

Crossref

12.

Yang

Y-B

,

Jiang

Y-L

.

Unconditional optimal error estimates of linearized second-order BDF Galerkin FEMs for the Landau-Lifshitz equation

.

Appl Numer Math

.

2021

;

159

:

21

-

45

.

Google Scholar

Crossref

13

Brezis

H

.

Functional analysis, Sobolev spaces and partial differential equations

.

NY

:

Springer New York

;

2010

.

Google Scholar

14

Boffi

D

,

Brezzi

F

,

Fortin

M

.

Mixed finite element methods and applications

.

Berlin, Heidelberg

:

Springer-Verlag

;

2013

.

Google Scholar

Crossref

15.

Lions

JL

.

Quelques Méthodes de Résolution des Problmes aux Limites Non Linéaires

.

Paris

:

Dunod

;

1969

.

Google Scholar

16.

Barrett

JW

,

Liu

W

.

Finite element approximation of the p-Laplacian

.

Math Comput

.

1993

;

61

(

204

):

523

-

37

.

Google Scholar

17.

Heywood

JG.

,

Rannacher

R

,

Turek

S

.

Artificial boundaries and flux and pressure conditions for the incompressible Navier-Stokes equations

.

Int J Numer Methods Fluids

.

1996

;

22

(

5

):

325

-

52

.

Google Scholar

Crossref

18.

Chen

Y-Z

,

Wu

L-C

.

Second order elliptic equations and elliptic systems

. Vol.

174

.

American Mathematical Soc.

,

Providence, Rhode Island

;

1998

.

Google Scholar

Crossref

19

Thomée

V

.

Galerkin finite element methods for parabolic problems

.

Berlin, Heidelberg

:

Springer

;

1984

.

1054

.

Google Scholar

20.

Koko

J

.

A MATLAB mesh generator for the two-dimensional finite element method

.

Appl Math Comput

.

2015

;

250

:

650

-

64

.

Google Scholar

21.

Koko

J

.

Efficient MATLAB codes for the 2d/3d stokes equation with the mini-element

.

Informatica

.

2019

;

30

(

2

):

243

-

68

.

Google Scholar

Crossref

2023

M.S. Daoussa Haggar and M. Mbehou

Published in the Arab Journal of Mathematical Sciences. Published by Emerald Publishing Limited. This article is published under the Creative Commons Attribution (CC BY 4.0) license. Anyone may reproduce, distribute, translate and create derivative works of this article (for both commercial and non-commercial purposes), subject to full attribution to the original publication and authors. The full terms of this license may be seen at http://creativecommons.org/licences/by/4.0/legalcode

Optimal error estimates of a linearized second-order BDF scheme for a nonlocal parabolic problem

1. Introduction

2. Preliminaries and main results

3. Error estimates for the semi-discrete problem

4. Error estimates for the fully discrete problem

5. Numerical results

6. Conclusion

References

Email Alerts

Cited By

Optimal error estimates of a linearized second-order BDF scheme for a nonlocal parabolic problem Open Access

1. Introduction

2. Preliminaries and main results

3. Error estimates for the semi-discrete problem

4. Error estimates for the fully discrete problem

5. Numerical results

6. Conclusion

References

Email Alerts

Suggested Reading

Related Chapters

Recommended for you

Cited By

Sharing Unavailable

Optimal error estimates of a linearized second-order BDF scheme for a nonlocal parabolic problem