Interaction Effect between Hull and Accommodation on Wind Drag Acting on a Container Ship

In this paper, we present our research on applying the commercial Computational Fluid Dynamics (CFD) code to investigate interaction effect between hull and accommodation on wind drag acting above the water hull surface of a full scale 1200 TEU container ship. With this purpose, aerodynamic performances and wind drag acting on the ship hull with and without accommodations have been computed. Analyzing the obtained CFD results, the interaction effect between hull and accommodation on aerodynamic performances and wind drag acting on the ship have been found. Various new accommodation shapes have been proposed for the original ship to reduce the interaction effect on wind drag. A drastic reduction in the interaction effect between hull and accommodation on wind drag acting on the ship has been achieved and the obtained results have been shown in this paper.


Introduction
Nowadays, research on reducing wind resistance acting on ship hull to save fuel consumption and improve economic efficiency is still an attractive topic in the field of marine transportation. This is especially true for container ships with a large windward area, due to their high above-water superstructure and multi-external forms, resulting from the position of containers on the deck, where wind drag accounts for a large percentage of the total resistance. Therefore, researchers around the world are of the opinion that research on reducing wind drag acting on container ships as well as the ships with high superstructure and or large windward areas is more and more important.
It goes without saying that in recent years, a large number of studies on the reduction in wind drag acting on ships' hulls as well as on optimal aerodynamic hull forms, the effect of side covers and domes on wind drag, the interaction effect between hull and accommodation, the effect of strong wind on resistance, the safety of the ship in strong wind, and so on, have been published. A comprehensive review is given as follows: There are many studies on applying commercial Computational Fluid Dynamics (CFD) to solve the aerodynamic performances of the ships. The most important point of the research is that: the CFD has been a popular and useful tool to solve ships' aerodynamic problems with fairly good accuracy. Using CFD modelling, various types of modified hull shapes with reduced wind drag acting on the ship have been developed [1][2][3][4][5][6][7].
Janssen, W.D. et al. (2017) presented a study on using both the commercial CFD code and wind tunnel model test to study aerodynamic performances and wind drag acting on a container ship hull. The 3D steady Reynolds Averaged Navier Stokes (RANS) CFD simulation method was used to solve the problem. The author concluded that CFD results were in good agreement with the experimental ones and that the CFD could be used to develop a new hull shape with reduced wind drag acting on the ship. The average absolute difference of wind drag obtained in CFD simulation and tunnel test ranged from 37.9% for a box shaped representation of the ship to only 5.9% for a more detailed model. Modelling the spaces in between container stacks decreased the average total wind load considerably. The average absolute difference of total wind load was 10.4%. Using a slender ship hull instead of a blunt ship one decreased total wind load by up to 5.9%. Taking into account wind tunnel blockage following the approach of the engineering science data showed an underestimation of up to 17.5% for the lateral wind load, as shown by comparing the CFD results obtained in the narrow domain with those of the wider domain [1].
Andersen, I.M.V. (2013) presented a study that investigates the influence of the container configuration on the deck of a 900TEU container ship on wind drag by wind tunnel test with scale model of 1:450. The author concluded that the results serve as an indication that the magnitude of wind force acting on a large container ship depends on the container configuration on deck rather than the purpose of assessing the full scale wind resistance of a given container ship. For the reduced longitudinal force in relative wind, it was advantageous to make the container configuration as smooth as possible and streamlining could reduce the force in the head wind. However, streamlining of the configuration on the aft deck was a trade-off as it increased the yaw moment compared to full load on the aft deck. The high container stacks in the configuration appeared to increase the longitudinal force more than the low or empty container bays in the configuration. The transverse force of concern in beam wind depends largely on the side area of the ship, and for a fully loaded ship it could be reduced by reducing the stacking height of the outermost stacks. The yaw moment in relative wind could be reduced by achieving full load on the aft of the ship, and it is possible to reduce the resistance through a combination of drift and increased rudder angle. Moreover, the author gave a general recommendation that the external form of the ship should be made as smooth as possible. In addition, it is important that the center of gravity of the side area was as far from the aft as possible [2]. Fujiwara, T. et al. (2009) presented an experimental investigation on wind force for a container ship with various external forms due to different positions of containers on the desk, in order to study aerodynamic specifications. Based on the obtained experimental results, an estimation method for wind load acting on the container ship has been proposed. In the tunnel experiment, a scale model of 1.5 m long, a mean wind velocity of 25 m/s and a Reynolds number of about 2.4 × 10 6 have been selected. Under these conditions, the flow field was turbulent and the drag coefficients were independent of the Reynolds number. The authors have shown that a good agreement between their proposed estimation method the tunnel experimental test has been obtained. They also showed that the new method is important for the calculation of speed as well as other characteristics in the operation stage of container ships [3].
In a study presented by Kim, Y. et al. (2015), several design concepts and devices on the superstructure of a container ship have been suggested and tested in a wind tunnel to estimate the wind drag reduction. The authors have also used CFD with the RANS simulation method to estimate wind drag acting on the ship. The results show that the gap protectors between container stacks and visors in front of upper deck have been found to be the most effective means for reducing wind drag acting on the ship. The CFD results agreed with the experimental measurements in the wind tunnel, and the wind drag acting on the modified ships could be reduced by up to 56% in the wind direction angle from 0 to 50 degrees [4]. Other researchers [5][6][7] also presented results on using CFD and experimental tests to develop a modified hull shape with reduced wind drag acting on a container ship. The authors have proposed modified hull shapes with attached side covers, a center wall, a "T" center wall and a dome at the bow deck of the container ship to reduce wind drag. By using side covers and a center wall, the container ship could reduce wind drag by up to 40% in the head wind direction. A dome at the bow of the ship could reduce the total wind drag acting on the container ship by up to 30% at wind direction angles of less than 30 degrees.
Other available studies focused on the aerodynamic performances and wind drag acting on different types of ships and offshore using both the wind tunnel measurement and CFD simulation. Jun, S. et al. (2020) investigated the resistance performance of ships, using the air resistance correction method [8]. Wnęk, A.D. et al. (2011Wnęk, A.D. et al. ( , 2015 focused on the wind load acting on a LNG carrier with a specific geometrical hull shape [9,10]. Saydam, A.Z. et al. (2018) conducted an evaluation of wind load acting on ships by CFD analysis [11]. He, N.V. et al. (2013He, N.V. et al. ( , 2016He, N.V. et al. ( , 2019 presented research on the interaction effect between hull and accommodation on aerodynamic performances and wind drag acting on a wood chip carrier hull [12][13][14]. Sugata, K. et al. (2010) studied the reduction in wind drag acting on the hull of a Non Ballast water Ship (NBS) [15], and so on [16][17][18][19].
The mesh of computed domain effects on the CFD results of a ship has been reported in many previous published works. Viola, I.M. et al. (2009) tested two different mesh numbers of 1 million and 6.5 million tetrahedral elements. The obtained CFD results, i.e., the wind forces, have been compared with those of the experimental ones. It has been observed that the wind forces computed by the turbulent viscous model Realizable k-ε was in good agreement with the experimental results and the differences between computed results in the two mesh numbers were lower than 5% [20]. Wnęk, A.D. et al. (2011Wnęk, A.D. et al. ( , 2015 used two models of a floating LNG platform and an LNG carrier in the 1:400 scale to conduct an experimental tunnel test and CFD computation to investigate the wind forces acting on the models. The turbulent viscous model k-ω SST and three different mesh generation techniques, which were CFD hexa (y+ = 5), CFD tetra (no prism) and CFD tetra y+ (y+ = 0.1), were used for CFD computation. The CFD results of CFD tetra (y+ = 0.1) showed the best agreement with the experimental data obtained in the wind tunnel test and the tetrahedral mesh was shown to be of good quality for the meshing of complex geometric shapes [9,10]. Saydam, A.Z. et al. (2017) used the mesh numbers from 1 million to 5 million elements of unstructured tetrahedral mesh to test CFD results independent of mesh number for computed wind forces acting on the ships. It has been shown that the CFD results were in good agreement with the experimental data obtained in wind tunnel test, and mesh independence of wind forces was up to more than 2 million elements [11]. Watanabe, I. et al. (2016) and Trieu, N.V. et al. (2017a, 2017b) used unstructured tetrahedral mesh with three mesh numbers of 2.2 million, 2.6 million and 3.8 million elements, with y+ of less than 25, to test CFD results independent of mesh. They showed that the CFD results given were in good agreement with the experimental data, and about 2% of the wind drag coefficient was different from the mesh number [5,7,21]. In our previous paper [12][13][14], we used the same unstructured tetrahedral mesh to investigate the wind drag acting on a wood chip carrier. It has been shown that the CFD results were in good agreement with those of the experimental data obtained in the towing tank test conducted at Osaka Prefecture University, Japan.
In this paper, a commercial CFD code ANSYS-Fluent has been used to investigate the interaction effect between the hull and accommodation of a container ship on wind drag. By applying the CFD, aerodynamic performances and wind drag acting on a 1200 TEU container ship have been computed for two different cases, namely the hull with and without accommodation on its deck, to determine the interaction effect between the hull and accommodation of the ship. From the obtained results, several new hull shapes and frontal accommodation shapes, with reduced interaction effects between hull and accommodation on wind drag acting on the ship, have been proposed in this paper.

Model Ship Used for Computation
The full scale 1200 TEU container ship, which is designed with an accommodation located at the aft, has been used as the reference model for computation. The principal particulars of the ship are given in Table 1, and Figure 1 shows the designed hull form of the original model.

Computational Domain and Boundary Conditions
In this research, all steps for the designed simulation domain, meshing and setting of the boundary conditions have been done as per the user guide published by the International Towing Tank Conference (ITTC) manual for using CFD, [22][23][24][25][26] and our previous publications [12][13][14]. In this

Computational Domain and Boundary Conditions
In this research, all steps for the designed simulation domain, meshing and setting of the boundary conditions have been done as per the user guide published by the International Towing Tank Conference (ITTC) manual for using CFD, [22][23][24][25][26] and our previous publications [12][13][14]. In this study, the computational domain has been designed with 6.5 L, 2.2 L and 0.75 L instead of 1200 m length, 450 m breadth and 150 m height. Mesh of the computed domain has been used as unstructured mesh with a triple prism layer and y+ has been taken to be less than 5. The popular turbulent viscous model k-ε for unsteady flow has been used. The velocity inlet is given at the inlet as the wind velocity, the pressure outlet is given at the outlet of the computed domain. The bottom, top and sides are given at the walls. To meet realistic case pressure, the outlet makes the open air condition possible, the limited dimension of the computed domain must be designed to be large enough. However, the dimensions of the computed domain has an effect on the mesh number and also an effect on running time. If the top and side of the domain is far enough, the effects on the CFD result become small. Therefore, we must design an appropriate domain and boundary condition. In computation, the time size step is of 0.005 s and 6000 time steps. All cases have been done by an Intel Xeon Gold 6138 2.0 GHz computer, with 3.7 GHz Turbo, 20 C with 64 GB of RAM. Figure 2 shows computed domain and mesh used for simulation. result become small. Therefore, we must design an appropriate domain and boundary condition. In computation, the time size step is of 0.005 s and 6000 time steps. All cases have been done by an Intel Xeon Gold 6138 2.0 GHz computer, with 3.7 GHz Turbo, 20 C with 64 GB of RAM. Figure 2 shows computed domain and mesh used for simulation.

Covegence of Mesh Using for Computation
In this research, the effect of mesh numbers generated in the computed domain on aerodynamic performance and wind drag acting on the original container ship has been investigated. The computed domain is shown in Figure 2, where twelve different mesh numbers have been used. The mesh numbers are from 0.05 to 21.52 million elements and value of y+ is taken from 1.08 to 502.49, respectively. The meshes have been used to compute the aerodynamic performances of the ship under the same conditions. Analyzing the obtained CFD results, the effect of mesh numbers on aerodynamic performances and wind drag as well as the state of independence of mesh numbers have been found. The computation was done in head wind condition at a wind velocity of 14 knot instead of a Reynolds number of 6.73 × 10 7 . Table 2 shows a detailed mesh of the computed domain.

Covegence of Mesh Using for Computation
In this research, the effect of mesh numbers generated in the computed domain on aerodynamic performance and wind drag acting on the original container ship has been investigated. The computed domain is shown in Figure 2, where twelve different mesh numbers have been used. The mesh numbers are from 0.05 to 21.52 million elements and value of y+ is taken from 1.08 to 502.49, respectively. The meshes have been used to compute the aerodynamic performances of the ship under the same conditions. Analyzing the obtained CFD results, the effect of mesh numbers on aerodynamic performances and wind drag as well as the state of independence of mesh numbers have been found. The computation was done in head wind condition at a wind velocity of 14 knot instead of a Reynolds number of 6.73 × 10 7 . Table 2 shows a detailed mesh of the computed domain.  Tables 3 and 4 show detailed results of wind drag, and the interaction effect of the mesh number on wind drag acting on the ship. Figure 3 shows the curve of the effect of the mesh number on wind drag, as well as the state of the mesh number becoming independent of the wind drag acting on the ship obtained in the computation. around 5% when the y+ is less than 50 and the same value when y+ is less than 7. When the y+ increases to over 50 but continues to be less than 500, the difference in wind drag increases by up to 42%, as shown. This is important in studying aerodynamic performances using the CFD tool. Apart from paying attention to calculating the wind drag acting on the ship, it is very useful to get a clear pattern of pressure distribution around the hull to understand how the flow has an effect on the hull, as well as factors affecting the wind drag, and so on.  Figure 3 shows the effect of mesh number on wind drag acting on the ship as well as the mesh convergence curve obtained in computation. From the results shown, it can be seen that when the mesh number increases to over two million instead of a y+ of less than 50, the effect of the mesh number on total wind drag acting on the ship drastically reduces and becomes zero when the y+ is less than 7, as shown. The obtained result is useful in the application of CFD to investigate aerodynamic performances of the ship.

Interaction Effect between Hull and Accommodation
In this section, aerodynamic performances of the ship have been computed in two different cases, namely the hull with and without accommodation on its deck. Analyzing the obtained CFD Tables 3 and 4 clearly show wind drag acting on the original ship in the different computed mesh numbers, where the wind drag acting on the ship is defined by the following equation [27,28]. where: C d is the wind drag coefficient. R x is the wind drag acting on the hull, N. A F(x) is the frontal projected area, m 2 . V is the ship velocity, m/s.
The results given in the Tables 3 and 4 show that the wind drag acting on the ship decreases slowly with increasing mesh number, although the wind drag is slightly different when the y+ is less than 50 and retains the same value when the y+ less than 7. The effect of mesh number on wind drag acting on the ship decreases with the increasing of the mesh number. The difference of wind drag is around 5% when the y+ is less than 50 and the same value when y+ is less than 7. When the y+ increases to over 50 but continues to be less than 500, the difference in wind drag increases by up to 42%, as shown. This is important in studying aerodynamic performances using the CFD tool. Apart from paying attention to calculating the wind drag acting on the ship, it is very useful to get a clear pattern of pressure distribution around the hull to understand how the flow has an effect on the hull, as well as factors affecting the wind drag, and so on. Figure 3 shows the effect of mesh number on wind drag acting on the ship as well as the mesh convergence curve obtained in computation. From the results shown, it can be seen that when the mesh number increases to over two million instead of a y+ of less than 50, the effect of the mesh number on total wind drag acting on the ship drastically reduces and becomes zero when the y+ is less than 7, as shown. The obtained result is useful in the application of CFD to investigate aerodynamic performances of the ship.

Interaction Effect between Hull and Accommodation
In this section, aerodynamic performances of the ship have been computed in two different cases, namely the hull with and without accommodation on its deck. Analyzing the obtained CFD results for the two cases, the interaction effect on aerodynamic performances of the hull and accommodation in head wind has been obtained. Figure 4 shows the pressure distribution around the ship in the two computed cases at a Reynolds number of 6.73 × 10 7 in head wind. results for the two cases, the interaction effect on aerodynamic performances of the hull and accommodation in head wind has been obtained. Figure 4 shows the pressure distribution around the ship in the two computed cases at a Reynolds number of 6.73 × 10 7 in head wind. Analyzing the results as shown in Figure 4, the red and yellow regions show high pressure and the blue region shows lower pressure acting on the model. Hence, the effect of hull and accommodation on pressure distribution over the hull surface of the ship could be evaluated. The interaction effect between hull and accommodation may also be determined by comparison of the wind drag acting on the hull and accommodation in both of the computed cases, as shown.
In this research, interaction effect between hull and accommodation on wind drag acting on the ship has been determined by the following Equation (2) where: Cd, % is the interaction effect between hull and accommodation. Cd (Hull with Acc) is wind drag coefficient acting on the hull with accommodation. Cd (Independent Hull and Acc) is total wind drag coefficient acting on hull without accommodation and independent accommodation. Tables 5-7 show, in detail, the wind drag acting on the hull, the accommodation, and the interaction effect between the hull and accommodation on the wind drag of the ship in head wind.  Analyzing the results as shown in Figure 4, the red and yellow regions show high pressure and the blue region shows lower pressure acting on the model. Hence, the effect of hull and accommodation on pressure distribution over the hull surface of the ship could be evaluated. The interaction effect between hull and accommodation may also be determined by comparison of the wind drag acting on the hull and accommodation in both of the computed cases, as shown.
In this research, interaction effect between hull and accommodation on wind drag acting on the ship has been determined by the following Equation (2) [12][13][14]: ∆C d , % = C d (Hull with Acc) − C d (Independent Hull and Acc) C d (Independent Hull and Acc) 100% (2) where: ∆C d , % is the interaction effect between hull and accommodation.
C d (Hull with Acc) is wind drag coefficient acting on the hull with accommodation. C d (Independent Hull and Acc) is total wind drag coefficient acting on hull without accommodation and independent accommodation. Tables 5-7 show, in detail, the wind drag acting on the hull, the accommodation, and the interaction effect between the hull and accommodation on the wind drag of the ship in head wind.
Analyzing the results shown in the Tables above, the wind drag acting on hull part and accommodation part in the case of the hull with accommodation is less than that of the case of accommodation independent of hull, by up to 42%. The total wind drag acting on the ship in the case of hull with accommodation is also less than that of the independent hull by up to 10%. For the wind drag acting on the hull part and the accommodation part, the interaction effect is about 42% for the hull part and 4% for the accommodation part. The interaction effect between the hull and accommodation on wind drag is around 10%, as shown.

Proposed Accommodation Shape for the Ship
The original container ship has been designed with an accommodation located at the aft of the ship. In this section, three new hull forms with an accommodation located at the frontal part of ship have been proposed. Three models with frontal accommodation, named N1, N2 and N3, have been used for computation to determine the interaction effect on aerodynamic performances of the proposed streamlined hull shapes. The dimensions of the models are the same ones as those of the original ship shown in Table 1. The computation for all ships have also been done under the same conditions as shown in the Section 2. Figure 5 shows the proposed new models N1, N2 and N3. In the proposed models, the N1 has the same accommodation shape of the original ship but is located at frontal hull, the N2 has a streamlined bow cover, and the N3 has a streamlined accommodation shape and it is located at frontal hull. Computation has been done for all the models under the same conditions to investigate the aerodynamic performances of the ships. In the proposed models, the N1 has the same accommodation shape of the original ship but is located at frontal hull, the N2 has a streamlined bow cover, and the N3 has a streamlined accommodation shape and it is located at frontal hull. Computation has been done for all the models under the same conditions to investigate the aerodynamic performances of the ships.

Reduced Interaction Effect of the New Hull Shapes on Pressure Distribution
In this section, the effect of the streamlined hull shapes on aerodynamic performances of the proposed container ships have been investigated by the CFD. Figures 6-11 show pressure distribution around and over the hull surface of the ships. Analyzing the obtained CFD results of pressure distribution around and over hull surfaces of the ships, the reasons for the increasing interaction effect between hull and accommodation have been clearly found.  In the above results, the blue region shows low dynamic pressure and, consequently, high static pressure acting on the ship. The results, as shown in Figure 6, clearly show reduced separation flow at the lower dynamic pressure region (blue color region) around the streamlined hull shape N3. Figure 7 shows dynamic pressure distribution around the ships at several horizontal planes of the computed domain.
The results, as shown in Figure 7, clearly show the effects of the accommodation shapes on the dynamic pressure distribution around the ships. The separation flow regions at the low dynamic pressure regions (blue color regions) at the aft of the accommodation N2 and N3 have been drastically reduced, as shown. Figure 8 shows the dynamic pressure distribution around several cross sections of the computed domain of the ships.
The results obtained in the above figures clearly show differences of dynamic pressure distribution around the ships. A drastic reduction in low dynamic pressure region (blue color) could be seen around the hull shapes of N2 and N3, as shown. These results clearly show the reduced interaction effect on ship aerodynamic performance of streamlined accommodation shapes which are located at frontal hull of the ship. There is a reason why a suitable frontal accommodation should be installed to reduce wind drag acting on the ship hull.  The results, as shown in Figure 7, clearly show the effects of the accommodation shapes on the dynamic pressure distribution around the ships. The separation flow regions at the low dynamic pressure regions (blue color regions) at the aft of the accommodation N2 and N3 have been drastically reduced, as shown. Figure 8 shows the dynamic pressure distribution around several cross sections of the computed domain of the ships. The results obtained in the above figures clearly show differences of dynamic pressure distribution around the ships. A drastic reduction in low dynamic pressure region (blue color) could be seen around the hull shapes of N2 and N3, as shown. These results clearly show the reduced interaction effect on ship aerodynamic performance of streamlined accommodation shapes which are located at frontal hull of the ship. There is a reason why a suitable frontal accommodation should be installed to reduce wind drag acting on the ship hull.

Interaction Effect between Hull and Accommodation on Wind Drag
In this section, wind drag acting on hull of the ship with and without accommodation has been investigated. By comparison, the obtained results of wind drag acting on hull of the ship in two different computed cases, namely hull with and without accommodation, interaction effect between hull and accommodation on wind drag acting on hull, have been found. Moreover, by comparison with the obtained results of wind drag acting on the ship with different accommodation shapes on the deck, the effect of accommodation shapes on wind drag has also been found. Figure 9 shows wind drag acting on the ships with different accommodation shapes on their decks in head wind. A detailed comparison of wind drag acting on the ships is shown in Table 8.

Interaction Effect between Hull and Accommodation on Wind Drag
In this section, wind drag acting on hull of the ship with and without accommodation has been investigated. By comparison, the obtained results of wind drag acting on hull of the ship in two different computed cases, namely hull with and without accommodation, interaction effect between hull and accommodation on wind drag acting on hull, have been found. Moreover, by comparison with the obtained results of wind drag acting on the ship with different accommodation shapes on the deck, the effect of accommodation shapes on wind drag has also been found. Figure 9 shows wind drag acting on the ships with different accommodation shapes on their decks in head wind. A detailed comparison of wind drag acting on the ships is shown in Table 8. The results, as shown in Figure 9, show wind drag coefficients acting on the ship with different frontal accommodation shapes. The results show drastically reduced wind drag acting on the ship by developed frontal accommodations N2 and N3. Table 9 and Figure 10 show detailed differences of wind drag acting on the ship in comparison with those of the original one. The results given in Table 9 and Figure 10 show drastically reduced wind drag acting on the ships; a reduction of up to 66.88% of the total wind drag has been obtained by developed frontal accommodation shape N3. The frontal box shape accommodation N1 and N2 are also smaller wind drag models, with a reduction of up to 33% of wind drag, as shown. The results of wind drag acting on the ships are in good agreement with the obtained pattern of pressure distribution around and over the hull surface of the ship.
Analyzing the pressure distribution around the ships is a useful tool for finding the reason behind the effects on the wind drag acting on the ship hull. Moreover, the results of the interaction effect between the hull and accommodation on wind drag are highly important from the marine container transportation point of view. Tables 10-12 show the results of the interaction effect between hull and accommodation on wind drag acting on the ship with the different frontal accommodation shapes. The results shown in Tables 10-12 clearly show wind drag acting on the ship with the different accommodation shapes in two computed cases, namely the hull with and without accommodation in head wind. The interaction effect between hull and accommodation on wind drag acting on the ship is drastically reduced due to the developed frontal accommodation shapes, as shown. The interaction effect on wind drag acting on the accommodation part is less than 10% and up to 65% of the hull part. The interaction effect on the total wind drag is about 10% of the frontal box shape accommodations N1 and N2, and up to 27% for the streamlined accommodation N3. Figure 11 shows a comparison of the interaction effects on wind drag acting on the ship with different frontal accommodation shapes.
The results show a drastically reduced interaction effect between hull and accommodation on wind drag of the model N3. The interaction effect between hull and accommodation of the original ship is higher than that of the model N3, by up to 70%, as shown. The obtained results are very useful to develop a new accommodation shape with reduced interaction effect and wind drag acting on the ship hull as little as possible.
For a streamlined accommodation shape like the proposed model N3, the interaction effect between hull and accommodation on wind drag has been drastically reduced, as shown. Most researchers are of the opinion that the streamlined accommodation is the best shape to reduce wind drag acting on the ship when in operation. However, the huge streamlined accommodation shape is difficult to be made and will also require a very high cost. Therefore, an accommodation shape which is simply made of only plat plates but still capable of producing a reduced wind drag could be a good choice for the ship building industry.

Conclusions
In this paper, aerodynamic performances of a container ship in head wind have been simulated by CFD. From the obtained results, the interaction effect between hull and accommodation has been investigated. Furthermore, the obtained results in this paper have shown that the newly proposed frontal accommodation shapes can be used to considerably reduce wind drag acting on the container ships, leading to more economic fuel consumption in marine transportation. The main conclusions are summarized as follows: • The effect of mesh number on aerodynamic performances of the reference container ship has been determined. The obtained results show that the effect of mesh number on wind drag acting on the container ship is about 10% when the value of y+ is less than 50 and the same value when y+ is less than 7.

•
The CFD results for the two computed cases, namely the original ship hull with and without accommodation, have been obtained. The effect of the hull on the wind drag is as high as 42%, while the interaction effect between hull and accommodation of the container ship on the total wind drag has been evaluated and is only about 10%.

•
The newly proposed frontal accommodation shapes' effect on aerodynamic performance as well as wind drag acting on hull of the ship. It has been shown that the developed frontal accommodation shapes affect both the pressure distribution around the ship hull and the separation regions around the accommodation. Hence, a suitable shape of the frontal accommodation could be installed to greatly reduce the total wind drag.

•
The effect of accommodation shapes on wind drag acting on hull of the ship has been investigated. Using the proposed frontal accommodation shapes, up to 30% of the total wind drag in head wind (N1), and up to 66.88% of the total wind drag acting on the hull at a Reynolds number of 8.65 × 10 7 in head wind (N3), could be reduced.

•
The interaction effect between the hull and accommodation of the ship, with different accommodation shapes, has been investigated. Drastically reduced interaction effects on wind drag acting on the hull, using the accommodation shape N2 (up to 49%) and using the accommodation shape N3 (up to 70%), have been shown.