Bankine源高阶边界元法求解势流问题.pdfVIP

　 辑第 14 卷第 1 期　　　　　 水 动 力 学 研 究 与 进 展 　　　 . , . 14, . 1 A Ser A V o l N o 　　1999 年 3 月　　　　　　　 　　　　　 . , 1999 JOU RN A L O F H YD ROD YN AM IC S M ar B ank in e 源高阶边界元法 求解势流问题 李谊乐　　刘应中　　缪国平 (上海交通大学船舶及海洋工程系, 上海 200030) 　　摘　要　本文以简单 Green 函数 1r 为基本解, 采用高阶边界元求解势流问题。利用 8 结点等 参元来逼近流域边界的几何形状和未知量, 利用变形的边界积分方程和坐标变换技术有效地处理 了在场点的立体角计算和数值积分中所谓的“奇异性”问题, 得到足够精度的结果。对零航速的半 球振荡问题的数值解表明, 1r 源高阶边界元法能给出满意的结果, 且消除了“不规则频率”现象。 关键词　高阶边界元法, 数值积分, R ank in e 源, 奇异性, 不规则频率 分类号　 353. 2 O 1　引言 自从H e ss Sm ith ( 1964) 首先利用面源法求解三维物体的无升力绕流问题, 三十多年 ( ) 来边界元方法 包括面源法和直接法 在流体力学、固体力学、传热等方面得到了广泛的应用。 常数面元法(CPM ) 应用最广, 但发现存在以下不足: 物理量分布间断, 需消耗较大的内存, 在 ( ) 需要边界线处数值的问题中 物体水波相互作用问题中水线处速度、裂缝处的应力等 无法得 到精确的数值。为克服CPM 的缺点, 又引入高阶边界元(HOB EM ) 的概念, 在单元上假定物理 量以二次或更高阶方式分布, 从而更精确地描述了边界几何形状和未知量分布, 并能有效地求 得边界线处的物理量。在波浪结构相互作用计算中, 刘勇辉用HOB EM 处理了流固耦合问题 ( ) ( ) ( ) 1988a , 并与CPM 进行了系统的比较 199 1 ; X u Yu e 1992 用九结点HOB EM 模拟了 三维全非线性波的翻卷; Boo K im ( 1995) 计算了非线性波与结构的相互作用。在耐波性计算 中,N ako s Sclavouno s ( 1990, 1993) 利用B 样条HOB EM 计算了 0. 25 时的频域和时域 解; H on g Ch io ( 1994) 利用B 样条八结点HOB EM 计算了 0. 25 和 0. 25 的频域解; Iw a sh ita ( 1996) 利用三次样条HOB EM 计算了双体船的水动力系数。 对于带自由表面的线性势流问题, 可以选取R ank in e Sou rce 或满足 自由面条件的Green ( ) 函数 刘应中, 1987 作为基本解。后者自动满足线性的 自由面条件和远方 N eumm an K elv in 辐射条件, 但其数值计算复杂, 且仅适用于简单的自由面条件, 应用上有很大的局限性, 刘勇辉 和 Iw a sh ita 等采用的是这种做法。R ank in e Sou rce 分布在所有边界, 能适合任何复杂的线性 或非线性边界条件, 计算简单。