ght=19 src="http://www.studa.net/Newspic/2006331/1040388725.gif" width=25 v:shapes="_x0000_i1037"> 的外法向单位向量。边界 由四条线段组成，上式左端第2项可以写成

                                           (4)

其中 为边长， 为第k边外法向通量。

    设向量U在单元内部保持不变，进一步离散(3)式，得到有限体积基本方程

                                        (5)

其中 满足

                                       (6)

    利用F(U)和G(U)的旋转不变性，引入旋转变换矩阵 及其逆矩阵 可改写（5）式为

                                (7)

其中

            ( 8)

    记(7)式右端为 ，则有

                                   (9)

采用两步Runge-Kutta法离散上式，时间精度也可以达到二阶，得到

                        (10)

    在每一单元的每一边，其通量按组合型TVD格式给出(例如对于单元i的第1边)

                                     (11)

式中 为单元i和卫星单元1的第l个右特征向量的算术平均。采用组合型TVD格式进行插值，

                         (12)

其中 代表单元i和卫星单元1的第l个特征速度的算术平均， 为二单元特征变量差分的算术平均， 为限量函数，采用文献[8]优选得到的MUSCL型限量函数，

             (13)

其中0.138正是计算溃坝波时，介于缓流（亚临界流）与急流（超临界流）之间的临界水深比。当初始下游与上游的水深比较大（湿底）时，宜采用单参数限量函数 ，否则，在干底时宜采用双参数限量函数 以避免在坝址附近的非物理扭曲现象。

                (14)

                  (15)         

其中

本文引用地址：http://www.uscc.com.cn/Discourse/Demonstrate/Science/Engineering/200609/Discourse_5099.htm

[1] [2] [3] [4] [5]

责任编辑：