Optimal ship tracking on a navigation route between two ports: a hydrodynamics approach

The optimal trajectory from Calcutta port to Mumbai port is charted for a tanker transshipping from the East coast to the West coast of India during rough weather. Rough weather is simulated over Indian seas using the state-of-the-art WAM numerical wave model (WAMDI Group in J Phys Oceanogr 18:1775–1810, 1988), assimilating satellite (IRS-P4) wind fields. These simulated wave fields and two-dimensional (2D) directional wave spectrum are an absolute representation of the irregular seaway. Hence, the same for the monsoon month of August 2000 formed the input basis for this study. Loss of ship speed due to the wave field (i.e., nonlinear motion of the tanker in waves) and associated sea-keeping characteristics in the seaway are estimated (Bhattacharya in Dynamics of marine vehicles, Wiley, New York, 1978). The approach adopted in this paper is unique in that it takes into account both voluntary and involuntary speed reductions of the ship. It helps in ship tracking by the optimum route using inverse velocity as the weight function for the path in an efficient way. Dijkstra’s algorithm [Numer Math 1(3):269–271, 1959] is applied in an iterative manner for determining the optimum track. The optimum track information has broad scope for use in modern shipping industry for obtaining safe and least-time routing by avoiding schedule delays with economic fuel consumption.     

Steady sailing performance of a hybrid-sail assisted bulk carrier

The steady sailing performance of a sail-assisted bulk carrier is investigated utilising towing-tank derived hydrodynamic derivatives and wind tunnel measured aerodynamic properties of the sails and the ship. The aerodynamic characteristics investigated include the ship hull at the fully-loaded draught, the sail–sail interaction effects for two sets of four identical hybrid-sails, and the sail–hull interaction effects for the same two sets of identical sails in the presence of the selected bulk carrier hull-form. This is in addition to lift–drag measurements of single isolated sails of each shape. The form of the two sets of soft sails was rectangular and triangular. This paper is concerned with assessing the benefits of a sail-assisted ship operation, and hence a steady-state rather than complete time-domain integrations of the governing equations are reported. The results of the completed analysis suggest that the benefits of the derived sail generated driving force are greater than the overhead of equipping the ship with a selected system of hybrid-sails. Sail-assisted ships could represent an important contribution to an improving global environment by reducing the demands for a driving force through the propeller.     

二维水翼非定常空泡流数值模拟

对非稳态空泡流现象的研究已成为当今一个研究重点。空泡非稳态脱落及其内部流动结构越来越受到研究人员的关注,但对其研究目前仍主要依靠实验手段。发展有效的数值模拟方法研究相关的机理具有重要的意义。本文从RANS方程入手,采用状态方程空化模型和计及可压缩性的SIMPLEC算法,实现了二维水翼CAV2003非定常空泡流动的数值模拟。研究了空泡从初生、发展、断裂,以及最终空泡在下游高压区溃灭的整个过程。计算表明,回射流并非完全液相的,而是汽液两相的。近壁面的回射流是导致空泡脱落的主要原因。最后,通过与第五届国际空泡会议上发表的7篇文献的对比,发现本文的计算结果是合理的。     

CFD simulation of 3-dimensional motion of a ship in waves: application to an advancing ship in regular heading waves

A new computational fluid dynamics simulation method has been developed for the unsteady motion of a ship advancing in waves. The objective is to evaluate the added resistance and predict the performance of a ship in waves. In this study, a finite volume method, in the framework of a boundary-fitted grid system, is employed. The motion of the ship is solved with six degrees of freedom by using the hydrodynamic forces and moments obtained from the solution of the simulation method. The marker–density–function method is employed to calculate the nonlinear free surface. This method is applied to the coupled motion problem of heaving and pitching. Received for publication on Nov. 15, 1999; accepted on Nov. 18, 1999     

粘性对方尾船波浪力的影响分析

航行船舶绕射问题的准确求解是预报船舶波浪中运动的基础.文章分别采用基于三维移动脉动源的边界元法和基于RANS方程的CFD方法求解了规则波中航行的方尾船模型的波浪力,并通过与试验结果比较验证了方法的可靠性.CFD计算时通过从波浪力中分离出粘性力,研究了粘性对方尾船波浪力的影响.结果表明,基于RANS方程的波浪力计算结果与试验结果和势流理论结果吻合良好,在长波中与试验值更接近.粘性对船模斜浪航行时的横摇力矩有较大影响.     

基于CFD的钻井船月池流场及附加阻力研究

研究了基于CFD方法的钻井船月池流场及阻力性能。针对某阶梯型月池的钻井船,采用定常方法分析其阻力性能,并与模型试验结果进行比较,分析钻井船月池阻力数值计算方法的可靠性;采用非定常方法模拟月池内流体的非定常流动,分析月池附加阻力产生的原因。计算方法及所得出的结论可为具有月池的船舶阻力性能研究提供参考。     

迎浪规则波中波浪增阻和船体垂向运动的数值预报(英文)

The numerical prediction of added resistance and vertical ship motions of one ITTC (International Towing Tank Conference) S-175 containership in regular head waves by our own in-house unsteady RANS solver naoe-FOAM-SL JTU is presented in this paper. The development of the solver naoe-FOAM-SJTU is based on the open source CFD tool, OpenFOAM. Numerical analysis is focused on the added resistance and vertical ship motions (heave and pitch motions) with four very different wavelengths (0.8Lpp≤λ≤1.5Lpp) in regular head waves. Once the wavelength is near the length of the ship model, the responses of the resistance and ship motions become strongly influenced by nonlinear factors, as a result difficulties within simulations occur. In the paper, a comparison of the experimental results and the nonlinear strip theory was reviewed and based on the findings, the RANS simulations by the solver naoe-FOAM-SJTU were considered competent with the prediction of added resistance and vertical ship motions in a wide range of wave lengths.     

The effect of sloshing on the sway motions of 2D rectangular cylinders in regular waves

In this paper, we investigated the effect of sloshing on the sway motions of two-dimensional rectangular cylinders in regular waves, bearing in mind possible applications for LNG-FPSO and LNG-FSRU. First, we carried out experiments for two models with different drafts, or the same draft but different filling ratios, in which the models were firmly connected to each other. The sway motion was measured with a noncontact video camera. This is an extension of Rognebakke and Faltinsen’s work for a single model (J Ship Res 47(3):208–221, 2003). It was found that the sway motion became small when the incident wave frequency was close to the lowest natural frequency of each model. The sway motion greatly increased when the wave frequency was higher than this frequency. The measured data were compared with numerical results obtained by a single-dominant multi-modal method; relatively good agreement was noted. However, the numerical results deviated from the experimental results near the lowest natural frequency of the smaller model, which was believed to be due to overturning waves, as observed during the experiment. Since this is out of the valid range for the single-dominant multi-modal method, other, more appropriate, methods such as the multi-dominant modal method must be applied instead.     

河口潮流垂线流速分布与床面阻力研究

本文从液体不稳定流微分方程出发,推导了河口均质潮流垂线流速分布公式,水流切力垂线分布公式以及床面阻力公式,并用实测资料作了验证,二者基本吻合。在此基础上,应用大量的实测资料统计分析,得到了反映潮汐河口床面阻力系数变化的表达式。     

桨舵干扰系统非定常水动力性能预估

基于对桨舵干扰的系统研究，本文提供了一个桨舵非定常干扰问题的理论计算方法，桨采用非定常涡格法，舵采用基于速度势的非定常面元法。在时域范围内通过迭代计算考虑螺旋桨和舵的相互影响。文中提供了一些算例考核，通过与实验结果对比，证明本文提供的方法是可靠的，可供实际应用。     

Time-domain TEBEM method for mean drift force and moment of ships with forward speed under the oblique seas

A numerical method for solving 3D unsteady potential flow problem of ship advancing in waves is put forward. The flow field is divided into an inner and an outer domain by introducing an artificial matching surface. The inner domain is surrounded by ship wetted surface and matching surface as well as part of the free surface. The free surface condition for the inner domain is formulated by perturbation about the double-body flow or uniform incoming flow assumption. The outer domain is surrounded by matching surface and the rest free surface as well as infinite far-field radiation boundary. The free surface condition for the outer domain is formulated by perturbation about uniform incoming flow. The simple Green function and transient free surface Green function are used to form the boundary integral equation (BIE) for the inner and outer domains, respectively. Taylor Expansion Boundary Element Method (TEBEM) is utilized to solve the double-body flow and inner domain and outer domain unsteady flow BIE. Matching conditions for the inner domain flow and outer domain flow are enforced by the continuity of velocity potential and normal velocity on the matching surface. Direct pressure integration on ship wetted surface is used to obtain the first-order and second-order wave forces (moments). The numerical predictions on the displacement, added resistance, sway mean drift force and yaw mean drift moment of the modified KVLCC2 ship at different forward speeds are investigated by the proposed TEBEM method. It is also compared with the other numerical results. The physical tank experiment results are also developed to validate the accuracy of numerical tank results. Compared with the experiment solutions, a good agreement can be obtained by TEBEM method.

     

Locomotion by mechanical pectoral fins

The objectives of this study were the development of a new device for maneuvering an underwater vehicle using the mechanism of a fish swimming, an experimental and theoretical analysis of the hydrodynamic characteristics of the device, and its application to maneuvering a fish robot. Observations and experimental analysis of the pectoral fins of a black bass (Micropterus salmoides) revealed that the locomotion of the fish, such as swimming forward at low speed, swimming backward, and turning in a horizontal plane is generated by using a combination of a feathering motion and a lead-lag motion of the pectoral fins. A mechanical pectoral fin making a feathering motion and a lead-lag motion generates a thrust force in a range of phase differences between both motions. The unsteady vortex lattice method, including the effect of viscosity, can express fairly well the unsteady forces acting on the mechanical pectoral fin in the range of phase differences where it exerts the thrust force. The fish robot, consisting of a model fish body and a pair of mechanical pectoral fins, can not only swim forward and turn in almost the same position, but can also swim in a lateral direction without changing the yaw angle. Translation of an article that appeared in the Journal of The Society of Naval Architects of Japan, vol. 182 (1997): The original article won the SNAJ prize, which is awarded annually to the best papers selected from the SNAJ Journal, JMST, or other quality journals in the field of naval architecture and ocean engineering.     

湍流入流下泵喷推进器力谱特性研究

[目的]旨在解决泵喷推进器在湍流入流作用下非定常力谱不明的问题.[方法]通过湍流生成栅格和频谱合成法相结合的方法,产生具有时间脉动和空间相干结构的湍流,然后将此种方法与模拟相结合,获得泵喷推进器非定常力的宽带谱.[结果]结果显示,所采用的计算方法能够获得泵喷推进器导管、定子和转子合理的非定常力宽带谱;泵喷推进器非定常力...     

潜伏式高速船型变航态自航数值模拟

针对一种具有半潜-水面航态变换功能的新型潜伏式高速船型,基于非稳态RANS方程组与VOF模型求解粘性流场,结合动态流域边界与滑移网格技术,对目标船在垂直面内三自由度运动及其喷水推进器叶轮随船定轴高速旋转的耦合运动进行数值模拟。结果表明：在相同主机平均转速条件下,所预报自航点体积傅氏数与自航模型试验结果的相对误差在7%以内;从半潜到水面的航态变换过程中,船体有显著的纵倾与升沉,船体主要克服压差阻力,船体兴波由相互叠加的首波、肩波、尾波逐渐发展为首尾分明的两个波系,在喷水推进器射流区内捕捉到生成和发展的高强度涡系。     

URANS and DES analysis for a Wigley hull at extreme drift angles

Vortical structures and associated instabilities for flows around the Wigley hull for a wide range of drift angles (10° ≤ α ≤ 60°) with free surface are identified and analyzed. Quantitative verification and validation are conducted on three systematically refined grids with comparison to the experimental data for α = 10°. Analysis of the flow pattern shows a strong correlation between the vortical structures and free-surface wave elevation. For α = 10° and 30°, the flows remain steady and vortices are generated at the keel and fore and aft perpendiculars of the hull. The strength and complexity of these vortices increase with increasing α. At α = 45°, flow becomes unsteady without any significant change in the main flow pattern. At α = 60°, a complex and unsteady flow field on the leeward side of the hull is formed with a large recirculation region from the aft to the fore end, which prevents the flow coming from below the keel from moving up and generating the keel vortices observed at lower drift angles. Karman-like and helical instabilities are analyzed. The effect of Froude number is more apparent for large than for small drift angles.     

A numerical investigation of the squat and resistance of ships advancing through a canal using CFD

As a ship approaches shallow water, a number of changes arise owing to the hydrodynamic interaction between the bottom of the ship’s hull and the seafloor. The flow velocity between the bottom of the hull and the seafloor increases, which leads to an increase in sinkage, trim and resistance. As the ship travels forward, squat of the ship may occur, stemming from this increase in sinkage and trim. Knowledge of a ship’s squat is necessary when navigating vessels through shallow water regions, such as rivers, channels and harbours. Accurate prediction of a ship’s squat is therefore essential, to minimize the risk of grounding for ships. Similarly, predicting a ship’s resistance in shallow water is equally important, to be able to calculate its power requirements. The key objective of this study was to perform fully nonlinear unsteady RANS simulations to predict the squat and resistance of a model-scale Duisburg Test Case container ship advancing in a canal. The analyses were carried out in different ship drafts at various speeds, utilizing a commercial CFD software package. The squat results obtained by CFD were then compared with available experimental data.     

非恒定流垂线流速分布规律的初探

从雷诺方程出发 ,通过对非恒定流的沿程阻力损失在垂线上分布规律的修正 ,可获得非恒定流和恒定非均匀流垂线流速分布公式。若定义均匀流垂线流速分布为A型 ,则非恒定流或恒定非均匀流垂线流速分布一般为B型和C型。B型分布在主流区 ,分布梯度大于A型 ;C型分布与A型相反 ,分布梯度一般为负     

用基于速度势的面元法预估导管桨的非定常性能

本文提供了一个计算导管螺旋桨的定常及非定常性能的基于速度势的非定常面元法,桨和导管的相互影响是通过在时域范围内迭代计算处理.关于桨泄出涡与导管泄出涡之间的非定常干扰在时域内匹配,本文推荐采用考虑它们非定常影响的简化处理方法,以减少计算时间,提高计算效率.算例计算结果表明,本方法可以有效地应用于导管桨的定常及非定常性能的预估.     

Numerical solutions for a two-dimensional airfoil undergoing unsteady motion

Continuous vorticity panels are used to model general unsteady inviscid, incompressible, and two-dimensional flows. The geometry of the airfoil is approximated by series of short straight segments having endpoints that lie on the actual surface. A piecewise linear, continuous distribution of vorticity over the airfoil surface is used to generate disturbance flow. The no-penetration condition is imposed at the midpoint of each segment and at discrete times. The wake is simulated by a system of point vortices, which move at local fluid velocity. At each time step, a new wake panel with uniform vorticity distribution is attached to the trailing edge, and the condition of constant circulation around the airfoil and wake is imposed. A new expression for Kutta condition is developed to study (i) the effect of thickness on the lift build-up of an impulsively started airfoil, (ii) the effects of reduced frequency and heave amplitude on the thrust production of flapping airfoils, and (iii) the vortex-airfoil interaction. This work presents some hydrodynamic results for tidalstream turbine.     

Effect of WO<Subscript>3</Subscript> nanoparticle loading on the microstructural,mechanical and corrosion resistance of Zn matrix/TiO<Subscript>2</Subscript>-WO<Subscript>3</Subscript> nanocomposite coatings for marine application

In this study, for marine application purposes, we evaluated the effect of process parameter and particle loading on the microstructure, mechanical reinforcement and corrosion resistance properties of a Zn-TiO₂-WO₃ nanocomposite produced via electrodeposition. We characterized the morphological properties of the composite coatings with a Scanning Electron Microscope (SEM) equipped with an Energy Dispersive Spectrometer (EDS). We carried out mechanical examination using a Dura Scan hardness tester and a CERT UMT-2 multi-functional tribological tester. We evaluated the corrosion properties by linear polarization in 3.5% NaCl. The results show that the coatings exhibited good stability and the quantitative particle loading greatly enhanced the structural and morphological properties, hardness behavior and corrosion resistance of the coatings. We observed the precipitation of this alloy on steel is greatly influenced by the composite characteristics.