求解河网一维N—S非线性方程组的一种简化延拓法.pdfVIP

第8期 广东水利水电 NJ【．8 2015年8月 CUANGD0NC WAIIFR RES0URCESAND HYDROPOWER Aug．2015 求解河网一维 N—S非线性方程组的一种简化延拓法 陈雪菲，周 斌 (汕尾市水利水电规划设计院，广东 汕尾 516600) 摘 要：随着计算能力的进步，河网中一维N—s方程应用越来越多，但数学上常用于求解N—s非线性方程组的拟牛顿 法收敛域小，不能保证迭代收敛，对工程应用带来了一定的困扰。该文提出了一种仪进行边界控制的简化延拓法，可大 幅拓展非线性方程组迭代求解的收敛域，同时能保持算法的简洁和较低的计算开销，经实例验证，方法可行。 关键词：河网；N—S方程求解；非线性方程组；延拓法 中图分类号：TV131．4 文献标识码 ：A 文章编号：1008一O112(2015)08—0019—03 在水利工程设计 中，常需进行河网的水力计算， 水力坡降。 河网的水力计算大多可归结为一维非恒定流问题，需 一 维 Ⅳ_一S方程常采用差分法求解。由于显式差 要建立一维 N—S方程组求解 。一维 N—S方程组具有 分格式的稳定性较差，一般工程常采用隐式差分格式 非线性的特点，数学上常将其拟牛顿化为线性方程组 求解 ，常见 的一 维 ，v—s方程 隐式差分格 式有 逐步逼近迭代求解 ¨I2J。由于迭代求解的收敛域有限， Preissmann格式和 66o一，o cM格式 。 而边界条件 的变化可能过于剧烈 ，导致迭代发散现象 2 Preissmann和 Abbott--Ionescu差分格式及其求解 时有发生，迫使计算 中断。鉴于河网N—s方程组发 Preissmann差分格式为隐式差分法，具有时间步 散主要发生在边界处，笔者在河网计算 中尝试仅在边 长大、数值稳定性好的优点，该格式的表达式为 ： 界处引入延拓法 来扩 大收敛 域，取得 了较好 的 一 Of 二 ： 二 二 r3、 效果 。 1 一维N—S方程及求解 (4) I=)t 2At N—S方程是流体力学 中描述粘性牛顿流体的基本 其 中系数和阻力项表达式为： 方程 ，形成于 19世纪。由于N—S方程是数学中最为 难解 的一类非线性方程 ，寻求精确解十分困难 ，致使 f(x，，t)： n~l+f )，+了 Ijn+l+f) (5) N—s方程产生初期工程运用不多。2O世纪60年代 以 式中 j为计算节点号；n为计算时段；0为变权系 来，随着计算机技术的飞速发展，N—S方程开始逐步 数；Ax为分段长度 ；At为分段时长；f为函数。 应用于工程建设 中，计算水力学也开始成为一 门独立 Abbott--Ionescu差分格式也是隐式格式，它将水位 的学科得到了飞速 的发展 。在我 国很多地区存在着河 和流量分别错开布置在不同的网格节点上，对于连续 网，河道相互交错，常规水力学较难求解，目前工程 方程(1)式和运动方程 (2)式采用了不同的近似格式， 界运用一维N—s方程研究河网问题越来越普遍 。 其 中连续方程使用 ¨： 河网一维 N—S方程…为 ： _y上— (6) + ：0 (11) Ot b