均匀流中潜艇水下运动表面尾迹的数值模拟

本文通过基于RANS方程和VOF方法的非定常粘性数值方法对潜艇在均匀流中水下运动的表面尾迹特征进行了探索性研究,较好地解决了用粘性流方法来模拟远场波浪中的数值耗散问题,数值计算得到的潜艇运动尾迹特征与诸多文献和理论分析结果基本一致.本文均匀流中的数值模拟将成为潜艇在分层流运动数值模拟的基础.     

参数化船型的阻力计算

以KCS船型为研究对象,用FRIENDSHIP软件创建CFD数值计算所需的几何模型；将船舶阻力分为兴波阻力和粘性阻力两部分计算,兴波阻力通过基于势流理论方法求解Euler方程得到,粘性阻力根据对RANS方程的求解来获取；将计算结果与实验数据进行对比分析.结果表明该参数化模型适用于CFD数值模拟.     

小水线面双体船粘性流数值模拟

利用商业软件FLUENT对一小水线面双体船(SWATH)的粘性绕流进行数值模拟,得到了不同航速下的三维粘性流场和粘性阻力,通过对计算结果的分析、比较,验证了FLUENT用于预报小水线面双体船粘性阻力和伴流分布的有效性和实用性。     

二维浮式结构物与完全非线性波相互作用数值分析

基于粘流理论、流固耦合理论和自由液面追踪技术——VOF(Volume of Fluid)法,以不可压缩流体的连续方程和N-S方程及结构动力学方程为基本控制方程,利用商用软件ANSYS Workbench的二次开发功能,建立了粘性数值波浪水池,数值模拟了完全非线性波的生成,计算分析了二维浮式结构物遭遇波浪的完全非线性现象,得到了甲板上浪水位高度和结构物的纵摇、垂荡、横摇等运动响应。结果表明:文中的方法可以用于波浪与浮式结构物相互作用的完全非线性数值模拟.     

全附体潜艇三维粘性流场数值计算方法研究

采用不同的湍流模型和网格形式对SUBOFF全附体模型进行三维粘性流场数值模拟,将计算所得的艉部伴流场同有关文献的试验结果进行比较,并以伴流分数对周向角的曲线形式展现出来,分析了湍流模型、艇体周向网格节点间距和y+值对全附体潜艇三维粘性流场数值计算结果的影响.     

基于CFD的规则波顺浪数值水池模拟

基于CFD（ Computational Fluid Dynamics）和粘性数值波浪水池技术,对顺浪航态的规则波波浪环境进行数值模拟,采用边界造波法进行波浪生成,VOF（Volume of Fluid）法模拟自由面。本文改进了传统的基于规则波的顶浪数值模拟方法,在特定波浪、不同航速条件下进行了顺浪航态的数值水池模拟。根据模拟呈现的动态效果,以及与Matlab仿真结果的对比,对试验的准确性进行验证。     

基于OpenFOAM的船舶与液舱流体晃荡在波浪中时域耦合运动的数值模拟

波浪中载液船舶运动激励舱内液体的晃荡,舱内液体晃荡产生的冲击力同时作用在舱壁上,进而影响船舶的运动姿态。波浪中船体水动力和时延函数是在势流理论范畴下采用切片法和脉冲响应函数方法计算获得的,液舱内液体非线性晃荡是基于粘性流理论实时计算模拟,两者耦合建立了波浪中载液船舶与液舱流体晃荡耦合的运动方程。论文基于开源CFD开发平台OpenFOAM,自主开发实现了船体运动与液舱晃荡的耦合计算程序,并进行了相应的数值模拟计算和验证工作。该方法完整地考虑了波浪、船体和液舱晃荡之间的耦合作用,并结合船体内外流场特点分别采用了势流和粘性流理论,具有较高的计算效率。通过数值模拟计算和模型实验研究表明,数值模拟计算能够清晰显现出液舱晃荡对船体全局运动影响,船体运动计算结果与模型实验结果吻合良好。     

基于CFD的规则波顺浪数值水池模拟

基于CFD(Computational Fluid Dynamics)和粘性数值波浪水池技术,对顺浪航态的规则波波浪环境进行数值模拟,采用边界造波法进行波浪生成,VOF(Volume of Fluid)法模拟自由面。本文改进了传统的基于规则波的顶浪数值模拟方法,在特定波浪、不同航速条件下进行了顺浪航态的数值水池模拟。根据模拟呈现的动态效果,以及与Matlab仿真结果的对比,对试验的准确性进行验证。     

用基于SST模型的DES方法数值模拟圆柱绕流

分离涡数值模拟方法(DES)在物面附近求解雷诺平均Navier-Stokes方程,在其他区域采用大涡模拟方法.兼具前者计算量小的优点和后者能模拟大分离湍流流动、计算精度高的优势.采用基于剪切应力传输(SST)两方程湍流模型的DES方法对粘性不可压缩流体的圆柱绕流问题进行数值模拟.通过对所得的速度场、压力场、阻力系数Cd、升力系数Cl、斯特罗哈数St等结果的分析及与文献上的实验和计算数据的比较,说明DES方法对于层流及高亚临界雷诺数的圆柱绕流模拟是合理的.     

剪切流下海洋立管涡激振动的三维数值模拟

对1根长为9.68 m的三维海洋立管涡激振动现象进行流固耦合模拟,立管的长径比为482,来流为剪切流。通过在立管长度方向上施加一固定的顶部预紧力,来更好地控制立管变形。三维粘性、不可压缩流体场采用非稳态的N-S方程和k-ω湍流模型进行三维CFD数值模拟,固体场采用基于三维实体单元的有限元方法进行模拟,通过一种新的方法 System Coupling实现流-固耦合交界面的数据交换。模拟结果与均匀流模拟结果相似,反映了多模态的振动特性,说明这种方法有一定的可行性。立管脱落呈现多种涡结构模式,稳定后呈"2S"。立管出现大的非对称弯曲变形现象,介于1阶模态和2阶模态之间转换,最终锁定在2阶模态。