On the numerical accuracy of the wave load distribution on a ship advancing in short and steep waves

Defining numerical uncertainty is an important part of the practical application of a numerical method. In the case of a ship advancing in short and steep waves, little knowledge exists on the solution behaviour as a function of discretisation resolution. This paper studies an interface-capturing (VOF) solution for a passenger ship advancing in steep (kA = 0.24) and short waves (L _w /L _pp = 0.16). The focus is to estimate quantitative uncertainties for the longitudinal distributions of the first–third harmonic wave loads in the ship bow area. These estimates are derived from the results of three systematically refined discretisation resolutions. The obtained uncertainty distributions reveal that even the uncertainty of the first harmonic wave load varies significantly along the ship bow area. It is shown that the largest local uncertainties of the first harmonic wave load relate to the differences in the local details of the propagating and deforming encountered waves along the hull. This paper also discusses the challenges that were encountered in the quantification of the uncertainties for this complex flow case.     

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

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.     

URANS simulations of catamaran interference in shallow water

This paper investigates the interference effects of wave systems on a multi-hull vessel in shallow water. A numerical analysis is made using the URANS code CFDSHIP-Iowa V.4 on the DELFT Catamaran model 372. The test matrix for numerical computations includes two separation distances (s = 0.17; 0.23) and the depth values of h/T = 8.2, 2.5 and 2, at several speeds ranging within Fr _H = 0.775–1.739. Numerical results are compared with the experimental data of the Bulgarian Ship Hydrodynamic Center, and verification and validation for resistance, sinkage and trim are also performed. Results show that, at critical speed (Fr _H ≈ 1), the presence of a finite depth significantly affects the catamaran total resistance, which, in shallower water, increases considerably with respect to deep water. At low h/T, small effects of the water depth on resistance occur at subcritical and supercritical speeds. The interference effects seem to be more relevant in shallow, rather than in deep water, with maximum IF values registered at critical speeds (Fr _H ≈ 1). Similarly to deep water, the lower the separation distance the greater the interference value. Moreover, in shallow water some negative interference is observed at Fr > 0.5. Wave patterns and wave profiles are analyzed and a comparison is made between several configurations of catamaran and a mono-hull vessel, in order to analyze how water depth and separation distance determine resistance and interference. Finally, a vortex instability study is also included.     

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.     

Prediction of wave parameters by using fuzzy inference system and the parametric models along the south coasts of the Black Sea

Forecasting of sea-state characteristics has a great importance in coastal and ocean engineering studies. Therefore, the purpose of this study was to investigate performances of Adaptive-Network-Based Fuzzy Inference System (ANFIS) and several parametric methods in the Black Sea. For this purpose, different fuzzy models with different input combinations were developed for two different wind data sources (TSMS and ECMWF) at two offshore buoy stations. It also aimed to apply several approaches to event-based data sets for wave predictions. Generally, in literature the tendency is to use time series data for wave predictions. In this kind of prediction approach, lagged time series data are taken as inputs and current or future variables are taken as output. In this study, event-based data for each independent storm were extracted from time series data. Simultaneous or concurrent data of wind speed, blowing duration, fetch length and wave characteristics were detected for each single storm. These event data were then used to set up models. The hindcast results were validated with significant wave height and mean wave period data recorded in Hopa and Sinop buoy stations. The performance of developed fuzzy models were also compared with that of four different parametric methods (Wilson, SPM, Jonswap, and CEM methods) applied for two wind data sources at both buoy stations. Finally, it was determined that in the prediction of both wave parameters (H _s and T _z) the ANFIS models (R = 0.66, squared correlation coefficient, and MAE = 0.37 m, mean absolute error, for the best model in prediction of H _s) were more accurate than the parametric methods (R = 0.63 and MAE = 0.75 m for the best model in prediction of H _s).     

Optimization of bow shape for a non ballast water ship

In this research,a commercial CFD code ＂Fluent＂ was applied to optimization of bulbous bow shape for a non ballast water ships（NBS）.The ship was developed at the Laboratory of the authors in Osaka Prefecture University,Japan.At first,accuracy of the CFD code was validated by comparing the CFD results with experimental results at towing tank of Osaka Prefecture University.In the optimizing process,the resistances acting on ships in calm water and in regular head waves were defined as the object function.Following features of bulbous bow shapes were considered as design parameters： volume of bulbous bow,height of its volume center,angle of bow bottom,and length of bulbous bow.When referring to the computed results given by the CFD like resistance,pressure and wave pattern made by ships in calm water and in waves,an optimal bow shape for ships was discovered by comparing the results in the series of bow shapes.In the computation on waves,the ship is in fully captured condition because shorter waves,λ/Lpp 0.6,are assumed.     

Resistance analysis of unsymmetrical trimaran model with outboard sidehulls configuration

The application of multi-hull ship or trimaran vessel as a mode of transports in both river and sea environments have grown rapidly in recent years.Trimaran vessels are currently of interest for many new high speed ship projects due to the high levels of hydrodynamic efficiency that can be achieved,compared to the mono-hull and catamaran hull forms.The purpose of this study is to identify the possible effects of using an unsymmetrical trimaran ship model with configuration（S/L） 0.1-0.3 and R/L=0.1-0.2.Unsymmetrical trimaran ship model with main dimensions： L=2000mm,B=200 mm and T=45 mm.Experimental methods（towing tank） were performed in the study using speed variations at Froude number 0.1-0.6.The ship model was pulled by an electric motor whose speed could be varied and adjusted.The ship model resistance was measured precisely by using a load cell transducer.The comparison of ship resistance for each configuration with mono-hull was shown on the graph as a function of the total resistance coefficient and Froude number.The test results found that the effective drag reduction could be achieved up to 17% at Fr=0.35 with configuration S/L=0.1.     

Comparison of CFD and EFD for the Series 60 C B= 0.6 in steady drift motion

 This paper presents comparisons of computational and experimental fluid dynamics results for boundary layers, wakes, and wave fields for the Series 60 C _B= 0.6 ship model in steady drift motion. The numerical method solves the unsteady Reynolds-averaged Navier–Stokes and continuity equations with the Baldwin–Lomax turbulence model, exact nonlinear kinematic and approximate dynamic free-surface boundary conditions, and a body/free-surface conforming grid. The experimental and computational conditions, i.e., Froude numbers of 0.16 and 0.316 for the experiments, and Froude numbers of 0 and 0.316 for the computations, allow comparisons of low and high Froude number results, respectively, which allows an evaluation of Froude number effects and validation of the computational fluid dynamics at both low and high Froude numbers. This article gives an overview of this numerical approach, and the computational conditions and uncertainty analysis are described. Results are presented for the wave and flow fields, with emphasis on the important flow features of drift- and wave-induced effects in comparison with the experiments. Finally, conclusions from the present study are given, together with recommendations for future work. Received: August 31, 2001 / Accepted: March 25, 2002     

Global wave loads on intact and damaged Ro-Ro ships in regular oblique waves

A nonlinear time-domain simulation method is presented for the prediction of dynamic global wave loads on a Ro-Ro ship at zero speed in regular oblique waves in an intact and a damaged condition. Numerical computations and model tests have been carried to investigate the structural responses of Ro-Ro ship Dextra to various wave amplitudes at three different wave headings (DTR-4.1-NEW-12.98, DEXTREMEL project BE97-4375, 1998; DTR-4.2-NEW-11.99, DEXTREMEL project BE97-4375, 1999). The results of numerical and experimental investigations for stern quartering waves are reviewed. Comparisons between predictions and measurements for global wave loads at the midship section of the intact and the damaged Ro-Ro ship show that the agreement between the theory and experiment for dynamic horizontal and vertical bending moments is excellent. On the other hand, correlation between the predictions and measurements for dynamic vertical shear force is better than that for dynamic horizontal shear force. Nevertheless, the calculated torsion moment values are higher than the measured values. As the wave amplitude is not small, the positive and negative peaks of global wave loads are no longer equal to each other as found in both the calculations and experiments. The dynamic vertical global wave loads in the damaged condition are larger than that in the intact condition.     

Comparison of Microbubble and Air Layer Injection with Porous Media for Drag Reduction on a Self-propelled Barge Ship Model

Ship resistance issues are related to fuel economy, speed, and cost efficiency. Air lubrication is a promising technique for lowering hull frictional resistance as it is supposed to modify the energy in the turbulent boundary layer and thereby reduce hull friction. In this paper, the objective is to identify the optimum type of air lubrication using microbubble drag reduction (MBDR) and air layer drag reduction (ALDR) techniques to reduce the resistance of a 56-m Indonesian self-propelled barge (SPB). A model with the following dimensions was constructed: length L?=?2000 mm, breadth B?=?521.60 mm, and draft T?=?52.50 mm. The ship model was towed using standard towing tank experimental parameters. The speed was varied over the Froude number range 0.11–0.31. The air layer flow rate was varied at 80, 85, and 90 standard liters per minute (SLPM) and the microbubble injection coefficient over the range 0.20–0.60. The results show that the ship model using the air layer had the highest drag reduction up to a maximum of 90%. Based on the characteristics of the SPB, which operates at low speed, the optimum air lubrication type to reduce resistance in this instance is ALDR.     

Experimental investigation on scour development and scour protection for offshore converter platform

Local scour around the offshore converter platform caused by current and wave will lead to the instability of foundation. A series of experiments are performed to investigate the scour development and scour protection of the offshore converter platform. The development and evolution characteristics of the local scour around the gravity based structure foundation under the actions of current flow and wave are analyzed. The results show that the edge scour mainly occurs at the lateral sides of the platform and the scour pits are symmetrical about the centerline of the platform in the streamwise flow direction. The scour depth around the gravity based structure increases with the increase of flow velocity. The inclination and collapse of the platform is observed at α = 90° and U₀ = 0.25 m/s. The sand ripples are observed around the gravity based structure under the wave action, and the scour depth increases with the increase of wave height and wave period. Based on the comparisons of different scour protection methods, the combined bionic grass-geotextile-riprap protection has the best scour protection efficiency for the offshore converter platform.     

Breaking wave at the bow of a fast displacement ship model

We investigated the flow structures under the bow wave generated by a fast displacement ship model (INSEAN model 2340) in the presence of wave breaking. The data acquired were also used for a detailed database for CFD validation. The mean and r.m.s. point-wise values of the wave height were measured by means of a finger probe. The intensity of the breaking wave was taken as the r.m.s value of the wave height. The mean velocity field under the free surface was measured at 0.15 L _PP and 0.2 L _PP downstream of the fore perpendicular by means of a 5-hole Pitot probe. Uncertainty assessment of the wave height and velocity field results was performed following the AIAA Standards S-071-1995. Preliminary CFD results from a RANSE code with a breaking model are shown in comparison with the measured data.     

Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids

Unsteady Reynolds averaged Navier–Stokes (URANS) computations of standard maneuvers are performed for a surface combatant at model and full scale. The computations are performed using CFDShip-Iowa v4, a free surface solver designed for 6DOF motions in free and semi-captive problems. Overset grids and a hierarchy of bodies allow the deflection of the rudders while the ship undergoes 6DOF motions. Two types of maneuvers are simulated: steady turn and zigzag. Simulations of steady turn at 35° rudder deflection and zigzag 20/20 maneuvers for Fr = 0.25 and 0.41 using constant RPM propulsion are benchmarked against experimental time histories of yaw, yaw rate and roll, and trajectories, and also compared against available integral variables. Differences between CFD and experiments are mostly within 10 % for both maneuvers, highly satisfactory given the degree of complexity of these computations. Simulations are performed also with waves, and with propulsion at either constant RPM or torque. 20/20 zigzag maneuvers are simulated at model and full scale for Fr = 0.41. The full scale case produces a thinner boundary layer profile compared to the model scale with different reaction times and handling needed for maneuvering. Results indicate that URANS computations of maneuvers are feasible, though issues regarding adequate modeling of propellers remain to be solved.     

Added resistance acting on hull of a non ballast water ship

In this paper, added resistances acting on a hull of non ballast water ship(NBS) in high waves is discussed. The non ballast water ships were developed at the laboratory of the authors at Osaka Prefecture University, Japan. In the present paper, the performances of three kinds of bow shapes developed for the NBS were theoretically and experimentally investigated to find the best one in high waves. In previous papers, an optimum bow shape for the NBS was developed in calm water and in moderated waves. For a 2 m model for experiments and computations, the wave height is 0.02 m. This means that the wave height is 15% of the draft of the ship in full load conditions. In this paper, added resistances in high waves up to 0.07 m for a 2 m model or 53% of the full load draft are investigated. In such high waves linear wave theories which have been used in the design stage of a ship for a long time may not work well anymore, and experiments are the only effective tool to predict the added resistance in high waves. With the computations for waves, the ship is in a fully captured condition because shorter waves, λ/Lpp0.6, are assumed.     

破碎波对直立堤作用的实验研究(英文)

Breaking waves can have tremendous destructive impact on vertical walls,yet they are poorly understood.By using particle imaging velocimetry(PIV) technology and high-precision pressure transducers,actual breaking wave loads on vertical walls were studied.By simultaneously comparing the flow field structure and wave pressure,the mechanisms of breaking wave pressure could be analyzed.The probability distribution of the peak value of the first impact of a breaking wave was investigated.The results showed that the impact pressure p is mainly distributed in the range of 0.25-2.75 ρν 2,with the greatest possible probability at p/ρ v2= 0.75.     

Research on the dynamic buckling of a typical deck grillage structure subjected to in-plane impact Load

The dynamic buckling of the main deck grillage would result in the total collapse of the ship hull subjected to a far-filed underwater explosion. This dynamic buckling is mainly due to the dynamic moment of the ship hull when the ship hull experiences a sudden movement under impact load from the explosion. In order to investigate the ultimate strength of a typical deck grillage under quasi-static and dynamic in-plane compressive load, a structure model, in which the real constrained condition of the deck grillage was taken into consideration, was designed and manufactured. The quasi-static ultimate strength and damage mode of the deck grillage under in-plane compressive load was experimentally investigated. The Finite Element Method (FEM) was employed to predict the ultimate strength of the deck grillage subjected to quasi-static in-plane compressive load, and was validated by comparing the results from experimental tests and numerical simulations. In addition, the numerical simulations of dynamic buckling of the same model under in-plane impact load was performed, in which the influences of the load amplitude and the frequency of dynamic impact load, as well as the initial stress and deflection induced by wave load on the ultimate strength and failure mode were investigated. The results show that the dynamic buckling mode is quite different from the failure mode of the structure subjected to quasi-static in-plane compressive load. The displacements of deck edge in the vertical direction and the axial displacements are getting larger with the decrease of impact frequency. Besides, it is found that the dynamic buckling strength roughly linearly decreased with the increase of initial proportion of the static ultimate strength P₀. The conclusions drawn from the researches of this paper would help better designing of the ship structure under impact loads.     

Shallow water effects on high order statistics and probability distributions of wave run-ups along FPSO broadside

Ship hydrodynamics in shallow water becomes especially complicated since the nonlinearities in both the incident waves and the wave–hull interactions will be affected by the water depth. For a ship-shaped Floating Production, Storage and Offloading unit (FPSO) operating in shallow water, the broadside often suffers from the wave run-up and green water incidents in non-collinear harsh ocean environments. By applying the methods of ordinary moments and L-moments and the empirical Weibull distribution on the data measured in a series of model experiments, the high order statistics and the exceedance probability distribution of the run-ups along the FPSO broadside are evaluated and the effects of the shallow water depth and the incident environments are analyzed in this paper. It is seen that both the incident waves and the wave run-ups are non-Gaussian in shallow water and that the wave run-up characteristics are significantly influenced by the water depth and the incident environments, while the contribution due to the vessel vertical motions is negligible for the FPSO used in this study. The exceedance probabilities of the wave run-ups show that the broadside will be more likely to suffer from serious wave run-up and green water incidents in shallower water, in a higher incident wave and a non-collinear environment, especially so at locations around the FPSO midship within a range of 3/8L_pp ∼ 5/8L_pp. The dependency of the shape and scale parameters of the wave run-up probability distributions on the locations and the environment is quantified by model tests. The present study leads to the conclusion that the wave run-up characteristics and the shallow water effects should be considered carefully in determining the wave loads and the freeboard of a large FPSO in non-collinear environment conditions.     

Large-eddy simulation of cross-flow around ship sections

The present work is motivated by phenomena occurring in the flow field around ship-like bodies with an incoming lateral flow (cross-flow, 90 ° drift angle). Three-dimensional unsteady flows around different ship sections are investigated by using computational fluid dynamic (CFD) tools with large-eddy simulation (LES) subgrid-scale turbulence model. The simulation results are compared to measurements at several Reynolds numbers in the 90–200,000 range. Focus in our investigation is on the characterization of the motion of vortex structures generated by the separated flow. Another target in the study is to obtain better knowledge of the hydrodynamic forces acting on the sections. Computed pressure and drag coefficients are compared with experimental measurements. The comparison between simulations and measurements shows that an LES model can predict the flow field around ship sections in detail and the hydrodynamic forces acting on the sections.     

Computational predictions of ship-speed performance

This paper examines ship-speed performance based on acomputational method. The computations are carried out under identical model conditions, i.e., resistance and self-propulsion tests, to predict the speed-power relationship. The self-propulsion point is obtained from the self-propulsive computational results of two propeller rotative speeds. The speed-power relationship in full scale is obtained through analyzing the computational results in model scale according to the model-ship performance analysis method of ITTC’78. The object ship is a VLCC. The limiting streamlines and the distribution of the pressure coefficient on the hull, the wake characteristics on the propeller plane, and the wave characteristics around a model ship are also investigated. After completing the computations, a series of model tests are conducted to evaluate the accuracy of the predictions by comparing the computational results with the experimental results.     

An insight into the slamming behaviour of large high-speed catamarans through full-scale measurements

The slamming behaviour of a large high-speed catamaran has been investigated through the analysis of full-scale trials data. The US Navy conducted the trials in the North Sea and North Atlantic region on a 98 m wave piercer catamaran, HSV-2 Swift, designed by Revolution Design Pty Ltd and built by Incat Tasmania. For varying wave headings, vessel speeds and sea states the data records were interrogated to identify slam events. An automatic slam identification algorithm was developed, considering the measured rate of change of stress in the ship’s structure coupled with the vessel’s pitch motion. This has allowed the slam occurrence rates to be found for a range of conditions and the influence of vessel speed, wave environment and heading to be determined. The slam events have been further characterised by assessing the relative vertical velocity at impact between the vessel and the wave. Since the ship was equipped with a ride control system, its influence on the slam occurrence rates has also been assessed.