A Finite Element Variational Multiscale Method Based on the Backward Euler Scheme for the Time-Dependent Navier-Stokes Equations

定义1  设Ω是在${{\mathbb{R}}^{2}} $上具有利普希茨连续边界的有界区域，那么有下面的Navier-Stokes方程：

其中u：$ \mathit{\Omega} \to {{\mathbb{R}}^{2}}$表示速度矢量，p：$ \mathit{\Omega} \to {{\mathbb{R}}}$是压力，f：$ \mathit{\Omega} \to {{\mathbb{R}}^{2}}$是流体驱动的体积力，ν＞0为流体粘性系数，u₀是使得$\nabla \cdot {\mathit{\boldsymbol{u}}_0} = 0 $的初始速度，并且${\mathit{\boldsymbol{u}}_t} = \frac{{\partial \mathit{\boldsymbol{u}}}}{{\partial t}} $.

定义2  对于定义1的数学问题，我们引进下面的希尔伯特空间：

其中：(·，·)表示空间L²(Ω)²或L²(Ω)的标准内积，($ \nabla \mathit{\boldsymbol{u}}$，$ \nabla \mathit{\boldsymbol{v}}$)和‖$ \nabla \mathit{\boldsymbol{u}}$‖₀为空间X上的一般标量积和范数.用字母c表示一个与时间步长和网格参数无关的正数而且可能在每个式子中代表的数值都不相同.

定义3^[5]   三线性项b(·，·，·)的定义为

它有如下的性质：

定义4   方程(1)-(4)的变分形式为：对于任意的t∈(0，T]，存在(u，p)∈X×M，使得

定义5^[7]   对于方程(8)的有限元离散，我们假设${T^\mu }\left( \mathit{\Omega} \right) = \left\{ K \right\}\left( {\mu = H, {\rm{ }}h, H > h} \right) $是准均匀的三角形网格剖分并且网格尺寸0＜μ＜1.细网格T^h(Ω)可以被认为是由粗网格加密而产生的.协调有限元(X_μ，M_μ)满足下面的inf-sup条件：存在常数β＞0使得

定义6^[8-11]   设速度空间X_μ满足$ \forall K \in {T^\mu }\left( \mathit{\Omega} \right), {\rm{ }}{v_\mu } \in {X_\mu }{\rm{, }}(\nabla {v_\mu }){|_K}$是线性的，那么在本文中提到的方法仅适用于速度限制在(P₂)²上的有限元对.如：Taylor-Hood元，P₂-P₀元和Scott-Vogelius(P₂-P₁^disc)元等.

定义7^[12]  我们定义$ V = \left\{ {v \in X:{\rm{ }}\left( {\nabla \cdot{\rm{ }}\mathit{\boldsymbol{v}}, {\rm{ }}\mathit{\boldsymbol{q}}} \right) = 0, {\rm{ }}\forall \mathit{\boldsymbol{q}} \in \mathit{\boldsymbol{M}}} \right\}{\rm{ }}{V_\mu } = \{ {\mathit{\boldsymbol{v}}_\mu } \in {X_\mu }:(\nabla \cdot{\rm{ }}{\mathit{\boldsymbol{v}}_\mu }, {\mathit{\boldsymbol{q}}_\mu }) = 0, \forall {\mathit{\boldsymbol{q}}_\mu } \in {M_\mu }\} $，则有如下的估计：

定义8^[12]   设$ {P_{{V_\mu }}}:Y \to V$是L²到V_μ上的正交投影，则满足$ (\mathit{\boldsymbol{\xi }} - {P_{{V_\mu }}}\mathit{\boldsymbol{\xi }}, {\mathit{\boldsymbol{v}}_\mu }) = 0, {\rm{ }}\forall \mathit{\boldsymbol{\xi }} \in Y, {\mathit{\boldsymbol{v}}_\mu } \in {V_\mu }$.

引理1^[13]  离散Gronwall引理：对于任意整数n≥0，令Δt，H和a_n，b_n，c_n，d_n是非负数，满足

和Δtd_n＜$1\forall n $.有

定义9^[5]   设数值格式中出现的变分多尺度稳定项为

这里$ {\smallint _{K, {\rm{ }}s}}\left( \cdot \right){\rm{d}}x$表示K上适当的高斯积分，该积分对于次数不超过s(s=m，1，m≥2)的多项式是准确的. α＞0是一个自定义的稳定项参数.

定义

其中P₀是常量元素K的空间.那么标准的L²-正交投影Π_μ：$ L \to {L_\mu }$有下面的性质：

注1  稳定项(11)还可以表示为：

根据定义9，我们给出Navier-Stokes方程的标准的有限元变分多尺度方法.

方法1  标准的有限元变分多尺度方法^[5].

给定u_μ⁰，存在${(\mathit{\boldsymbol{u}}_\mu ^{n + 1}, \mathit{\boldsymbol{p}}_\mu ^{n + 1})_{n \ge 0}} \in ({X_\mu }, {\rm{ }}{M_\mu }) $使得：

令时间步长的尺寸Δt满足0＜Δt＜1，${t_n} = n\Delta t, {\rm{ }}n = 0, {\rm{ }}1, {\rm{ }} \ldots , {\rm{ }}N - 1 $，和$ N = \frac{T}{{\Delta t}}$. φ¹表示函数φ在时间t₁时的值，并且${\varphi ^{\frac{1}{2}}} = \frac{1}{2}({\varphi ^1} + {\varphi ^0}) $.初始速度$\mathit{\boldsymbol{u}}_\mu ^0 = {P_{{V_\mu }}}{\mathit{\boldsymbol{u}}_0} $.

引理2^[5]  令f∈L²(0，T；H^-1(Ω)²)和u₀∈L²(Ω)².则格式(15)的解是稳定的且满足任意的0＜l≤N

引理3^[5]   Navier-Stokes方程的精确解(u， p)满足u∈L^∞(0，T；H¹(Ω)²)，u_tt∈L^∞(0，T；H¹(Ω)²)，和u_ttt∈L^∞(0，T；L²(Ω)²).那么由格式(15)计算的全离散解有如下估计：

两水平有限元变分多尺度方法如下：

方法2  两水平有限元变分多尺度方法.

给定u_H⁰，u_h⁰，存在(u_h¹，p_h¹)_n≥0∈X_h×M_h.

1) 寻找粗网格上的一个解(u_H¹，p_H¹)_n≥0∈(X_H，M_H)使得

2) 寻找细网格上的一个解(u_h¹，p_h¹)_n≥0∈(X_h，M_h)使得

注2

定理1   Navier-Stokes方程的精确解(u，p)满足u∈L^∞(0，T；H¹(Ω)²)，u_tt∈L^∞(0，T；L²(Ω)²)，那么由式(19)-(20)得到的全离散解有下面的估计式：

证  令$ (\mathit{\boldsymbol{u}}({t_{n + 1}}), \mathit{\boldsymbol{p}}({t_{n + 1}})) = ({\mathit{\boldsymbol{u}}^{n + 1}}, {\mathit{\boldsymbol{p}}^{n + 1}}), {\rm{ }}({\mathit{\boldsymbol{e}}^{n + 1}}, {\mathit{\boldsymbol{\eta }}^{n + 1}}) = ({\mathit{\boldsymbol{u}}^{n + 1}} - {\mathit{\boldsymbol{u}}_h}^{n + 1}, {\mathit{\boldsymbol{p}}^{n + 1}} - {\mathit{\boldsymbol{p}}_h}^{n + 1}), n = 1, 2, \ldots , N - 1$.当时间t=t_n+1时，式(8)减去式(20)得：$\forall ({\mathit{\boldsymbol{v}}_h}, {\rm{ }}{\mathit{\boldsymbol{q}}_h}) \in ({X_h}, {\rm{ }}{M_h}) $

由于

所以

从而令

且

将v_h=φ_hⁿ⁺¹代入式(22)利用式(5)和$ 2\left( {a - b, {\rm{ }}a} \right) = {a^2} - {b^2} + {\left( {a - b} \right)^2}$得：

其中：λ_h∈M_h是pⁿ⁺¹的近似值.

根据定义7中V_h的定义，有

由于χⁿ⁺¹-χⁿ⊥V_h且φ_hⁿ⁺¹∈V_h，利用投影算子我们有

所以由式(23)可得：

现在利用施瓦兹不等式，Young不等式和式(13)对式(24)的右边进行估计：

利用式(5)-(7)和Young不等式估计下面的三线性项

对式(24)最后一项$ \left( {{\mathit{\boldsymbol{u}}_t}({t_{n + 1}}) - \frac{{{\mathit{\boldsymbol{u}}^{n + 1}} - {\mathit{\boldsymbol{u}}^n}}}{{\Delta t}}, {\varphi _h}^{n + 1}} \right)$利用泰勒展式，施瓦兹不等式和Young不等式，得：

将式(25)-(33)代入式(24)，有

将式(34)乘以2Δt，从n=1加到N-1且φ_h⁰=0.得：

最后，当Δt足够小时，应用离散Gronwall引理和式(10)有

利用三角不等式，可得式(21).

在本节中，我们利用FreeFem++软件^[14]进行一些实验验证理论预测的正确性. Taylor-Hood元用于空间离散化.粗网格上的非线性迭代的迭代限差为10^-6，并且非线性系统由牛顿迭代法求解.值得注意的是，对于非线性迭代，标准的变分多尺度方法和两水平变分多尺度方法的稳定项可以近似为

其中j表示非线性迭代的次数.

选择Navier-Stokes方程的精确解为：

其中解的区域Ω=[0, 1]×[0, 1] $\subset {{\mathbb{R}}^{2}} $，且ν=1.0×10^-7，T=0.01.

定理1的误差估计在理论上预测了能量范数对于$\mathcal{O} $ (h²)的空间收敛速度.由此可设α=0.1 h²，$ H={{h}^{\frac{1}{2}}}$.表 1给出了数值结果.从表 1可知：本文的方法对空间和时间离散是一阶收敛的，同时也表明我们的理论预测是正确的.

为了对比本文的方法和文献[5]中的方法，令网格尺寸$ H=\frac{1}{8}, h=\frac{1}{64}$，时间步长$\Delta t=\frac{1}{800} $，但是粘性系数分别为ν=0.01，0.001，0.000 1，0.000 01和0.000 001的情况下求解.数值结果在表 2中给出.由表 2可知：本文的两水平变分多尺度方法得到的解精确度和标准的网格变分多尺度方法大体一致，但是本文的方法可以节约一半以上的计算时间.

本文给出了全离散速度的误差估计.对比标准的变分多尺度方法，本文的方法可以节约很多计算时间.