# Finite difference/element method for time-fractional Navier-Stokes equations

###### Abstract.

We apply a composite idea of semi-discrete finite difference approximation in time and Galerkin finite element method in space to solve the Navier-Stokes equations with Caputo derivative of order . The stability properties and convergence error estimates for both the semi-discrete and fully discrete schemes are obtained. Numerical example is provided to illustrate the validity of theoretical results.

Keywords: Time-fractional Navier-Stokes equations; finite difference approximation; finite element method; error estimates; numerical examples

AMS Subject Classification: 76D05, 65N30, 65N12

Faculty of Mathematics and Computational Science, Xiangtan University

Xiangtan, Hunan 411105, P.R. China

Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Faculty of Science

King Abdulaziz University, Jeddah 21589, Saudi Arabia

e-mail: zouguangan@henu.edu.cn, yzhou@xtu.edu.cn, bashirahmad_qau@yahoo.com, aalsaedi@hotmail.com

Manuscript received xx

## 1. Introduction

In this paper, we study the following Navier-Stokes equations with time-fractional derivative in a bounded subset of with a smooth boundary :

(1.1) |

where represents the Caputo-type fractional derivative of order , denotes the velocity field at a point and , is viscosity coefficient, represents the pressure field, is the external force and is the initial velocity.

Notice that the problem (1.1) reduces to the classical Navier-Stokes equations (NSEs) for . The existence and non-existence of solutions for the NSEs have been discussed in [1]. Chemin et al.[2] studied the global regularity for the large solutions to the NSEs. Miura [3] focused on the uniqueness of mild solutions to the NSEs. Germain [4] presented the uniqueness criteria for the solutions of the Cauchy problem associated to the NSEs. The existence of global weak solutions for supercritical NSEs was discussed [5]. The lower bounds on blow up solutions for the NSEs in homogeneous Sobolev spaces were studied in [6]. The numerical methods for solving the NSEs have been investigated by many authors [7,8,9,10,11,12,23]. The study of time-fractional Navier-Stokes equations (TFNSEs) has become a hot topic of research due to its significant role in simulating the anomalous diffusion in fractal media. There are also some analytical methods available for solving the TFNSEs. Momani and Odibat [13] applied Adomian decomposition method to obtain the analytical solution of the TFNSEs. In [14, 15], the homotopy perturbation (transform) method was used to find the analytical solution of the TFNSEs. Wang and Liu [16] solved TFNSEs by applying the transform methods. Concerning the existence of global and local mild solutions to TFNSEs, see Carvalho-Neto and Gabriela [17], Zhou and Peng [18]. Moreover, Zhou and Peng [19] investigated the existence of weak solutions and optimal control for TFNSEs, while Peng et al.[20] presented the rigorous exposition of local solutions of TFNSEs in Sobolev space. However, one can notice that there are only a few works related to the numerical solution of the TFNSEs. The details of meshless local Petrov-Galerkin method based on moving Kriging interpolation for solving the TFNSEs can be found in the literature [21]. The purpose of this paper is to present finite difference/element method to obtain the numerical solution of TFNSEs.

The rest of the paper is arranged as follows. In Section 2, we give some notations and preliminaries. Section 3 deals with a semi-discrete scheme for the TFNSEs, which is based on a mixed finite element method in space. We also discuss the stability and error estimates of this semi-discrete scheme. In Section 4, we use a finite difference approximation to discrete time direction to obtain the fully discrete scheme. The stability and error estimates for the discrete schemes are also found. In Section 5, numerical results are discussed to confirm our theoretical analysis. Conclusions are given in the final section.

## 2. Notations and preliminaries

In this section, we present some preliminary concepts of the functional spaces. Firstly, we introduce the following Hilbert spaces:

where the space is associated with the usual inner product and the norm . The space is associated with the following inner product and equivalent norm:

Denote by and the closed subsets of and respectively, which are given by

We denote the Stokes operator by , in which is the -orthogonal projection of onto . The domain of is and let with the norm . Observe that , and .

Next, we define the Riemann-Liouville fractional integral operator of order () as (see [22])

(2.1) |

with .

The Caputo-type derivative of order , in (1.1) is defined by

(2.2) |

Further, the operator is defined as

(2.3) |

where stands for the gamma function .

Next we introduce the following continuous bilinear forms and on and respectively as follows:

and the trilinear form on is given by

It is well-known that the trilinear form has the following properties:

In terms of the above notations, the weak formulation of problem (1.1) is as follows: find for all such that for all :

(2.4) |

In [19], Zhou and Peng discussed the existence and uniqueness of weak solutions for the problem (2.4). The objective of the present work is to obtain the numerical solution of the problem at hand.

## 3. Finite element method for space discretization

Let be a mesh of with a mesh size function , which is the diameter of element containing . Assuming be the largest mesh size of , we introduce the mixed finite element subspace of and define the subspace of as

Let denote the -orthogonal projection defined by

With the above notations, we need some further basic assumptions on the mixed finite element spaces (Refs.[8,10,12]).

(A1) Approximation. For each , there exist approximations such that

(3.1) |

(A2) Inverse estimate. For any , the following relations hold:

(3.2) |

(A3) Stability property. For any , the well-known inf-sup condition holds:

(3.3) |

where is a constant.

Further, the following classical properties hold:

(3.4) | |||

(3.5) |

The standard finite element Galerkin approximation for (2.4) holds as follows: Find for all such that for all , we have

(3.6) |

With the above semi-discrete approximation, a discrete analogue of the Stokes operator is defined as via the condition for all . The trilinear form satisfies the following properties:

(3.7) | |||

(3.8) |

Theorem 3.1. For any and , let be the solution of equation (3.6). Then there exists a positive constant such that

where is constant.

Proof. Taking , in (3.6) and using the Young’s inequality, we get

(3.9) |

where denotes the dual operator in .

Applying the integral operator (2.3) to both sides of (3.9) and using the Young’s inequality, we obtain

where with .

In view of the inequality

it follows that

This completes the proof.

Theorem 3.2. For any , , let and be the solutions of equations (2.4) and (3.6) respectively. Then there exists a positive constant such that

(3.10) |

Proof. Setting , we deduce from (2.4) and (3.6) that

(3.11) |

Taking and in (3.11), we get

Using the properties of together with Young’s inequality, we obtain

(3.12) |

Applying (2.3) to both sides of (3.12), we have

(3.13) |

By means of the generalized integral version of Gronwall’s lemma [24], we get

(3.14) |

Furthermore, setting and in (3.11), and using inf-sup condition (3.3), combining (3.1)-(3.5),(3.8),(3.14) and using the integral operator (2.3) in (3.15), we conclude that

(3.15) |

This completes the proof.

## 4. Finite difference method for time discretization

The discretization of time-fractional derivative can be found in [25-32] and references therein. Here, we will introduce a uniform grid by discretizing the temporal domain given by the points: for , with the time-step size . Hence, the Riemann-Liouville fractional integral operator of order can be discretized as follows:

(4.1) |

where and the truncation error is given by

Therefore, we have

Lemma 4.1. (see [25]) If , then

(4.2) |

where .

Lemma 4.2. (see [25]) For and , let the coefficient be given by (4.1). Then

Applying the integral operator (2.3) to both sides of (3.6), we obtain

(4.3) |

Let and be the numerical solutions of and at respectively. By (4.2) and (4.3), our full discrete scheme of equation (2.4) can be defined by seeking such that for all :

(4.4) |

where .

Theorem 4.1. For any , the full discrete scheme (4.4) is unconditionally stable, and that

Proof. Setting in (4.4), we get

(4.5) |

Taking and in (4.5), we have

Making use of Young’s inequality, we obtain

that is,

Assuming and , the following inequality holds

Setting and in (4.4), we get

By the elementary identity and the Young’s inequality, we have

Together with Lemma 4.2 (ii) (that is, ) and , we obtain

(4.6) |

The proof is completed.

Lemma 4.3. Let be the viscosity coefficient and that

(4.7) |

Then

(4.8) |

where is defined by (3.8) and is constant.

Proof. Making use of (4.6) and (4.7), we get

that is,

By the property of , we get

The proof of the lemma is completed.

Theorem 4.2. For , let and be the solutions of equations (3.6) and (4.4) respectively. There exists a constant such that

(4.8) |

Proof. Let and . Using (4.2)-(4.4) and noting , we deduce

(4.9) |

For , taking and in (4.9), we have

By Cauchy-Schwarz inequality and Lemma 4.3, we get

(4.10) |

Let us assume that for . In order to show that the first inequality in (4.8) holds for