System-based investigation on 4-DOF ship maneuvering with hydrodynamic derivatives determined by RANS simulation of captive model tests

For the non-negligible roll-coupling effect on ship maneuvering motion, a system-based method is used to investigate 4-DOF ship maneuvering motion in calm water for the ONR tumblehome model. A 4-DOF MMG model is employed to describe ship maneuvering motion including surge, sway, roll, and yaw. Simulations of circular motion test, static drift and heel tests are performed by solving the Reynolds-averaged Navier-Stokes (RANS) equations, after a convergence study quantifying the necessary grid spacing and time step to resolve the flow field adequately. The local flow field is analyzed for the selected cases, and the global hydrodynamic forces acting on the ship model are compared with the available experiment data. Hydrodynamic derivatives relating to sway velocity, yaw rate, and heel angle are computed from the computed force/moment data using least square method, showing good agreement with those obtained from EFD data overall. In order to investigate further the validity of these derivatives, turning circle and zigzag tests are simulated by using the 4-DOF MMG model with these derivatives. The trajectories and the time histories of the kinematic variables show satisfactory agreement with the data of free-running model tests, indicating that the system-based method coupled with CFD simulation has promising capability to predict the 4-DOF ship maneuvering motion for the unconventional vessel.     

Numerical Study on Air Turbines with Enhanced Techniques for OWC Wave Energy Conversion

In recent years, the oscillating water column (OWC) wave energy converter, which can capture wave energy from the ocean, has been widely applied all over the world. As the essential part of the OWC system, the impulse and Wells turbines are capable of converting the low pressure pneumatic energy into the mechanical shaft power. As an enhanced technique, the design of endplate or ring attached to the blade tip is investigated numerically in this paper. 3D numerical models based on a CFD-software FLUENT 12.0 are established and validated by the corresponding experimental results from the reports of Setoguchi et al. (2004) and Takao et al. (2001). Then the flow fields and non-dimensional evaluating coefficients are calculated and analyzed under steady conditions. Results show that the efficiency of impulse turbine with ring can reach up to 0.49 when φ=1, which is 4% higher than that in the cases for the endplate-type and the original one. And the ring-type Wells turbine with fixed guide vanes shows the best performance with the maximal efficiency of 0.55, which is 22% higher than that of the original one. In addition, the quasi-steady analysis is used to calculate the mean efficiency and output-work of a wave cycle under sinusoidal flow condition. Taking all together, this study provides support for structural optimization of impulse turbine and Wells turbine in the future.     

CFD simulation of airflow over a regular array of cubes. Part I: Three-dimensional simulation of the flow and validation with wind-tunnel measurements

Air flow inside an array of cubes is simulated. Cubes (edge length 0.15 m) are arranged in a regular array, separated by 0.15 m in the streamwise and spanwise directions. Numerical simulations are performed based on Reynolds-averaged Navier–Stokes equations (RANS), solved in a computational fluid dynamics model (CFD), with standard k–ε turbulent closure (two prognostic equations are solved for the turbulent kinetic energy k and its dissipation ε, respectively). Simulations are validated against wind-tunnel data using a technique based on hit-rate calculations, and calculated statistical parameters. The results show that the horizontal velocity is very well modelled, and despite some discrepancies, the model that fulfils the hit-rate test criteria gives useful results that are used to investigate three-dimensional (3-D) flow structures. The 3-D analysis of the flow shows interesting patterns: the centre of the canyon vortex is at 3/4 of the canyon height, and stronger downward than upward motions are present within the canyon. Such behaviour is explained by the presence of a compensation flow through the side of the canyon, which enters the canyon from the upper part and exits from the lower part. This complex 3-D structure affects the tracer dispersion, and is responsible for pollutant transport and diffusion.     

平台障碍物对海气通量观测影响的数值分析

基于计算流体动力学模型（CFD）模拟了背景风风速为2 m·s-1、4 m·s-1、6 m·s-1下通量观测平台周围的风场分布,确定了不同风速下平台障碍物对观测影响最小的风向范围。发现风速不同,平台对背景风的影响范围不同,风速越大,平台对背景风的影响越大。运用涡动相关法计算不同风向条件下的湍流通量,定量分析由于平台干扰对海气通量计算的影响,发现选取迎风条件下的数据计算得到的海气通量会比实际偏低。     

大神堂人工鱼礁礁区布局的合理性研究(英文)

以CFD(计算流体动力学)技术为手段,研究了不同配置组合人工鱼礁区周围的流场效应。研究中采用了四种配置组合,分别为"444"三列,"464"三列,"12345"五列三角形,"34543"五列正六边形进行数值模拟,结果说明大神堂人工鱼礁礁区的礁群配置组合能够较好地发挥礁区的调控作用。     

Simplified CFD modeling for bilge keel force and hull pressure distribution on a rotating cylinder

A horizontal, circular cylinder fitted with one bilge keel is forced to rotate harmonically around its axis. The bilge keel load and hull pressure distribution are investigated. A fully submerged condition (infinite fluid), and three partly-submerged conditions are considered. A two-dimensional numerical study is performed, and the results are validated against recently published experimental data by van’t Veer et al. [30]. In addition, comparisons for mass and drag coefficients are also made with experimental data for plate in infinite fluid (Keulegan and Carpenter [8]), and wall-mounted plate (Sarpkaya and O’Keefe [9]) in oscillatory flow.A Navier–Stokes solver based on the Finite Volume Method is adopted for solving laminar flow of incompressible water. The free-surface condition is linearized by neglecting the nonlinear free-surface terms and the influence of viscous stresses in the free surface zone, while the body-boundary condition is exact. This simplified modeling of the problem required the mesh to be fine only around the bilge keels, leading to a total number of cells around N ≃ 1 ×10⁴, which reduced computational cost significantly.The influence of draft and amplitude of oscillations on the bilge keel force and hull pressure distribution are considered. The bilge keel force is presented in terms of non-dimensional drag and mass coefficients including higher harmonic components. The numerical results are also compared with the industry standard empirical method for calculation of roll damping proposed by Ikeda et al. [4]. In general, a good agreement between the results of the present numerical method and the experimental data is obtained and the differences with those predicted by the empirical method are addressed.     

Numerical modelling of flow structures over idealized transverse aeolian dunes of varying geometry

A Computational Fluid Dynamics (CFD) model (PHOENICS™ 3.5) previously validated for wind tunnel measurements is used to simulate the streamwise and vertical velocity flow fields over idealized transverse dunes of varying height (h) and stoss slope basal length (L). The model accurately reproduced patterns of: flow deceleration at the dune toe; stoss flow acceleration; vertical lift in the crest region; lee-side flow separation, re-attachment and reversal; and flow recovery distance. Results indicate that the flow field over transverse dunes is particularly sensitive to changes in dune height, with an increase in height resulting in flow deceleration at the toe, streamwise acceleration and vertical lift at the crest, and an increase in the extent of, and strength of reversed flows within, the lee-side separation cell. In general, the length of the separation zone varied from 3 to 15 h from the crest and increased over taller, steeper dunes. Similarly, the flow recovery distance ranged from 45 to >75 h and was more sensitive to changes in dune height. For the range of dune shapes investigated in this study, the differing effects of height and stoss slope length raise questions regarding the applicability of dune aspect ratio as a parameter for explaining airflow over transverse dunes. Evidence is also provided to support existing research on: streamline curvature and the maintenance of sand transport in the toe region; vertical lift in the crest region and its effect on grainfall delivery; relations between the turbulent shear layer and downward forcing of flow re-attachment; and extended flow recovery distances beyond the separation cell. Field validation is required to test these findings in natural settings. Future applications of the model will characterize turbulence and shear stress fields, examine the effects of more complex isolated dune forms and investigate flow over multiple dunes.     

基于WRF和CFD软件结合的风能资源数值模拟试验研究

运用中尺度数值模式WRF与法国CFD软件MeteodynwT相结合的方法（WRF／WT）,进行了广东省海陵岛地区的水平分辨率100m×100m的风能资源数值模拟试验,采用海陵岛上7座测风塔观测资料对WRF／WT模式的模拟风场进行误差检验,并与WRF／WAsP模式系统对单点风能参数模拟误差进行对比,研究WRF／WT模式系统在风电场微观选址和分散式风电开发利用中应用的可行性。结果表明：中尺度模式与CFD软件结合的数值模拟方法对区域风能资源分布趋势的模拟比单纯应用CFD软件更准确;WRF／WT模式系统应用于复杂地形风能资源数值模拟评估是可行的,其对区域风能资源参数分布模拟的准确率与WRF／WAsP对2km范围内风能资源参数模拟的准确率相当;WRF／WT模式系统在风速频率分布不满足Weibull分布的情况下和陡峭地形条件下有较好的模拟效果,相对WRF／wAsP有明显优势。今后需进一步研究中尺度模式与CFD软件的衔接方法,以及对中尺度模式模拟结果的误差订正。     

多级可分离式AUV系统设计与仿真分析

多级可分离式自主水下航行器(MS-AUV)可搭载多种不同功能的载荷舱,并在目标海域完成载荷舱布放。其在海洋探测,海防等领域具有较高的应用价值。对MS-AUV进行了主体结构和分离机构的设计和研究,其中功能载荷舱与航行器本体采用真空连接方式,实现多级连接,降低了机构的复杂性,提高了载荷舱分离运动的稳定性和安全性。为了研究载荷舱与航行器分离的安全性,采用CFD和六自由度(6-DOF)刚体运动学的耦合方法来模拟其在水下的分离运动。在仿真过程中,采用弹簧光顺和局部重构混合网格生成方法的非结构化动态网格可以很好地适应大距离多体分离运动。仿真结果表明,在一定工况条件下,提高MS-AUV初始航行速度V0和分离机构弹簧刚度K有利于载荷与航行器本体分离。通过上述研究工作,验证了载荷分离方法的可行性,缩短研制周期,降低设计成本,对多级可分离式AUV样机的设计制造具有重要的参考价值。