Shape and Spacing Effects on Curvy Twin Sail for Autonomous Sailing Drone

There is an accelerating requirement for ocean sensing where autonomous vehicles can play an essential role in assisting engineers, researchers, and scientists with environmental monitoring and collecting oceanographic data. This paper is performed to develop a rigid sail for the autonomous sailing drone. Our study aims to numerically analyze the aerodynamic characteristics of curvy twin sail and compare it with wing sail. Because racing regulations limit the sail shape, only the two-dimensional geometry was open for an optimization. Therefore, this study’s first objective was to identify the aerodynamic performance of such curvy twin sails. Simultaneously, a secondary objective was to estimate the effect of the sail’s spacing and shapes. A viscous Navier-Stokes flow solver is used for the numerical aerodynamic analysis. The 2D aerodynamic investigation is a preliminary evaluation. The results have shown that the curvy twin sail designs have improved lift, drag, and driving force coefficient compared to the wing sails. The spacing between the port and starboard sails of curvy twin sail is an important parameter. The spacing is 0.035L, 0.07L, and 0.14L shows the lift coefficient reduction because of dramatically stall effect, while flow separation is improved with spacing is 0.21L, 0.28L, and 0.35L. Significantly, the spacing 0.28L shows the maximum high pressure at the lower area and the small low pressure area at leading edges, so the highest lift is generated.


Introduction
The marine environment in the sea is facing problems such as eutrophication, acidification, and oxygen deficiency. The need is growing for ocean sensing to benefit marine safety, environmental monitoring and climate research. An unmanned autonomous vehicle that can assist surveying the sea by acquiring long-term placed bound measurements, reporting pollution, hydrographic and meteorological data could be a helpful and cheap alternative to a commercial research vessel.
The researchers are now developing an autonomous sailing drone, which will be used as an experimental platform for autonomous data acquisition in the sea. The application drives requirements on the vehicle, where it needs to have robustness, reliability, and energy efficiency. The requirement is to develop a vehicle that can survive in any weather conditions all year under free conditions.
Utilization of wing sails in autonomous sailing drones has been an object for research in maritime robotics for over 15 years. Researches on wing sails have been conducted by many researchers. Several attempts to develop such systems have been made, for instance, by (Elkaim, 2001) and (Neal & Sauze, 2008) [1,2]. Elkaim has completed experimental aerodynamic performance of self-trimming wing sail for autonomous surface vehicles [3]. In 2008, Aberystwyth University designed a sailing robot using double wing sail as an aero-propulsion system [4]. A vehicle called "Sail drone" has been launched (Jekins, Meinig, Lawrenece-Salavas, & Tabisola, 2015) [5]. BMW Oracle was updated by a 60-meter high, multi-element wing sail, which won the 33 rd America's Cup in 2010. The rigid wing sail's flexible size and excellent performance have been paid attention by the shipping industry [6]. Recently, Sweden Research Group have developed an autonomous sailboat with free-rotating wing sail [7]. The research in the field shows that wing sails have good aerodynamic performance, good control capabilities, precision, and low power consumption. In order to improve the aerodynamic performance characteristic, different methods are used, such as trailing flap [8,9], leading slat [10], deformed flaps [11], and leading-edge turbulences [12], etc. However, difficulties in designing wing sails for autonomous sailing drone has been identified. Because of working in the long term autonomous oceanography, the rigid sail's simple design and good aerodynamic performance are required. Therefore, the workability improvement of the sail in aerodynamic performance must be considered. A two dimensional curvy twin sail which installed on a Sailing drone is proposed and preliminary investigated.

The Geometry of the Curvy twin sail
The general specifications of Sailing Drone are presented in Figure 1 and Table 1.

Material CFRP
The curvy twin sail's geometry is a simplified 2D configuration of two-element rigid sail as Figure 2 and Table 2. The non-dimensional number for characterizing the flow of fluid, Reynolds number is defined by Equation (1): Where µ is the viscosity coefficient of air, no experimental data of optimized sail was available to validate the numerical results. Therefore, a proper numerical setting was performed to produce reliable results. The aerodynamic of NACA 0018 based on experiments were simulated with Re = 6x10 5 (the wind speed is assumed 8 m/s). Moreover, NACA 0018 was used to develop a wing sail by The Royal Institute Technology. The NACA 0018 wing sail simulation is compared with the curvy twin sail.  The total aerodynamic force can be decomposed in two ways. The first method considers the driving (FR) and heeling (FH) components of the total force. This is used when studying the balance of the aero-hydrodynamic forces. Nevertheless, when the aerodynamic behavior is analyzed on its own, it is common practice to decompose the aerodynamic force in its lift (L) and drag (D) components, as plotted in Figure 3. = cos + D sin It is obvious that the goal is to maximize the driving force and simultaneously minimize the heeling force. When beating against the wind, the drag's contribution is negative because it increases the heeling force and reduces the driving forces. Therefore, in this point of sail, the objective is to maximize the lift and minimize the drag. As the apparent wind angle increases through the reaching point of sail, the drag's contribution to the driving force becomes more important. At the last stage, when downwind, the maximum driving force will equivalent to maximum drag. This work's primary purpose is to study the sail's shape and spacing, which are considered optimized parameters. By two dimensional simulation, the aerodynamic performance of 5 different sail shape was carried out to find the most appropriate curvy shape and spacing. The aerodynamic of NACA 0018 based on experiments were simulated. Moreover, NACA 0018 was used to develop a wing sail by The Royal Institute Technology [7]. The NACA 0018 simulation is compared with the two dimensional curvy twin sail.

Computational Domain
The calculation domain for computational simulation is presented in Figure 4. The sail is located at the origin, while a semi-circle of radius 15L forms the inlet. The cross-stream domain width is 30L, while the down-stream length is 40L. These domain dimensions ensure domain-independent results [13, 14 and 15]. AOA is adjusted by rotating the direction of airflow. The boundary conditions are described in Table 3. The boundary conditions on domain sides are symmetry, while the downstream boundary is outlet. The velocity at the inlet is uniform, and the value is the same as the free flow velocity.

Mesh details and mesh independence
The quality of mesh determines the numerical result's accuracy. Different numerical mesh parameters were examined, and a mesh independency test was done by observing the value of the computed lift coefficient at an angle of attack of 5 degrees. Figure 5 shows the adopted mesh shapes where structured mesh was used [16,17]. The mesh distribution following the hexahedral was also refined to obtain accurate results. Table 4 shows the obtained lift coefficient of the NACA 0018 at an angle of attack of 5 for different mesh sizes included in the mesh independency test. The computed lift coefficient varies from 0.45888 for mesh 1 to 0.4825 for mesh 4. The difference between the lift coefficient obtained by the coarse mesh and fine mesh (mesh 1 and mesh 4) is 5%. However, the difference between mesh 4 results and mesh 3 results is only 0.18%. Therefore, mesh 3 is adopted throughout the present work.

Computational Approach
To investigate the flow phenomena, an appropriate CFD flow solver must be selected. A CFD commercial code software is used to simulate in a parallel computation environment, along with a module for partitioning the volumetric grid as a pre-processor in parallel flow computations. CFD software solver can provide direct flow solutions and ad joint solutions for potential, Euler, Navier-Stokes, and Reynolds Averaged Navier-Stokes (RANS) governing equations. It uses a Finite Volume Method (FVM) for spatial discretization. Both explicit and implicit methods are available for time integration, and central difference or upwind methods can also be used for spatial discretization. The advanced numerical techniques of residual smoothing and structured mesh methods are also available to improve the robustness and convergence of the flow solution. For the numerical analysis, the solver for the two-dimensional k-epsilon realizable governing equations which describe the mass conservation, momentum, and energy in a viscous fluid was used [18,19]. The flow around the sail has a very low speed. It is important to consider the viscous effect for a more accurate flow solution. The setup for the computational model is illustrated in Table 5 [20, 21, and 22].

Simulation type
Single phase flow

Realizable -epsilon
Acting fluid Air at 25 o C

Computational model validation
In order to verify the numerical results, the lift coefficient of the two dimensional NACA 0018 wing sail model in a free flow environment is compared with the experimental results. Figure 6 shows the predicted variation of CL with AOA for the NACA 0018 using the realizable -epsilon turbulence models. A comparison with the experimental results [23], also included in Figure 6, shows that the model is in close agreement with the experiments for AOA varying from 0 to 14 degrees. Therefore, the realizable -epsilon turbulence model is selected for further simulations.

Shape Optimizations
The shape optimization was illustrated in Figure 7. The curvy twin sail's shapes was chosen randomly. Firstly, the lift and drag coefficients are investigated. Figure 8 illustrates the predicted lift and drag coefficients with different sail shapes. To evaluate the performance of curvy twin sail designs, they have been compared with the NACA 0018 wing sail. The curvy twin sail's lift and drag coefficients are superior to the NACA 0018 wing sail. When the angle of attack is from 0 o to 75 o , the lift coefficient of curvy twin sails is higher than NACA wing sail. The maximum lift coefficient of the NACA 0018 wing sail is 1.5 with AOA = 15 o while the curvy twin sails are about 2 with AOA around 35 o . However, due to the curvy shape, the comprehensive performance of curvy twin sail in case 1-1 is higher than others. Therefore, it is necessary to consider sail shape when choosing the shape parameters. Because the curvy twin sail's lift is generated by lift of sail's two components consisting of starboard and port sails, the lift coefficients need to be investigated separately, as shown in Figure 9. Although the drag coefficient in curvy twin sail designs is higher than NACA 0018 wingsail, the drag may also form part of the propulsion force for the sailing drone in downwind conditions. So the aerodynamic performance for the curvy twin sail propose in this paper is more excellent.   The port sail mainly creates the favorable lift force while the starboard sail makes an adverse lift in starboard wind conditions. In contrast, the function of port and starboard sails are interchangeable in port wind conditions. It is convenient only to consider working in starboard wind conditions. The lift coefficient of port sail reaches a peak when the angle of attack is around 15 o , but the total lift coefficient is not maximum. It is attributed to the starboard sail generated dramatic minor lift. When the angle of attack is over 55 o , the starboard sail is not effected on the lift coefficient. Although there are unequal lifts at the port and starboard sail in curvy twin sail' case optimization, the difference in drag among them is unremarkable. The flow phenomenon between the port and starboard components of curvy twin sail is different. For port sail, at angle of attack from 5 o to 15 o , the improvement of the lower area smooth flow and high velocity flow area between two sails took place due to the formation of high pressure region at the lower area and low pressure region at the area between two sails. The difference between these high and low pressure region cause a lift on the port sail. As the angle of attack increases, the flow separation appears at both lower and area between two sails. It leads to a reduction in flow velocity, which hinders the port sail's ability to create lift. For starboard sail, the stall is caused by flow separation at almost angle of attack, resulting in raising pressure at the upper area, so the adverse lift is created. Consequently, there is a reduction in total lift generated by curvy twin sail. Therefore, the airflow separation at angle of attack 35 o , the total lift is able to be maximum. When angle of attack increases from 15 o to 90 o , the pressure at the spacing area gradually balances with pressure at the upper area, so the adverse lift reduces to 0. Figure 10 shows the velocity streamline and contours on sails at AOA = 15 o . With the curvy twin sails, the flow separation line appears at the upper area. It explains that the pressure at upper area of twin curvy sails is higher than NACA 0018 wing sails, as shown in Figure 11. Besides, the spacing between the port and starboard sails generates the Laval nozzle effect; hence the increase in velocity make a reduction in pressure at this area. Because of curvy shape of port sails, the pressure distribution of twin curvy sails at the lower area is different and higher than NACA 0018 wing sails.
The high pressure region at the lower area in case 1-1 is more extended than others, so generated lift is better than other cases. It can be seen from Figure 12 that the upper area of curvy twin sails has a large flow separation at AOA = 35 o , and the opening of this area is unnoticeable in comparison with the NACA 0018 wing sail. This phenomenon leads to similar pressure at the upper area, as shown in Figure 13. Nevertheless, the Laval effect still occurs at the area between two sails, the low pressure region maintained. The low velocity at the lower area of curvy twin sails created a higher pressure than the NACA 0018 wing sail. Consequently, the lift generated by curvy twin sails is more effective than NACA 0018 wing sails. The flow separation appears at the starboard sail's upper area for curvy twin sails at almost angle of attack. With the increase in angle of attack, the flow separation of starboard sail of the curvy twin sails expands and becomes larger. It is explained that raising in pressure at the upper area of the starboard sail creates adverse lift. The lift generated in case 1-1 reaches at peak at AOA = 35 o and a little higher than other cases. Because of the curvy shape, the lower area and area between the two sails flow characteristics are different. The smoother flow at the area between two sails in case 1-1 causes lower pressure, as shown in figure 13. This phenomenon leads to little lift improvement in case 1-1 in comparison with other cases.

Driving force coefficient comparison
When the true wind direction (VT) is constant and sailing drone permit sailing at constant speed (VS), the relative wind acting on sails is called apparent wind (VA) as illustrated in figure 14.

Figure 14 The velocity triangle
The relationship among true wind, sailing drone speed and apparent wind is give equations (14,15): To investigate the effect of the flow interaction, which has a significant effect on sail's aerodynamic performance, a range of numerical aerodynamic analyses were conducted for sails by varying the true wind direction from 0° to 180° in intervals of 5°. It is supposed that sailing drone permit sailing at course sailed direction and VS = 0 m/s. Therefore, the apparent wind is equal to the true wind. The requirement in generating driving force at port side and starboard side working conditions is similar, so the symmetric design of curvy twin sail is applied. The apparent wind direction is determined from 0° to 180°; the lift and drag generated by sail are calculated. These forces are resolved into driving force (FR) and the heeling force (FH). The driving force is maximized, corresponding with each certain apparent wind angle. For certain apparent wind speed and angle, the angle of attack adjusted to generate different driving force values, and the maximum driving force is determined. Suppose that the heeling force (FH) is equal to horizontal side force (FS) [16,17]. For apparent wind angles aligned with the entry point of the sail, the sail lift is the predominant component of propulsion. For apparent wind angles behind the sail, lift diminishes, and drag increases as the predominant component of propulsion. The driving force coefficient in all sail designs is illustrated in Figure 15. It is shown that the curvy twin sail is more effective than NACA 0018 wing sail in terms of the driving force.

Figure 15 Driving force coefficient
The driving force is totally calculated by drag and lift. At the most angle of attack, the drag and lift of curvy twin sail are higher than NACA 0018 sail, so the generated driving force is improved. The case 1-1 curvy twin sail is found that the driving force was more effective than NACA 0018 sail at the apparent wind angle from 20 o to 180 o . Moreover, the case 1-1 curvy twin sail is able to generate a greater driving force than other cases. The reason for this phenomenon is attributed to the higher lift and drag coefficient at the different angles of attack. Although the lift and drag ratio is not improved, the lift coefficient development of curvy twin sail can make a driving force improvement.

Spacing optimizations
The two dimensional flow around the curvy twin sail at different spacing S values is studied in detail. In order to better study the effect of the spacing on aerodynamic characteristics and established stall of the curvy twin sail, we analyzed the lift and drag characteristics when the spacing is changed, as shown in Table 6.   Figures 17, 18, 19, and 20. As the velocity streamlines shown, although the velocity flow at the area between two sails increases when the spacing S reduces and leading to a lower pressure region, the upper area separation point of flow is more significant. As a result of the curvy twin sails' upper area flow is more stall, the higher pressure region is formed above the sail. These higher pressure regions reduce the differential pressure causing a lower lift force.   With the increase in spacing S from 0.035L to 0.14L, the upper area's flow separation narrows considerably. However, the upper area's flow separation is unchanged when the spacing S is over 0.14L. It explains that the pressure at the upper area in case 2-1; 2-2; 2-3 is higher than in case 2-4; 2-5; 2-6, as shown in Figure 18. Although the smaller spacing between the port and starboard sails is chosen, the higher velocity flow is generated. The vortex appeared in case 2-1; 2-2; 2-3 leads to a reduction in lift coefficient. It can be seen from Figure 19 that the upper area of sails in case 2-1 has the largest flow separation at AOA = 25 o in comparison with other cases. Although the twin curvy sails shape in spacing optimization is similar, the spacing effects on the flow at the lower area of the curvy twin sail. The spacing 0.28L show the maximum high pressure region at lower area and the small low pressure area at leading edges, so the higher lift force is generated in comparison with other spacing value, as shown in Figure 20.

Driving force coefficient comparison
To further verify the driving force performance of curvy twin sail, the lift and drag coefficients are converted to the driving force coefficients, as illustrated in Figure 21. When the apparent wind angle is less the 20 o , the spacing 0.035L and 0.07L twin curvy sails and NACA 0018 wing sail are unable to generated driving force. In contrast, the spacing 0.28L and 0.35L twin curvy sail can generate driving forces from 10 o to 180 o . Moreover, the case 2-5 corresponding with spacing 0.28L is found that the driving force is most effective.

Conclusions
This study investigated the proposed curvy twin sail and compared it with the NACA 0018 wing sail on the aerodynamic performance. Since sail shape and the spacing between two sail components were used to introduce as optimized parameters. The CFD two dimensional model was generated by using CFD commercial code and validated by published experimental data. The lift, drag, and driving force coefficients were obtained using the same Reynolds number.
The shape optimization results show that the port sail mainly creates the favorable lift force while the starboard sail makes an adverse lift in starboard wind conditions. The port and starboard sails' function is interchangeable in port wind conditions, so the symmetric design is a requirement for the curvy twin sail. When the sail shape changed, the flow separation phenomenon appeared around the sail is considerable. The curvy twin sail showed the performance improvement in lift, drag, and driving force coefficient. This improvement can be attributed to the airstream appearance between two sail components and their curvy shape.
The spacing optimization mainly studies the influence of the spacing on the aerodynamic characteristics of curvy twin sail. When the spacing is 0.035L, 0.07L and 0.14L, the stall has been dramatically affected the lift coefficient reduction, while flow separation improvement with spacing is 0.21L, 0.28L and 0.35L. Especially, the spacing 0.28L shows the maximum high pressure at lower area and the small low pressure area at leading edges, so the higher lift is generated compared to other spacing values.