[计算机]经验正交函数在气候数值模拟中的应用.pdfVIP

目 种 墨 1991年4月 B 辑 第 4期 经验正交函数在气候数值模拟中的应用 张邦林 丑纪范 兰州大学大气科学系，兰州 730000) 摘 要 本文给出了一个缩减大气环流模式自由度的新方法，理论基础是强迫耗散非线 性系统的高维相空间将演化到一个低维的吸引子上．具体方法是对模武的一个现实 作经验正交函数(EOF)分解，从而找出支撑吸引子的少数 自由度，并 以此为基底获 得简化模型．用一个理论模型进行了数值模拟试验，试验证实了此方法的可行性和 有效性． 美键词：气候模拟，数值计算方法，经验正交函数 早在 50年代初，Kos~oropoB就指出：一般说来，在为了考虑某一现象而采取 的数学 模型中，研究者的技巧在于寻找一个非常简单的位相空间Q(即系统的各种可能状态 ∞的集 合)，使得当我们把实际过程换成点 ∞在这个空间中的因果式的变迁过程时，仍能抓住实际过 程的各 ／卜主要方面 u． 大气环流为一组偏微分方程所描述，系统具有无穷多个 自由度，为了进行数值计算，必须 进行离散化，用有限个参数(设为)来表征大气的状态，现在的大气全球模式，r．／已达 500000, 需要用现代世界上最大的计算机，而且n的数 目仍在继续增长．但是，值得注意的是 ，不论 n 多大，仍然是某种近似，在同等的近似程度下，如何选取合适的参数来描述系统的状态，使r．／最 大可能地减小，这样不仅节省计算量，还降低了对计算机的要求，使本来受机器能力限制而不 能做的问题也变得可以计算，显然是值得探讨的工作．本文从理论和数值试验两个方面对此 问题作了研究． 一 、 理 论 基 础 大气作为气候系统中的一个成员，是一个强迫耗散的非线性系统，郝柏林指出：“由于耗散 系统由高维相空间收缩到低维吸引子的演化，实际上是一个归并 自由度的过程，耗散消耗掉大 量小尺度的较快的运动模式，使决定系统长期行为的有效 自由度数 目减少，许多 自由度在演化 过程中成为 无‘关变量 ，最终剩下支撑起吸引子的少数 自由度．如果描述系统状态所选取的 宏观变量集合中，恰好包括了t一 ∞时起作用的自由度，那就会有一个比较成功的宏观描述 ． 本文 1990年 1月 17日收到，1990年 10月 1931收到修改稿 国家教委高等学拉博士点专硕科研基金资助课题 第4期 张邦林等：经验正交函数在气候数值模拟中的应用 443 如所知 ，大气动力 热‘力的偏微分方程组，可以写成 Hilbert空间中的算子方程 ，n +Ⅳ(og) +Pog： ， (1) f 。 这里 BP为正定 自伴算子，Ⅳ为反伴算子．利用算子的特征可以证明整体吸收集及其中不变 点集的存在性．随着时间f的不断增长，系统将愈来愈靠近某个不变点集，该点集反映了系统 的终态．从物理上讲，也就是系统向外源的适应，在相当广泛的条件下，可以证明该点集是有 限维的(与Folas和 TemanI证明Navier-Stokes方程的吸引子是有限维的，Ghidagliat证明粘 性不可压缩流动的各种方程的吸引子均具有有限维数的方法类似)；可以给出吸引子维数的 估计(用类似于Constanin等 估计非对流型湍流，Folas等 估计对流型方程的吸引子维数 的方法)．任何数值模式的状态变量都是有限的(设为月)，已经把问题 由 空阃转到 了 空 间了，转换中应保持算子的性质不变，即仍可以写为 空问的形如(1)式的算子方程．可以证 明在 中同样存在一个整体吸收的不变点集，且其测度为零，有理 由认为 中