自由流线问题有限元分析.PDFVIP

　 辑第 14 卷第 3 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 14, . 3 A Ser A V o l N o 　　1999 年 9 月　　　　　　　JOU RNAL O F H YDROD YNAM ICS　　　　　Sep. , 1999 ① 自由流线问题有限元分析 胡玉棠 (浙江省水利水电科学研究院, 杭州 310004) 林炳尧 (浙江省河口海岸研究所, 杭州 310016) 　　摘　要　在文 [7 ] [8 ] [9 ] 的基础上, 研究自由流线问题有限元解的通用计算方法。将各结点 的坐标, 代替流函数值, 定义为有限元方程的广义坐标, 使自由流线理论中对应于正、逆问题以及 空腔流动问题的有限元方程有相同形式。用 Galerk in 方法构成关于流速的变分方程, 解决了弱解 求导的困难。给出二种类型的空腔流动及机翼逆问题三个典型算例, 表明方法有较好的适应能力。 关键词　自由流线问题, 有限元, 通用算法, 溢流面通气槽, 叶栅, 空腔流 分类号　 351 O 1　引言 自由流线问题是流体力学一个重要课题。近年, 这一经典问题的研究获得很大进展。这问 题的主要困难在于: 一是自由流线的位置未知, 因此, 求解域待定; 二是自由流线上有压力及流 线两个边界条件, 在重力流动中, 问题总是非线性的。 自由流线问题有限元分析方法可以归纳为两种类型: ( 1) 　以自由流线上两个条件中的一个, 作为边界条件, 在初始域上, 求解固定域的 L ap lace 方程; 再根据求出的流动参量分布和另一个条件修正自由流线, 以此构成迭代关系; (2) 　利用可变域泛函表达式, 构造既能概括 L ap lace 方程, 又能满足两个定解条件的统 一变分原理。计算时, 将自由流线坐标和流动参量处理成广义坐标, 一并求解。 文 [7 ] 讨论了流函数模式统一变分原理。因为沿着自由流线流函数法向导数不为零, 引 用不同的贯截条件, 会得到不同形式的变分原理。事实上, 存在着四种流函数模式的变分原理 ( ) 表 1 。 各类变分原理有相同的基本泛函, 文 [8 ] 给出基本泛函的 H u ssu in 阵元素估计式 2 V 1, r i, rj 　≈ [ r g rad f ] 式中, f 为流动参数, 或者 ; r 为结点座标变分的允许方向。参数 [ r g rad f ] 越大, 有限 元方程解收敛越快。因为流速势函数在自由表面上的梯度平行于自由表面, 将自由表面的切向 取作为结点坐标的允许变分方向又没有意义, 因此, L uke [1 ] 给出的流速势函数模式变分原理 ① 本文于 1998 年 6 月 22 日收到。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 352 水　动　力　学　研　究　与　进　展　　　　　　　　　1999 年第 3 期 很难用于有限元计算[8 ]。 表 1　自由面有势流的各种变分原理 模式 类别 变分方程