浅水环流问题的有限元分析

考虑涡流粘性,底部摩擦,Coriolis力和重力作用,以单层模型的竖直向积分建立二维流动的环流控制方程。用Galerkin加权余量法建立环流边值问题的弱形式。用有限元法和分裂时间法分别离散空间和时间进行数值近似计算。采用线性内插的“人工光滑”处理以消除短波长噪音干扰,而避免已有的大阻尼系数方法等缺陷。为节省计算内存容量,程序中建立稀疏矩阵的一维紧凑存贮格式,排除了全部“零”元素的容量要求。