适用于非静力大气模式的完全平衡多矩约束有限体积方法

大气数值模式离散垂直动量方程时,一般而言气压梯度力和重力不易维持严格的静力平衡关系。为精确平衡数值离散的垂直气压梯度力和重力,基于高精度多矩约束有限体积方法引入完全平衡数值公式,即以满足静力平衡关系的热力学参考态对重力源项进行数值离散构造,发展了适用于非静力大气的完全平衡多矩约束有限体积方法。一维标准数值试验表明,完全平衡多矩约束有限体积方法能在较粗糙的计算网格点上保持静力平衡参考态的数值计算误差在计算机的单精度（10-6）和双精度（10-14）水平,在具有小量级扰动的初始条件下,完全平衡多矩约束有限体积方法能较好地模拟扰动的传播,二维非静力热泡试验进一步验证了完全平衡多矩约束有限体积方法对非静力大气运动的模拟能力。数值试验结果验证了所发展方法的完全平衡属性和适用性,这为非静力大气模式发展提供了良好参考价值。