大气强迫耗散的非线性动力学.pptVIP

ψ*的表达式是类似的，即 (7.25) (7.26) 值得注意的是，上述选择与y＝0，πL处的刚壁条件?Fi/?x＝0有关。矩形区域x是从零到2π，y是从零到π，所以对整个区域积分，这样可得常微分方程组如下： (7.27) (7.32) 其中， 7.3.2.1 仅考虑地形影响 考虑ψ*和K等于零的守恒情况，如果ψC，ψM和ψN初始为零，以后它们将一直保持为零[由(7.29)、(7.30)和(7.31)式可以看出]，这样方程(7.27)-(7.32)可以化为 (7.33) 它们的定常解(即平衡态)，为 ψL＝0，ψA＝任意值，ψK＝hn1ψA/bn1 其中 方程组(7.33)的线性化小扰动方程为 (7.34) 对(7.34)第三个方程求微分，并将(7.34)第一、二个方程式代入，即得 (7.35) 其特征方程为 (7.36) (7.37) 如果 和 ，当 趋向于零时，λ→∞，其解eλt也趋于无穷大， 变成共振增幅。 7.3.2.2 强迫耗散系统中的多平衡态(只考虑模态1) 引进纬向驱动和耗散，保留ψA*，其它项中的ψ*为零。ψA*是驱动Rossby数 ，U0为纬向风的振幅。同样考虑ψC，ψM，ψN＝0(初始为零，以后保持为零)，这样原方程组(7.27)-(7.32)改写为 (7.38) 平衡态情况下，有 (7.39) (7.39)第一个方程式和第二个方程式相等，即可求解三次方： (ψA?ψA*)[(αmψA?βn1)2+K2]+h01hn1ψA=0 (7.40) 在图7.3中，在ψA*，ψL/hn1平面上的曲线表示(7.39)第二个方程式和直线表示(7．39)第一个方程式，其中K＝10-2，L/α＝l/4，n＝2，ψA*＝0.05，0.1，0.2和h0/H＝0.05，0.1，0.2。直线与曲线的交点表示定常解，即平衡态。 由图易见，对于合适的参数值，ψA有三个根，两个在靠近共振值ψA=βn1/αn1附近，第三个靠近ψA*。这种平衡态(即ψC=ψM=ψN=0)的稳定性，可由方程组(7.27)-(7.32)的一次近似方程级的特征根的实部的符号情况来判断，即如下六阶行列式为 这一六阶的特征方程可以分为二个三次方程，正好对应于第一和第二个y模态。这是二个σ的三次方程，分别称为第一波型和第二波型，可用数值方法算出其特征根实部的最大值，如为负，则所有特征根的实部均为负，解是稳定的；如为正，则是不稳定的。表7.1是取K＝10-2，L/α＝1/4，n＝2和h0/H＝0.2所算得的结果。 由表可见，当只有一个平衡态，一般是稳定的；当有三个平衡态时，ψA的中间那个值对第一波型的扰动是不稳定的，其余二个是稳定的；但ψA*充分大时，ψA的最小的值对第二波型的扰动变得不稳定。于是存在这样一个区域，当ψA*取的值在这个区域中时，存在两个稳定的平衡态。就表7.1所取的K、L/α、n和h0/H的值而言，ψA*＝0.2就是这个区域中的值。这两个稳定的平衡态中，一个是“低指数”，有着较强的波动分量和较弱的接近线性共振的纬圈气流；另一个是“高指数”，有着较弱的波动分量和相当强的远离线性共振的纬圈气流。 为进一步证实这个结论， Charney对方程(7.21)用差分方法进行数值计算，对同一区域取了16×16个格点对ψA*＝0.2，h0/H＝0.2，L/α＝1/4，n＝2，K＝10-2的高截断谱模式有两个稳定的平衡态，(ψA，ψK，ψL)＝(0.154，0.0377，-0.0019)和(0.0294，-0.0309，-0.0071)，分别用它作为初值计算，计算了750时步(约一个月)发现都已充分接近定常状态，再将格点上的值分解成谱分量得(ψA，ψK，ψL)为(0.153，0.0399，-0.0026)和(0.0290，-0.0312，-0.0088)，这与谱模式相当一致。由此Charney认为阻塞形势是一种准稳定的低指数接近于共振特性的平衡态，是热力、地形强迫、耗散和非线性三者共同作用所产生的。 7.3.2.3 考虑完整的模态 对给定的初始状态，方程组(7.27)-(7.32)规定了其随时间演变的情况，不过我们的兴趣现在并不在这种演变的详细细节，而是要知道当时间足够长后所达到的与强迫的外源相一致的动力学状态是怎样的；我们想要知道的还不只是某些初值出发的这种状态，而是所有可能的初值出发的这种状态。我们建立一个六维的相空间(ψA，ψK，ψL，ψC，ψM，ψN)，要确定的是方程组(7.27)-(7.32)在相空间中的所有吸引子及其吸引域，如果这个数学问题得到解决，