新型深吃水多立柱FDPSO变截面立柱涡激运动特性研究

文章采用FLUENT软件结合分离涡法对某新型深吃水多立柱FDPSO变截面立柱涡激运动特性开展了研究。将变截面立柱涡激振动系统简化为质量-弹簧-阻尼模型,引入雷诺平均应力模型求解不可压缩粘性Navier-Stokes方程,通过计算出流场作用于柱体的瞬时升力与阻力,并基于UDF程序求解运动微分方程同时运用动网格技术实现流场更新,实现了对变截面立柱涡激运动的数值模拟。对不同来流速度下变截面柱绕流和涡激运动特性进行了研究,分析了变截面柱阻力系数、升力系数、运动幅值和运动轨迹等,研究发现：在高雷诺数单柱绕流模拟中采用三维数值模型更具优势,变截面柱后方同一时间内可能有多个旋涡产生,旋涡之间相互影响,抑制了旋涡对柱体的升阻力作用;在折合速度6~9范围内,变截面柱未出现明显“锁定”现象;变截面柱涡激运动的横荡振幅峰值在Ur=4.5,约为1.2D;变截面柱流向平衡位置随着折合速度增加而增加;变截面柱运动轨迹出现了典型的“8”字形,但变截面柱的轨迹相对较紊乱。     

Numerical analysis on free surface waves and stern viscous flow of a ship model

This paper deals with numerical techniques for computing the viscous flow past a ship hull with and without a free surface using a Reynolds-averaged Navier-Stokes solver with global conservation. In the first technique, a coarse grid is used to find an approximate solution to the free surface problem. Interpolation of a fine grid is subsequently carried out, and a more exact solution, particularly in the boundary layer and wake, is obtained. In the second technique, a modified Baldwin-Lomax model is introduced to compute the viscous flow with and without a free surface. These numerical techniques are applied to simulations of the flow around a Series 60 and an SR196C ship model. The results are compared with measurement data, and the usefulness of the numerical techniques is demonstrated.     

Numerical analysis of ducted propeller performance under open water test condition

This paper investigates the open water performance of the Ka-series propellers at various pitch and expanded area ratios in combination with the 19A duct by employing the panel method panMARE and the RANSE code ANSYS-CFX. An efficient method, Caly, is developed in order to generate the 3D-geometry and the surface numerical grid of the ducted propellers. Caly can be coupled with ANSYS-TurboGrid to automatically produce 3D-grids for the RANSE solver. The numerical results are compared with published experimental data and the flow details are concluded and compared. The influences that the grid resolution, the panel arrangements of duct and blade, and the flow in gap between inside wall of the duct and blade tip on the numerical results are studied. Grids verification, turbulence model dependency analysis and Reynolds number scantling are also discussed.     

用naoe-FOAM-SJTU求解器数值模拟船舶在波浪上的操纵运动问题（英文）

Ship maneuvering in waves includes the performance of ship resistance, seakeeping, propulsion, and maneuverability. It is a complex hydrodynamic problem with the interaction of many factors. With the purpose of directly predicting the behavior of ship maneuvering in waves, a CFD solver named naoe-FOAM-SJTU is developed by the Computational Marine Hydrodynamics Lab(CMHL) in Shanghai Jiao Tong University. The solver is based on open source platform OpenFOAM and has introduced dynamic overset grid technology to handle complex ship hull-propeller-rudder motion system. Maneuvering control module based on feedback control mechanism is also developed to accurately simulate corresponding motion behavior of free running ship maneuver. Inlet boundary wavemaker and relaxation zone technique is used to generate desired waves. Based on the developed modules, unsteady Reynolds-averaged Navier-Stokes(RANS) computations are carried out for several validation cases of free running ship maneuver in waves including zigzag, turning circle, and course keeping maneuvers. The simulation results are compared with available benchmark data. Ship motions, trajectories, and other maneuvering parameters are consistent with available experimental data, which indicate that the present solver can be suitable and reliable in predicting the performance of ship maneuvering in waves. Flow visualizations, such as free surface elevation, wake flow, vortical structures, are presented to explain the hydrodynamic performance of ship maneuvering in waves. Large flow separation can be observed around propellers and rudders. It is concluded that RANS approach is not accurate enough for predicting ship maneuvering in waves with large flow separations and detached eddy simulation(DES) or large eddy simulation(LES) computations are required to improve the prediction accuracy.     

浅水中低速斜航船舶水动力预报及验证与确认分析

浅水中的斜航船舶受到浅水阻塞效应和不对称流的综合影响。为预报该运动中的船舶水动力,文章采用基于定常雷诺平均纳维—斯托克斯方程的计算流体动力学方法,对浅水中做斜航运动的船舶粘性绕流场进行数值模拟。考虑低航速运动的特点,忽略航速影响下的自由面兴波,由数值计算得到水动力系数在漂角影响下的变化规律。针对计算精度问题,在数值模拟中从验证和确认角度分析和评估计算结果：通过网格收敛性分析分析数值误差与不确定度;结合试验数据考察计算模型的误差。此外,从计算区域尺度、湍流模型、边界条件、船体下沉和纵倾作用方面对模型误差的影响因素进行探讨,可为改进计算模型、提高数值模拟精度提供参考依据。     

Validation of a dynamic mooring model coupled with a RANS solver

Standard design procedures and simulation tools for marine structures are aimed primarily for use by the offshore oil and gas. Mooring system restoring forces acting on floating offshore structures are obtained from a quasi-static mooring model alone or from a coupled analysis based on potential flow solvers that do not always consider nonlinear mooring-induced restoring forces, fluid structure interactions, and associated hydrodynamic damping effects. This paper presents the validation of a dynamic mooring system analysis technique that couples the dynamic mooring model with a Reynolds-averaged Navier-Stokes (RANS) equations solver. We coupled a dynamic mooring model with a RANS equations solver, and analyzed a moored floating buoy in calm water, regular and irregular waves and validated our motion and mooring force predictions against experimental measurements. The mooring system consisted of three catenary chains. The analyzed response comprised decaying oscillating buoy motions, linear and quadratic damping characteristics, and tensile forces in mooring lines. The generally favorable comparison of predicted buoy motions and mooring forces to experimental data confirmed the reliability of our implemented coupling technique to predict system response. Additional comparative results from a potential flow solver demonstrated the benefits of the coupled dynamic mooring model with RANS equations. The successful validated tool of coupling the dynamic mooring model with the RANS solver is available as open source, and it shows the potential of the coupled methodology to be used for analyzing the moored offshore structures.     

散货船阻力预报中湍流模型的应用

采用直接求解RANS方程的方法对20艘散货船的绕流场和阻力进行了数值模拟,比较了标准模型、RNG模型及SST模型对数值计算结果的影响.通过与试验结果的详细对比,发现只要选取合适的湍流模型,数值模拟结果能够达到较高的精度,在1% -2%之间.这表明该数值计算方法能够满足工程的实际使用要求,可以在工程中推广应用.     

风机塔架与叶片相互作用的数值模拟(英文)

Numerical simulations of wind turbine blade-tower interaction by using the open source OpenFOAM tools coupled with arbitrary mesh interface(AMI) method were presented.The governing equations were the unsteady Reynolds-averaged Navier-Stokes(RANS) which were solved by the pimpleDyMFoam solver,and the AMI method was employed to handle mesh movements.The National Renewable Energy Laboratory(NREL) phase VI wind turbine in upwind configuration was selected for numerical tests with different incoming wind speeds(5,10,15,and 25 m/s) at a fixed blade pitch and constant rotational speed.Detailed numerical results of vortex structure,time histories of thrust,and pressure distribution on the blade and tower were presented.The findings show that the wind turbine tower has little effect on the whole aerodynamic performance of an upwind wind turbine,while the rotating rotor will induce an obvious cyclic drop in the front pressure of the tower.Also,strong interaction of blade tip vortices with separation from the tower was observed.     

Finite-volume simulation of flows about a ship in maneuvering motion

A finite-volume method of computing the viscous flow field about a ship in maneuvering motion was developed. The time-dependent Navier-Stokes equation discretized in the generalized boundary-fitted curvilinear coordinate system is solved numerically. A third-order upwind differencing scheme, a marker and cell (MAC)-type explicit time marching solution algorithm and a simplified subgrid scale (SGS) turbulence model are adopted. The simulation method is formulated, including the movement of a computational grid fitted to the body boundary that allows computation of the flow field around a body under unsteady motion. To estimate the maneuvering ability of a ship, the accurate prediction of the hydrodynamic forces and moments of the hull is important. Therefore, experimental methods of finding the hydrodynamic forces of a ship in maneuvering motion, such as the oblique towing test, the circular motion test (CMT) and planar motion mechanism (PMM) test, were established. Numerical simulation methods for those captive model experiments were developed introducing computational fluid dynamics (CFD). First, numerical methods for steady oblique tow and steady turn simulation were developed and then extended to unsteady forced motion. Simulations were conducted about several realistic hulls, and the results were verified by comparisons with measured results obtained in model experiments. Hydrodynamic forces and the moment, the longitudinal distribution of the hydrodynamic lateral force, and the pressure distribution on the hull surface showed good agreement.     

Cavitating Propeller Performance in Inclined Shaft Conditions with OpenFOAM: PPTC 2015 Test Case

In this paper, we present our analysis of the non-cavitating and cavitating unsteady performances of the Potsdam Propeller Test Case (PPTC) in oblique flow. For our calculations, we used the Reynolds-averaged Navier-Stokes equation (RANSE) solver from the open-source OpenFOAM libraries. We selected the homogeneous mixture approach to solve for multiphase flow with phase change, using the volume of fluid (VoF) approach to solve the multiphase flow and modeling the mass transfer between vapor and water with the Schnerr-Sauer model. Comparing the model results with the experimental measurements collected during the Second Workshop on Cavitation and Propeller Performance – SMP’15 enabled our assessment of the reliability of the open-source calculations. Comparisons with the numerical data collected during the workshop enabled further analysis of the reliability of different flow solvers from which we produced an overview of recommended guidelines (mesh arrangements and solver setups) for accurate numerical prediction even in off-design conditions. Lastly, we propose a number of calculations using the boundary element method developed at the University of Genoa for assessing the reliability of this dated but still widely adopted approach for design and optimization in the preliminary stages of very demanding test cases.     

带自由面肥大船粘性绕流场的数值模拟

采用直接求解RANS方程的方法对带自由面的肥大型船舶的流场和阻力进行了数值计算,湍流模型选用RNGk-ε两方程模式,自由面的处理采用VOF方法.数值计算结果与已公开试验数据的定量比较相当吻合.计算结果的分析表明:文中发展的CFD计算方法可以满足于带自由面肥大型船舶快速性工程预报的要求.     

基于OpenFOAM的螺旋桨敞水性能预报方法

为深入研究螺旋桨周围的粘流问题,基于面向对象的开源CFD计算平台OpenFOAM,选择DTMBP4119桨作为对象,利用RANS方程计算了桨的敞水性能,并分别考察了网格依赖性和离散格式的影响。同时分析了不同半径处叶剖面的压力分布及桨前后的流场速度分布。为保证对流项离散的稳定性和高精度性,采用了二阶NVD格式Gamma混合差分格式,选取k-wSST模型作为湍流模型,压力速度耦合采用SIMPLE松弛算法。计算结果与试验数据吻合较好。研究建立了基于OpenFOAM的螺旋桨敞水性能预报方法。     

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

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

舵桨相互作用下现代船舶绕流数值研究（英文）

Reducing the fuel consumption of ships presents both economic and environmental gains. Although in the past decades,extensive studies were carried out on the flow around ship hull, it is still difficult to calculate the flow around the hull while considering propeller interaction. In this paper, the viscous flow around modern ship hulls is computed considering propeller action. In this analysis, the numerical investigation of flow around the ship is combined with propeller theory to simulate the hull-propeller interaction. Various longitudinal positions of the rudder are also analyzed to determine the effect of rudder position on propeller efficiency. First, a numerical study was performed around a bare hull using Shipflow computational fluid dynamics(CFD) code to determine free-surface wave elevation and resistance components.A zonal approach was applied to successively incorporate Bpotential flow solver^ in the region outside the boundary layer and wake, Bboundary layer solver^ in the thin boundary layer region near the ship hull, and BNavier-Stokes solver^in the wake region. Propeller open water characteristics were determined using an open-source MATLAB code Open Prop, which is based on the lifting line theory, for the moderately loaded propeller. The obtained open water test results were specified in the flow module of Shipflow for self-propulsion tests. The velocity field behind the ship was recalculated into an effective wake and given to the propeller code that calculates the propeller load. Once the load was known, it was transferred to the Reynolds-averaged Navier-Stokes(RANS) solver to simulate the propeller action. The interaction between the hull and propeller with different rudder positions was then predicted to improve the propulsive efficiency.     

Numerical analysis of surface piercing propeller in unsteady conditions and cupped effect on ventilation pattern of blade cross-section

The aim of this study is to calculate hydrodynamic performance and ventilation flow around wedge, 2D blade and 3D surface piercing propeller (SPP), using computational fluid dynamic based on Reynolds-averaged Navier–Stokes method. First, numerical analyses for two-phase fluid flow around the wedge and 2D blade section (cupped and non-cupped) are presented. Flow ventilation, pressure distribution and forces are determined and compared with experimental data. Then, the method is extended to predict the hydrodynamic performance of propeller SPP-841B. The propeller exhibits a cupped blade. In the simulated configuration, SPP is one-third submerged (I = h/D = 0.33) and is working at various loadings with full ventilation occurring at low advance coefficient (J). The open water performance, pressure distribution, forces/moments and ventilation pattern on the SPP-841B model are obtained and compared with experimental data. The numerical results are in good agreement with experimental measurements, especially at high advance coefficient.     

Single- and multiobjective design optimization of a fast multihull ship: numerical and experimental results

Numerical optimization of the initial design of a fast catamaran (high-speed sealift research model B, HSSL-B) has been carried out through a simulation-based design (SBD) framework, based on an advanced free-surface unsteady Reynolds-averaged Navier–Stokes (URANS) solver and a potential flow solver, and global optimization (GO) algorithms. The potential flow computational fluid dynamics (CFD) SBD was used to guide the more expensive URANS CFD SBD. The fluid-dynamic analysis of the flow past the catamaran proved that the use of the URANS solver was fundamental in dealing with the multihull interference problem. In the case investigated, the separation distance was small and the viscous flow quite distorted by the proximity of the hulls, so that only viscous solvers could correctly capture the flow details. Sinkage and trim effects, due to the high speed range and again to the small separation distance investigated, are also relevant. The initial HSSL-B geometry and three optimization problems, including single- and multiobjective optimization problems, proposed by designers from Bath Iron Works, were successfully optimized/solved, and finally an experimental campaign was carried out to validate the optimal design. A new verification and validation methodology for assessing uncertainties and errors in simulation-based optimization was used based on the trends, i.e., the differences between the numerically predicted improvement of the objective function and the actual improvement measured in a dedicated experimental campaign, including consideration of numerical and experimental uncertainties. Finally, the success of the optimization processes was confirmed by the experimental measurements, and trends for total resistance, sinkage, and trim between the original and optimal designs were numerically and experimentally verified and validated.     

复合材料螺旋桨水动力特性的流固耦合数值模拟

为了分析复合材料螺旋桨变形对其水动力性能的影响,利用FLUENT和ANSYS结构模块建立了一种流固耦合的方法。基于此方法分析了桨叶变形对敞水曲线,桨叶附近流场及其表面压力的影响。研究结果表明,初始几何为无侧斜无纵倾的螺旋桨变形后纵倾发生改变,使推力和扭矩系数变大,且推力系数和扭矩系数的增值随进速系数减小而增大。变形后的螺旋桨桨叶表面压力增大,压力系数变化最大值可达40%,螺旋桨轴向诱导速度变化最大值可达18.7%。     

带前置定子的导管桨片空泡数值预报及梢涡空泡起始判别研究

对水下航行体后带前置定子的导管桨片空泡和梢涡空泡起始进行研究。首先,采用基于混合网格的RANS求解器结合Singal空泡模型,数值模拟带前置定子导管桨片空泡形态,并与试验结果进行对比,预报的片空泡形态与试验结果吻合良好,表明用该方法预报水下航行体后带前置定子导管桨片空泡可行;其次,对水下航行体后带前置定子的导管桨多个方案全湿流动进行了数值模拟,通过梢部端面的最小压力值判断空泡起始,与试验结果进行比对,其吻合良好,表明可用该方法判别梢涡空泡起始。     

A porous cell method for the simulation of fluid-solid interactions

There are many ways of describing a solid,porous or fluid region of the computational domain when solving the Navier-Stokes equations(NSE)for flow motions.Amongst these the porous cell method is one of the most flexible approaches.In this method,a parameter is defined as a ratio of the volume open to water and air in a calculation cell to its cell volume.In the calculation,the same numerical procedure is applied to every cell and no explicit boundary conditions are needed at solid boundaries.The method is used to simulate flow through porous media,around solid bodies and over a moving seabed.The results compare well with experimental data and other numerical results.In our future work the porous cell method will be applied to more complex fluid-solid interaction situations.     

CFD在潜艇外形方案比较中的应用

以水滴型潜艇为母型进行改造,提出了新型椭圆截面和碟形截面的潜艇外形设计.水滴型潜艇的型值选用"大青花鱼"号的简化型值,在船长与排水量不变的条件下通过改变截面形状给出了后两种新艇型的型值,建立了计算模型.采用有限体积法求解不可压缩的雷诺平均(RANS)方程数值模拟了三种艇体周围的粘性流场,通过对三者的粘性绕流和阻力性能的比较,给出了三种方案各自的优缺点.