在立方球体坐标系下推导出浅水波方程的形式,使用M DQ方法和4阶R ung-K u tta方法求解该方程,计算了余弦钟绕球平流算例的三种情况,并且给出了该算例在立方球体坐标系下6个区间的速度矢量的表达式.为了解决由于边界插值所引起的数值不...在立方球体坐标系下推导出浅水波方程的形式,使用M DQ方法和4阶R ung-K u tta方法求解该方程,计算了余弦钟绕球平流算例的三种情况,并且给出了该算例在立方球体坐标系下6个区间的速度矢量的表达式.为了解决由于边界插值所引起的数值不稳定性,引入了4阶的人工粘性项(-4).通过对测试方案的求解和分析,可以看出该方法的优点:①完全显式格式;②时间上的4阶精度;③仅仅通过改变一个参数就可以获得空间的高阶精度;④该数值方案有较强的适用性.展开更多
文摘在立方球体坐标系下推导出浅水波方程的形式,使用M DQ方法和4阶R ung-K u tta方法求解该方程,计算了余弦钟绕球平流算例的三种情况,并且给出了该算例在立方球体坐标系下6个区间的速度矢量的表达式.为了解决由于边界插值所引起的数值不稳定性,引入了4阶的人工粘性项(-4).通过对测试方案的求解和分析,可以看出该方法的优点:①完全显式格式;②时间上的4阶精度;③仅仅通过改变一个参数就可以获得空间的高阶精度;④该数值方案有较强的适用性.
