The Optimal Error Estimate of the Fully Discrete Locally Stabilized Finite Volume Method for the Non-Stationary Navier-Stokes Problem

Abstract

This paper proves the optimal estimations of a low-order spatial-temporal fully discrete method for the non-stationary Navier-Stokes Problem. In this paper, the semi-implicit scheme based on Euler method is adopted for time discretization, while the special finite volume scheme is adopted for space discretization. Specifically, the spatial discretization adopts the traditional triangle $P_1-P_0$ trial function pair, combined with macro element form to ensure local stability. The theoretical analysis results show that under certain conditions, the full discretization proposed here has the characteristics of local stability, and we can indeed obtain the optimal theoretic and numerical order error estimation of velocity and pressure. This helps to enrich the corresponding theoretical results.

1. Introduction

Recently, due to the characteristics of simple implementation, the finite volume method has been widely used in many scientific research and engineering fields. It has obtained many ideal numerical simulation and calculation results, and often used to solve the complex engineering calculation problems well. Nevertheless, compared with a wide range of application scenarios, its theoretical analysis, such as stability and convergence analysis, is far behind, which inevitably shadow benefits of the finite volume method, so it needs to be studied continuously. Among them, the theoretical analysis of Navier Stokes process is one of the important field.

Finite volume method is an effective method to solve differential equations. In the past several decades, the calculation methods used to solve Navier-Stokes problems have developed rapidly. The results are richer and richer, but to our dismay, the theoretical analysis of the algorithm is still insufficient [1,2]. As we all know, based on the current advanced computing equipment, simple numerical methods are easy to distribute and suitable for large-scale computing, which makes them the hope of solving complex problems, such as incompressible Navier-Stokes equations. Among them, using simple coordinated low-order elements and local stabilization method is a good choice [3,4].

However, it is well known that low order coordinated finite volume function pairs, such as $P_1-P_0$, are unstable for numerical solution of Navier-Stokes equations. The common way to overcome this shortage is to use local stabilization technique, that is, add “macro element condition” to improve the stability of the algorithm. This kind of low order method has been widely analyzed and applied, and has been proved to be effective in practice [5,6,7,8,9]. The basic idea of this method was first proposed by Boland and Nicolaides, and has been vigorously developed since then. The recent work of Wen [10] and Li [11] paves the way for the numerical analysis of Stokes and Navier Stokes problems. In addition, He [12] and Li [13] have given the locally stable finite volume method for partial spatial discretization of Navier-Stokes problems, which has a good effect except some fully discrete results.

This paper continues to analyze the convergence results of using FVM to solve two-dimensional time-dependent Navier-Stokes equations, so as to enrich the relevant theories. Here we will still focus on $P_1-P_0$ pairs. For this purpose, let’s assume $T_h$ is the uniform and regular triangulation of $\overline{\mathsf{\Omega }}$. It should be reminded that the finite element space here does not have the inf-sup condition for $(x_h,m_h)$, so a similar skill in the paper [11,12,13] is required. Because these papers mainly discuss the spatial discrete case, this paper studies a approximation based on time-discretization is Euler semi-implicit and space-discretization is $P_1-P_0$ Locally stable FVM.

For brevity, this paper assumes $T_h$ is a regular triangulation that satisfies the general regular condition [14,15]. Let the mesh size is h and the time step is $0<\Delta t<1$, the theoretical results of the optimal order convergence of the fully discrete FVM based on the low order coordinated finite element local stabilization are as follows:

For such a finite volume solution $(u_h^n,p_h^n)$ obtained by the fully discrete locally stabilized FVM, we derive in this paper the following order error estimates:

$$\Delta t\sum _{k=1}^N\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_k\right)-u_h^k){\parallel }_0^2+\tau \left(t_n\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h^n)\parallel }_0^2\le \kappa (h^2+\Delta t),$$

(1)

$$\begin{array}{c}\hfill \Delta t\sum _{k=1}^N\tau \left(t_k\right)\parallel p\left(t_k\right)-p_h^k{\parallel }_0^2+{\tau }^2\left(t_n\right){\parallel p\left(t_n\right)-p_h^n\parallel }_0^2\le \kappa (h^2+\Delta t),\end{array}$$

(2)

where $N=\frac{T}{\Delta t},t_n=n\Delta t\in (0,T],\tau \left(t\right)=min\{1,t\}$.

The rest of the paper is organized as follows. Section 2 introduces some basic concepts and function definitions related to Navier-Stokes problem and the locally stable FVM. Some basic results are prepared in Section 3. Section 4 mainly analyzes the error estimation of time semi-discretization based on Euler semi-implicit scheme. Section 5 proves the error estimations of time and spatial fully discretization. Section 6 contains some numerical results and a summary of the article is included in Section 7.

2. Foundation of Finite Volume Method for the Navier-Stokes Problem

This article consider the follow non-stationary Navier-Stokes equations

$$\left\{\begin{array}{c}{u}_t-\nu \Delta u+(u·\nabla )u+\nabla p=f,divu=0\forall (x,t)\in \mathsf{\Omega }\times (0,T],\hfill \\ u(x,0)=u_0{\left(x\right)\forall x\in \mathsf{\Omega };u(x,t)|}_{\partial \mathsf{\Omega }}=0\forall t\in (0,T],\hfill \end{array}\right.$$

(3)

where $\mathsf{\Omega }$ be a bounded domain in $R^2$ assumed to have a common continuous boundary $\partial \mathsf{\Omega }$ (stronger than Lipschitz continuity) [14,16]; $u=u(x,t)=(u_1(x,t),u_2(x,t))$ are the velocity vector and $p=p(x,t)$ is the pressure, $f=f(x,t)$ is the body force, $u_0\left(x\right)$ is the initial velocity, and $\nu >0$ is the viscosity.

For the convenience of analysis of problem (3), we introduce the following common abbreviations

$$X=H_0^1{\left(\mathsf{\Omega }\right)}^2,Y=L^2{\left(\mathsf{\Omega }\right)}^2,M=L_0^2\left(\mathsf{\Omega }\right)=\{q\in L^2\left(\mathsf{\Omega }\right);{\int }_{\mathsf{\Omega }}qdx=0\},$$

$$H=\{v\in Y;\mathrm{div}v=0,v·n{|}_{\Gamma }=0\}.V=\{v\in X;divv=0\}.$$

where V is the closed subset of X and H is denote the closed subset of Y, respectively. The spaces $L^2{\left(\mathsf{\Omega }\right)}^m,m=1,2,4$, are endowed with the common norm denoted by $(·,·)$ and ${\parallel ·\parallel }_0$, respectively. The norms of the Hilbert space $H_0^1\left(\mathsf{\Omega }\right)$ and X are

$$\left((u,v)\right)=(\nabla u,\nabla v),\parallel u\parallel ={\left((u,u)\right)}^{1/2}.$$

For more information about the above marks, we refer the reader to [14,16,17]. We also need to denote $A=-\Delta$ as the Laplace operator and $\overline{A}=-P\Delta$ as the Stokes operator, where P is the $L^2$-orthogonal projection of Y onto H.

It is known [17] that

$${\parallel Av\parallel }_0^2\le {\parallel v\parallel }_2^2\le c{\parallel \overline{A}v\parallel }_0^2\forall v\in H^2{\left(\mathsf{\Omega }\right)}^2\cap V,$$

(4)

$${\parallel v\parallel }_0^2\le {\gamma }_0{\parallel v\parallel }^2\forall v\in {X,\parallel v\parallel }^2\le {\gamma }_0{\parallel v\parallel }_2^2\le c{\parallel Av\parallel }_0^2\forall v\in D\left(A\right),$$

(5)

where $D\left(A\right)=H^2{\left(\mathsf{\Omega }\right)}^2\cap X$, ${\gamma }_0$ is positive constant depending only on $\mathsf{\Omega }$. C, like the quantities $C_i,i=1,2,\cdots$, appear subsequently, is a positive constant depending on $\mathsf{\Omega }$.

This paper uses following kind of continuous bilinear forms $a(·,·)$ and $d(·,·)$ on $X\times X$ and $X\times M$, respectively,

$$a(u,v)=\nu \left(\right(u,v\left)\right)\forall u,v\in X,$$

$$d(v,q)=-(v,\nabla q)=(q,\mathrm{div}v)\forall v\in X,q\in M.$$

and a trilinear form on $X\times X\times X$

$$\begin{array}{cc}\hfill b(u,v,w)& =\frac{1}{2}((u·\nabla )v,w)-\frac{1}{2}((u·\nabla )w,v)\forall u,v,w\in X.\hfill \end{array}$$

With the above notations, the general variational formulation of problem (3) is: get $(u,p)\in L^{\infty }(0,T;H)\cap L^2(0,T;X)\times L^2(0,T;M)$ satisfy:

$$\left\{\begin{array}{c}(u_t,v)+a(u,v)-d(v,p)+d(u,q)+b(u,u,v)=(f,v),t\in (0,T],\hfill \\ u\left(0\right)=u_0.\hfill \end{array}\right.$$

(6)

for all $(v,q)\in (X,M)$.

In order to get the fully dispersed error estimates, we need the following smoothness.

Theorem 1.

Assume some continuity of $u_0$ and $f\left(x\right)$ are valid [14,17]. Then problem (6) admits a unique solution $(u,p)$ satisfying the following estimates:

$$\begin{array}{c}\hfill {\parallel u\left(t\right)\parallel }^2+{\int }_0^t\left(\parallel u_t{\parallel }_0^2+{\parallel Au\parallel }_0^2+{\parallel p\parallel }_1^2\right)ds\le \kappa ,\\ \hfill \tau \left(t\right)\left({\parallel Au\left(t\right)\parallel }_0^2+{\parallel p\left(t\right)\parallel }_1^2+{\parallel u_t\left(t\right)\parallel }_0^2\right)+{\int }_0^t\tau \left(s\right){\parallel u_t\parallel }^{2}ds\le \kappa ,\\ \hfill {\tau }^2\left(t\right)\parallel u_t{\left(t\right)\parallel }^2+{\int }_0^t{\tau }^2\left(s\right)(\parallel u_{tt}{\parallel }_0^2+\parallel A{u}_t{\parallel }_0^2+\parallel p_t{\parallel }_1^2)ds\le \kappa ,\end{array}$$

(7)

for all $t\in [0,T]$, where κ is a positive constant.

Now, we consider the fully discrete locally stabilized FVM for two-dimensional time-dependent incompressible Navier-Stokes Equation (3). For convenience, let ${\tau }_h={\tau }_h\left(\mathsf{\Omega }\right)$ a partitioning of $\overline{\mathsf{\Omega }}$ into triangles satisfied the regular in the usual sense (see [5,14]). $(X_h,M_h)$ are the corresponding finite element subspace of $(X,M)$. ${\mathsf{\Gamma }}_h$ is the set of all interelement boundaries.

In order to define the finite volume method for Equation (3), We need introduce a popular configuration of dual partition $T_h^{*}$ for $T_h$: the interior point $p_i$ is chosen to the barycenter of element $K_i\in T_h$, and the midpoint $m_{ij}$ on side of $\overline{v_i{v}_j}$. See Figure 1. This type of dual partition is locally regular if $T_h$ is locally regular.

Corresponding to Hilbert space $V,H$, we define the following finite element velocity subspaces

$$X_h=\{v\in X:v_i{|}_K\in P_1\left(K\right),i=1,2\forall K\in {\tau }_h\},$$

and the pressure subspace

$$M_h=\{q\in M:q{|}_K\in P_0\left(K\right)\forall K\in {\tau }_h\}.$$

The finite volume dual of velocity is $X_h^{*}$

$$X_h^{*}=\{v\in L^2{\left(\mathsf{\Omega }\right)}^2:v_i{|}_{K^{*}}\in P_0\left(K^{*}\right),i=1,2\forall K^{*}\in T_h^{*}\}.$$

Let interpolation operator $I_h^{*}$:$X_h\to X_h^{*}$.

$$I_h^{*}{u}_h=\sum _{x_i\in N_h}{u}_h\left(x_i\right){\chi }_i\left(x\right).$$

The $L^2$ dimensional reduction $P_h:X\to X_h$ and $J_h:M\to M_h$ are defined as follows:

$$(P_{h}v,v_h)=(v,v_h)\forall v\in Y,v_h\in X_h,(J_{h}q,q_h)=(q,q_h)\forall q\in M,q_h\in M_h.$$

(8)

The finite volume forms of velocity $\tilde{a}(·,·)$ on $X_h\times X_h$ is,

$$\tilde{a}(u_h,I_h^{*}{v}_h)=\nu \left((u_h,I_h^{*}{v}_h)\right)=-\nu \sum _{K_i^{*}\in T_h^{*}}{\int }_{\partial K_i^{*}}{v}_h\left(x_i\right)\frac{\partial u_h}{\partial \mathbf{n}}ds\forall u_h,v_h\in X_h,$$

where $\mathbf{n}$ is the unit outnormal vector. The finite volume form $\tilde{d}(·,·)$ of pressure on $X_h\times M_h$ is defined as

$$\tilde{d}(I_h^{*}{v}_h,p_h)=-(I_h^{*}{v}_h,\nabla p_h)=\sum _{K_i^{*}\in T_h^{*}}{\int }_{\partial K_i^{*}}{p}_h{v}_h\left(x_i\right)·\mathbf{n}ds\forall u_h\in X_h,p_h\in M_h,$$

To facilitate the analysis, we need the following two trilinear forms.

$$\begin{array}{cc}\hfill \tilde{b}(u_h,v_h,I_h^{*}{w}_h)=((u_h·\nabla )v_h,I_h^{*}{w}_h)& +\frac{1}{2}(\left(\mathrm{div}{u}_h\right)v_h,I_h^{*}{w}_h),\hfill \\ \hfill \overline{b}(u_h,v_h,w_h-I_h^{*}{w}_h)=((u_h·\nabla )v_h,I_h^{*}{w}_h-I_h^{*}{w}_h)& +\frac{1}{2}(\left(\mathrm{div}{u}_h\right)v_h,w_h-I_h^{*}{w}_h)\hfill \\ \hfill & \forall u_h,v_h,w_h\in X_h,\hfill \end{array}$$

The last time difference part is

$$(u_{ht},I_h^{*}{v}_h)=\sum _{K_i^{*}\in T_h^{*}}{\int }_{K_i^{*}}{v}_h\left(x_i\right)u_{ht}dx\forall v_h\in X_h.$$

The finite volume form of the right side is

$$(f,I_h^{*}{v}_h)=\sum _{K_i^{*}\in T_h^{*}}{\int }_{K_i^{*}}{v}_h\left(x_i\right)fdx\forall v_h\in X_h.$$

For the convenience of reading, we introduce the following generalized form

$$\tilde{\mathcal{B}}((u_h,p_h);(I_h^{*}{v}_h,q_h))=\tilde{a}(u_h,I_h^{*}{v}_h)-\tilde{d}(I_h^{*}{v}_h,p_h)+d(u_h,q_h).$$

In this paper, the norms are defined as following:

$$\begin{array}{cc}\hfill \parallel u_h{\parallel }_{0,h,K}& =\begin{array}{c}{\left[{\displaystyle \frac{S_Q}{3}(u_{P_i}^2+u_{P_j}^2+u_{P_k}^2)}\right]}^{1/2},\hfill \end{array}\hfill \\ \hfill \parallel {\tilde{A}}_h^{1/2}{u}_h{\parallel }_{0,h,K}& =\begin{array}{c}{\left\{\left[{\left(\frac{\partial u_h\left(Q\right)}{\partial x}\right)}^2+{\left(\frac{\partial u_h\left(Q\right)}{\partial y}\right)}^2\right]S_Q\right\}}^{1/2},\hfill \end{array}\hfill \end{array}$$

$$\begin{array}{cc}& \parallel u_h{\parallel }_{0,h}={\left({\displaystyle \sum _{K\in T_h}{\parallel u_h\parallel }_{0,h,K}^2}\right)}^{1/2},\hfill \\ & \parallel {\tilde{A}}_h^{1/2}{u}_h{\parallel }_0={\left({\displaystyle \sum _{K\in T_h}{\parallel {\tilde{A}}_h^{1/2}{u}_h\parallel }_{0,h,K}^2}\right)}^{1/2},\hfill \\ & \parallel u_h{\parallel }_{1,h}={\left(\parallel u_h{\parallel }_{0,h}^2+{\parallel {\tilde{A}}_h^{1/2}{u}_h\parallel }_0^2\right)}^{1/2},\hfill \end{array}$$

where $S_Q$ is the area of $▵P_i{P}_j{P}_k$ (see Figure 2).

To describe the locally stabilized formulation of the non-stationary Navier-Stokes problem, we use the classic not overlap macroelement partitioning ${\Lambda }_h$ [18]. For every macroelement $\mathcal{K}$ in ${\Lambda }_h$, the set of interelement(small finite element) edges is denoted by ${\mathsf{\Gamma }}_{\mathcal{K}}$, and the length of an edge $e\in {\mathsf{\Gamma }}_{\mathcal{K}}$ is denoted by $h_e$.

With the above definitions, a locally stabilized formulation of the non-stationary Navier-Stokes problem (3) can be stated as follows.

Definition 1.

Locally stabilized finite volume formulation for non-stationary Navier-Stokes: Find $(u_h,p_h)\in (X_h,M_h),$ such that for all $(v,q)\in (X_h,M_h)$

$$\left\{\begin{array}{c}(u_{ht},I_h^{*}v)+{\tilde{\mathcal{B}}}_h((u_h,p_h);(I_h^{*}v,q))+\tilde{b}(u_h,u_h,I_h^{*}{w}_h)=(f,I_h^{*}v),\hfill \\ u_h\left(0\right)=u_{0h}.\hfill \end{array}\right.$$

(9)

where

$${\tilde{\mathcal{B}}}_h((u_h,p_h);(I_h^{*}v,q))=\tilde{\mathcal{B}}((u_h,p_h);(I_h^{*}v,q))+\beta {\mathcal{C}}_h(p_h,q)\forall (u_h,p_h),(v,q)\in (X_h,M_h),$$

$${\mathcal{C}}_h(p,q)=\sum _{\mathcal{K}\in {\Lambda }_h}\sum _{e\in {\mathsf{\Gamma }}_{\mathcal{K}}}{h}_e\underset{e}{\int }{\left[p\right]}_e{\left[q\right]}_{e}ds,$$

for all $p,q$ in the algebraic sum $H^1\left(\mathsf{\Omega }\right)+M_h$, and ${[·]}_e$ is the jump operator across $e\in {\mathsf{\Gamma }}_{\mathcal{K}}$ and $\beta >0$ is the local stabilization parameter.

In order get the regularity of the above definitions, we need the following stability results [5,8,9,12].

Theorem 2.

For any two neighboring macroelements ${\mathcal{K}}_1$ and ${\mathcal{K}}_2$ with ${\int }_{{\mathcal{K}}_1\cap {\mathcal{K}}_2}ds\ne 0$, if there exists $v\in X_h$ such that

$$suppv\subset {\mathcal{K}}_1\cup {\mathcal{K}}_2\mathrm{and}{\int }_{{\mathcal{K}}_1\cap {\mathcal{K}}_2}v·nds\ne 0.$$

(10)

Then

$$\begin{array}{cc}\hfill |{\tilde{\mathcal{B}}}_h((u,p);(I_h^{*}v,q))|& \le c(\parallel u\parallel +{\parallel p\parallel }_0\left)\right(\parallel v\parallel +{\parallel q\parallel }_0)\forall (u,p),(v,q)\in (X,M),\hfill \\ \hfill \alpha (\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0+\parallel p_h{\parallel }_0)& \le \underset{(v_h,q_h)\in (X_h,M_h)}{sup}\frac{{\tilde{\mathcal{B}}}_h((u_h,p_h);(I_h^{*}{v}_h,q_h))}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0+{\parallel q_h\parallel }_0}\hfill \end{array}$$

(11)

for all $(u_h,p_h)\in (X_h,M_h)$, and

$${\mathcal{C}}_h(p,q_h)=0,{\mathcal{C}}_h(p_h,q)=0,{\mathcal{C}}_h(p,q)=0\forall p,q\in H^1\left(\mathsf{\Omega }\right),p_h,q_h\in M_h,$$

(12)

where $\beta \ge {\beta }_0>0$ and $\alpha >0$ are two constant.

3. Technical Preliminaries

The main task of this section is to prepare many basic estimates which will help the error analyses for the finite volume solution $(u_h,p_h)$.

Since the bilinear forms $\left((u_h,v_h)\right)$ and $\left((u_h,I_h^{*}{v}_h)\right)$ are coercive on $X_h\times X_h$, they generate invertible operators $A_h:X_h\to X_h$ and ${\tilde{A}}_h:X_h\to X_h$ respectively through the condition:

$$\begin{array}{c}\hfill (A_h{u}_h,v_h)=(A_h^{1/2}{u}_h,A_h^{1/2}{v}_h)=\left((u_h,v_h)\right)\forall u_h,v_h\in X_h,\\ \hfill ({\tilde{A}}_h{u}_h,I_h^{*}{v}_h)=\left((u_h,I_h^{*}{v}_h)\right)\forall u_h,v_h\in X_h.\end{array}$$

Moreover, we also need the discrete gradient operators:

$$(v_h,{\nabla }_h{q}_h)=-d(v_h,q_h)\forall (v_h,q_h)\in (X_h,M_h),$$

Firstly, we have the following classical properties (see [8,19])

$$\begin{array}{cc}\hfill \parallel v_h\parallel \le c{h}^{-1}\parallel v_h{\parallel }_0,{\parallel v_h\parallel }_{\infty }& \le {c|lnh|}^{1/2}\parallel v_h\parallel \forall v_h\in X_h,\hfill \\ \hfill \parallel P_{h}v\parallel & \le c\parallel v\parallel \forall v\in X.\hfill \\ \hfill \parallel v-I_h{v\parallel }_0+h\parallel v-I_{h}v\parallel & \le c{h}^2{\parallel v\parallel }_2\forall v\in D\left(A\right),\hfill \\ \hfill \parallel v-I_h^{*}{v\parallel }_0& \le ch\parallel v\parallel \forall v\in X,\hfill \end{array}$$

(13)

where ${\parallel v\parallel }_{\infty }={\parallel v\parallel }_{L^{\infty }}={\mathrm{ess}.\sup }_{x\in \mathsf{\Omega }}{\parallel v\left(x\right)\parallel }_0$. For the $P_1-P_0$ triangular element, It follows from (4), (5) and (13) that [20,21]

$$\begin{array}{c}\hfill \parallel v_h{\parallel }_0\le {\gamma }_0\parallel v_h\parallel ,\parallel v_h\parallel \le {\gamma }_0\parallel A_h{v}_h{\parallel }_0,\parallel A_h{v}_h{\parallel }_0\le c{h}^{-1}\parallel v_h\parallel \forall v_h\in X_h,\\ \hfill \parallel v_h{\parallel }_{0,h}\le {\gamma }_0^{\prime }\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0,\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0\le {\gamma }_0^{\prime }{\parallel {\tilde{A}}_h{v}_h\parallel }_0\forall v_h\in X_h,\\ \hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0\le c{h}^{-1}\parallel v_h{\parallel }_{0,h},\parallel {\tilde{A}}_h{v}_h{\parallel }_0\le c{h}^{-1}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h\parallel }_0\forall v_h\in X_h,\\ \hfill c_1\parallel v_h{\parallel }_0\le \parallel v_h{\parallel }_{0,h}\le c_1^{\prime }\parallel v_h{\parallel }_0,c_1\parallel A_h{v}_h{\parallel }_0\le \parallel {\tilde{A}}_h{v}_h{\parallel }_0\le c_1^{\prime }{\parallel A_h{v}_h\parallel }_0\forall v_h\in X_h.\end{array}$$

(14)

As for the trilinear forms $\tilde{b}$ and $\overline{b}$,we can deduce the following results.

Lemma 1.

If $u_h,v_h,w_h\in X_h$, we have

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& + |\overline{b}(u_h,v_h,w_h-I_h^{*}{w}_h)|\le c_2\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0{\parallel {\tilde{A}}_h^{\frac{1}{2}}{w}_h\parallel }_0,\hfill \\ \hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& \le c_2\parallel u_h{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{w}_h{\parallel }_0^{1/2}{\parallel {\tilde{A}}_h{w}_h\parallel }_0^{\frac{1}{2}},\hfill \\ \hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& \le c_2\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0^{1/2}\parallel {\tilde{A}}_h{u}_h{\parallel }_0^{\frac{1}{2}}\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0{\parallel w_h\parallel }_0.\hfill \\ \hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& \le c_2{|lnh|}^{\frac{1}{2}}\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0{\parallel w_h\parallel }_0,\hfill \end{array}\end{array}$$

(15)

Proof. 

Since the following discrete analogue of the Sobolev inequality holds [17]

$$\parallel {\varphi }_h{\parallel }_{L^4}\le c\parallel {\varphi }_h{\parallel }_0^{1/2}{\parallel A_h^{1/2}{\varphi }_h\parallel }_0^{1/2}\forall {\varphi }_h\in X_h.$$

(16)

If $u_h,v_h,w_h\in X_h,$, we have

$$\left|\tilde{b}(u_h,v_h,I_h^{*}{w}_h)\right|\le c\parallel u_h{\parallel }_{L^4}\parallel A_h^{1/2}{v}_h{\parallel }_0{\parallel w_h\parallel }_{L^4},$$

combining the above formula with (14) and (16), we can get the first formula in (15).

To prove the others in (15), we need the follow discrete results [17,22], namely for any $h>0$,

$$\begin{array}{cc}\hfill \parallel {\varphi }_h{\parallel }_{L^{\infty }}& \le c_3{|lnh|}^{1/2}\parallel v_h\parallel ,\hfill \\ \hfill \parallel {\varphi }_h{\parallel }_{L^6}& \le c_3\parallel {\varphi }_h\parallel \forall {\varphi }_h\in X_h,\hfill \\ \hfill \parallel {\varphi }_h{\parallel }_{L^{\infty }}+{\parallel \nabla {\varphi }_h\parallel }_{L^3}& \le c_3\parallel {\varphi }_h{\parallel }^{1/2}{\parallel A_h{\varphi }_h\parallel }_0^{1/2}\forall {\varphi }_h\in X_h.\hfill \end{array}$$

(17)

For any $u_h,v_h,w_h\in X_h$, we have [9]

$$\begin{array}{ccc}\hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& \le \hfill & \hfill c_4\parallel u_h{\parallel }_{L^{\infty }}\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0\parallel w_h{\parallel }_0+c_4\parallel \nabla u_h{\parallel }_{L^3}\parallel v_h{\parallel }_{L^6}{\parallel w_h\parallel }_0,\\ \hfill |\tilde{b}(v_h,u_h,I_h^{*}{w}_h)|& \le \hfill & \hfill c_4\parallel v_h{\parallel }_{L^6}\parallel \nabla u_h{\parallel }_{L^3}\parallel w_h{\parallel }_0+c_4\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0\parallel u_h{\parallel }_{L^{\infty }}{\parallel w_h\parallel }_0,\\ \hfill |\tilde{b}(u_h,v_h,I_h^{*}{w}_h)|& \le \hfill & \hfill c_4\parallel u_h{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h\parallel \parallel w_h{\parallel }_{L^{\infty }}+c_4\parallel u_h{\parallel }_0\parallel \nabla v_h{\parallel }_{L^3}{\parallel w_h\parallel }_{L^6},\end{array}$$

which together with (14), (17) imply (15).  □

Similar to the results in [12], we also need to define the projection operator $({\tilde{R}}_h,{\tilde{Q}}_h):(X,Y)\to (X_h,M_h)$ as

$$\begin{array}{c}\hfill {\tilde{\mathcal{B}}}_h(({\tilde{R}}_h(u,p),{\tilde{Q}}_h(u,p));(I_h^{*}{v}_h,q_h))={\tilde{\mathcal{B}}}_h((u,p);(I_h^{*}{v}_h,q_h))\\ \hfill \forall (u,p)\in (X,M),(v_h,q_h)\in (X_h,M_h).\end{array}$$

Due to Theorem 2, we know that $({\tilde{R}}_h,{\tilde{Q}}_h)$ is well defined and have the properties [12]:

$$\parallel {\tilde{A}}_h^{\frac{1}{2}}({\tilde{R}}_h(u,p)-u){\parallel }_0+\parallel {\tilde{Q}}_h{(u,p)-p\parallel }_0\le {ch(\parallel u\parallel }_2+{\parallel p\parallel }_1),$$

(18)

for all $(u,p)\in (H^2{\left(\mathsf{\Omega }\right)}^2\cap X,H^1\left(\mathsf{\Omega }\right)\cap M)$.

Beside, we need the specific result in He et al. [12].

Theorem 3.

Under the assumptions of Theorems 1 and 2, $(u_h,p_h)$ satisfies

$${\int }_0^t\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h){\parallel }_0^{2}ds+\tau \left(t\right)\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t\right)-u_h\left(t\right)){\parallel }_0^2+{\tau }^2\left(t\right){\parallel p\left(t\right)-p_h\left(t\right)\parallel }_0^2\le \kappa h^2,$$

(19)

for all $t\in [0,T]$.

Since our error analysis for the time discretization depends heavily on these regularity estimates, we then provide some smoothness estimates of $(u_h,p_h)$. The main idea is similar to the work in He et al. [9,23].

Theorem 4.

Under the assumptions of Theorems 3, the finite volume solution $(u_h,p_h)$ satisfies

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\left(t\right)\parallel }_0^2& +{\int }_0^t(\parallel {\tilde{A}}_h^{-1}(u_{htt}+{\nabla }_h{p}_{ht}){\parallel }_0^2+\parallel u_{ht}{\parallel }_0^2+\parallel {\tilde{A}}_h{u}_h{\parallel }_0^2)ds\le \kappa ,\hfill \\ \hfill {\tau }^2\left(t\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_{ht}\left(t\right)\parallel }^2& +{\int }_0^{\tau }{\tau }^2\left(s\right)(\parallel u_{htt}+{\nabla }_h{p}_{ht}{\parallel }_0^2+\parallel {\tilde{A}}_h{u}_{ht}{\parallel }_0^2)ds\le \kappa ,\hfill \\ \hfill {\tau }^2\left(t\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_{ht}\left(t\right)\parallel }^2& +{\int }_0^{\tau }{\tau }^2\left(s\right)(\parallel u_{htt}+{\nabla }_h{p}_{ht}{\parallel }_0^2+\parallel {\tilde{A}}_h{u}_{ht}{\parallel }_0^2)ds\le \kappa ,\hfill \end{array}\end{array}$$

(20)

for all $t\in [0,T]$.

Proof. 

Note from (13) and (14) we can deduce the following estimates:

$$\begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\left(t\right)\parallel }_0& \le \parallel {\tilde{A}}_h^{\frac{1}{2}}(u_h\left(t\right)-A_h^{1/2}{P}_{h}u\left(t\right)){\parallel }_0+c\parallel u\left(t\right)\parallel \hfill \\ \hfill & \le c{h}^{-1}\parallel u_h{\left(t\right)-u\left(t\right)\parallel }_0+c\parallel u\left(t\right)\parallel ,\hfill \\ \hfill \parallel {\tilde{A}}_h{u}_h{\parallel }_0& =\underset{v_h\in X_h}{sup}\frac{\left|\right({\tilde{A}}_h{u}_h,I_h^{*}{v}_h\left)\right|}{\parallel I_h^{*}{v}_h{\parallel }_0}\le (c{h}^{-1}\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h){\parallel }_0+{\parallel Au\parallel }_0).\hfill \end{array}$$

and from Theorems 1 and 3, we have

$$\begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\left(t\right)\parallel }_0^2& +{\int }_0^t{\parallel {\tilde{A}}_h{u}_h\parallel }_0^{2}ds\le c\left({\parallel u\left(t\right)\parallel }^2+{\int }_0^t{\parallel Au\parallel }_0^{2}ds\right)\hfill \\ \hfill & +c{h}^{-2}\left(\parallel u\left(t\right)-u_h{\left(t\right)\parallel }_{0,h}^2+{\int }_0^t{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h)\parallel }_0^{2}ds\right)\le \kappa ,\hfill \\ \hfill \tau \left(t\right)\parallel {\tilde{A}}_h{u}_h{\left(t\right)\parallel }_0^2& \le c\tau \left(t\right)(h^{-2}\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t\right)-u_h\left(t\right)){\parallel }_0^2+\parallel Au\left(t\right){\parallel }_0^2)\le \kappa ,\hfill \end{array}$$

for all $t\in [0,T]$. Then using the similar method in [9], we can get these estimates (20).  □

Finally, in order to get the upper bounders of velocity and pressure in the time related case, we state the classical Gronwall lemma used in [24].

Lemma 2.

Let $C_0$ and $a_k,b_k,c_k,d_k$, for integers $k\ge 0$, be nonnegative numbers such that

$$a_n+\Delta t\sum _{k=0}^n{b}_k\le \Delta t\sum _{k=0}^{n-1}{d}_k{a}_k+\Delta t\sum _{k=0}^{n-1}{c}_k+C_0\forall n\ge 1.$$

(21)

Then,

$$a_n+\Delta t\sum _{k=0}^n{b}_k\le exp\left(\Delta t\sum _{k=0}^{n-1}{d}_k\right)\left(\Delta t\sum _{k=0}^{n-1}{c}_k+C_0\right)\forall n\ge 1.$$

(22)

The following is dual Gronwall lemma.

Lemma 3.

Given integer $m>0$ and let C and $a_k,b_k,c_k,d_k$, for integers $0\le k\le m$, be nonnegative numbers such that

$$a_n+\Delta t\sum _{k=n+1}^m{b}_k\le \Delta t\sum _{k=n+1}^m{d}_k{a}_k+\Delta t\sum _{k=n+1}^m{c}_k+C,0\le n\le m.$$

(23)

Then,

$$a_n+\Delta t\sum _{k=n+1}^m{b}_k\le exp\left(\Delta t\sum _{k=n+1}^m{d}_k\right)\left(\Delta t\sum _{k=n+1}^m{c}_k+C\right),0\le n\le m.$$

(24)

4. Error Estimates for Semi-Discrete for Time Depended Navier-Stokes Equations

In this section we consider the time discretization of the locally stabilized and get some useful estimates. Let T time to stop calculation and N the time corresponding step. So we have

$$\Delta t=\frac{T}{N},t_n=n\Delta t,n=0,1,\cdots ,N.$$

For the first part, we need to analysis the errors of finite element Original case. It’s well known that the common Euler semi-implicit scheme applied to the spatially discrete problem (9) can be described as:

$$\begin{array}{c}\hfill (d_t{u}_h^n,I_h^{*}{v}_h)+{\tilde{\mathcal{B}}}_h((u_h^n,p_h^n);(I_h^{*}{v}_h,q_h))+\tilde{b}(u_h^{n-1},u_h^n,I_h^{*}{v}_h)=(f^n,I_h^{*}{v}_h),\end{array}$$

(25)

for all $(v_h,q_h)\in (X_h,M_h)$, where $u_h^0=u_{0h}$ is starting value and

$$d_t{u}_h^n=\frac{1}{\Delta t}(u_h^n-u_h^{n-1}),f^n=\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}f\left(t\right)dt,$$

$${\overline{p}}_h\left(t_n\right)=\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}{p}_h\left(t\right)dt.$$

To deduce the discretization error $(e_h^n,{\mu }_h^n)=(u_h\left(t_n\right)-u_h^n,{\overline{p}}_h\left(t_n\right)-p_h^n)$, we integrate and differentiate (9), respectively, to get

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \frac{1}{\Delta t}(u_h\left(t_n\right)-u_h\left(t_{n-1}\right),I_h^{*}{v}_h)& +\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}{\tilde{\mathcal{B}}}_h((u_h\left(t\right),p_h\left(t\right));(I_h^{*}{v}_h,q_h)))dt\hfill \\ & +\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}\tilde{b}(u_h\left(t\right),u_h\left(t\right),I_h^{*}{v}_h)dt=(f^n,I_h^{*}{v}_h),\hfill \end{array}\end{array}$$

(26)

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill (u_{htt}+{\nabla }_h{p}_{ht},I_h^{*}{v}_h)+\tilde{a}(u_{ht},I_h^{*}{v}_h)& -\tilde{d}(I_h^{*}{v}_h,u_{ht})+d(u_{ht},q_h){+}_h(p_{ht},q_h)\hfill \\ & +\tilde{b}(u_{ht},u_h,I_h^{*}{v}_h)+\tilde{b}(u_h,u_{ht},I_h^{*}{v}_h)=(f_t,I_h^{*}{v}_h),\hfill \end{array}\end{array}$$

(27)

for all $(v_h,q_h)\in (X_h,M_h)$.

Subtracting (25) from (26) and using (27) and the relation:

$$\begin{array}{c}\hfill \varphi \left(t_n\right)-\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}\varphi \left(t\right)dt=\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1}){\varphi }_t\left(t\right)dt\end{array}$$

for all $\varphi \in H^1(t_{n-1},t_n;F)$ for some Hilbert space F, we have

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill (d_t{e}_h^n,v_h)& +{\tilde{\mathcal{B}}}_h((e_h^n,{\mu }_h^n);(I_h^{*}{v}_h,q_h))+\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{v}_h)\hfill \\ & +\tilde{b}(u_h^{n-1},e_h^n,I_h^{*}{v}_h)=(E_n,I_h^{*}{v}_h),\hfill \end{array}\end{array}$$

(28)

for all $(v_h,q_h)\in (X_h,M_h)$, where $(e_h^0,{\mu }_h^0)=(0,0)$ and

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill (E_n,I_h^{*}{v}_h)=& -\frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1})(u_{htt}+{\nabla }_h{p}_{ht},I_h^{*}{v}_h)dt\hfill \\ & +\tilde{b}\left({\int }_{t_{n-1}}^{t_n}{u}_{ht}dt,u_h\left(t_n\right),I_h^{*}{v}_h\right).\hfill \end{array}\end{array}$$

(29)

Lemma 4.

Under the assumptions of Theorem 3 the error $E_n$ satisfies the following bounds:

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Delta t\sum _{n=1}^m{\parallel {\tilde{A}}_h^{-1}{E}_n\parallel }_0^2& \le \kappa \Delta t^2,1\le m\le N,\hfill \\ \hfill \Delta t\sum _{n=1}^m{\tau }^i\left(t_n\right){\parallel {\tilde{A}}_h^{-1/2}{E}_n\parallel }_0^2& \le \kappa \Delta t^{i+1},1\le m\le N,i=0,1,\hfill \\ \hfill \Delta t\sum _{n=1}^m{\tau }^i\left(t_n\right){\parallel E_n\parallel }_{0,h}^2& \le \kappa \Delta t^i,1\le m\le N,i=0,1,2.\hfill \end{array}\end{array}$$

(30)

Proof. 

By using (14), (15) and (29), we derive

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{-1}{E}_n{\parallel }_0& \le \frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1})\parallel {\tilde{A}}_h^{-1}(u_{htt}+{\nabla }_h{p}_{ht}){\parallel }_{0}dt+c\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h\left(t_n\right){\parallel }_0{\int }_{t_{n-1}}^{t_n}{\parallel u_{ht}\parallel }_{0}dt\hfill \\ & \le c\Delta t^{1/2}{\left({\int }_{t_{n-1}}^{t_n}\left(\parallel {\tilde{A}}_h^{-1}(u_{htt}+{\nabla }_h{p}_{ht}){\parallel }_0^2+\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h\left(t_n\right){\parallel }_0^2{\parallel u_{ht}\parallel }_0^2\right)dt\right)}^{1/2}.\hfill \end{array}\end{array}$$

(31)

Applying Theorem 4 in (39), we obtain

$$\begin{array}{cc}\hfill \Delta t\parallel {\tilde{A}}_h^{-1}{E}_n{\parallel }_0^2& \le c\Delta t^2{\int }_{t_{n-1}}^{t_n}\left(\parallel {\tilde{A}}_h^{-1}(u_{htt}+{\nabla }_h{p}_{ht}){\parallel }_0^2+{\parallel u_{ht}\parallel }_0^2\right)dt.\end{array}$$

(32)

Utilizing Theorem 4, if we sum (32) from $n=1$ to $n=m$, we can derive the first inequality in (30) directly.

Next, we deduce from (19) and (29) that

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{-1/2}{E}_n{\parallel }_0& \le \frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1}){\parallel {\tilde{A}}_h^{-1/2}(u_{htt}+{\nabla }_h{p}_{ht})\parallel }_{0}dt\hfill \\ & +c{\int }_{t_{n-1}}^{t_n}\parallel {\tilde{A}}_h{u}_h\left(t_n\right){\parallel }_0{\parallel u_{htt}\parallel }_{0}dt\hfill \\ & \le c({\int }_{t_{n-1}}^{t_n}\left((t-t_{n-1}){\parallel {\tilde{A}}_h^{-1/2}(u_{htt}+{\nabla }_h{p}_{ht})\parallel }_0^2\right.\hfill \\ & +\Delta t\left.\parallel {\tilde{A}}_h{u}_h\left(t_n\right){\parallel }_0^2{\parallel u_{ht}\parallel }_0^2\right)dt{)}^{1/2}.\hfill \end{array}\end{array}$$

(33)

Because

$$\begin{array}{c}\hfill \tau \left(t_n\right)\le \tau \left(t\right)+\Delta t,\Delta t\le \tau \left(t_n\right),t-t_{n-1}\le \tau \left(t\right),\forall t\in [t_{n-1},t_n],\end{array}$$

(34)

We can deduce from (33) and Theorem 4 to get

$$\begin{array}{c}\hfill {\tau }^i\left(t_n\right)\parallel {\tilde{A}}_h^{-1/2}{E}_n{\parallel }_0^2\Delta t\le c\Delta t^{1+i}{\int }_{t_{n-1}}^{t_n}(\tau \left(t\right)\parallel {\tilde{A}}_h^{-1/2}(u_{htt}+{\nabla }_h{p}_{ht}){\parallel }_0^2+\parallel u_{ht}{\parallel }_0^2)dt,\end{array}$$

(35)

for $i=0,1$. Summing (35) from $n=1$ to $n=m$, we derive the second inequality in (30).

As for the last one, Deriving from (14) and (29), we have

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \parallel E_n{\parallel }_{0,h}& \le \frac{1}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1}){\parallel u_{htt}+{\nabla }_h{p}_{ht}\parallel }_{0}dt\hfill \\ & +c\parallel {\tilde{A}}_h{u}_h\left(t_n\right){\parallel }_0{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u_h\left(t_n\right)-u_h\left(t_{n-1}\right))\parallel }_0.\hfill \end{array}\end{array}$$

(36)

Hence, Formula (36) and Theorem 4 imply that

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\tau }^i\left(t_n\right){\parallel E_n\parallel }_{0,h}^2\Delta t& \le c{\tau }^i\left(t_n\right){\int }_{t_{n-1}}^{t_n}{(t-t_{n-1})}^2{\parallel u_{htt}+{\nabla }_h{p}_{ht}\parallel }_0^{2}dt\hfill \\ & +c{\tau }^i\left(t_n\right)\parallel {\tilde{A}}_h{u}_h\left(t_n\right){\parallel }_0^2{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u_h\left(t_n\right)-u_h\left(t_{n-1}\right))\parallel }_0^2\Delta t.\hfill \end{array}\end{array}$$

(37)

Similarly, summing (37) from $n=1$ to $n=m$ and using Theorem 4 deduces

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Delta t\sum _{n=1}^m{\tau }^i\left(t_n\right){\parallel E_n\parallel }_{0,h}^2& \le c\sum _{n=1}^m{\tau }^i\left(t_n\right){\int }_{t_{n-1}}^{t_n}{(t-t_{n-1})}^2{\parallel u_{htt}+{\nabla }_h{p}_{ht}\parallel }_0^{2}dt\hfill \\ \hfill & +c{\tau }^i\left(t_1\right)\parallel {\tilde{A}}_h{u}_h\left(t_1\right){\parallel }_0^2{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u_h\left(t_1\right)-u_h\left(t_0\right))\parallel }_0^2\Delta t\hfill \\ & +c\Delta t^2\sum _{n=2}^m{\tau }^i\left(t_n\right)\parallel {\tilde{A}}_h{u}_h\left(t_n\right){\parallel }_0^2{\int }_{t_{n-1}}^{t_n}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_{ht}\parallel }_0^{2}dt\hfill \\ & \le c\Delta t^i\sum _{n=1}^m{\int }_{t_{n-1}}^{t_n}{\tau }^2\left(t\right){\parallel u_{htt}+{\nabla }_h{p}_{ht}\parallel }_0^{2}dt\hfill \\ & +c\parallel {\tilde{A}}_h^{\frac{1}{2}}(u_h\left(t_1\right)-u_h\left(t_0\right)){\parallel }_0^2\Delta t^i+c\Delta t^i\sum _{n=2}^m{\int }_{t_{n-1}}^{t_n}\tau \left(t\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_{ht}\parallel }_0^{2}dt\hfill \\ & \le \kappa \Delta t^i,\hfill \end{array}\end{array}$$

for i = 0,1,2, which yields the third one in (30).  □

Now, let’s discuss the second part: the error of time discrete duality argument corresponding to (25). Firstly, the dual problem corresponding to (25) usually describes as: find $({\mathsf{\Phi }}_h^{n-1},{\mathsf{\Psi }}_h^{n-1})\in (X_h,M_h)$ such that, for all $(v_h,q_h)\in (X_h,M_h),$

$$\begin{array}{cc}\hfill (v_h,I_h^{*}\left(d_t{\mathsf{\Phi }}_h^n\right))& -{\tilde{\mathcal{B}}}_h((v_h,q_h);(I_h^{*}{\mathsf{\Phi }}_h^{n-1},{\mathsf{\Psi }}_h^{n-1}))\hfill \\ \hfill & -\tilde{b}(v_h,u_h\left(t_n\right),I_h^{*}{\mathsf{\Phi }}_h^{n-1})-\tilde{b}(u_h\left(t_{n-1}\right),v_h,I_h^{*}{\mathsf{\Phi }}_h^{n-1})=(v_h,I_h^{*}{z}^n),\hfill \end{array}$$

(38)

where ${\mathsf{\Phi }}_h^m=0$.

For the dual part, we also need the existence, uniqueness and regularity of problem (38), so we introduce the following results:

$$\begin{array}{cc}\hfill \parallel P_{h}B(u_h\left(t_n\right),·){\parallel }_{*}& =\underset{v_h\in X_h}{sup}\frac{\parallel P_{h}B(u_h\left(t_n\right),v_h){\parallel }_0}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0},\hfill \\ \hfill \parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}& =\underset{v_h\in X_h}{sup}\frac{\parallel P_{h}B(v_h,u_h\left(t_n\right)){\parallel }_0}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{v}_h{\parallel }_0}.\hfill \end{array}$$

With the similar method in He [9], we have the following two lemmas.

Lemma 5.

Assume that the assumptions of Theorem 3 is valid and $\Delta t$ satisfies

$$\frac{8}{\nu }{c}_2^2|lnh|\underset{t\in [0,T]}{max}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h\left(t\right)\parallel }_0^2\Delta t\le 1.$$

(39)

Then,

$$\Delta t\sum _{n=1}^N\left(\parallel P_{h}B(u_h\left(t_{n-1}\right),·){\parallel }_{*}^2+\parallel P_{h}B{(·,u_h\left(t_n\right)\parallel }_{*}^2\right)\le \kappa .$$

(40)

Lemma 6.

Assume that the assumptions of Theorem 3 are valid and that $\Delta t$ satisfies (39). Then, problem (38) admits a unique solution $({\mathsf{\Phi }}_h^{n-1},{\mathsf{\Psi }}_h^{n-1})\in (X_h,M_h)$ for a given ${\mathsf{\Phi }}_h^n$. Furthermore, the solution sequence ${\{{\mathsf{\Phi }}_h^n,{\mathsf{\Psi }}_h^n\}}_{n=0}^{m-1}$ of problem (38) satisfies the following bound

$$\underset{0\le r\le m}{sup}\parallel {\tilde{A}}_h^{\frac{1}{2}}{\mathsf{\Phi }}_h^r{\parallel }_0^2+\Delta t\sum _{n=1}^{m-1}\parallel d_t{\mathsf{\Phi }}_h^n{\parallel }_0^2+\Delta t\sum _{n=0}^{m-1}\parallel {\tilde{A}}_h{\mathsf{\Phi }}_h^n{\parallel }_0^2\le \kappa \Delta t\sum _{n=1}^m{\parallel z^n\parallel }_0^2.$$

(41)

5. Error Analysis for Time and Spacial Discrete

In this section we proof the upper bounds for the error $(e_h^n,{\mu }_h^n)=(u_h\left(t_n\right)-u_h^n,{\overline{p}}_h\left(t_n\right)-p_h^n)$ in $L_2$ and $H_1$ norms, and deduce the last optimal order estimates. Firstly, we have

Lemma 7.

Assume that the assumptions of Theorem 3 are valid and $\Delta t$ satisfies (39). Then, the error $(e_h^n,{\mu }_h^n),1\le n\le N,$ satisfies the following bound:

$$\parallel e_h^m{\parallel }_{0,h}^2+\Delta t\sum _{n=1}^m\left(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)\right)\le \kappa \Delta t,1\le m\le N.$$

(42)

Proof. 

Setting $v_h=e_h^n,q_h={\mu }_h^n$ in (28) and using (14) we obtain

$$\begin{array}{cc}\hfill \frac{c_1}{2\Delta t}(\parallel e_h^n{\parallel }_0^2& +\parallel e_h^n-e_h^{n-1}{\parallel }_0^2-\parallel e_h^{n-1}{\parallel }_0^2)+\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)\hfill \\ \hfill & +\tilde{b}(e_h^{n-1},u\left(t_n\right),I_h^{*}{e}_h^n)+\tilde{b}(u\left(t_n\right),e_h^n,I_h^{*}{e}_h^n)\hfill \\ \hfill & =(E_n,e_h^n)\le \frac{\nu }{4}\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+{\nu }^{-1}{\parallel A_h^{-1/2}{E}_n\parallel }_0^2.\hfill \end{array}$$

(43)

From (14) and (15), we deduce that

$$\begin{array}{cc}\hfill |\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{e}_h^n)|& \le c\parallel e_h^{n-1}{\parallel }_0\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n\parallel }_0\hfill \\ \hfill & \le \frac{\nu }{8}\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\kappa \parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2{\parallel e_h^{n-1}\parallel }_0^2,\hfill \\ \hfill |\tilde{b}(u\left(t_n\right),e_h^n,I_h^{*}{e}_h^n)|& \le c\parallel e_h^n{\parallel }_0\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n\parallel }_0\hfill \\ \hfill & \le \frac{\nu }{8}\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\kappa \parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2{\parallel e_h^n\parallel }_0^2.\hfill \end{array}$$

Combining (43) with the above estimate yields

$$\begin{array}{cc}\hfill \parallel e_h^n{\parallel }_0^2& -\parallel e_h^{n-1}{\parallel }_0^2+(\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n))\Delta t\hfill \\ \hfill & \le \kappa \left(\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2(\parallel e_h^{n-1}{\parallel }_{0,h}^2+\parallel e_h^n{\parallel }_{0,h}^2)+\parallel {\tilde{A}}_h^{-1/2}{E}_n{\parallel }_0^2\right)\Delta t.\hfill \end{array}$$

(44)

Summing (44) from $n=1$ to $n=m$, we obtain

$$\begin{array}{cc}\hfill \parallel e^m{\parallel }_{0,h}^2& +\Delta t\sum _{n=1}^m(\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n))\hfill \\ \hfill & \le \kappa \Delta t\left(\sum _{n=1}^m\parallel P_{h}B(·,u_h\left(t_{n+1}\right)){\parallel }_{*}^2(\parallel e_h^{n-1}{\parallel }_{0,h}^2+\parallel e_h^n{\parallel }_{0,h}^2)+\sum _{n=1}^m{\parallel {\tilde{A}}_h^{-1/2}{E}_n\parallel }_0^2\right).\hfill \end{array}$$

(45)

Applying Lemmas 3 and 6 in (45) and by (14) yields (43).  □

Lemma 8.

Under the assumptions of Lemma 7, we have

$$\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^m{\parallel }^2+\Delta t\sum _{n=1}^m{\parallel d_t{e}_h^n\parallel }_{0,h}^2\le \kappa ,1\le m\le N.$$

(46)

Proof. 

We derive from (28) that

$$\begin{array}{cc}\hfill (d_t{e}_h^n,I_h^{*}{v}_h)& +\tilde{a}(e_h^n,I_h^{*}{v}_h)-\tilde{d}(I_h^{*}{v}_h,{\mu }_h^n)+d(d_t{e}_h^n,q_h)+\beta {\mathcal{C}}_h(d_t{\mu }_h^n,q_h)\hfill \\ \hfill & +\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{v}_h)+\tilde{b}(u_h^{n-1},e_h^n,I_h^{*}{v}_h)=(E_n,I_h^{*}{v}_h),\hfill \end{array}$$

(47)

for all $(v_h,q_h)\in (X_h,M_h)$. Setting $(v_h,q_h)=(d_t{e}_h^n,{\mu }_h^n)\Delta t$ in (47), we obtain

$$\begin{array}{cc}\hfill \Delta t\parallel d_t{e}_h^n{\parallel }_{0,h}^2& +\frac{\nu }{2}(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2-\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2)+\frac{\beta }{2}({\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)-{\mathcal{C}}_h({\mu }_h^{n-1},{\mu }_h^{n-1}))\hfill \\ \hfill & +\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{d}_t{e}_h^n)\Delta t+\tilde{b}(u_h\left(t_{n-1}\right),e_h^{n-1},I_h^{*}{d}_t{e}_h^n)\Delta t\hfill \\ \hfill & +\tilde{b}(e_h^{n-1},e_h^n,I_h^{*}{d}_t{e}_h^n)\Delta t=(E_n,I_h^{*}{d}_t{e}_h^n)\Delta t\hfill \end{array}$$

(48)

From (5) and Lemma 1, we get that

$$\begin{array}{cc}\hfill |\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{d}_t{e}_h^n)|& +|\tilde{b}(u_h\left(t_{n-1}\right),e_h^{n-1},I_h^{*}{d}_t{e}_h^n)|\hfill \\ \hfill \le \frac{1}{8}{\parallel d_t{e}_h^n\parallel }_{0,h}^2& +\kappa \left(\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2+{\parallel P_{h}B(u_h\left(t_{n-1}\right),·)\parallel }_{*}^2\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0^2,\hfill \\ \hfill \left|b\right(e_h^{n-1},e_h^n,I_h^{*}{d}_t{e}_h^n\left)\right|\Delta t& \le c\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2+\kappa \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0,\hfill \\ \hfill \left|\right(E_n,I_h^{*}{d}_t{e}_h^n\left)\right|& \le \frac{1}{4}\parallel d_t{e}_h^n{\parallel }_{0,h}^2+c{\parallel E_n\parallel }_0^2.\hfill \end{array}$$

Combining (48) with the above estimates and using Lemma 7 yields

$$\begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2-{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0^2& +{\nu }^{-1}\beta ({\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)-{\mathcal{C}}_h({\mu }_h^{n-1},{\mu }_h^{n-1}))+{\nu }^{-1}{\parallel d_t{e}_h^n\parallel }_{0,h}^2\Delta t\hfill \\ \hfill & \le \kappa \left(\right(\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2+\parallel P_{h}B(u_h\left(t_{n-1}\right),·){\parallel }_{*}^2)\hfill \\ \hfill & +\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2+c\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0+c{\parallel E_n\parallel }_{0,h}^2)\Delta t.\hfill \end{array}$$

(49)

Summing (48) from $n=1$ to $n=m$ and applying Lemmas 5 and 7, we obtain

$$\begin{array}{cc}\hfill \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}^m{\parallel }_0^2& +{\nu }^{-1}\Delta t\sum _{n=1}^m(\parallel d_t{e}_h^n{\parallel }_{0,h}^2\hfill \\ \hfill & \le \kappa \sum _{n=1}^m(\Delta t(\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2+\parallel P_{h}B(u_h\left(t_{n-1}\right),·){\parallel }_{*}^2)\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2\hfill \\ \hfill & \left(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2+\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0\right)+\Delta t{\parallel E_n\parallel }_{0,h}^2)\le \kappa .\hfill \end{array}$$

Combining the above estimate with (14) which yields(46).  □

Lemma 9.

Under the assumptions of Lemma 7, we have that

$$\tau \left(t_m\right){\parallel e_h^m\parallel }_{0,h}^2+\Delta t\sum _{n=1}^m\left(\parallel e_h^n{\parallel }_{0,h}^2+\tau \left(t_n\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n\parallel }_0^2+\tau \left(t_n\right){\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)\right)\le \kappa \Delta t^2.$$

(50)

for all $1\le m\le N$.

Proof. 

Setting $v=e_h^n,q={\mu }_h^n,z^n=e_h^n$ in (38) and $v_h={\mathsf{\Phi }}_h^{n-1},q_h={\mathsf{\Psi }}_h^{n-1}$ in (28) we obtain

$$\begin{array}{cc}\hfill (e_h^n,I_h^{*}{d}_t{\mathsf{\Phi }}_h^n)& -{\tilde{\mathcal{B}}}_h((e_h^n,{\mu }_h^n);(I_h^{*}{\mathsf{\Phi }}_h^{n-1},{\mathsf{\Psi }}_h^{n-1}))-\tilde{b}(e_h^n,u_h\left(t_n\right),I_h^{*}{\mathsf{\Phi }}_h^{n-1})\hfill \\ \hfill & -\tilde{b}(u_h\left(t_{n-1}\right),e_h^n,I_h^{*}{\mathsf{\Phi }}_h^{n-1})=(e_h^n,I_h^{*}{e}_h^n),\hfill \end{array}$$

(51)

$$\begin{array}{cc}\hfill (d_t{e}^n,I_h^{*}{\mathsf{\Phi }}_h^{n-1})& +{\tilde{\mathcal{B}}}_h((e_h^n,{\mu }_h^n);(I_h^{*}{\mathsf{\Phi }}_h^{n-1},{\mathsf{\Psi }}_h^{n-1}))+\tilde{b}(u_h^{n-1},e_h^n,I_h^{*}{\mathsf{\Phi }}_h^{n-1})\hfill \\ \hfill & +\tilde{b}(e_h^{n-1},u_h\left(t_n\right),I_h^{*}{\mathsf{\Phi }}^{n-1})=(E_n,I_h^{*}{\mathsf{\Phi }}_h^{n-1}).\hfill \end{array}$$

(52)

Adding (52) to (51), we obtain

$$\begin{array}{cc}\hfill \parallel e_h^n{\parallel }_{0,h}^2& =\frac{1}{\Delta t}((e_h^n,I_h^{*}{\mathsf{\Phi }}^n)-(e_h^{n-1},I_h^{*}{\mathsf{\Phi }}^{n-1}))-\tilde{b}(d_t{e}_h^n,u_h\left(t_n\right),I_h^{*}{\mathsf{\Phi }}^{n-1})\Delta t\hfill \\ \hfill & -\tilde{b}(e^{n-1},e^n,I_h^{*}{\mathsf{\Phi }}_h^{n-1})-(E_n,I_h^{*}{\mathsf{\Phi }}_h^{n-1}).\hfill \end{array}$$

(53)

It follows from (14) and (15) that

$$\begin{array}{cc}\hfill |\tilde{b}(e_h^{n-1},e_h^n,I_h^{*}{\mathsf{\Phi }}_h^{n-1})|& \le c\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0{\parallel {\tilde{A}}_h^{\frac{1}{2}}{\mathsf{\Phi }}_h^{n-1}\parallel }_0,\hfill \\ \hfill |\tilde{b}(d_t{e}_h^n,u_h\left(t_n\right),I_h^{*}{\mathsf{\Phi }}_h^{n-1})|& \le c\parallel d_t{e}_h^n{\parallel }_{0,h}\parallel P_h(·,u_h\left(t_n\right)){\parallel }_{*}{\parallel {\tilde{A}}_h^{\frac{1}{2}}{\mathsf{\Phi }}_h^{n-1}\parallel }_0,\hfill \\ \hfill \parallel (E_n,I_h^{*}{\mathsf{\Phi }}_h^{n-1}){\parallel }_0& \le \parallel {\tilde{A}}_h^{-1}{E}_n{\parallel }_0{\parallel {\tilde{A}}_h{\mathsf{\Phi }}_h^{n-1}\parallel }_0.\hfill \end{array}$$

Combining (53) with the above estimates yields

$$\begin{array}{cc}\hfill \parallel e_h^n{\parallel }_{0,h}^2\Delta t& =(e_h^n,I_h^{*}{\mathsf{\Phi }}_h^n)-(e_h^{n-1},I_h^{*}{\mathsf{\Phi }}_h^{n-1})+c\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0{\parallel {\tilde{A}}_h^{\frac{1}{2}}{\mathsf{\Phi }}^{n-1}\parallel }_0\Delta t\hfill \\ \hfill & +c\Delta t^2\parallel d_t{e}_h^n{\parallel }_{0,h}\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}\parallel {\tilde{A}}_h^{\frac{1}{2}}{\mathsf{\Phi }}_h^{n-1}{\parallel }_0+\parallel {\tilde{A}}_h^{-1}{E}_n{\parallel }_0{\parallel {\tilde{A}}_h{\mathsf{\Phi }}_h^{n-1}\parallel }_0\Delta t,\hfill \end{array}$$

(54)

with $e_h^0={\mathsf{\Phi }}^m=0$. Summing (54) from $n=1$ to $n=m$ and applying Lemma 6, we obtain

$$\begin{array}{cc}\hfill \Delta t\sum _{n=1}^m{\parallel e_h^n\parallel }_{0,h}^2& \le \kappa \Delta t\left(\sum _{n=1}^m\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\sum _{n=1}^m{\parallel {\tilde{A}}_h^{-1}{E}_n\parallel }_0^2\right)\hfill \\ \hfill & +c\Delta t^2\left(\sum _{n=1}^m{\parallel d_t{e}_h^n\parallel }_{0,h}^2\Delta t\right)\left(\sum _{n=1}^m{\parallel P_{h}B(·,u_h\left(t_n\right))\parallel }_{*}^2\Delta t\right).\hfill \end{array}$$

(55)

Applying Lemmas 4, 5, 7 and 8 in (55), we get and

$$\Delta t\sum _{n=1}^m{\parallel e_h^n\parallel }_{0,h}^2\le \kappa \Delta t^2,\forall 1\le m\le N.$$

(56)

Now, multiplying (55) by $\tau \left(t_n\right)$ and noting

$$\tau \left(t_n\right)\le \Delta t+\tau \left(t_{n-1}\right),\Delta t\le \tau \left(t_{n-1}\right),2\le n\le N,e_h^{n-1}=0,forn=1,$$

we get

$$\begin{array}{cc}\hfill \tau \left(t_n\right){\parallel e_h^n\parallel }_0^2& -\tau \left(t_{n-1}\right)\parallel e_h^{n-1}{\parallel }_0^2+\tau \left(t_n\right)(\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n))\Delta t\hfill \\ \hfill & \le c\parallel P_{h}B(·,u_h\left(t_n\right)){\parallel }_{*}^2\tau \left(t_{n-1}\right){\parallel e_h^{n-1}\parallel }_0^2\Delta t\hfill \\ \hfill & +c\left(\parallel e_h^{n-1}{\parallel }_0^2+\tau \left(t_n\right){\parallel {\tilde{A}}_h^{-1/2}{E}_n\parallel }_0^2\right)\Delta t.\hfill \end{array}$$

(57)

Summing (57) from $n=1$ to $n=m$ and applying Lemmas 3 and 5, we deduce that

$$\begin{array}{cc}\hfill \tau \left(t_m\right){\parallel e_h^m\parallel }_{0,h}^2& +\Delta t\sum _{n=1}^m\tau \left(t_n\right)(\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n))\hfill \\ \hfill & \le \kappa \Delta t\sum _{n=1}^m\left(\parallel e_h^{n-1}{\parallel }_{0,h}^2+\tau \left(t_n\right){\parallel A_h^{-1/2}{E}_n\parallel }_0^2\right).\hfill \end{array}$$

(58)

Applying Lemma 4 and (56) in (58), we have

$$\tau \left(t_m\right){\parallel e_h^m\parallel }_{0,h}^2+\Delta t\sum _{n=0}^m\tau \left(t_n\right)\left(\nu \parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n)\right)\le \kappa \Delta t^2.$$

(59)

this and (56) yield (50).  □

Lemma 10.

Under the assumptions of Lemma 7, the error $(e_h^n,{\mu }_h^n)$ satisfies the following bound:

$${\tau }^2\left(t_m\right)\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^m{\parallel }_0^2+\Delta t\sum _{n=1}^m{\tau }^2\left(t_n\right){\parallel d_t{e}_h^n\parallel }_{0,h}^2\le \kappa \Delta t^2,1\le m\le N.$$

(60)

Proof. 

Multiplying (49) by ${\tau }^2\left(t_n\right),$ noting ${\tau }^2\left(t_n\right)\le {\tau }^2\left(t_{n-1}\right)+3\tau \left(t_{n-1}\right)\Delta t$ and using Theorem 4, we have that

$$\begin{array}{cc}\hfill {\tau }^2\left(t_n\right)(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2& +{\nu }^{-1}\beta {\mathcal{C}}_h({\mu }_h^n,{\mu }_h^n))-{\tau }^2\left(t_{n-1}\right)(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2+{\nu }^{-1}\beta {\mathcal{C}}_h({\mu }_h^{n-1},{\mu }_h^{n-1}))\hfill \\ \hfill & +{\nu }^{-1}{\tau }^2\left(t_n\right)\parallel d_t{e}_h^n{\parallel }_{0,h}^2\Delta t\le c{\tau }^2\left(t_n\right){\parallel E_n\parallel }_{0,h}^2\Delta t\hfill \\ \hfill & +c\tau \left(t_{n-1}\right)(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0^2+{\nu }^{-1}\beta {\mathcal{C}}_h({\mu }_h^{n-1},{\mu }_h^{n-1}))\Delta t.\hfill \end{array}$$

(61)

Summing (61) from $n=1$ to $n=m$ and applying Lemmas 4 and 9, we have

$${\tau }^2\left(t_m\right)\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^m{\parallel }_0^2+\Delta t\sum _{n=1}^m{\tau }^2\left(t_n\right){\parallel d_t{e}_h^n\parallel }_{0,h}^2\le \kappa \Delta t^2.$$

(62)

 □

Lemma 11.

Under the assumptions of Lemma 7, the error ${\mu }_h^n={\overline{p}}_h\left(t_n\right)-p_h^n$ satisfies the following bound:

$${\tau }^2\left(t_n\right){\parallel {\overline{p}}_h\left(t_n\right)-p_h^n\parallel }_0^2\le \kappa \Delta t,1\le m\le N.$$

(63)

Proof. 

From Theorem 2, (5), (14) and (28), we deduce that

$$\begin{array}{cc}\hfill \parallel {\mu }_h^n{\parallel }_0\le c{\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n\parallel }_0& +c\parallel d_t{e}_h^n{\parallel }_{0,h}+c\left(\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0+{\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h\left(t_{n-1}\right)\parallel }_0\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n\parallel }_0\hfill \\ \hfill & +c\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h\left(t_n\right){\parallel }_0\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}{\parallel }_0+c{\parallel E_n\parallel }_{0,h}.\hfill \end{array}$$

Hence

$$\begin{array}{cc}\hfill {\tau }^2\left(t_n\right){\parallel {\mu }_h^n\parallel }_0^2\Delta t& \le \kappa \left(\tau \left(t_n\right)\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^n{\parallel }_0^2+\tau \left(t_{n-1}\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0^2\right)\Delta t+\kappa {\parallel {\tilde{A}}_h^{\frac{1}{2}}{e}_h^{n-1}\parallel }_0^2\Delta t^2\hfill \\ \hfill & +c{\tau }^2\left(t_n\right)\parallel d_t{e}_h^n{\parallel }_{0,h}^2\Delta t+c{\tau }^2\left(t_n\right){\parallel E_n\parallel }_{0,h}^2\Delta t.\hfill \end{array}$$

(64)

Summing (64) for n from $n=1$ to $n=m$ and applying Lemmas 4, 9 and 10, we deduce that

$$\Delta t\sum _{n=1}^m{\tau }^2\left(t_n\right){\parallel {\mu }_h^n\parallel }_0^2\le \kappa \Delta t^2,0\le m\le N,$$

(65)

which yields (63).  □

Theorem 5.

Under the assumptions of Lemma 7, the error $(u\left(t_n\right)-u_h^n,p\left(t_n\right)-p_h^n)$ satisfies the following optimal bounds:

$$\begin{array}{c}\hfill \Delta t\sum _{n=1}^N\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h^n){\parallel }_0^2+\tau \left(t_m\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_m\right)-u_h^m)\parallel }_0^2\le \kappa (h^2+\Delta t),\\ \hfill \Delta t\sum _{n=1}^N\tau \left(t_n\right)\parallel p\left(t_n\right)-p_h^n{\parallel }_0^2+{\tau }^2\left(t_m\right){\parallel p\left(t_m\right)-p_h^m\parallel }_0^2\le \kappa (h^2+\Delta t),\end{array}$$

(66)

for all $t_m\in (0,T]$.

Proof. 

Integration by parts directly can show

$$\begin{array}{cc}\hfill \Delta t\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h\left(t_n\right)){\parallel }_0^2& \le 2{\int }_{t_{n-1}}^{t_n}{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h)\parallel }_0^{2}dt\hfill \\ \hfill & +{\int }_{t_{n-1}}^{t_n}{(t-t_{n-1})}^2{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u_t-u_{ht})\parallel }_0^{2}dt.\hfill \end{array}$$

(67)

Summing (67) from $n=1$ to $n=N$, using (19) and noting $t-t_{n-1}\le \tau \left(t\right)$, we obtain

$$\Delta t\sum _{n=1}^N\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h\left(t_n\right)){\parallel }_0^2\le 2{\int }_0^T{\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h)\parallel }_0^{2}dt+\kappa h^2,$$

(68)

Combining (68) with (14), Theorem 3, Lemmas 7 and 9 yields the first estimate in (66).

$$\begin{array}{c}\hfill \Delta t\sum _{n=1}^N\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h^n){\parallel }_0^2+\tau \left(t_m\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_m\right)-u_h^m)\parallel }_0^2\le \kappa (h^2+\Delta t),\end{array}$$

(69)

Moreover, by (5), (9), (11), (14) and (20), we deduce that

$${\int }_0^t{\parallel p_h\parallel }_0^{2}ds\le c{\int }_0^t\left(\parallel u_{ht}{\parallel }_0^2+\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0^4+{\parallel f\parallel }_0^2\right)ds\le \kappa ,0\le t\le T.$$

(70)

Hence, we obtain from Theorems 1 and 3 that

$$\begin{array}{cc}\hfill \tau \left(t_1\right)\parallel p\left(t_1\right)& -{\overline{p}}_h\left(t_1\right){\parallel }_0\le \Delta t^{1/2}\left({\tau }^{1/2}\left(t_1\right){\parallel p\left(t_1\right)\parallel }_0+{\left({\int }_0^{t_1}{\parallel p_h\parallel }_0^{2}ds\right)}^{1/2}\right)\hfill \\ \hfill & \le \kappa \Delta t^{1/2},\hfill \end{array}$$

(71)

$$\begin{array}{cc}\hfill \tau \left(t_n\right)\parallel p\left(t_n\right)& -{\overline{p}}_h\left(t_n\right){\parallel }_0\le \tau \left(t_n\right)\parallel p\left(t_n\right)-\overline{p}\left(t_n\right){\parallel }_0+\tau \left(t_n\right){\parallel \overline{p}\left(t_n\right)-{\overline{p}}_h\left(t_n\right)\parallel }_0\hfill \\ \hfill & \le \frac{2}{\Delta t}{\int }_{t_{n-1}}^{t_n}(t-t_{n-1})\tau \left(t\right)\parallel p_t{\parallel }_{0}dt+\frac{2}{\Delta t}{\int }_{t_{n-1}}^{t_n}\tau \left(t\right){\parallel p-p_h\parallel }_{0}dt\hfill \\ \hfill & \le c\Delta t^{1/2}{\left({\int }_{t_{n-1}}^{t_n}{\tau }^2\left(t\right){\parallel p_t\parallel }_0^{2}dt\right)}^{1/2}+\kappa h\le \kappa \left(h+\Delta t^{1/2}\right).\hfill \end{array}$$

(72)

Combining (63) with (65)–(72) yields

$${\tau }^2\left(t_m\right){\parallel p\left(t_m\right)-p_h^m\parallel }_0^2\le \kappa \left(h^2+\Delta t\right),$$

(73)

for all $t_m\in (0,T]$.

Besides, noting that

$$\begin{array}{cc}\hfill \Delta t\sum _{n=2}^N\tau \left(t_n\right){\parallel p\left(t_n\right)-{\overline{p}}_h\left(t_n\right)\parallel }_0^2& \le \sum _{n=2}^N\frac{2\tau \left(t_n\right)}{\Delta t}({\parallel {\int }_{t_{n-1}}^{t_n}(p\left(t_n\right)-p\left(t\right))dt\parallel }_0^2\hfill \\ \hfill +\parallel {\int }_{t_{n-1}}^{t_n}(p-p_h){dt\parallel }_0^2)& \le 4\Delta t{\int }_{t_1}^T{\tau }^2\left(t\right)\parallel p_t{\parallel }_0^{2}dt+4{\int }_{t_1}^T\tau \left(t\right){\parallel p-p_h\parallel }_0^{2}dt.\hfill \end{array}$$

(74)

Recalling (Lemma 4.3) in [12], we have

$$\begin{array}{cc}\hfill (u_t-u_{ht},I_h^{*}v)& +{\tilde{\mathcal{B}}}_h((e_h,{\mu }_h);(I_h^{*}v,q))+\tilde{b}(u,u-u_h,I_h^{*}v)+\tilde{b}(u-u_h,u,I_h^{*}v)\hfill \\ \hfill & -\tilde{b}(u-u_h,u-u_h,I_h^{*}v)=0\forall (v,q)\in (X_h,M_h),\hfill \end{array}$$

(75)

where $e_h={\tilde{R}}_h(u,p)-u_h,{\mu }_h={\tilde{Q}}_h(u,p)-p_h.$ Using (9), (11), (14) and (75), we obtain

$${\int }_0^T\tau \left(t\right){\parallel {\mu }_h\parallel }_0^{2}dt\le c{\int }_0^T\left(\tau \left(t\right)\parallel u_t-u_{ht}{\parallel }_0^2+\parallel {\tilde{A}}_h^{\frac{1}{2}}(u-u_h){\parallel }_0^2{(\parallel u\parallel }^2+\parallel {\tilde{A}}_h^{\frac{1}{2}}{u}_h{\parallel }_0^2)\right)dt.$$

It follows from (5), (7), (18)–(20), and the above estimate that

$$\begin{array}{cc}\hfill {\int }_0^T\tau \left(t\right){\parallel p-p_h\parallel }_0^{2}dt& \le 2{\int }_0^T\tau \left(t\right)(\parallel p-{\tilde{Q}}_h{(u,p)\parallel }_0^2+\parallel {\mu }_h{\parallel }_0^2)dt\hfill \\ \hfill & \le c{h}^2{\int }_0^T{(\parallel Au\parallel }_0^2+{\parallel p\parallel }_1^2)dt+\kappa h^2\le \kappa h^2.\hfill \end{array}$$

(76)

Substituting (76) into (74) and using (73) and Theorem 1, we get

$$\Delta t\sum _{n=1}^N\tau \left(t_n\right){\parallel p\left(t_n\right)-{\overline{p}}_h\left(t_n\right)\parallel }_0^2\le \kappa (h^2+\Delta t).$$

(77)

Combining (77) with (65), we have

$$\Delta t\sum _{n=1}^N\tau \left(t_n\right){\parallel p\left(t_n\right)-p_h^n\parallel }_0^2\le \kappa (h^2+\Delta t).$$

(78)

This and (69) yield (66).  □

6. Single Numerical Example

In this section, some numerical results are computed to test the rationality of the theoretical analysis. Because it is difficult to obtain the analytical solution of the general problem governed by the Navier-Stokes equation, we show the relevant numerical results through an example with analytical solution for simplicity. So we consider the following model problem in the unit square area $[0,1]\times [0,1]$. Here this example might as well takes $\nu$ as 0.01. Only the velocities and pressure are given here. The right term f of the equations can be obtained by bringing the relationship between $p\left(x\right)$ and $u(x,y)=(u_1(x,y),u_2(x,y))$ into the NS equations, and the initial values of $u_1(x,y),u_2(x,y)$ and $p(x,y)$ can be obtained by bringing $t=0$ into the calculation.

Now, consider a unit square domain with an exact solution given by

$$\begin{array}{cc}\hfill & u(x,y)=(u_1(x,y),u_2(x,y)),\hfill \\ \hfill & u_1(x,y)=e^{-8{\pi }^2\nu }10x^2{(x-1)}^{2}y(y-1)(2y-1),\hfill \\ \hfill & u_2(x,y)=-e^{-8{\pi }^2\nu }10x(x-1)(2x-1)y^2{(y-1)}^2,\hfill \\ \hfill & p(x,y)=e^{-16{\pi }^2\nu }10(2x-1)(2y-1).\hfill \end{array}$$

(79)

f is determined by (3). It can be verified that such $u_1(x,y),u_2(x,y)$ satisfy the non divergence condition.

For simplicity, we can record the time and spatial discretization of the problem as follows:

$$\left\{\begin{array}{c}A{u}^{n+1}+N\left(u^n\right)u^{n+1}+B{p}^{n+1}=f^{n+1},\hfill \\ -B{u}^{n+1}+\beta C=0.\hfill \end{array}\right.$$

(80)

where the matrices in (80) correspond to the differential operators: $A\backsim -diag(v\Delta +1/\Delta t),N\left(u^n\right)\sim u^n·\nabla ,B\sim \nabla ,C\sim J_{h}div,C\sim {\mathcal{C}}_h(·,·)$, and I is the identity matrix. The right-hand side $f^{n+1}$ contains the source term.

To make the next iterations less complex, here are a few new notations. Let $W=u^{n+1}$, $q=p^{n+1}$, then we can further record the above equation as follows:

$$vA+N\left(w\right)w+Bq=f-B^{\mathrm{T}}w+\beta C=0.$$

Besides, in order to improve the calculation efficiency, we can generally adopt Newton iterative method to solve the above nonlinear problems. The typical calculation steps are as follows:

$$\begin{array}{cc}\hfill & \left(1\right)R=f-\left(vA{v}^{\mathrm{old}}+N\left(v^{\mathrm{old}}\right)\right)v^{\mathrm{old}}-B{q}^{\mathrm{old},}r=-B^T{v}^{\mathrm{old}};\hfill \\ \hfill & \left(2\right)\left(vA{v}^{\mathrm{mid}}+N\left(v^{\mathrm{old}}\right)\right)v^{\mathrm{mid}}+N\left(v^{\mathrm{mid}}\right)v^{\mathrm{old}}+B{q}^{\mathrm{mid}}=R,B^T{v}^{\mathrm{mid}}=r;\hfill \\ \hfill & \left(3\right)v^{\mathrm{new}}=v^{\mathrm{old}}+v^{\mathrm{mid}},q^{\mathrm{new}}=q^{\mathrm{old}}+q^{\mathrm{mid}}.\hfill \end{array}$$

(81)

It is worth noting that since Newton iterative method requires high initial values, we need to use the following Picard method to obtain the initial values of Newton iterative method:

$$\begin{array}{cc}\hfill & \left(1\right)R=f-\left(vA{v}^{\mathrm{old}}+N\left(v^{\mathrm{old}}\right)\right)v^{\mathrm{old}}-B{q}^{\mathrm{old}},r=-B^T{v}^{\mathrm{old}};\hfill \\ \hfill & \left(2\right)\left(vA{v}^{\mathrm{mid}}+N\left(v^{\mathrm{old}}\right)\right)v^{\mathrm{mid}}+B{q}^{\mathrm{mid}}=R,B^T{v}^{\mathrm{mid}}=r;\hfill \\ \hfill & \left(3\right)v^{\mathrm{new}}=v^{\mathrm{old}}+v^{\mathrm{mid}},q^{\mathrm{new}}=q^{\mathrm{old}}+q^{\mathrm{mid}}.\hfill \end{array}$$

(82)

In this way, we can finally transform the time-dependent Navier Stokes problem into a large-scale linear system of equations and solve it through such links: Euler time discretization → finite volume space discretization → Newton iterative transformation → Picard format transformation → large linear equations → solving equations.

Based on the above description, we can get some simulation results ($\beta =6$). The following Figure 3 is the result of one iteration based on the initial value. The time step here is $\Delta t=0.001$, and the spatial grid is divided into two congruent triangles on the basis of $100\times 100$ equidistant rectangular grid. We only pay attention to that only $1/5$ of the data in both X and Y directions are selected.

The following Figure 4 is the result after 1000 iterations with time step $\Delta t=0.001$. Here, the reference value of streamline is the same as that when $t=0.001$.

Compared Figure 2 with Figure 3, it can be seen from Figure 3 that the streamline and flow field are weakened accordingly. It is easy to understand that since the interpretation constructed here is decaying, both velocity and pressure show a decaying trend. This also shows that our numerical method maintains strong robustness, and the calculation results are more intuitive.

The following table shows the $L^2$ error order of velocity and pressure after 1000 iterations with $20\times 20,40\times 40,80\times 80$ grids selected under the condition of $\Delta t=0.001$.

From the above

Table 1. Convergence order of spatial solution the FVM ($\nu =0.01$, $\beta =6$).

h	$\frac{\parallel {\tilde{\mathit{A}}}_{\mathit{h}}^{\frac{1}{2}}(\mathit{u}-{\mathit{u}}_{\mathit{h}}){\parallel }_0}{\parallel {\tilde{\mathit{A}}}_{\mathit{h}}^{\frac{1}{2}}{\mathit{u}\parallel }_0}$	$\frac{\parallel \mathit{u}-{\mathit{u}}_{\mathit{h}}{\parallel }_0}{{\parallel \mathit{u}\parallel }_0}$	$\frac{\parallel \mathit{p}-{\mathit{p}}_{\mathit{h}}{\parallel }_0}{{\parallel \mathit{p}\parallel }_0}$
1/20	0.0699111	0.0035712	0.0771281
1/40	0.0394281	0.0009511	0.0440669
$\frac{1/20}{1/40}$	1.773	3.754	1.750
1/80	0.0209534	0.0002638	0.0241342
$\frac{1/40}{1/80}$	1.718	3.605	1.673

, we can see the optimal $L^2$ convergence rate, almost 2 for velocities and 1 for pressure are really obtained, which confirm the numerical analysis above.

The following Figure 5 shows the error curve calculated based on $100\times 100$ spatial grid with different time steps: $\Delta t=0.01,0.005,0.0025,0.00125$ (`dt’ in Figure 5 is $\Delta t$). It can be seen from here that the initial error tends to increase, but the error also decreases with the decrease of flow field energy, which shows that the time iteration is stable.

The following

Table 2. Numerical results of the FVM ($\nu =0.005$, $\beta =10$).

t	1	2	3	4	5	6	7	8
$\frac{\Delta t=0.01}{\Delta t=0.005}$	1.9513	1.9513	1.9456	1.9424	1.9404	1.9390	1.9380	1.9373
$\frac{\Delta t=0.05}{\Delta t=0.0025}$	1.9253	1.9253	1.9225	1.9210	1.9200	1.9193	1.9188	1.9184
$\frac{\Delta t=0.0025}{\Delta t=0.00125}$	1.9126	1.9126	1.9112	1.9104	1.9099	1.9096	1.9093	1.9092

shows the error ratio of different time steps: $dt$ = 0.01, 0.005, 0.0025, 0.00125:

The error ratio curve in the Figure 6 below shows the convergence of one section of time and is relatively stable.

Table 2 and Figure 6 tell us the optimal convergence rate of time is 1 which is consistent with the theoretical analysis.

Due to time constraints, our numerical results only show these. It is also worth noting that if the solution does not decay but increases, and if the growth rate is fast, the error of numerical results will increase with the increase of numerical calculation time. At that time, the method may not converge or inefficient. This requires a little attention in specific applications.

7. Conclusions

After detailed theoretical analysis, this article finally proves that if we use the finite volume method based on $P_1-P_0$ element to approximate the non-stationary Navier-Stokes equation, we can achieve the follow optimal numerical error estimation:

$$\begin{array}{c}\hfill \Delta t\sum _{n=1}^N\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_n\right)-u_h^n){\parallel }_0^2+\tau \left(t_m\right){\parallel {\tilde{A}}_h^{\frac{1}{2}}(u\left(t_m\right)-u_h^m)\parallel }_0^2\le \kappa (h^2+\Delta t),\\ \hfill \Delta t\sum _{n=1}^N\tau \left(t_n\right)\parallel p\left(t_n\right)-p_h^n{\parallel }_0^2+{\tau }^2\left(t_m\right){\parallel p\left(t_m\right)-p_h^m\parallel }_0^2\le \kappa (h^2+\Delta t).\end{array}$$

(83)

The optimal error estimate (83) shows that the time discretization of Euler method is 1 order and the space discretization is 2 order in this space-time full discretization finite volume method, which is consistent with the theoretical optimal order error estimation of $P_1-P_0$ element in solving differential equations. Although the proof process is challenging and cumbersome, the optimal result is also obvious and certain. However, this work is still helpful to reveal some special aspects of the finite volume method which is different from finite element and other methods in solving complex differential equations. Therefore, it is helpful to improve the corresponding numerical analysis theory.

In addition, with the continuous research and exploration of using neural network to solve differential equations recently [25,26,27], although very gratifying numerical results have been obtained, but the effectiveness and convergence of neural network-based methods for solving differential equations are still unclear, and there is a lack of theory. However, we know gradually that the key point of neural network in the calculation of differential equations is to use low-order functions for numerical approximation, which is similar to the basic principle of using low-order continuous functions to discretize and approximate differential equations, except that the former is global approximation while the latter is piecewise approximation. Therefore, improving and enriching the low-order function approximation theory also provides some important reference in further understanding the neural network in solving differential equations efficiently, It is worth studying.