计算复杂边界渗流的蒙特卡罗方法.pdf

第 17 卷第 4 期 计 算 力 学 学 报 V o l. 17 N o. 4 2000 年 11 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS N ovem ber 2000 文章编号:(2000) 计算复杂边界渗流的蒙特卡罗方法 1 1 2 吉庆丰 , 　　刘　超 , 　　郑邦民 ( 1 扬州大学, 江苏 扬州 225009; 2 武汉水利电力大学, 湖北 武汉 430072) 摘　要: 本文运用不规则游动网格的蒙特卡罗方法求解复杂边界渗流问题。与以往各种数值方法相 比, 蒙特卡罗方法在计算渗流问题时要更灵活、更方便, 可以根据需要独立地计算出渗流区域内任意 一点的渗透压力、渗流流速。 关键词: 渗流; 数值计算; 蒙特卡罗方法; 随机游动模型 中图分类号: O 242 2; O 3573　　文献标识码: A 1　前　言 渗流问题是水利工程常见的一种问题, 实际工程中, 人们常常要对渗流进行计算, 以求得 需要部位的渗透压力、渗流流速。渗流流场计算的方法很多, 简单问题常用经验法、流网法计 算, 较复杂问题常用模拟试验和数值计算来进行。数值计算的具体算法也很多, 用得较多的是 有限元法。有限元法可以较好地解决复杂边界问题、非均匀渗流问题, 但在网格结点很多的情 况下, 其计算工作量还是很大的。 蒙特卡罗方法是一种具有独特风格的数值计算方法, 它的基本思想是: 为了求解问题, 首 先建立一个随机游动模型, 使它的参数等于问题的解, 然后通过对模型的观察或抽样试验来计 算所求参数的统计特征, 最后给出所求解的近似值[ 1, 2 ]。与其它数值方法相比较, 蒙特卡罗方法 ( ) 具有如下一些显著的优点: 1 对线性问题, 可以单独就计算区域内任意一点作计算, 其值无需 ( ) ( ) 通过联解其它结点值求出, 计算单点的工作量大大减小。 2 计算结果形象直观。 3 在有自由 ( ) 面的流动问题计算中, 自由面的调整确定变得更简单更灵活。 4 由于随机游动在平面上进行 和在空间上进行无根本差别, 方法易于推广到三维问题。 许多情况下, 渗流问题只需计算局部区域的物理量, 比如闸基渗流问题, 人们常常只对闸 底板上的渗透压力、闸基下游出口边界处的渗流流速感兴趣。蒙特卡罗方法很灵活、很方便, 可 以根据需要独立地计算出渗流区域内任意一点的物理量。因此, 运用蒙特卡罗方法求解渗流问 题应该是非常适合的。本文运用不规则游动网格的蒙特卡罗方法, 讨论复杂边界渗流问题的计 算, 给出了具体计算的随机游动模型及算例。 2　渗流问题的数学模型 不可压缩稳定渗流的微分方程式是 H H H (kx ) + (ky ) + (k z ) = 0 ( 1) x x x y z z 收稿日期: 　　基金项目: 扬州大学水利动力工程省重点实验室开放课题资助 ( ) 　　作者简介: 吉庆丰 1964～ , 男, 博士, 副教授 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All r