贴体人坐标变换技术在渗流计算中的应用.pdfVIP

2007 2007年5月 西安石油大学学报(自然科学版) May ofXi’an Science Vd．22No．sl 第22卷第s1期 Journal ShiyouUniversity(NaturalEdition) 文章编号：1673—0“X(2007)s1—0145—03 贴体坐标变换技术在渗流计算中的应用 赵鸿，刘玉玲 (西安理工大学水利水电学院，陕西西安710048) 摘要：将贴体坐标变换技术应用于闸坝地基渗流数值计算中，建立了曲线坐标系下闸坝地基渗流的 数学模型，给出了边界条件和求解方法．计算实例表明，基于贴体坐标变换的有限差分法的计算结 果精度与有限元法相当，但计算效率较有限元方法有很大的提高．该方法具有简单、快捷的优点， 是求解复杂边界渗流问题的一种有效方法． 关键词：曲线贴体坐标；有限差分法；渗流 中图分类号：们玛19 文献标识码：A 地基渗流控制不当是导致闸坝发生渗流事故和 破坏的重要原因，对闸坝地基进行渗流计算与分析 k磬坞势=0， (1) 式中，H为水头；也，忌。分别是z，Y方向的渗透系 具有十分重要的实际意义．近年来关于闸坝的渗流 数．水头H与流速的关系为 计算用得较多的是有限元法，如：黄蔚等[1]对求解 ，aH ．aH ，．、 含自由面的三维无压渗流场的边界不确定的非线性 2一定y Lz， “2一意z百j，u ji， 问题，在传导矩阵的基础上，提出了单元的有限元算 式中，U，础分别为z，Y方向的流速． 法．黄春娥等[2]用条分法与有限元法相结合分析边 方程(1)中只有一个未知数，给定适当的边界 坡的稳定性，能适应边界条件较为复杂的渗流场．本 条件就有定解． 文采用贴体坐标变换技术模拟复杂的渗流区域，将 I．2 贴体坐标变换技术【4] 复杂的渗流区域变换成规则的矩形区域，在矩形区 在笛氏坐标系下使用矩形网格划分计算区域会 域上采用有限差分法数值离散求解变换后的控制方 将复杂边界概化为锯齿形边界，若使用差分法计算 程，运用追赶法求解方程，最终实现对整个渗流场的 时有可能出现计算奇点，边界附近数值解误差较大． 流态计算．本文的方法具有有限元法处理复杂计算 因此本文采用贴体坐标变换技术，将物理计算域 区域方便的优点，同时又克服了有限元法计算工作 (z，Y)变换成矩形的计算区域(拿，刁)．在矩形计算 量大的缺陷，是求解闸坝地基渗流问题较为理想的 区域上求解变换后的控制方程，这样可以避免对边 一种方法． 界条件的插值，使计算精度提高．坐标变换方程为 1数学模型