二维物体入水数值模拟.pdfVIP

华 中科 技 大学 学 报 自然科 学 版. .第 卷 第 期 正 .. . . . 月 二维物体入水数值模拟 李云波 李亚军 胡始宏 哈尔滨工程大学 多体船技术实验室; 船舶工程学院,黑龙江哈尔滨 摘要 基于势流理论和边界元法模拟了二维剖面以恒定速度入水问题,采用双点模型处理自由面物面交点、 阶龙格一库塔法对自由面进行迭代求解,并利用自由面平滑技术和网格重新划分技术阻止自由面不稳定现象, 根据射流区截断模型处理射流区,采用辅助函数法求解物面压强分布以保证数值解的准确性和稳定性.探讨 了网格划分对数值结果的影响以及剖面斜升角对最大压强值以及垂向力的影响.预报了不同斜升角的 型 剖面的自由面形状和剖面压强分布、最大压强值以及垂向力,并与相似解进行了比较.比较结果验证了本文方 法的有效性和准确性. 关键词 势流理论;边界元法;二维剖面;剖面入水;压强分布;垂向力 中图分类号 . 文献标志码文章编号? ? : , ,,? .? , . ? ? .? .~. ? . . ?,, ?..;;? ; ;; 物体入水砰击问题有着广泛且重要的应用, 算剖面压强,边界元法获得的结果和相似解符合 其对高速船结构设计和水动力性能都有显著的影 良好.文献 基于柯西定理模拟该问题,所得结 响.二维剖面入水模型在模拟甲板砰击以及甲板 果除射流区外与相似解符合良好.文献 对二维 上浪中有着广泛的应用,研究二维剖面入水数值 对称 型剖面自由入水问题进行了实验研究,探 模拟有着重要意义.文献应用边界元方法模拟 讨了物体质量、初始落水高度以及斜升角对压力 该问题,用截断模型来处理射流区并用差分法计 系数的影响.文献 ? 数值模拟了二维 型剖 收稿日期? . 作者简介 李云波 一 ,女,教授,: .. 基金项目 国家自然科学基金资助项目 .第 期 李云波,等:二维物体入水数值模拟 面自由入水,其采用相似解作为初始条件,所得数 为计算域边界.当场点在 上时,则式为 值结果与实验结果符合较好.文献测量了二维 叩型剖面恒定速度入水时的自由面升高以及物体 令 ? ,则有 / , 由下式求得: 所受压力.本研究基于势流理论应用边界元法模 一 拟了二维剖面恒定入水问题,采用双点模型、 阶 龙格一库塔法、射流区截断模型、自由面平滑技术 将物面 自由面 、远方边界面 以及 和网格重新划分技术对问题进行研究,以保证数 底部 分别划分为直线单元并分布线性源,则 值解的准确性和稳定性.预报了剖面自由面形状 方程可离散为线性方程组,解之可求得物面上 和压强分布以及最大压强和垂向力,并与相似解 的速度势及自由面上的速度势导数.求得物面速 进行了比较.数值模拟结果表明数值结果与相似 度势后,物面压强 由伯努利方程可得 解符合较好. ? 一一 ÷, \ 一 , 数值模拟 式中: 为大气压强; 为水的密度. . 数值方法 . 定解条件及积分方程的建立 在自由面与物面交点等拐点处,由于法向量 考虑一个二维、左右对称剖面以恒定速度 的不连续性导致速度势在拐点处连续而速度势导 垂直入水,初始时刻水面静止,并且此过程中没有 数不连续,若处理不当则会导致数值发散.用双点 空泡产生.建立笛卡尔坐标系如图 所示:图中 模型 处理该点可以得到比较精确的结果.在每 时刻求解方程 得到自由面上的速度势后,可用 阶龙格一库塔法对式中自由面条件式进行迭 / 代得到下一时刻的自由面和自由面速度势,由于 ? / \ 式数值解不够精确会导致出现自由面锯齿现 ? 象,用 点 阶平滑公式 可移除锯齿现象.在迭/ 代过程中,自由面上的节点可能靠得太近或太远, , , , 在每一时刻采用 次样条对自由面网格进行重新 图 坐标系定义和问题描述 划分可保证方程数值解的精确性.斜升角较小 轴和静止自由面重合, 轴和剖面对称轴重合且 的物体砰击过程中会形成非常薄的射流区,自由 垂直向上.假设流体不可压缩且不考虑黏性,流动 面和物面贴得太近会导致射流区尤其是自由面与 无旋,则存在速度势 在流域内满足的基本方程、 物面交点附近区域的物理量数值解有较大误差, 边界条件和初始条件为 甚至导致计算发散,需要特殊处理.采用截断模 / /一 流域内 ; 型。。 截去射流区可保证数值解的稳定性,由于射/ 一一 物面; 流区附近的物面压强接近于大气压,因此采用截/ ? / ,/ / 断模型不会对物体受力产生显著影响.在求得每 自由面; 时刻的物面速度势后,可用差分法直接求解式 / ? 。/ ? 自由面; 右端括号里的第 项进而求得物面压强,但数值 计算表明差分法计算物面压强容易带来误差,采/ 一 底部; 用辅助函数法口可以得到更为精确的结果./ 一 远方边界面, 式