半拉格朗日完全能量守恒计算格式地研究.pdfVIP

半拉格朗日完全能量守恒计算格式的研究 陈嘉滨 (中国科学院大气物理研究所．北京100029) 详细摘要 一、前言 目前半拉格朗日计算格式广泛应用于数值天气预报业务，提高了计算时效和计算精度。 但它存在的问题是原来偏微分方程的整体守恒性得不到保证．长时间数值积分会带来计算误 差。世界上使用的守恒格式，是指在时问保持连续情况下的守恒，典型格式是^r8kwa的守恒 格式。在时间离散情况下．这样的守恒格式却不能保证物理量的守恒。以上守恒性问题是局 限于欧拉空同。国外研究半拉格朗日空问的守恒格式尚不多见。本研究是在欧拉空问显式完 全平方守恒格式的基础上推广和发展。研究半拉格朗日空间守恒性问题。 二、一维浅水波方程的分析 欧拉空间完全平方守恒格式的关键．是方程的右端项的微分算子的离散空闯算子，具 有反对称性质．他们的平方变量即能量可写为辐散形式。在拉格朗日空间．方程也应具有 此性质才能推广。对一维重力波方程 芈d一絮 、…’ d* 警=一尹瓦3u (2．2) u是速度，母是位势。其能量方程为 舞(％z)；丢(血乒)：一麦(譬)一她≯3．。_u，u (：．3) 显然，上式右端项不可能构造出具有反对称性质的离散空间算子。 引入辐散的拉格郎日表示[3】 嚣=夏i差(如) (2．4) 其中缸表示随流体运动的无限小的长度同距。於是，在拉格朗日空间一维重力液方程 (2．1)一(2．2)可改写为 -著(UDx)—．m{乏(譬) (z．5) 署(纰)=一Dx号螽(u) (2．6) 其中 ≠=而，u=缸，m=瓶 方程(2．5)一(2．6)的总能量方程为 ·25· 墨(掣抵)=一溉￡(u譬)：．a(u￡2) (：．，， 方程(2．7)的左端项表示在无限小的空间缸内的能量。重力渡方程经过引入散度的拉格 ， 朗日表示式，其形状改变为方程(2．5)和(2．6)。他们的总能量方程(2．7)的右端项是辐散 形式。速表明方程组(2．5)一(2．6)右端的微分算子的离散空间算子是可构造成具有反对 ’ 称性质的。因此．我们就可以在拉梏朗日空闯应用酞拉空间完垒平方平方守恒格式。 兰、在拉格朗日空间的完全平方守恒格式构造 类似於欧拉空间显式完全平方平方守恒格式．在拉格朗日空间中．发展方程数值求解 问题可归结为如下形式 警+胛=G，1imF=一01(x) (3．1) o“ v．O 其中dF／是随俸微商；日sH(F，*．t)是一个非线性箅子。G是源汇项。在(，．f)空 闯．离敬化并布网(M，w)．其中h是空间步长，r为时间步长。于是可对方程(3．I)构造 相应的半拉格朗日平方守恒计算格式如下 竺∑旦+A一：∞ (3．2) 其中角码n和n+1分别表示现在和未来时伺层；下角码j表示标准两格点坐标，而*则 表示流体质点经过r时间到达标准隔格点j的出发点坐标。这个出发点坐标*一般不与 标准网格j重合。流体质点在出发点*的变量，是由标准