1. 引言
阻塞形势可视为一种生命史相对较长而结构持续稳定的大振幅孤立系统 [1] [2],为此,Malguzzi和Malamotte-Rizzoli [3] [4] 采用KDV (孤立波)动力学方法来阐述阻塞高压的形成、维持及崩溃。然而,在大气中观测到的阻塞高压和切断低压都不是波状流动,而是孤立涡动 [5]。为此,Flierl [6] 和Mc Willianms等 [7] [8] 在正压流体中,利用强非线性作用,求出了孤立涡解。这种孤立涡解的流线分布,形成反对称的一对偶极子,两个极子相互非线性作用,导致整个偶极子向东或向西移动。孤立涡具有非频散特征,能较好地揭示阻塞的长时间维持。然而,只有当第一类Bessel函数的内解波数在很窄的范围内,才会有这类偶极子型的孤立涡解。
罗德海 [1] 提出的非线性Shrӧdinger型包络Rossby孤立波理论,既能反映偶极子阻塞的孤立波特征,又能体现偶极子阻塞的衰减机制,较好地揭示了大气中偶极子阻塞的形成过程和衰退机制。为此,本研究通过对具有耗散效应的正压准地转涡度方程的简化,获得带有强迫耗散的非线性Shrӧdinger方程。由于本研究在非线性Shrӧdinger方程中保留了非齐次的耗散强迫项,即扰动项,要严格求解一般是不可能的,而将扰动项视为小量处理,采用Yan [9] [10] [11] 所建立的基于分离变量法的直接微扰理论对非齐次的非线性Shrӧdinger方程进行分析。Liu [12] 等采用直接微扰法,讨论了北极变暖引起中高纬度阻塞加强、维持时间变长的一种动力机制,解决了“北极变暖为何会容易使得偶极型阻塞维持和加强,从而导致冬季北半球的极端冷事件 [13] - [19] ”这一困扰气象界多年的难题。本研究将采用类似的方法,结合孤立子的拓扑结构,讨论强迫耗散对阻塞高压的作用。
本研究首先通过多尺度变换和摄动法将正压模式简化为带有扰动项的非线性Shrӧdinger方程;其次,分析不同基本纬向速度下的定常Shrӧdinger孤立子的形态;最后,通过孤立子直接微扰理论,研究耗散对阻塞形势的影响。
2. 具有强迫耗散准地转模式的Shrӧdinger方程
考虑具有地形效应的正压准地转涡度方程:
(1)
式中,
为地转流函数;
为科里奥利力变率;
是地廓线;
,为二维Jacobi算子。
侧边界条件为刚壁条件的无量纲形式:
(2)
取流函数为
(3)
其中,
是常数,它是切变气流中线性长波的相速度;
是无量纲Rossby数,表征非线性的强弱,
为扰动函数。将(3)式带入方程(1),得到扰动流函数
的方程
(4)
同样,将(3)式带入边界条件(2)中,得到关于扰动
的边界条件
(5)
令
(6)
将(6)式带入(4)式,略去撇号“
”,得到
(7)
下面,我们将以无量纲的扰动方程(7),讨论强迫耗散对偶极子型阻塞的影响。
3. 具有强迫耗散准地转模式的Shrӧdinger方程
根据(7)式,要使耗散与非线性之间达到平衡,需满足以下条件:
![](//html.hanspub.org/file/2-2810120x26_hanspub.png)
其中,
。而在大气运动中,耗散往往比非线性效应更弱,为此,我们设
(8)
其中,
。
由于大气运动中,存在多种时空尺度的运动。根据多尺度变换方法 [20] [21],引入缓变坐标:
![](//html.hanspub.org/file/2-2810120x30_hanspub.png)
由于
,则
为缓变量。根据叶笃正和巢纪平 [22] 的研究结果,
为大尺度特征时间,T为慢变外源强迫特征时间,X为长波特征尺度,
为超长波特征尺度。于是,作如下变换(Pedlosky 1970):
(9)
根据摄动法,将对扰动流函数
,按WKB (小参数)方法 [23] 展开:
(10)
考虑到扰动流的边界条件(5),有如下的扰动方程边界条件:
(11)
将(8)式、(9)式和(10)式带入到扰动方程(7)中,合并
项,可以得到关于
的各阶问题。对于
问题,满足如下方程:
(12)
其中
![](//html.hanspub.org/file/2-2810120x43_hanspub.png)
在边界条件(11)下,方程(12)有如下的解形式:
(13)
其中,A为复振幅,k为纬向波数,
为常数,cc表示其前项的共轭。关于
的特征问题如下:
![](//html.hanspub.org/file/2-2810120x47_hanspub.png)
对于
问题,满足如下方程:
(14)
方程(14)是关于
的,将(13)式带入到(14)式中,得到
(15)
由于算子
具有一个形如
的解,而方程的(15)具有
的项,从而将导致共振。为此,要消除共振项,即使
![](//html.hanspub.org/file/2-2810120x55_hanspub.png)
成立,其中,
。这时方程(15)可以化为
(16)
根据边界条件(11),可知方程(16)的解具有如下形式:
(17)
比较(16)式和(17)式,有
(18)
其中,
![](//html.hanspub.org/file/2-2810120x60_hanspub.png)
为此,根据(18)式,可将B取为如下形式:
(19)
将(19)式带入到方程(18)中,得到关于
的特征方程如下:
![](//html.hanspub.org/file/2-2810120x63_hanspub.png)
对于
问题,满足如下方程:
(20)
与
问题一样,通过消除共振项,得到
问题(20)的消除久期项的条件,即关于振幅A方程
(21)
其中
![](//html.hanspub.org/file/2-2810120x69_hanspub.png)
下面,将以准地转流的扰动振幅方程(21)为基础,讨论阻塞流场的几何结构,以及强迫耗散对阻塞结构的影响。
4. 无强迫耗散的阻塞流场
在讨论地形强迫对阻塞流场的影响之前,需要对阻塞流场的拓扑结构作清晰的了解。为此,令
,则(21)式化为:
(22)
对于没有地形强迫的情况,即
,则(22)式为标准的非线性Schrӧdinger方程,此时,其有如下形式的单孤立子解:
(23)
其中,
,
,
和
是四个实参数。
表征孤立子的高度(波幅),
表征孤立子的速度,
和
表征孤立子的初始位置和初始位相。此时,对应的流函数为:
(24)
由于初始位相和初始位置仅能改变孤立子的中心位置,而不改变其振幅和结构,因此,我们考虑初始位置
,初始位相
时的流场,利用欧拉公式,重写流函数。注意到cc为共轭,为此,(24)式化为
(25)
其中,
,
。利用准地转流
与纬向速度u和经向速度v的关系,(25)式可化为如下的向量方程组:
(26)
为了获得方程组(26)的平衡点,设(26)式右边为零,则可解得
(27)
在波数m和纬向平均速度
和相速度
一定的情况下,方程(28)有两种情况:
(I) ![](//html.hanspub.org/file/2-2810120x93_hanspub.png)
此时,方程(27)有解,向量方程(26)有平衡点
,其中,
。对固定的k而言,其对应的Jacobian矩阵具有如下的形式:
![](//html.hanspub.org/file/2-2810120x96_hanspub.png)
根据阻塞形势的尺度以及平均西风纬向风速(罗德海1989),取
,
,此时,
,其对应的流场为如图1的偶极子型阻塞流场。
![](//html.hanspub.org/file/2-2810120x101_hanspub.png)
Figure 1. The blocking flow field of unheated single soliton solution
图1. 无加热单孤立子解的阻塞型流场
(II) ![](//html.hanspub.org/file/2-2810120x102_hanspub.png)
此时,方程(27)无解,向量方程(26)没有平衡点,但(26)式的第二式右边在
条件下恒为零,其在
处有退化的中心结构,对应的流场为如图2的近似纬向平均流场。
![](//html.hanspub.org/file/2-2810120x106_hanspub.png)
Figure 2. The approximate latitudinal flow field of unheated single soliton solution
图2. 无加热单孤立子解的近似纬向流场
从方程(27)可以看出,在波数m和纬向平均速度
一定的情况下,孤立子的波幅大小
决定着流场类型,而与孤立子的速度
没有关系。上面的分析可以看出,没有强迫耗散,仅在大气的非线性作用,会形成阻塞流场。由于大气运动是耗散过程,其耗散效应相对于大尺度运动是小项,为此,下面将耗散效应视为微扰项,分析微扰对孤立子波幅的影响,进而讨论非耗散效应对阻塞流场的影响。
5. 强迫耗散下Schrӧdinger孤立子的变化情况
由于
,故地形强迫在慢时空尺度
中仍可视为小项。从而,本节将以地形强迫项为微扰项,分析地形强迫对孤立子波幅的影响,为此,考虑NLSE的近似解:
(28)
其初始条件为:
(29)
为了消除久期项,引入慢时间变量变换
(30)
在此变换下,有
![](//html.hanspub.org/file/2-2810120x115_hanspub.png)
将(28)、(29)和(30)式带入(22)式,并线性化,得到如下具有初始扰动的线性方程:
(31)
其中,
为S的共轭。零阶近似解
满足标准NLS方程,即(23)式。
根据零阶近似解(23)式,设满足方程(22)的
具有以下形式的解
(32)
由(31)可知,
(33)
是关于慢时间变量
的函数,而
与时间t无关。为此,根据(33)式,可知
(34)
其中
![](//html.hanspub.org/file/2-2810120x126_hanspub.png)
由于方程(31)中含有S及其共轭
,为了更方便的运用变量z,将考虑关于S及其共轭
的方程组。为此,将(34)式带入到方程(31),得到
(35)
令
![](//html.hanspub.org/file/2-2810120x130_hanspub.png)
则方程(35)可以写成矩阵形式:
(36)
其中,
,称为第三Pauli (泡利)矩阵。而W和
定义如下:
![](//html.hanspub.org/file/2-2810120x134_hanspub.png)
![](//html.hanspub.org/file/2-2810120x135_hanspub.png)
为了求解方程(36),需要分析线性算子
及其伴随算子
的本征值问题,以构造正交完备化的本征函数集,即
(37)
其中
![](//html.hanspub.org/file/2-2810120x139_hanspub.png)
由于孤立子解
具有零边值条件
![](//html.hanspub.org/file/2-2810120x141_hanspub.png)
所以,可以用类似于反散射变换过程中求Jost解的方式来获得本征值问题(37)的本征值:(具体过程见 [24] [25] )
![](//html.hanspub.org/file/2-2810120x142_hanspub.png)
其中,k为实数,即算子方程(37)具有连续谱。对应的正交完备基
和
分别为:
![](//html.hanspub.org/file/2-2810120x145_hanspub.png)
其中
![](//html.hanspub.org/file/2-2810120x146_hanspub.png)
![](//html.hanspub.org/file/2-2810120x147_hanspub.png)
利用基函数
将W展开,得到方程:
(38)
其中,
,
和
为展开系数。
将展开式(38)带入到方程(37)中,由于积分后各展开系数只是时间t的函数,为此,得到关于展开系数
,
和
的常微分方程。通过展开系数的解,并考虑到(32)和(33)式,可得关于孤立子波幅M及速度
随慢变时间
变化的公式:
(39)
其中
![](//html.hanspub.org/file/2-2810120x159_hanspub.png)
在得到扰动对孤立子影响的关系式(39)后,下面将利用这一关系式来讨论地形强迫对偶极子型阻塞形势的影响。
6. 强迫耗散对阻塞形势的影响
为了进一步讨论强迫耗散对阻塞形势的影响,将耗散项项,即扰动表达式
带入到(39)式,利用留数定理所求的公式 [26]:
![](//html.hanspub.org/file/2-2810120x161_hanspub.png)
并考虑到矩阵次幂公式:
![](//html.hanspub.org/file/2-2810120x162_hanspub.png)
求(39)式,得到孤立子移动速度
对于慢变时间
的变化:
![](//html.hanspub.org/file/2-2810120x165_hanspub.png)
设
的初始值为
,注意到
关于慢变时间
是定常的,则可得到孤立子波幅M对于慢变时间
的变化:
. (40)
根据阻塞形势的尺度以及平均西风纬向风速的特征,取
,
,
,
,
[1],得到耗散作用下的孤立子波幅增长趋势图像(图3):
![](//html.hanspub.org/file/2-2810120x177_hanspub.png)
Figure 3. Suppresive effect of forced disspation on soliton wave amplitude
图3. 强迫耗散对孤立子型波幅的抑制效应
上面的分析表明:耗散作用对阻塞形势具有抑制作用,这种抑制作用随时间的负指数函数变化。
7. 结论与讨论
本研究运用多尺度变换和摄动法简化具有耗散效应的正压准地转涡度方程,得到带有扰动项的非线性Shrӧdinger方程,分析了非线性Shrӧdinger孤立子的拓扑结构,并采用孤立子直接微扰理论研究强迫耗散对阻塞结构的作用。结果表明:1) 在没有耗散的情况下,当基本纬向速度和孤立子波幅满足不同的条件时,定常Shrӧdinger孤立子具有两种形态:当基本纬向速度较小而孤立子波幅较大时,流场具有偶极子型的孤立波,表征高指数的经向环流;当基本纬向速度较大而孤立子波幅较小时,流场是退化的中心结构,表征低指数的纬向环流。2) 耗散效应对阻塞发展具有抑制作用,这种抑制作用随时间的负指数函数变化。
基金项目
国家自然科学资金项目41675095、41730959。