非线性兴波阻力计算.pdfVIP

／ 一一 第 20卷 第 5期 武 汉 交 通 科 技 走 学 学 报 Vo1．20 No 5 l996年 10月 Jaurna10fW uhanTransportationUniversity October1996 非线性兴波阻力计算 ． 连 迎 高 (武汉变通科技大学船舶及海洋丁程系 武汉 430063) 摘要 ：叙述 计算船舶非线性兴 波阻力的一种方法 ，采用面元 法满足完全非线性边条 的势流场 ．特点是采用三角形单元 ，蔼【面 自自由面上置和 自湿水面下漕 ，每戎迭代均适 应新的波面重新作网格划分 ，所采用的选代方式未发生不收敛同题 ，只需少数迭代戎数 即达到收敛 ．对Todd60系列 一0．6～0．75及不同L／B，B／丁作 计算 ，所得计算结 笪 美键词：l巳圭呈生卷丝}墼焦地 桶1；橙＼们轴Ⅱ 中圈法分类号 ：U661．311 甩理论计算决定船舶兴渡阻力大小，讨论船型参数对其影响 ，一直是造船工作者梦寐 以求 的 ．近 20a计算机迅猛发展推动了计算流体力学的发展，也促进 了兴波阻力计算 日趋成熟 ， 其问起了重要作用的首推 Dawson方法_1]，Dawson方法的重要性在于 ：其一是采用 了Rank iile源取代了满足 自由面线性边条的Kelvin源 ，这一基本解的变动，使得 向非线性推广更为可 能；其二是采取单边差分来满足辐射条件，这一变化也是摒弃了Kelvin源后不得 已采取的措 施，但也正是 由于这一措施并取得成功，才得 以摆脱线性满足辐射条件的约束，此外 ，甩简单的 Rankine源无需作复杂的无穷振荡积分，计算工作量大减，也是其一大优点． 受到Dawson方法计算结果令人满意 的鼓励，各 目竞相编制商用程序，如荷兰的 DAW一 ， 已工作了近 10a，据报导每年应用于约 60条商船，取得 了大量经验 ．1986年以来各目研 究者又开始推广到非线性问题 ，首先是SsPA一系列博士论文 ，如Xia口，Ni口j，Kiml{j，然后是 Jensencs]，Musker[6j，Kin和 Lucc~[ ，Raven[等人．非线性 问题 需要在事先未知波面 (实际上 是解的一部分)上满足 自由面边条 ，因此H能用迭代方法 ．而收敛 问题则是迭代方法的焦点 ． 为了解决解的收敛和稳定问题， 采用了高阶面元 J，Ki1311和Lueos提出了一参数族差分格 式 j，在试算大方形系数 ( =O8)时又提出内外层迭代方式 ，Jensen则用网格错位法满足 辐射条件 ，并甩Taylor展开近似计算新的 自由面位置 j．Reven采用了 自由面上置源并详细 讨论 了解法各步骤is：．总之 ．在各方面做了大量工作 ，不过看来 目前还没有达到共识 ，有一些做 法还需要深人研究认识． 非线性兴波阻力计算 的基本要求是 ：1)自由面边界条件要包括完整的非线性项 ；2) 自由 收稿 日期 ：1996—08—27 李世谟 ：男．72岁 ，教授 *国家 自然科学蕈盒资助 ，批准号 ·520 · 武汉交通科技大学学报 1996年 第 20卷 面边界条件要在真实的波面上被满足；3)物面应当是真实波面下的湿水面． 本文是在满足此基本要求下进行计算的 ．其它细节有 ：采取边界积分方法；以Rankine源 为基本解；源面从 自由面上置一定距离；物面源分布 自物面下潜于物面之 内；所有源面取三角 形单元 ，避免单元之间有阊隙，源强作线性分布 ；采用网格错位法以满足辐射条件 ，不借助于差 分格式；影响函数用解析表达，不 人【人为因素． 阻力系数计算结果与 1o、dd6O系列试验值 比较，方形系数 Cb一0．6～O．75，弗氏数n = 0．18--0．34范围内均 令人满意 ．限于硬件条件 ，网格数较少，也未作网格数影响的讨论 -文 末指 出前人未及