蒙特卡罗法在势流计算中的应用研究 - 水科学进展.pdfVIP

第 15 卷 第 4 期 水 科 学 进 展 Vol15 ,No4 2004 年 7 月 　　 ADVANCES IN WATER SCIENCE 　 Jul , 2004 　 蒙特卡罗法在势流计算中的应用研究 1 2 吉庆丰 , 郑邦民 ( 1 扬州大学水利科学与工程学院 , 江苏 扬州 　225009 ; 2 武汉大学水利水电学院 , 湖北 武汉 　430072) 摘要 : 针对水动力学实际问题多存在复杂几何边界的状况 , 提出了用不规则游动网格求解偏微分方程的蒙特卡罗 法 , 建立了相应的随机游动模型 。选择具有复杂自由面的堰闸流动问题作为算例 , 验证了新方法的正确性 。与有限 元法相比 , 蒙特卡罗法解势流等线性问题时更灵活 , 可以根据需要 , 单独计算流动区域内任意一点的未知物理量 , 且所用计算容量较少 。 关 　键 　词 : 蒙特卡罗法 ; 不规则游动网格 ; 随机游动模型 ; 复杂边界 ; 势流 ( ) 中图分类号 : TV1312 　　　文献标识码 : A 　　　文章编号 : 2004 蒙特卡罗法不仅可以求解随机性问题 , 而且可以求解诸如非线性方程 、多重积分 、常微分方程和偏微分方 程等多种确定性问题[ 1 ,2 ] 。与有限差分法 、有限元法等常用数值方法相比, 蒙特卡罗法在求解确定性问题时虽 显得有些 “笨”, 但仍有显著的优点 , 比如对线性问题 , 可以单独就计算区域内任意一点作计算 , 其值无需通 过联解其它结点值求出; 计算稳定性较好 ; 便于并行计算等 。因此 , 对蒙特卡罗法进行研究具有一定的理论和 实际应用价值 。 以前用蒙特卡罗法求解水动力学问题时 , 采用的游动网格均为规则网格 , 如矩形网格 , 对不规则边界问 题 , 由于边界的近似处理 , 计算精度受到了很大的影响 。为了克服这一缺点 , 本文提出了不规则游动网格的蒙 特卡罗法 。这种方法改善了蒙特卡罗法求解复杂边界问题的精度 , 拓宽了蒙特卡罗法在水动力学中的应用范 围。本文以理论上较成熟 、计算亦较简单的势流问题为例 , 介绍不规则游动网格蒙特卡罗法的具体做法 , 以便 突出方法本身的优点 。 1 　随机游动模型的建立 蒙特卡罗法的关键在于建立一个正确可行的随机游动模型。为简便起见 , 本文以二维 Laplace 方程的边值 问题为例 , 给出不规则游动网格的随机游动模型 。方程及边界条件为 ΔΦ( )