明渠非恒定流问题的数值计算方法.DOCVIP

明渠非恒定流问题的数值计算方法 ???? 关键词：有限元；微分方程1 特征线法 特征线法由Aabott M．B．于 1966 年提出， 其基本思路是先将圣·维南基本方程组（2-1）（2-2）通过线性组合变换为沿特征线的常微分方程组，即顺、逆特征线方程和顺、逆特征方程；然后对常微分方程组结合初始条件和边界条件进行数值计算，而要求的圣·维南方程组的离散解就是求顺、逆特征线的交点。此法物理概念明确，数学分析严谨，计算精度较高。 2 直接差分法 在流体力学数值解法中，广泛应用的是有限差分法，其基本思路是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，以Taylor级数展开等方法，将控制方程中的微商用差商代替进行离散，再建立代数方程组来求解。即在离散的网格点上把各偏导数项化为差商来求数值解：选择合适的基本方程，确定相应的定解条件，然后将微分方程离散到差分网格上求解。差分方法的基本要求是差分方程能准确地代表微分方程。有限差分法的近似处理会引进误差，包括差分方程的截断误差、差分方程解的离散误差，以及计算过程中引入的舍入误差。这些误差对计算过程和计算结果的影响，反映在差分方程的数学特性上。与此相应的是差分方程的相容性、收敛性和稳定性。 由于采用不同的偏差商代替偏导数，基本微分方程可演变为多种形式的差分方程，显式中差分方程往往是一组线性代数方程。而隐式中差分方程往往是一组非线性代数方程，需迭代求解。常用的显式差分格式有扩散格式和菱形格式（蛙步法）（Leap-frog scheme）。常用的隐式差分格式有Crank-Nicholson 格式和 普莱斯曼（Preissmann ）四点偏心格式等。其中 Preissmann 四点隐式格式应是差分格式中较为广泛的。大量的非恒定流模拟软件（如RootCanal、CanalMan、CARIMA、MODIS2 和 DUFLOW，MYBC/Y1 等）以及河（渠）网非恒定流方程数值离散方法基本上均使用该格式。 3 有限元法 有限元法主要用来求解二元明渠非恒定流问题。主要针对边界条件复杂，水面宽阔的二元明渠水流如河口、海洋的潮汐等水力计算，由于河流或渠道中的非恒定流一般属一元流问题，采用有限元求解并无明显优点，故求解仍以特征线法、差分法为主。在这里不再具体讲述有限元法。 4 高阶浅水理论的数值计算方法 Boussinesq方程是描述浅水中具有色散性和非线性的波动的方程，其色散性介于无色散的浅水非线性长波和Stokes色散波之间。 非线性Boussinesq方程对波浪传播模拟的研究，最早是Peregrine(1967)采用有限差分法进行的，为以后应用Boussinesq方程开辟了道路。其后的研究大都是在改进模式的方程上进行的，smallmanandcooper(1989)采用预报一校正法模拟了Set一Down波的产生和传播用差分法求解微分方程时，差分方程的精度至少要达到微分方程中微商的最高阶数。在Boussinesq方程中最高阶导数为三阶，则该方程中各项差分的精度都必须在三阶以上，如果方程中某些项只达到二阶精度，则会产生数值色散，并将干扰原方程的物理色散性。因而在Boussinesq方程的求解中，设计高精度的格式十分重要，特别是影响解的非线性的对流项。 5 用紧致各式构造差分方程 紧致格式是由数学中的三次样条理论发展而来的,设u（x）为一个函数，其在网格点xi上的函数值为u（xi）= ui ,三次样条函数SΔ（u,x）= SΔ（x）在区间[a,b]上满足一阶导数和二阶导数连续，同时在每个小区间xi-1≤x≤xi中可以写成三次多项式，并满足SΔ（ui,xi）= ui,函数u（x）属于C4[a,b]。设mi=S?驻(xi)，mi=S?驻(xi)，通过推导得出离散点上一阶导数m和二阶导数M所满足的两个方程，也是构造紧致差分各式的两个基本关系式（取空间点为等距分布）： 式中:?驻x为结点的空间间距。利用上两式我们就可以构造高精度的差分格式。 6 预报-校正法构造差分方程 预报-校正格式(PrediCtor-CorrectorMethod)的基本思想是，从时刻tn推进到时刻tn 1分成两步：首先用稳定性较好的一阶精度格式计算出tn 1时刻的近似值，然后在时间区(tn，tn 1)上用二阶精度格式计算tn 1时刻的近似值。预报-校正格式可以采用全隐格式，非线性精度较高，且格式比较稳定。但该法要采用迭代求解，计算时间较长；当计算区域较大时，该法计算量过大。 7 ADI法构造差分方程 交替方向法是1955年PeacemanRachford提出的交替方向隐式法(AlternatingDireetionImplieitMethod)，简称ADI法。这是一种显式和隐式交替使用的方法，在求解多维抛物线