ON RESISTANCE REDUCTION OF A HULL BY TRIM OPTIMIZATION

Transcription

1 Brodogradnja/Shipbilding/Open access Volume 69 Number 1, 2018 Xujian Lyu Haiwen Tu De Xie Jianglong Sun ISSN X eissn ON RESISTANCE REDUCTION OF A HULL BY TRIM OPTIMIZATION Abstract UDC : Original scientific paper The paper aims at conducting trim optimization for a hull to reveal the influence of trim on wave resistance by a potential-based panel method coupled with a response surface method. First, a numerical program for solving the linear free-surface flow problem of a hull moving with a uniform speed in calm water is built by the panel method. The S60 hull model is used to validate the numerical procedure. Next, calculation for hull is performed with two different trims at a wide range of Froude number; resistance test is conducted to validate the numerical method in demonstrating the influence of trim on wave resistance. Finally, a response surface of wave resistance is constructed with respect to variations of trim and Froude number, using the database of wave resistance calculated by the surface method. In this way, a framework is developed to perform trim optimization. The optimum trim point for the present hull shows a significant improvement in both wave resistance and total resistance, compared with that of even keel and the worst trim point. The optimization framework is proved to be effective in energy saving due to resistance reduction. Key words: Trim optimization; Wave resistance; Panel method; Response surface 1. Introduction In the face of global climate change, the importance of addressing greenhouse gas emissions attracts worldwide attention gradually. For shipping industry, the problems of greenhouse gas emission and fuel efficiency in operation are becoming more and more important due to the requirement of Environmental Ship Index (ESI) from International Maritime Organization (IMO) [1, 2]. When it comes to green shipping, one of the main concerns about ship hydrodynamics is resistance reduction. As it is well known, employing excellent hull forms or optimizing hull lines with low resistance in ship design stage is an effective measure for this goal [3-6]. However, it usually requires modifications to hull lines, which may cost lots of money and time for ships in operation. Another option to optimize the resistance performance without modifications to hull is increasing the operational efficiency of ships [7]. Traditionally, ships have been optimized for minimum resistance at the designing trim condition, which is the even keel trim condition or trimming somewhat aft for most 1

2 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization ships. Thus, one might think that the optimum trim condition, when wave-making resistance is considered, should not always be the designing one. With typical wave-making resistance taken into account, the paper is therefore to explore the optimum trim condition for a certain hull at given Froude numbers and displacement. As one of the most widely used tools for wave-making calculation in hydrodynamics, the potential-based panel method deals with the Laplace boundary-value problem with solving integral equation only at boundaries of the flow field [8-11]. Therefore, it is usually more effective in wave resistance calculation than other tools such as viscous flow based Computational Fluid Dynamics (CFD). Ghassemi and Kohansal studied the generation of wave due to moving hydrofoil in steady streams close to free surface, employing a potential based panel method [12]. A modified Rankine panel method was used to improve the prediction of wave pattern and wave-making resistance for ships of full form by Peng et al. [13], and an alternative method was first developed to satisfy the free surface boundary conditions on arbitrary body-fitted lines instead of streamlines. Lv et al. employed a potentialbased boundary element method to determine a better trim condition with lower wave resistance of the hull. However, without an optimization model, the trim had not been proved to be optimum for the hull. Iakovatos et al. [14] experimentally studied the impact of trim angles on resistance. In addition, Sherbaz and Duan assessed the influence of trim on ship resistance with CFD as well [15]. More recently, Sun et al. [16] developed a response surface based program with database from experiment and CFD to perform trim optimization. The optimized results had been validated by real ship test as well. However, experiment and CFD require a lot of money or computer resource input to assess the influence of trim on ship resistance. In this study, an efficient and economical framework is developed, to get resistance reduction based on trim optimization for a certain hull. Other than using a double model based boundary element method for wave resistance calculation [17], a direct potential-based panel method is applied to build a database for wave resistance. The Series 60 hull is utilized as a verification example and model test for a certain hull is conducted to verify the influence of trim on wave resistance as well. Finally, the response surface method, with surface fitting by Lagrange interpolation function instead of polynomial function, is applied to perform trim optimization, and the influence of trim on wave resistance and total resistance is evaluated. 2. Statement of the Boundary Value Problem Let us consider a ship moving with a constant speed U in calm water. A ship-fixed Cartesian coordinate system O-xyz is defined with origin O located at the mid-ship on the calm water surface and the x-axis positive towards the ship stern. The y-axis positive extends to starboards and the z-axis is vertically towards (Fig. 1). It is assumed that the fluid is incompressible and inviscid, and the flow is irrotational. The total velocity potential can be expressed as a sum of the oncoming velocity potential and the perturbation potential, namely U x (1) Fig. 1 Definition sketch of the coordinate system 2

3 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, Conservation of the mass applied to the potential flow gives the Laplace equation as a governing equation 2 0 Then motion of the flow can be uniquely defined by imposing the boundary condition on the boundary surface as follows: 1) Hull boundary condition U n on SH (3) n where n is the unit vector normal to the boundary, defined positive when pointing into the fluid region. 2) Linearized free surface boundary condition can be simplified as 2 2 U z g x 2 0 on SF (4) where g is gravity acceleration. In addition, a four-point finite upstream scheme is used to ensure that the free-surface waves vanish upstream of the disturbance. By applying Green s theorem, the equation about the perturbation potential can be expressed as where where S H S F G - G d S 0 n n For the problem of flow around the hull. Equation (4) can be written as G ( q) G 4 E ( p) ( q) ds d ( ) d S G S q S H n SH n SF n 1 on S E 2 1 on S H F Green s function G is given by 1 1 G R ( p, q ) R '( p, q ) q q q R ( x ) ( y ) ( z ) is the distance between field point and the source point, and R' ( x ) ( y ) ( z ) is the distance between field point and the image of source point. The hull surface SH and the free surface SF are discretized into a number of NH and NF panels respectively to get an approximate solution of the problem [18]. The perturbation potential and its normal derivative are assumed constant at each panel and equation n about and can be obtained as n 3 (2) (5) (6) (7) (8)

4 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization NH ( ) NF q ( q) ( p) ( q) I I I HD HS FS j 1 nq j 1 nq where, and are influence coefficients proposed by Morino et al. [19] and the details can be referred to the literature [20]. As the free surface should be satisfied on the true surface, an iterative method is used to obtain the solution of equation (9). To get acceptable results, the iteration number is about 4~7 in this paper. The wave resistance coefficient I HD I HS I FS R 1 C C i n S N /2 w () 2 H w p xi W i (1/ 2) U SW SW i 1 S W where is the wetted surface area of the hull in calm water, is the water density, Wi is the area of the panel element, is x-component of the unit vector normal to the surface panel element, and C p n xi is the pressure coefficient determined by the Bernoulli s equation. The wave height can be calculated with the linearized kinematic free surface condition U ( x, y) on z 0 g x S (9) (10) (11) 3. Sample calculation and validation The Series 60 hull model (Cb = 0.6) is used as an example to check the reliability for the numerical procedure of present panel method [21]. The free surface is discretized (-1 x / Lpp 2, -1 y / Lpp 1) by quadrilateral panels and the hull by 1346 panels, 5 6 for the stern and for other parts of the hull respectively. Calculation is conducted for the S60 hull with fixed trim and sinkage and computed wave resistance coefficient is compared with that obtained from the literature by Tarafder & Suzuki [18], as shown in Fig. 2. Though there are some differences between the results of wave resistance from present study and that from literature at some Froude numbers, the overall agreement of the results is quite acceptable. In addition, numerical results on wave profiles along the body surface at Fn = 0.22 and Fn = 0.34 are shown in Figs 3 and 4 respectively. The variation trend of wave profiles from present study and literature matches well with each other. The present results show higher accuracy than that from literature when compared with experimental results, as shown in the figures. The comparisons show that the present numerical method is able to catch typical flow characters around the hull, and thus it can be concluded that the numerical procedure is feasible for wave resistance and wave profile prediction. Subsequently, the numerical method is applied for the calculation of a HUST (Huazhong University of Science & Technology) hull [17], of which the main parameters are summarized in table 1. Both the hull surface and the free surface domain (-1 x / Lpp 2, -1 y / Lpp 1) are discretized by quadrilateral panels, as shown in Figs 5 and 6. The displacement of the hull is kept constant during the calculation. To provide further validation of the numerical procedure, resistance tests of the HUST hull were performed at the towing tank laboratory of the School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan, P. R. China. The principle dimension of the towing tank is m [16]. According to the ITTC guidelines, the wave resistance coefficient of HUST hull from model test can be written as C (1 k) C +C (12) t f w 4

5 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, Fig. 2 Wave resistance of S60 hull Fig. 3 Wave profiles of the S60 hull at Fn = 0.22 Fig. 4 Wave profiles of the S60 hull at Fn = where Ct is the total resistance coefficient got from model test, Cf is the frictional resistance coefficient determined by the ITTC 57 correlation line and k is the form factor determined by Prohaska method adopted by ITTC in The trim condition of present hull is defined by the aft draft Ta and the fore draft Tf, that is Trim T T (13) f a where Trim > 0 indicates trim by bow and Trim < 0 indicates trim by stern, as shown in Fig. 7. 5

6 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization Table 1 Typical panel arrangement for the hull Descriptions Parameters Value Volume of displacement (m3) Wetted surface area Sw (m2) Length between perpendiculars LPP (m) Beam B (m) Draft D (m) Block coefficient CB Fig. 5 Typical panel arrangement for the HUST hull Fig. 6 Typical panel arrangement for the free surface Fig. 7 Sketch of trim conditions for the hull 6

7 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, Comparisons of wave resistance for the hull at two different trim conditions between numerical results and experiment results are shown in Fig. 8. The variation trend of calculated and experimental wave resistance with respect to trim conditions agrees quite well at a wide range of Froude number (i.e. Fn > 0.17), even at the inflection point (i.e. Fn 0.192, Fn 0.207) where the trim variations leads to few changes in wave resistance. During the experiment, the HUST hull is tested with free sinkage and it s found that sinkage of the hull is more significant at Fn=0.203 than that at other Froude numbers. Since the wave resistance calculation was conducted with fixed sinkage for the HUST hull, one possible reason for the peak difference between wave resistance curves from calculation and experiment (Fn = in Fig. 8) is that, the fixed sinkage during calculation covers up the influence of significant free sinkage to wave resistance. Typical wave profiles for the hull at Fn = and are shown in Figs 9 and 10 respectively. The variation trend of wave heights over trims from calculation shows a good agreement with that from experiment as well. In addition, calculated wave patterns provide further verification by detailing wave height changes with different trim conditions around the hull, as depicted in Figs 11 and 12. Thus, it is can be concluded that the present numerical procedure is able to properly demonstrate the wave resistance changes caused by trim adjustment for the HUST hull at Fn > Trim optimization To explore the relationship between a few explanatory variables (i.e. trim condition, Froude number) and a response variable (i.e. wave resistance), the response surface methodology is introduced to conduct trim optimization [16, 22]. Instead of using a seconddegree polynomial model to approximate the surface, the Lagrange interpolation function is applied to make sure that the effective points from calculation is included in the response sutrface. Fig. 8 Wave resistance of the HUST hull 7

8 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization Fig. 9 Wave profiles along the HUST hull at Fn=0.203 Fig. 10 Wave profiles along the HUST hull at Fn=0.217 Fig. 11 Wave patterns around the HUST hull at Fn=0.203 Fig. 12 Wave patterns around the HUST hull at Fn=

9 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, To create a database, the wave resistance of the HUST hull at different Froude numbers and trims is calculated with above numerical method and summarized in Table 2. The response surface of wave resistance with respect to Froude numbers and trims is constructed based on the wave resistance database, as depicted in Fig. 13. The trim optimization of the HUST hull is conducted by using the response surface, in which a response of wave resistance is influenced by trim and Froude number. In present trim optimization for the HUST hull with the displacement at a certain initial draft (i.e. D = mm), the design variables are trim and ship speed based Froude number. The objective is to find the optimal value of trim to minimize the wave resistance of the hull at a given Froude number. The optimization model is given as follows: C f ( Trim, Fn, Dis ) (14) w min f ( Dis, Fn, Trim) Dis Dis0 st. Fn Fn0 Trim1 Trim Trim2 where f(trim, Fn, Diss) is the response surface, Dis0 is present volume displacement of the hull, Fn0 is the Froude number based on given speed, and Trim1 and Trim2 define the variation range of the trim condition. The overall optimization procedure is demonstrated in Fig. 14. There exists an optimal trim curve (may be made up by several curve portions) on the response surface, and the corresponding wave resistance is the minimum. The intersection point of the optimal trim curve and the Froude number isoline represents the optimal point of the minimum wave resistance at present ship velocity, as shown in Fig. 13 and detailed in Figs To evaluate the trim optimization results, the wave resistance reduction ratio is defined as w ( Cw Optimum Cw) 100% C w and the total resistance reduction ratio is defined as t ( Ct Optimum Ct) 100% C t where CwOptimum and CtOptimum are wave resistance coefficient and total resistance coefficient at the optimum trim condition, respectively. Cw and Ct are wave resistance coefficient and total resistance coefficient at the reference trim condition point, respectively. The optimized trim condition is able to provide 9.7%-11.5% of wave resistance reduction and 3.5%-3.8% of total resistance reduction, compared with the even keel condition for the HUST hull. What s more, Table 2 Numerical results of wave resistance for the HUST hull Tirm (mm) Fn C w (15) (16) (17)

10 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization Fig. 13 Response surface of wave resistance Fig. 14 Optimization procedure based on response surface 10

11 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, Fig. 15 Froude number isoline of response surface at Fn = Fig. 16 Froude number isoline of response surface at Fn = Fig. 17 Froude number isoline of response surface at Fn = Table 3 Resistance reduction ratio at optimum trim condition point Fn C woptimum α w related to Even-keel point 11.5% 26.2% 9.7% α w related to Worst point 20.2% 35.8% 23.7% α t related to Even-keel point 3.5% 7.2% 3.8% α t related to Worst point 6.8% 11.3% 11.8% the reduction ratio can be as high as 20.2%-35.8% of wave resistance and 6.8%-11.3 of total resistance while the comparison is taken between the optimized trim condition and the worst trim condition, as shown in Table 3. Obviously, the contribution of trim optimization to resistance reduction for the HUST hull is significant. 11

12 Xujian Lyu, Haiwen Tu, On resistance reduction of a hull by trim optimization 5. Conclusion Traditionally, hull lines are designed or optimized for minimum resistance at the designing trim condition, which is the even keel trim condition or trimming somewhat aft for most ships. To achieve reduction of resistance in a different way, the paper in which typical wave making of a hull is taken into account therefore aims at exploring a trim for the HUST hull, with minimum wave resistance at given Froude numbers and displacement. After being validated by results from calculation of the S60 hull and resistance test of the HUST hull, a potential-based panel method coupled with a response surface method is used to develop a framework for trim optimization. The optimum trim condition for the hull shows a significant improvement in both wave resistance and total resistance. The optimization framework is effective in energy saving due to resistance reduction, not only for the present HUST hull, but also for general hull models. Acknowledgement: The authors gratefully acknowledge the financial support from the National Natural Science Foundation of China (NSFC) under Project Number & , and the financial support from the Opening Fund of the Hubei Key Laboratory of Naval Architecture & Ocean Engineering Hydrodynamics under Project Number as well. REFERENCES [1] Traut, M., Gilbert, P., Walsh, C., Bows, A., Filippone, A., Stansby, P., Wood, R.: Propulsive power contribution of a kite and a Flettner rotor on selected shipping routes, Applied Energy, 113: , [2] Gibbs, D., Rigot-Muller, P., Mangan, J., Lalwani, C.: The role of sea ports in end-to-end maritime transport chain emissions, Energy Policy, 64: , [3] Huang, F., Wang, L., Yang, C.: A new improved artificial bee colony algorithm for ship hull form optimization, Engineering Optimization, DOI: / X , [4] Ayob, A. F., Nik, W. W., Ray, T., Smith, W. F.: Hull surface information retrieval and optimization of high speed planing craft, IOP Conference Series: Materials Science and Engineering, [5] Tahara, Y., Peri, D., Campana, E. F., Stern, F.: Single- and multiobjective design optimization of a fast multihull ship: numerical and experimental results, Journal of Marine Science and Technology, 16(4): , [6] Sun, J., Lv, X., Liu, W., Ning, H., Chen, X.: Research on a method of hull form design based on wavemaking resistance optimization, Polish Maritime Research, 19(3): 16-25, [7] Banks, C., Turan, O., Incecik, A., Theotokatos, G., Izkan, S., Shewell, C., Tian, X.: Understanding ship operating profiles with an aim to improve energy efficient ship operations, Proceedings of the Low Carbon Shipping Conference, London, [8] Lee, C., Newman, J. N.: Computation of wave effects using the panel method, Numerical Models in Fluid Structure Interaction, 42: , [9] Yan, H., Liu, Y.: An efficient high-order boundary element method for nonlinear wave wave and wave-body interactions, Journal of Computational Physics, 230(2): , [10] Belibassakis, K. A., Gerostathis, T. P., Kostas, K. V., Politis, C. G., Kaklis, P. D., Ginnis, A. I., Feurer, C.: A BEM-isogeometric method for the ship wave-resistance problem, Ocean Engineering, 60: 53-67, [11] Liu, S. K., Papanikolaou, A., Zaraphonitis, G.: Prediction of added resistance of ships in waves, Ocean Engineering, 38(4): , [12] Ghassemi, H., Kohansal, A. R.: Wave generated by the NACA 4412 hydrofoil near free surface, Journal of Applied Fluid Mechanics, 6(1): 1-6,

13 On resistance reduction of a hull by trim optimization Xujian Lyu, Haiwen Tu, [13] Peng, H., Ni, S., Qiu, W.: Wave pattern and resistance prediction for ships of full form, Ocean Engineering, 87: , [14] Iakovatos, M. N., L Iarokapis, D. E., Andtzabiras, G. D.: Experimental investigation of the trim influence on the resistance characteristics of five ship models, Development in Marine Transportation and Exploitation of Sea Resources, 1: 23-32, [15] Sherbaz, S., Duan, W.: Ship Trim Optimization: Assessment of Influence of Trim on Resistance of MOERI Container Ship, The Scientific World Journal, 2014, [16] Sun, J., Tu, H., Chen, Y., Xie, D., Zhou, J.: A study on trim optimization for a container ship based on effects due to resistance, Journal of Ship Research, 60(1): 30-47, [17] Lv, X., Wu, X., Sun, J., Tu, H.: Trim optimization of ship by a potential-based panel method, Advances in Mechanical Engineering, 5: , [18] Tarafder, M. S., Suzuki, K.: Wave-making resistance of a catamaran hull in shallow water using a potentialbased panel method, Journal of Ship Research, 52(1): 16-29, [19] Morino, L., Chen, L. T., Suciu, E. O.: Steady and oscillatory subsonic and supersonic aerodynamics around complex configurations, AIAA Journal, 13(3): , [20] Tarafder, M. S., Suzuki, K.: Computation of wave-making resistance of a catamaran in deep water using a potential-based panel method, Ocean Engineering, 34(13): , [21] Toda, Y., Stern, F., Longo, J.: Mean-flow measurements in the boundary layer and wake and wave field of a Series 60 Cb=0.6 ship model - part 1: Froude numbers 0.16 and 0.316, Journal of ship Research, 36(4): , Submitted: Accepted: Xujian Lyu 1, 3 Haiwen Tu 2 De Xie 2, 3, 4, Jianglong Sun 2, 3, 4, * 1 School of Energy and Power Engineering, Nanjing University of Science and Technology, Nanjing , P. R. China 2 School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Hubei, P. R. China 3 Hubei Key Laboratory of Naval Architecture and Ocean Engineering Hydrodynamics, Huazhong University of Science and Technology, Hubei, P. R. China 4 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai, P. R. China * Corresponding author (J. Sun): sundapao@126.com ( ); (Fax); (Tel) 13