理论数学  >> Vol. 8 No. 4 (July 2018)

基于Navier-Stokes方程分析飞机抖振问题
A Class of Spectral Method for Solving Unsteady Navier-Stokes Equations

DOI: 10.12677/PM.2018.84060, PDF, HTML, XML, 下载: 665  浏览: 941  科研立项经费支持

作者: 李建森, 刘 俊, 张 姣, 范馨月:贵州大学数学与统计学院,贵州 贵阳

关键词: Navier-Stokes方程非线性Galerkin法抖振问题Legendre谱方法误差估计Navier-Stokes Equation Nonlinear Galerkin Method Buffet Problem Legendre Spectral Method Error Estimation

摘要: 针对不定常的Navier-Stokes方程,给出了一种非线性Galerkin-Legendre谱方法(GL方法),该方法是将勒让德谱方法和非线性Gelerkin方法结合起来,使得速度与压力的混合元空间不需要满足离散的Babuška-Brezzi不等式条件,给出了GL方法的谱格式,并给出该方法的稳定性证明和误差估计。关于飞机抖振问题已经由基于不可压缩N-S方程的方法的研究取代了原来粘、位流分开计算的研究,因此用谱方法研究飞机抖振问题,通过误差分析,减少阻力,达到节约能源的目的。
Abstract: This paper focuses on the issue about using numerical computation to solve unsteady Navier-Stokes equations. We present a scheme called nonlinear Galerkin-Legendre Spectral method, which combines the Legendre spectral method with nonlinear Galerkin-Legendre method. Therefore, it is not necessary to meet Babuska-Brezzi inequality condition in the velocity space and pressure space. In this paper, we give a spectral scheme of nonlinear Galerkin-Legendre methods for solving 2D N-S equations. What’s more, we also provide the error estimation and proof the stability of this scheme. The method based on the separate calculation of viscous flow and potential flow, was used to solve the problem of plane buffeting, and has been replaced by the method based on the incompressible N-S equations. Therefore, it is necessary to use the spectral method to study the chattering problem of aircraft and to reduce the resistance by error analysis and to save energy.

文章引用: 李建森, 刘俊, 张姣, 范馨月. 基于Navier-Stokes方程分析飞机抖振问题[J]. 理论数学, 2018, 8(4): 450-458. https://doi.org/10.12677/PM.2018.84060

2. 经典Galerkin方法下的结论

2.1. 预备知识

双线性 设 a ( u , v ) 是一个双线性泛函 [6] ,则

a ( u , v ) = Ω u x v x d x ( u , v H 0 1 (Ω))

三线性设 a ( u , v , w ) 是一个三线性泛函,则

a ( u , v , w ) = 1 2 Ω i , j = 1 n [ u i v j x i w j u i w j x i v j ] d x ( u , v , w H 0 1 (Ω))

对于 a ( u , v , w ) 有下列性质(参见 [7] ):

a ( u , v , w ) = a ( u , v , w ) ; a ( u , v , v ) = 0

a ( u , v , w ) = C u H 0 1 v H 0 1 w L 2 ( Ω ) ; u , v H 0 1 ( Ω ) , w L 2 (Ω)

2.2. Galerkin方法的弱形式

Ω 是一个有界区域。考虑下面定常N-S方程 [8] :

问题I求 ( u , p ) 满足

{ u t ν Δ u + u u + p = f ( Ω ) u = 0 ( Ω ) u = 0 ( Ω ) (1)

其中u表示速度,p表示压力,f表示外力, ν 是雷诺数的倒数,是常数。

如下记号:

X = H 0 1 ( Ω ) , M = L 0 ( Ω ) = { q L 2 ( Ω ) ; Ω q d x = 0 }

问题I的变分形式为:

问题I* ( u , p ) ( X , M ) 满足

( u t , v ) + a ( u , v ) + a 1 ( u , u , v ) b ( p , v ) + b ( q , u ) = ( f , v )

( v , q ) ( X , M )

其中

a ( u , v ) = ν Ω u v d x , b ( q , v ) = Ω q d i v v d x

,

( f , v ) = Ω f v d x .

X N M N 分别是X和M空间中次数不超过N的多项式空间,则(1)的谱格式为:

{ ( u t , v N ) + a ( u N , v N ) + a 1 ( u N , u N , v N ) b ( p N , v N ) = ( f , v N ) ( u N , q N ) = 0 , v N X N , q N M N (2)

从文献 [7] 可以得到以下结论。

定理1.1 如果Babuška不等式成立,那么

u u N 1 + β N p p N L 2 N 1 m u m + N s u s

其中N为插值多项式的个数, β N 满足Babuška不等式。

3. 非线性Galerkin-Legendre谱方法

3.1. 准备知识

先引进一些记号。设

S N = s p a n { L 0 ( x ) , L 1 ( x ) , , L N ( x ) } ,

Q 2 N = s p a n { L N 1 ( x ) , L N ( x ) , , L 2 N ( x ) } ,

W 2 N = { w Q 2 N / w ( ± 1 ) = 0 } ,

S 2 N = Q 2 N S N 2 = s p a n { L 0 ( x ) , L 1 ( x ) , , L 2 N ( x ) } ,

V 2 N = V N W 2 N = { v S 2 N / v ( ± 1 ) = 0 } .

注意到 V N 中的元素均在 x = ± 1 处为零,即在边界上为零。

P N L 2 ( Ω ) V N 的正交投影算子,即

( u P N u , v ) = 0 , v S N , u L 2 (Ω)

Π N H 0 1 ( Ω ) V N 的投影算子,即

a ( u Π N u , v ) = 0 , v S N , u H 0 1 (Ω)

3.2. 非线性Galerkin-Legendre方法

问题I的非线性Galerkin-Lengedre格式为 [9]

问题I** u N V N ,和 u 2 N W 2 N ,使得

{ d d t ( u N , ϕ ) + a ( u N , ϕ ) + a 1 ( u N + u 2 N , u N + u 2 N , ϕ ) b ( p N , ϕ ) = ( f , ϕ ) , ϕ V N a ( u 2 N , ω ) + a 1 ( u N + u 2 N , u N , ω ) b ( p 2 N , ω ) = ( f , ω ) ( q , v ) = 0 , ω W 2 N (3)

引理2.1根据三线性泛函的性质,有以下结论

u N 2 C ( u 1 2 + f 2 )

以上引理在 [10] 中已经得证。

定理2.1 在引理2.1成立下,对充分大的N,我们有

u u N u 2 N C N 1 s + K Δ t 2 , p p N p 2 N M N d 1 2 m

其中C为与 u s 有关的常数,K为与 u t 有关的常数,M为与 p m 有关的常数, d 1 ,N为插值多项式的个数。

证明:因为 Π 2 N 的定义满足

{ d d t ( u , ϕ + ω ) + ν ( ( 2 N u ) x , ( ϕ + ω ) x ) + a 1 ( u , u , ϕ + ω ) = ( f , ϕ + ω ) , ϕ V N , ω W 2 N b ( q , ( ϕ + ω ) ) = 0 , ϕ V N , ω W 2 N

e = 2 N u ( u N + u 2 N ) , ξ = N u u N , η = 2 N u u N u 2 N .

显然 ξ V N η W 2 N e = ξ + η ,特别地

e 1 2 = ξ 1 2 + η 1 2

由(1)式与(3)相减

d d t ( u , ϕ + ω ) d d t ( u N , ϕ ) + a ( e , ϕ + ω ) + a 1 ( u , u , ϕ ) a 1 ( u N + u 2 N , u N + u 2 N , ϕ ) + a 1 ( u , u , ω ) a 1 ( u N + u 2 N , u N , ω ) = 0

整理得到

a ( e , ϕ + ω ) + ( ξ t , ϕ ) = [ a 1 ( u N + u 2 N , u N + u 2 N , ϕ ) a 1 ( u , u , ϕ ) ] + [ a 1 ( u N + u 2 N , u N , ω ) a 1 ( u , u , ω ) ] + ( ( Π N u u ) t , ϕ ) ( u t , ω )

不妨取 ϕ = ξ ω = η ,于是

e 1 2 + 1 2 d d t ξ 2 = [ a 1 ( u N + u 2 N , u N + u 2 N , ξ ) a 1 ( u , u , ξ ) ] + [ a 1 ( u N + u 2 N , u N + u 2 N , η ) a 1 ( u , u , η ) ] + [ ( ( Π N u u ) t , ξ ) ( u t , η ) ] I 1 + I 2 + I 3

其中

I 1 = a 1 ( u N + u 2 N , u N + u 2 N , ξ ) a 1 ( u , u , ξ ) = a 1 ( u N + u 2 N , u 2 N u , ξ ) a 1 ( u N + u 2 N , e , ξ ) a 1 ( u 2 N u + e , u , ξ ) I 11 + I 12 + I 13

I 2 = a 1 ( u N + u 2 N , u N , η ) a 1 ( u , u , η ) = a 1 ( u N + u 2 N , u N u , η ) a 1 ( u N + u 2 N , ξ , η ) a 1 ( u 2 N u + e , u , η ) I 21 + I 22 + I 23

注意到 I 12 + I 22 = 0 。根据Young不等式

I 13 + I 23 = a 1 ( u 2 N u + e , u , e ) C ( u 2 N u 1 + e 1 ) u 1 e ν 8 e 1 2 + M ( u 2 N u 1 2 + e 2 )

因为 N 1 C u 1 ,应用Young不等式,Poincare不等式以及 η 1 e 1 ,对充分大的N,有

I 11 + I 21 = a 1 ( e 2 N u , u 2 N u , e ) + a 1 ( e N u , 2 N u N u , η ) ν 4 e 1 2 + M ( u 2 N u 1 2 + e 2 ) + M N 2 N u 2 N u 1 2

因为 η W 2 N

I 3 = [ ( ( Π N u u ) t , ξ ) ( ( I P N 2 ) u t , η ) ] ν 8 ( ξ 1 2 + η 1 2 ) + C ( ( Π N u u ) 1 2 + ( I P N 2 ) u t 1 2 ) ν 8 e 1 2 + C ( ( Π N u u ) 1 2 + ( I P N 2 ) u t 1 2 ) ν 8 e 1 2 + Δ t 2 u t 1 2

综上

e 1 2 + 1 2 d d t ξ 2 M e 2 + M N 2 N u 2 N u 1 2 + Δ t 2 u t 1 2 .

由Gronwall不等式

u u N u 2 N N 1 s u s + Δ t 2 u t 1 2

其中 s 1 ,N为基函数个数。

现在对压力作误差估计。

V N = { v N H 0 1 ( Ω ) / ( q N , v N ) = 0 , q N M N } 在文献 [11] 中有以下结果:

引理2.2 如果Babuška条件成立,则

β N p p N p 2 N inf u N V N u 2 N W 2 N | u u N u 2 N | + inf q N M N p q N

下面我们根据Babuška不等式性质以及速度误差的相关结论有

(1)式减去(3)式

e 1 2 ( p p N p 2 N , v N ) = 0

由Babuška不等式,找任意 q N M N q 2 N M 2 N

β N q p N p 2 N sup v N H 0 1 ( Ω ) ( q N p N p 2 N , v N ) | v N | = sup v N H 0 1 ( Ω ) α ( u u N u 2 N , v N ) + ν ( ( u u N u 2 N ) , v N ) ( p q N , v N ) | v N |

注意到 v H 0 1 ( Ω )

因此

β N q p N p 2 N β N inf { p q N + q N p N p 2 N } inf q N M N p q N + u u N u 2 N

根据 e 1 2 N 1 s u s ,考虑到 P N P N 2 正交,文献 [5] 有以下结论

β N C ( α , ν ) N d 1 2 , ( d = 2 3 )

因此

N d 1 2 p p N p 2 N N m p m

综上结论,我们有如下误差估计

| u u ˜ N | 1 + p p ˜ N L 2 C N 1 s + M N d 1 2 m + K Δ t 2

其中C是与 u s 有关的常数,K是与 u t 有关的常数,M是与 p m 有关的常数, d 2

4. 数值算例

我们考虑如下的N-S方程(飞机抖振数学模型):

{ u t ν Δ u + u u + p = 0 ( Ω ) d i v u = 0 ( Ω ) u = 0 ( Ω ) (4)

其中 Ω = [ 1 , 1 ] × [ 1 , 1 ] 0 t T ,速度空间是一个狄利克雷边界,时间上使用差分法(显示欧拉格式),压强空间上用次数不超过N的多项式空间近似逼近,得到如下格式:

( u N k + 1 u N k Δ t , v N ) + ν ( ( u N k + 1 ) x , ( v N ) x ) + ( u N k + 1 ( u N k ) x , v N ) + ( p N k + 1 p N k Δ t , v N ) = 0 (5)

现对u做Legendre多项式插值,即

u N k k = 0 N 2 u ˜ k ϕ k ( x ) ,

对p做Legendre多项式插值,即

p N k k = 0 N p ˜ k ϕ k

u ˜ k p ˜ k 是插值系数, ϕ k 是插值基函数。

将这些近似多项式和代入上述(5),则有

[ 1 Δ t M + K D T D 0 ] [ U k + 1 p k + 1 p k ] = [ 1 Δ t M U k D T p k 0 ]

这是一个线性方程组形式,算法的步骤如下:

{ Step 1 : ( 1 Δ t M + K ) U ˜ k + 1 = 1 Δ t M U k D T p k Step 2 : D M 1 D T Φ k + 1 = 1 Δ t D U ˜ k + 1 Step 3 : U k + 1 = U ˜ k + 1 Δ t M 1 D T Φ k + 1 Step 4 : p k + 1 = Φ k + 1 + p k

已知方程(4) ( u , p ) 的精确解

u ( x , y , t ) = π sin t ( sin 2 π y sin 2 π x , sin 2 π x sin 2 π y )

p ( x , y , t ) = sin t cos π x sin π y

处理不可压缩N-S方程的谱方法与传统的变分法以及有限元法处理N-S方程的区别很明显,从表1表2可以看出相同的时间间隔,谱方法得出来的结果比经典格式的结果更加精确;表3中的误差表

Table 1. The comparison of error estimates for velocity in two different formats with N = 16 and T = 1

表1. N = 16和T = 1两种不同格式下的速度的误差估计的比较

Table 2. The comparison of error estimates for pressure in two different formats with N = 16 and T = 1

表2. N = 16和T = 1两种不同格式下的压强的误差估计的比较

Table 3. The error estimates and convergence order for different time intervals and N with T = 1 and T = 2

表3. T = 1和T = 2时不同时间间隔和不同N的误差估计以及收敛阶

示在T = 1、T = 2时刻下的速度误差和压力误差总和,从表3中可以看出不同时间间隔下,时间间隔对方程的误差估计影响更大,该方程的收敛阶可以用 Δ t 2 来刻画,并且速度收敛阶是1;在相同时间间隔和不同N下,N对方程的误差估计影响更大,该方程的收敛阶可以用N来刻画,且随着基函数的个数N的增加,误差越来越少,所得的精度越高,收敛结果也符合谱方法的收敛阶,以N的负指数次幂速度收敛精确解。

5. 结语

研究提出了一类求解定常纳维斯托克斯方程的谱方法,基于二维不可压缩N-S方程给出了LG格式,并且得出相关误差结果,比较了该方法与经典方法和有限元方法的区别,从数值结果我们可以发现LG谱方法下速度逼近的收敛率是最优的,而压强逼近的收敛率也很好。从表格上可以看出当解 ( u , p ) 是光滑的,即使对较小的N也能得到很好的结果,而且计算量不大,这充分体现了LG谱方法的优越性。基于N-S方程的GL方法的数值分析,可以看出对于飞机抖振问题,可以通过减少与空气阻力,达到节约能源的目的。

基金项目

贵州省科学技术基金项目(黔科合J字[2013]2128),贵州省大数据重点实验室开放课题 (2017BDKFJJ012)。

NOTES

*通讯作者。

参考文献

[1] Temam, R. (2013) Navier-Stokes Equations, Theory and Numerical Analysis. AMS Chelsea Edition, American Mathematical Society, Rhode Island.
[2] Marion, M. and Temam, R. (1989) Nonlinear Galerkin Methods. SIAM Journal on Numerical Analysis, 2, 1139-1157.
https://doi.org/10.1137/0726063
[3] 罗振东, 朱江, 王会军, 等. 定常的Navier-Stokes方程的非线性Galerkin/Petrov最小二乘混元法[J]. 应用数学和力学, 2011, 7(7): 697-887.
[4] Shen, J. and Wang, L.-L. (2011) Spectral Methods, Springer, Heideberg, Dordrecht, London, New York.
[5] 何银年, 黄艾香. 非定常N-S方程的非线性Galerkin算法及其误差估计[J]. 系统科学与数学, 1991, 19(1): 116-122.
[6] 向新民. 谱方法的数值分析[M]. 北京: 科学出版社, 2000: 153-156.
[7] France, L.P. and Hughes, T.J. (1998) Two Classes of Mixed Finite Element Methods. Computer Methods in Applied Mechanics and Engineering, 69, 89-129.
https://doi.org/10.1016/0045-7825(88)90168-5
[8] 牟让科, 杨永年. 飞机抖振问题研究进展[J]. 应用力学学报, 2001, 18(s1): 143-149.
[9] Shen, J. and Yu, H. (2010) Efficient Spectral Sparse Grid Methods and Applications to High Dimensional Elliptic Problems. SIAM Journal on Scientific Computing, 32, 3228-3250.
https://doi.org/10.1137/100787842
[10] 胡园园, 谢江, 张武, 等. 二维不可压缩Navier-Stokes方程的并行谱有限元法求解[J]. 计算机应用, 2017(1): 42-47
[11] Liu, J. and Fan, X.Y. (2017) Legendre-Galerkin Methods for Nonlinear Partial Differential Equations. Advance in Intelligent Systems Research, 23, 1951-6851.
https://doi.org/10.2991/icmia-17.2017.71