绕船体自由表面流的数值模拟

本文采用RANS方程和VOF算法数值模拟了系列60粘性自由面流动，比较了幂指数律，二阶迎风，QUICK格式对计算结果的影响，详细给出了波形，阻力，伴流等的计算结果，与试验及其它计算结果比较显示本文的计算方法具有较好的计算精度，可用于船型性能分析和优化。     

应用分区有限体积法数值计算舵粘性水动力

通过数值求解不可压缩粘性流的RANS方程计算均匀来流中的船舵水动力。采用标准的κ-ε湍流模式封闭方程，并用二阶有限体积法在非交错分区贴体网格上对守恒方程进行离散；采用Rhie和Chow的插值方法确定单元面的质量通量，用SIMPLE法求解压力分布，用壁函数法处理壁面边界条件。对一种NACA冀型舵在两种不同雷诺数下大攻角范围的流动进行了计算，得到的舵水动力与试验结果和他人的计算结果进行比较，在舵未达到失速角前吻合良好，而所确定的舵失速角与最大升力值则均偏小。失速角后舵背流面的流动分离点也能捕捉到。最后分析了雷诺数对舵水动力的影响。     

带舵船体自航条件下周围流动的数值模拟

本文介绍了对一艘带舵油轮在自航条件下的周围流动所进行的数值模拟。船体和舵周围的粘性流场是用数值求解假设为不可压缩的纳维尔－斯托克斯方程得到的。利用等效体积力分布来考虑螺旋桨的影响。通过无限叶数的螺旋桨理论估算螺旋桨体积力。计算结果与试验数据进行的比较证实了本方法的有效性，表明该方法可作为船舶设计的有效的工具。     

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

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

绕Wigley船自由表面粘性流场计算

本文采用商用粘性流场求解软件COMET计算带自由液面绕Wigley船的粘性流动，网格数约为216000。计算中使用标准k-ε湍流模式和壁函数。用HRIC算法确定自由面。同时域步进法得到稳态解，计算的在傅氏数为0.30时的阻力结果与模型试验测量值吻合良好。预报的船体周围的波型及船体表面的波面形状也较合理。     

船舶在浅水域横向停靠运动粘性水动力计算

研究船舶在浅水域的低速横向漂移运动对研究船舶停靠运动特性具有重要意义.船舶停靠通常是一个低速慢漂的过程,在这个过程中,由于船舶特定的运动形式和几何形状,粘性水动力起主要作用,而兴波很小,其影响可以忽略,可以将自由面作为刚壁处理,从而使问题简化并大大减小数值计算的计算量.文中以标准Wigley船型为对象,对船舶在浅水中横向漂移运动的粘性流场和水动力进行了计算.为了得到准确、稳定的计算结果,分别选取了不同的计算区域、网格数量和湍流模型进行计算,并将计算结果与他人的计算结果和模型试验结果进行了比较,确定了模拟船舶低速停靠运动粘性流场的有效的数值方法.     

k—ω湍流模式在船舶粘性流场计算中的应用

湍流模式的选取对精确预报船舶艉部流场至关重要，本文利用SSTk-ω湍流模型计算了SR196A油船模型的粘性绕流场，船艉三个截面上的轴向速度等值线和横向速度矢量的计算结果和文献提供的试验值进行了对比，结果符合得很好，同时计算得到的阻力系数也与试验值进行了比较。结果同样令人满意，表明SSTk-ω湍流模型在船舶粘性流场的计算中有较强的应用能力。     

两种形式支柱小水线面双体船粘性阻力计算方法研究

对船舶粘性阻力工程及理论计算方法作了介绍,构建了一典型直支柱SWATH数学船型,变换得到斜支柱SWATH,运用数值与工程估算两种方法对SWATH粘性阻力进行计算,通过对比分析,验证了工程估算的适应性以及数值计算的准确性。因此,两种计算方法在SWATH设计初始可以作为其粘性阻力计算和船型优化分析的有效工具。     

串列螺旋桨水动力性能的数值预报

为了分析串列螺旋桨的水动力性能,本文运用计算流体动力学理论,结合雷诺时均 RANS方程和相对运动参考坐标系对其三维定常粘性流动进行数值模拟。应用 Fortran语言编制程序计算螺旋桨的型值点,并采用三次样条曲线拟合各点,建立串列桨三维模型。以某一串列螺旋桨作为研究对象,得到螺旋桨的推力系数、转矩系数以及流域内速度分布等水动力特性参数,并给出敞水性能曲线。计算结果与试验数据吻合较好,验证了数值方法的可行性和准确性。     

串列螺旋桨水动力性能的数值预报

为了分析串列螺旋桨的水动力性能,本文运用计算流体动力学理论,结合雷诺时均RANS方程和相对运动参考坐标系对其三维定常粘性流动进行数值模拟。应用Fortran语言编制程序计算螺旋桨的型值点,并采用三次样条曲线拟合各点,建立串列桨三维模型。以某一串列螺旋桨作为研究对象,得到螺旋桨的推力系数、转矩系数以及流域内速度分布等水动力特性参数,并给出敞水性能曲线。计算结果与试验数据吻合较好,验证了数值方法的可行性和准确性。     

绕船体自由面周围三维粘性流场的数值模拟

本文采用有限体积法通过求解不可压缩的雷诺平均(RANS)方程数值模拟了包括兴波的三维船体周围的粘性流场,湍流模式使用了子网格尺度模式(SGS)和Baldwin-Lomax模式相结合的混合模式.对于自由表面的处理,采用了任意拉格朗日-欧拉方法,网格为不仅与物体表面贴体,而且与自由表面贴体的动网格,即随着自由表面的变化要不断地重新划分网格.虽然此方法需要很长的计算时间,但能较好地描述船体的兴波情况.本文计算了系列60船模在Fn=0.316,Re=1.9×106时带自由面的粘性流场,计算结果与试验结果吻合较好.，An finite - volume method solving of incompressible RANS Equations is developed for the numeri-cal simulation of three - dimensional viscous flow with free - surface about a ship. A hybrid turbulence mod-el by combining the sub - grid scale (SGS) model and the Baldwin- Lomax model is used in this pa-per. The arbitrary - Lagrange - Euler formulation is devised for the treatment of free surface boundary condi-tion. The boundary - fitted coordinate system is fitted not only to the ship hull surface but also to the freesurface, so the computational grid is regenerated to follow the free surface deformation. Although this methodneeds spend much longer computational time, but it can fairly display wave pattern made by a ship. The nu-merical simulation has been carried out for viscous flow with free surface past Series 60 ship hull at Fn =0. 316, Re = 1.9 × 106. The computed results are in reasonable agreement with experimental measure-ments.     

潜艇水面与水下粘性绕流数值模拟

本文采用求解RANS方程的方法结合四种湍流模型,对于带有不同附体的SUBOFF模型尾流场进行了数值模拟。数值预报的桨盘面处不同半径上的轴向无量纲速度与试验结果进行了对比,结果表明湍流模型在数值模拟中起到重要作用。潜艇水面航行性能十分重要,因而对于潜艇自由液面绕流的数值模拟备受关注。本文采用VOF方法对于两条潜艇模型在不同傅汝德数下的自由液面绕流进行了数值模拟。计算得到的阻力、波形与试验结果吻合较好。文中也探讨了潜艇自由液面绕流的一般特性。并验证了用CFD手段预报潜艇粘性流场的能力。     

三体船阻力数值计算及方案优选

采用CFD技术模拟了船池试验三体船型在不同速度下的粘性流场,计算中考虑了自由表面效应影响,数值模拟结果与试验结果进行了比较,然后在此基础上计算了侧体对应于不同纵横向位置时的9种三体船方案的总阻力.分析比较各方案的计算结果后,得出三体船阻力性能较优时侧体布局的一般规律,即侧体纵向位置变化对三体船静水阻力性能影响较大,当侧体横向位置不变,在较高速度时侧体纵向向后布置对阻力性能更有利.最后按阻力性能,从9种计算方案中得出优选方案.     

基于RANS方程的船模兴波阻力分离计算方法及其适用范围研究

运用空气-水两相流RANS方程计算了两型船模的粘性兴波流场,提出了兴波阻力的分离计算方法,并研究了该方法的适用范围.结果表明,在长度傅汝德数低于0.4时,文中对船模兴波阻力及波高的计算结果都与模型试验基本相符,但当长度傅汝德数高于0.4时,计算结果与试验还存在一定的差距.     

重叠网格在船舶CFD中的应用研究

文章采用RANS方法和重叠网格计算了带自由液面的船舶绕流问题。计算网格采用重叠网格的型式,自由液面的模拟采用单相Level-Set方法,Reynolds应力采用k-ω模型,采用体单元有限差分方法和PISO算法求解RANS方程。文中简要描述了重叠网格和单相Level Set自由液面模拟方法的数学模型及求解。通过对S60单体船型约束模和自由模型采用重叠网格的数值求解及与试验结果的比较表明该重叠可较好地模拟带自由液面船舶自由态绕流问题。另采用重叠网格对一简单的双体船和三体船进行了数值求解,计算结果也表明重叠网格和单相Level-Set方法可较好地模拟带自由液面的船舶绕流问题。     

潜艇近水面航行兴波特征研究

潜艇作为一种水下航行器，其近水面运动时产生的绕流场远比无界直航运动时要复杂得多，其近水面航行时对自由液面的影响不容忽视。本文采用粘性流体力学计算方法，对SUBOFF潜艇近水面航行时的粘性绕流场进行模拟，对自由液面的兴波进行计算，并对波形特点进行深入分析。通过本文的研究，得到潜艇近水面航行时自由液面兴波的特点和相关规律。     

实尺度喷水推进船拖泵工况数值模拟与分析

对于四泵推进的喷水推进船,在巡航工况时中间加速泵通常处于锁轴状态,其拖曳阻力的大小对喷水推进器的选型以及船泵机的最优匹配有着重要影响。然而,拖泵阻力很难通过船模试验的方法获得。为此,该研究在验证均匀和非均匀条件喷水推进器数值模型的准确性基础上,采用数值试验的方法对18节航速下某双泵推进喷水推进船的实尺度"船体+两台喷水推进器"系统带自由液面的流场进行了数值模拟,计算此时喷泵拖曳阻力及其所占船体阻力的百分比。以此喷泵拖曳阻力作为参考,对尺寸与上述喷水推进泵相近的某四泵推进喷水推进船的喷泵进行了选型和设计,并对该船在18节航速下加速泵拖曳阻力的大小进行了计算,进一步验证选型时拖曳阻力取值的合理性。为消除尺度效应的影响采用实尺度模型对"船体+四台喷水推进器"系统带自由液面的非定常流场进行计算,并探索了大尺度条件下船泵系统考虑自由液面和重力影响的非定常计算方法。     

潜艇应急上浮稳性研究

通过对应急上浮运动及高压空气吹除主压载水舱规律的简化,给出了潜艇应急上浮运动方程.将上层建筑的背水量作为临时加载,推导出应急上浮稳性计算方法.将非定常势流理论用于应急上浮运动流场分析与水动力计算.考虑了非线性自由表面条件,在应急上浮运动中随时进行流场以及自由表面的运动计算,从而实现了应急上浮运动与流场的耦合求解.采用将艇体及自由表面划分成四边形单元的方法计算流场,用Adams方法计算自由表面运动和艇体上浮运动.文中给出了某潜艇应急上浮运动和稳性的算例,计算结果验证了方法的可靠性.     

船体艏部水动力性能优化

提出一种以势流兴波阻力理论Rankine源方法为基础,结合SHIPFLOW软件为计算工具,利用CAD-CFD集成平台FRIENDSHIP-Framework软件进行变形优化,研究船舶的最小兴波阻力型线优化设计的方法,并考察了兴波优化得到的船型总阻力变化情况。在型线优化过程中,以兴波阻力系数为目标函数,排水量变化范围为约束条件,在Wigley船体前端增加一个利用Feature建模技术参数化生成的球艏并调整艏部型线使得船体表面光顺。选取球鼻艏形状的各项参数作为基本设计变量,利用DOE方法对船艏进行优化,获得了设计航速下兴波阻力较小的船型,验证了所提方法进行船艏型线优化的有效性。相应的考察变形及优化前后总阻力变化情况表明:在高傅汝德数情况下增加球艏所带来的粘性阻力的增加小于兴波阻力的减小量,总阻力得到了改善,优化后得到的球艏能在进一步减小兴波的同时减小总阻力。此外,还运用所提方法对3100TEU船型的船艏,利用Delta Shift方法进行变形,在设计航速下,将变形的参数作为设计变量,利用DOE方法进行优化设计。结果显示:在排水量限制范围内当球鼻艏向上向前伸展一定长度时可以降低兴波阻力。与此同时,由于优化前后船体湿表面积变化很小,粘性阻力的变化并不明显,兴波的减小则使得总阻力得到了改善。     

舰船气泡尾流数值模拟(英文)

Based on a volume of fluid two-phase model imbedded in the general computational fluid dynamics code FLUENT6.3.26, the viscous flow with free surface around a model-scaled KRISO container ship (KCS) was first numerically simulated. Then with a rigid-lid-free-surface method, the underwater flow field was computed based on the mixture multiphase model to simulate the bubbly wake around the KCS hull. The realizable k-ε two-equation turbulence model and Reynolds stress model were used to analyze the effects of turbulence model on the ship bubbly wake. The air entrainment model, which is relative to the normal velocity gradient of the free surface, and the solving method were verified by the qualitatively reasonable computed results.