Advanced
A RANS modelling approach for predicting powering performance of ships in waves
A RANS modelling approach for predicting powering performance of ships in waves
International Journal of Naval Architecture and Ocean Engineering. 2014. Jun, 6(2): 418-430
Copyright © 2014, The Society of Naval Architects of Korea
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Published : June 30, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Björn Windén
Stephen Turnock
Dominic Hudson
Abstract
In this paper, a modelling technique for simulating self-propelled ships in waves is presented. The flow is modelled using a RANS solver coupled with an actuator disk model for the propeller. The motion of the ship is taken into consideration in the definition of the actuator disk region as well as the advance ratio of the propeller. The RPM of the propeller is controlled using a PID-controller with constraints added on the maximum permissible RPM increase rate. Results are presented for a freely surging model in regular waves with different constraints put on the PID-controller. The described method shows promising results and allows for the studying of several factors relating to self-propulsion. However, more validation data is needed to judge the accuracy of the model .
Keywords
INTRODUCTION
The performance of a ship when travelling in a seaway has always been a hot topic. This is true now more than ever with ever increasing requirements for environmental performance and rapid transit of goods. The introduction of the Energy Efficiency Design Index (EEDI) as a design constraint for new ships from 2013 ( IMO, 2011 ) is likely to lead to an increased need for design modifications that lead to improved efficiency without negatively impacting other factors such as transport capacity, safety, etc. Accurately predicting performance in waves and being able to judge the impact of subtle design changes is therefore identified as an area of high interest. From a hydrodynamical point of view, understanding the interaction between the ship, the sea state that it is in and the propeller that is powering it will give an answer to most questions regarding performance. A schematic of this interaction created by Windén (2012) is shown in Fig. 1 .
PPT Slide
Lager Image
Hydrodynamics related to performance in waves.
To properly model this interaction, strong linking is needed between flow conditions around the hull and the performance of the propeller. Currently, the only method that allows for detailed flow features around the hull and the propeller to be modelled is a Navier-Stokes (NS) based solver.
APPROACH
One reason why proper verification and validation of NS based modelling of ships in waves is hard to achieve is the complexity of the model itself. For self-propulsion, a way to model the rotating propeller has to be added on top of the already complex dynamical geometry. This leads to a difficulty in identifying the most important sources of error. For example incorrect phasing of forces and motions has been identified as a problem in recent workshops ( Larsson et al., 2014 ). If the motions are seen to be out of phase, it is hard to tell if this is due to incorrect forcing, a problem with the mesh deformation scheme or the coupling scheme between forces and motions. This means that the predicted propulsive performance of the ship has many associated error sources and it may be hard to clearly quantify the underlying reason for the error.
One of the most well adopted ways to address a problem with many unknowns is to start with something that is known or can be found out easily and gradually build understanding from this. Applying this reasoning to the problem of modelling ship performance in waves means removing many of the uncertainties mentioned above. RANS is chosen as a modelling technique due to its advantage in computational efficiency over other NS based methods. A validation study based on experiments by Journée (1992) of the forces and moments on a fixed Wigley hull in waves predicted using the open source CFD software OpenFOAM has been carried out by Windén et al. (2012) . Fig. 2 shows encouraging results for the prediction of the amplitude of the forces at different wavelengths but a problem with the phasing still remains even when the hull is fixed. The study also compares results with a nonlinear Boundary Element Method (BEM) by Kjellberg (2011) which seemed to predict the phasing slightly better while giving a greater error in predicting the amplitudes of the forces.
PPT Slide
Lager Image
Results of validation case.
The next step after the forces and moments are validated is either to study the motions or to study propulsion. In this paper, the second issue is addressed. The objective is to identify more about what factors should be considered when modelling selfpropelled ships in waves using RANS. In the simulations in this paper, the model is free to surge and the RPM of the propeller is controlled using a PID-controller. This means that the influence of the propeller thrust on the surge for different control conditions can be studied.
GEOMETRY AND MESHING
The setup of the numerical towing tank is shown in Fig. 3 . The size of the computational domain is set so that Lf = L and La = 4 L . Waves are generated and dissipated using the relaxation-based wave generation toolbox waves2Foam ( Jacobsen et al., 2012 ). The relaxation zones are designed so that Lg = Ld = 2 λ . Forward speed is modelled using a steady current which is set corresponding to Fn = 0.2 . The waves are chosen to have an amplitude of ζ =0.023m which corresponds to about 12% of the model draught. The wavelength is set to be equal to the model length so that λ = L .
PPT Slide
Lager Image
Geometry of the numerical towing tank.
The mesh is created using an automated refinement strategy described by Windén et al. (2012) . The final mesh size lies around 6M cells with a y + value of around 2. The mesh is designed to ensure
  • • an adequate distribution of cells in the free surface region
The particulars of the chosen hull are shown in Table 1 .
Particulars of tested hulls.
PPT Slide
Lager Image
Particulars of tested hulls.
where a 2 gives the shape of the hull as
PPT Slide
Lager Image
PROPELLER MODELLING
The Wigley hull is a purely academical hullform and has no consideration of propeller positioning in the definition of its shape. The propeller disc is therefore positioned half a radius behind the aft perpendicular, with a lower point flush with the keel and with a diameter of 85% of the draught. This is seen as representative of a typical positioning of the propeller on a ship. The positioning of the propeller disc relative to the rear profile of the used Wigley hull is shown in Fig. 4 . Furthermore a hub with a radius of 10% of the total disk radius is used as well as a shaft extending into the hull and ending in an ellipsoid cap to minimise separation.
PPT Slide
Lager Image
Positioning of propeller disc behind Wigley hull.
The propeller is modelled using a body force approach, meaning that an external force is added to the RANS momentum equation to represent the work done on the fluid by the propeller. The body force due to the propeller Fv is defined as a force per unit volume and is incorporated in the momentum equation as
PPT Slide
Lager Image
In Eq. (2), the left hand side represents the momentum change of the fluid. The average velocity components are denoted as ui with i =1,2,3. The right hand side of Eq. (2) represents stresses where the body force vector Fv is added to the pressure gradient, the viscous stresses and the Reynolds stress tensor τij which is modelled using the turbulence model k ω SST ( Menter et al., 2003 ).
- Active cell identification
The magnitude and direction of Fv in every cell in the domain is decided by the location of that cell in relation to the propeller to be modelled. Cells within the propeller disk will be active (i.e. be eligible to receive an Fv >0. To not perform unnecessary calculations, these cells must be identified before any body force is added. Firstly, the centroid of the propeller in the undisturbed state is denoted xp 0 , the initial centre of gravity of the hull CG 0 and the orientation of the propeller axis in relation to the initial state O 0 . An arbitrary rotation and translation of the initial state is shown if Fig. 5 . The rotation tensor due to the rotation of the hull in three degrees of freedom (pitch, roll and yaw) is denoted as Q so that
PPT Slide
Lager Image
and the planar shift of the centre of gravity due to motions is denoted P . With this notation, the location of the propeller centroid xp at an arbitrary position and orientation of the hull is
PPT Slide
Lager Image
and the orientation of the propeller axis O (normal to the propeller disk plane) is
PPT Slide
Lager Image
PPT Slide
Lager Image
Arbitrary rotation and translation of initial state.
A cell is given a non-zero value of Fv if it lies within the radius of the propeller disk R but outside of the hub radius RH and within the thickness of the disk d . A vector vI is defined as going between the current propeller centroid xp and the centre of cell I . The distance from the propeller centroid along the propeller normal axis (
PPT Slide
Lager Image
) is found from
PPT Slide
Lager Image
The local radius of cell I , RI is found by projecting the vector vI onto the propeller plane. Since the normal vector to the propeller plane is the vector O = ( O 1 , O 2 , O 3 ) the projection of vI onto the propeller plane is achieved by
PPT Slide
Lager Image
The definitions of RI , dI and vI are shown in Fig. 6
PPT Slide
Lager Image
Definitions of propeller disk extent.
From this definition it follows that
PPT Slide
Lager Image
- Thrust and torque distribution
The body force exerted on a certain cell within the propeller disk can be calculated in several ways. The most accurate representation is to model the lift and drag of the blades along radius by separating them into finite elements, each having a certain foil shape and attitude to the oncoming flow. This is known as Blade Element Momentum theory (BEMt.) Simplifications to this can be made by allowing the blades to sweep over a certain distance to give average values in larger sectors of the disk. A further step is to assume that the propeller rotates fast enough for the thrust and torque to be smeared over the entire disk at any given time step. The approach to modelling the blades depends on the amount of information is needed about the instantaneous effect of the blades on the surrounding fluid. This must be put in relation to the time scales of other phenomena around the propeller such as vortices separating from the stern area and how well the interaction of these with the propeller is to be modelled.
In this case the frequency of the waves passing and thus the variations in resistance and surge are deemed most important. The relation between the passing period of wave crests e T at different λ / L in relation to different rotational frequencies 1/ Tr for the propeller are shown in Fig. 7 . The rotational frequency is given as the full scale RPM found as RPMfull = α −0.5 RPMmod el where α is taken as an example scale factor of 40 ( Lpp of full scale ship = 120 m when using the model described in Table 1 ).
PPT Slide
Lager Image
Contours of the ratio Te / Tr .
As seen in Fig. 7 , there will be more than 4 rotations of the propeller for every wave encounter in most situations. For low values of λ / L , the wave encounter frequency is significant. However, as seen in Fig. 2 , the added forces due to the waves at these wavelengths are low. This means that the motions will be low and that the main simplification made by using a smeared distribution will be to ignore the periodic variation of inflow velocity due to the waves. This is something that will be true regardless of the wavelength. For the example used in this paper of λ / L = 1 , there will be 6-12 rotations for each wave encounter. This is considered enough for a smeared approach to be applied to the thrust and torque distributions. The thrust and torque distribution inside the propeller disk is determined by a radial shape function which is zero at the hub and tip and approximately follows the Goldstein (1929) optimum distribution. A non-dimensional radius rs is defined as
PPT Slide
Lager Image
A shape function describing the thrust distribution is given as
PPT Slide
Lager Image
and equally for the torque
PPT Slide
Lager Image
The shape functions fK and fQ are shown in Fig. 8 .
PPT Slide
Lager Image
Shape functions fK and fQ .
The shape functions are normalized so that the summation value over the entire disk corresponds to the desired thrust and torque
PPT Slide
Lager Image
PPT Slide
Lager Image
where FK is the thrust and FQ is the torque. γ is a unit vector describing the direction of the circumferential force in each cell and is defined as
PPT Slide
Lager Image
The direction of rotation can be altered by replacing Eq. (14) with
PPT Slide
Lager Image
The body force for a cell I is now given as
PPT Slide
Lager Image
The thrust and torque are decided by the coefficients of thrust and torque KT and KQ as well as the current working RPM.
PPT Slide
Lager Image
PPT Slide
Lager Image
Since the propeller blades are not explicitly modelled in this case, the variation of KT and KQ with the RPM is not known. For stock propellers such as the Wageningen B-series, these can be found by interpolation of tabulated values based on experimental data for example from Bernitsas et al. (1981) . Here, a 5th order polynomial fit to the open water characteristics of a MARIN 7967 propeller as described by Carrica et al. (2013) is used where
PPT Slide
Lager Image
PPT Slide
Lager Image
In this case, the inflow speed is taken as the base forward speed U 0 plus the velocity of the propeller centroid relative to the moving reference frame. This gives the advance coefficient J as
PPT Slide
Lager Image
In Eq. (21), care must be taken with the signs of U 0 and ∂xp / ∂t so that they, with the coordinate system used, have the same definition of positive speed.
Eqs. (17) and (18) as well as the shape functions (Eqs. (10) and (11)) assume that the flow into the propeller is uniform so that the local advance ratio is constant along the circumference. If the propeller is working behind a hull, the wake created by the ship will distort this condition. The thrust and torque varies periodically because of the surging motion through Eq. (21). However, the presence of the waves would also mean a periodic variation of thrust and torque due to an unsteady wake, something that is not considered here. By applying a BEMt method post run-time to estimate the thrust and torque variations behind the Wigley hull in this study, Windén et al. (2013) found a thrust variation of 16% around the mean value due to the unsteady wake.
PROPELLER CONTROL
Since the thrust produced by the propeller is a function of the RPM, it is possible vary the RPM based on given criteria for the thrust. In the interest of maintaining a constant service speed in waves, the RPM can be varied as to maintain the thrust needed to keep constant forward speed. Since the undisturbed forward speed is modelled with a current in this case, keeping constant speed is equivalent with keeping zero surge motion.
This can be achieved in two ways. The most direct approach is to iterate in each time step to find the thrust that will ensure force balance in the surge direction. The other alternative is to use a control function to steer the RPM towards the correct value based on information available from previous time steps. The second approach is preferable here, both in terms of computational time but, more importantly because it more closely represents a situation that could be applied on a real ship. Because of the nonlinear properties of marine propulsion systems, to achieve optimal propulsion performance it is necessary to apply some sort of control algorithm to govern the propeller RPM ( Xiros, 2002 ). This is particularly important when ship dynamics are also considered as a factor ( Kashima and Takata, 2002 ). A similar approach as the one described in this paper with a body force propeller and PID regulator for the RPM was applied by ( Carrica et al., 2008 ) for studying broaching events.
- PID-controller
The PID controller is defined by letting the RPM in a time step be a function of the RPM in the previous time step and a correction value based on an error e , the rate of change of the error and the integral value of the error since the start of the simulation.
PPT Slide
Lager Image
An appropriate definition of the error e must be chosen to link the thrust from the propeller to the surge movement of the hull. Care is taken here to only include variables that could be feasibly measured on a real ship. For example, the hydrodynamic drag as well as the thrust is available in the RANS-solver and the difference between them could make up the error e . However, neither the thrust nor the resistance can be feasibly measured on board a ship in real time. For this reason, the acceleration of the centre of gravity is chosen as the main control variable. The acceleration can be measured in real time by placing an accelerometer at an arbitrary location provided the position of that location relative to the centre of gravity is known. The error e is thus based on the acceleration in the direction of travel ( a x 1 ) .
To accentuate the connection with the resistance, the acceleration is multiplied by the mass of the ship to give an estimate of the current force acting on the hull. Furthermore, to avoid motion with constant velocity (in this case meaning that the hull moves with constant velocity other than the one specified) a penalty of 1000 is put on the integral value of acceleration. The error is now defined as
PPT Slide
Lager Image
- RPM limiting
In addition to the PID-controller governing the RPM, some restrictions must be applied to make the simulation more realistic. The RPM should not be allowed to exceed a certain value since that would mean a great risk of cavitation and damage to machinery in a real propulsion system. Furthermore, large acceleration of the propeller rotation should also be disallowed since the torque needed to achieve such acceleration lies beyond the capability of both engines and propeller shafts. Furthermore, RPM increases are usually limited in marine engines to avoid large heat gradients etc. A typical maximum value for how fast the propeller rotation is allowed to change is around 1 RPM per minute in full scale. If the example scale factor of 40 is applied in this case, this is equivalent to a maximum increase of 0.67 RPM/s in model scale. Furthermore, typical working shaft-RPM for most large merchant vessels lie between 40 and 100 RPM. This corresponds to between 250 and 630 RPM in the example model scale.
RESULTS
Simulations are started from an initialised solution where the steady wave pattern of the hull is allowed to develop. At the initial state, the resistance is recorded and the starting RPM is set as to match this resistance with a corresponding thrust. The hull is then subjected to regular waves with the first crest reaching the hull about t / Te = 3.
PID coefficients and initial RPM value.
PPT Slide
Lager Image
PID coefficients and initial RPM value.
All the information about the propeller extent and position are shown in Table 3
Positioning and extent of propeller.
PPT Slide
Lager Image
Positioning and extent of propeller.
Nine different values on the maximum permissible RPM change rate ∂RPM / ∂t are tested as well as a reference case where the acceleration is unlimited. The tested constraints are listed in Table 4 .
RPM constraints.
PPT Slide
Lager Image
RPM constraints.
Controllers 2-5 performs similarly in terms of surge motion compared to controller 1, significant differences are only found for controllers 6-10. In Figs. 9 , 10 and 11, the forward speed, surge and RPM for controllers 1, 7 and 10 are shown. Finally, the delivered power, which can be calculated as PD = 2 π FQRPM /60 (Molland et al., 2011) is shown in Fig. 12 .
PPT Slide
Lager Image
Changes in forward speed due to waves with different values of ∂RPM / ∂t allowed.
PPT Slide
Lager Image
Changes in surge due to waves with different values of ∂RPM / ∂t allowed.
PPT Slide
Lager Image
Development of propeller RPM with different values of ∂RPM / ∂t allowed.
PPT Slide
Lager Image
Changes in delivered power with different values of ∂RPM / ∂t allowed.
As seen in Fig. 9, the hull experiences a periodical oscillation in forward speed due to the waves. There is also a drop in the mean velocity. This is due to the mean increase in resistance due to waves RAW . It is clear that none of the controller constraints are generous enough to influence the periodical velocity notably. The RPM in all cases follows the maximum permitted increase. However, the higher the permitted value of ∂RPM / ∂t is, the faster the ship can overcome RAW with an increased thrust and return to the original forward speed. This however comes at the cost of a higher power delivered to the propeller as seen in Fig. 12 . Controller 8 is very successful in keeping constant forward speed, however the RPM increases (and subsequent power) required to do so are not in any way related to reality which is why controller 8 is excluded from the results figures.
As seen in Figure 10 , even though the hull returns to its initial forward speed, the integral value of the velocity means that it has been given a large distortion in surge. While this is not a problem in real applications, in this case it means that the mesh has been distorted. In the interest of accuracy future implementations should take this into consideration.
If the actual powering performance in waves is sought, it would not be suitable to construct a controller that puts a penalty on the surge drift. Instead it would be preferable to initialize the simulation from a fixed hull in waves and use the RPM needed to overcome RAW as the initial condition.
CONCLUSIONS
The controller function described in this paper is a simplified version which mostly serves to demonstrate the opportunities of using open source RANS modelling with external inputs to more closely represent real operating conditions. Much more advanced controlling algorithms could be employed here, see for example Xiros (2002) . Furthermore, the actuator disk model has many limitations. To properly include the propeller geometry and to draw more detailed conclusions about the effects of unsteady loading, a BEMt approach is preferred. The future aim for this model is to replace the actuator disk formulation with the BEMt model developed by Phillips et al. (2009 ; 2010) .
This study has shown that it is possible to model and discuss around several aspects of self-propulsion using the described method. There is however need for proper validation against experimental data to ensure the accuracy of the model. The benefit of using RANS CFD to do this rather than experimental methods is that the entire flow field can be accessed and analysed to study the influence of e.g. changes in bow and stern shape and the subsequent change in performance. The future goal is to apply the methodology described in this paper to more realistic hull forms to be able to conduct detailed studies on more of the factors affecting ship performance in waves.
Finally, the conclusions drawn regarding the controller constraints may become relevant with future engine development (with for example fuel cell powered electrical engines) where the RPM increase should be able to happen faster than with current marine diesels.
Acknowledgements
All RANS-CFD simulations conducted as a part of this paper uses the open source CFD package OpenFOAM. The research has been conducted within the framework of the research project entitledship design for enhanced sustainabilitywhich is sponsored by the Lloyd’s Register Educational Trust whose support is gratefully acknowledged.
This paper has been selected from the Proceedings of PRADS 2013, reviewed by referees and modified to meet guidelines for publication in IJNAOE.
References
Bernitsas M. , Ray D. , Kinley P. 1981 Kt, Kq and efficiency curves for the Wageningen B-series propellers, Technical report Department of Naval Architecture and Marine Engineering College, The University of Michigan Ann Arbor
Carrica P. , Ismail F. , Hyman M. , Bhushan S. , Stern F. 2013 Turn and zigzag maneuvers of a surface combatant using a URANS approach with dynamic overset grids Journal of Marine Science and Technology 18 166 - 181
Carrica P. , Paik K. , Hosseini H. , Stern F. 2008 URANS analysis of a broaching event in irregular quartering Seas Journal of Marine Science and Technology 13 395 - 407
Goldstein S. 1929 On the vortex theory of screw propellers Proceedings of the Royal Society of London 123 (792) 440 - 465
2011 MARPOL annex VI energy efficiency amendments, Technical report, Resolution MEPC IMO London 203 (62)
Jacobsen N.G. , Fuhrman D.R. , Fredse J. 2012 A wave generation toolbox for the open-source CFD library: OpenFoamⓒ International Journal Numerical Methods in Fluids 70 1073 - 1088
Journée J. 1992 Experiments and calculations on 4Wigley hull forms in head waves, Technical report Delft university of Technology Delft
Kashima T. , Takata J. 2002 An optimal control of marine propulsion system considering ship dynamics In Proceedings of the 2002 IEEE International Conference on Control Applications 1058 - 1063
Kjellberg M. 2011 Fully nonlinear unsteady three-dimensional boundary element method for force predictions on a restrained hull in waves. Thesis for the Degree of Licentiate of Engineering Chalmers University of Technology
Larsson L. , Stern F. , Visonneau M. 2014 Numerical ship hydrodynamics Springer New York An assessment of the Gothenburg 2010 Workshop 185 -
Larsson L. , Stern F. , Visonneau M. 2003 Ten years of industrial experience with the SST turbulence model Heat and Mass Transfer 4 625 - 632
Molland A. , Turnock S. , Hudson D. 2011 Ship resistance and propulsion: practical estimation of ship propulsive power Cambridge University Press New York
Phillips A. , Turnock S. , Furlong M. 2009 Evaluation of manoeuvring coefficients of a self-propelled ship using a blade element momentum propeller model coupled to a Reynolds averaged Navier Stokes flow solver Ocean Engineering 36 (15-16) 1217 - 1225
Phillips A. , Turnock S. , Furlong M. 2010 Accurate capture of rudder-propeller interaction using a coupled blade element momentum-RANS approach Ship Technology Research(Schiffstechnik) 57 (2) 128 - 139
Windén B. 2012 Powering performance of ships in waves. M.Phil/Ph.D Transfer thesis University of Southampton
Windén B. , Turnock S. , Hudson D. 2012 Validating force calculations using OpenFOAMⓒon a Fixed Wigley hull in waves In Proceedings of the 15th Numerical Towing Tank Symposium 170 - 175
Windén B. , Turnock S. , Hudson D. 2013 Predicting powering performance changes for ships in offshore conditions from small design modifications In Proceedings of the 23rd International Offshore and Polar Engineering Conference (ISOPE13) 773 - 780
Xiros N. 2002 Robust control of diesel ship propulsion Springer-Verlag London