Next Article in Journal
A Numerical Simulation of a Variable-Shape Buoy Wave Energy Converter
Next Article in Special Issue
Laser Powder Bed Fusion of a Topology Optimized and Surface Textured Rudder Bulb with Lightweight and Drag-Reducing Design
Previous Article in Journal
Effect of Seal Locations of Pump-Turbine on Axial Hydraulic Trust
Previous Article in Special Issue
VPP Coupling High-Fidelity Analyses and Analytical Formulations for Multihulls Sails and Appendages Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving the Downwind Sail Design Process by Means of a Novel FSI Approach

1
Dipartimento di Ingegneria, Università degli Studi di Palermo, Viale delle Scienze, 90127 Palermo, Italy
2
Facoltà di Ingegneria e Architettura, Università degli Studi di Enna Kore, Cittadella Universitaria, 94100 Enna, Italy
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2021, 9(6), 624; https://doi.org/10.3390/jmse9060624
Submission received: 29 April 2021 / Revised: 26 May 2021 / Accepted: 31 May 2021 / Published: 4 June 2021

Abstract

:
The process of designing a sail can be a challenging task because of the difficulties in predicting the real aerodynamic performance. This is especially true in the case of downwind sails, where the evaluation of the real shapes and aerodynamic forces can be very complex because of turbulent and detached flows and the high-deformable behavior of structures. Of course, numerical methods are very useful and reliable tools to investigate sail performances, and their use, also as a result of the exponential growth of computational resources at a very low cost, is spreading more and more, even in not highly competitive fields. This paper presents a new methodology to support sail designers in evaluating and optimizing downwind sail performance and manufacturing. A new weakly coupled fluid–structure interaction (FSI) procedure has been developed to study downwind sails. The proposed method is parametric and automated and allows for investigating multiple kinds of sails under different sailing conditions. The study of a gennaker of a small sailing yacht is presented as a case study. Based on the numerical results obtained, an analytical formulation for calculating the sail corner loads has been also proposed. The novel proposed methodology could represent a promising approach to allow for the widespread and effective use of numerical methods in the design and manufacturing of yacht sails.

1. Introduction

In recent years, the exponential growth of computational resources at a very low cost has allowed for the widespread use of numerical methods in the nautical field as well as at an industrial level, not only for scientific research or highly competitive purposes. Nevertheless, to date, advanced numerical methods are commonly used in the nautical field only for competitive applications, but it is conceivable that in the future these methods could also be used in other application fields. Recent racing yachts, as demonstrated by the last edition of the America’s Cup, have achieved very high performances thanks to the massive improvements in yacht and sail design, materials, and fabrication. The design of the hulls, sails, and rigging of competitive racing yachts needs more and more detailed research and development that can be obtained by combining advanced computational resources and experimental studies. In particular, with regard to the design of sails, based on traditional and empirical manufacturing processes, nowadays, the best sail designers aim to use more new research and development tools [1].
To improve sail design and to better understand sail aerodynamics, numerical simulations and experimental testing are the most used and reliable methods [2,3,4,5]. These methods have led to a better comprehension of the complex phenomena and the physics underlying sailing and, consequently, through them, a remarkable increase of the performances of yacht sails has been obtained in recent years [6,7]. Experimental full-scale testing is the most accurate and reliable method, but, because of cost constraints, its use is only recommended for a few exceptional cases [7,8] or to validate numerical methods [2,9,10]. Instead, numerical simulations are cheaper than experimental methods, but in order to obtain reliable and accurate results, they need long setup times by highly qualified and experienced users and they require high-performing computational resources.
Among the numerical methods, computational fluid dynamics (CFD) simulations are the most widely used in the design process of high-performance sailing yachts [8]. CFD analyses can be effectively used to design high-performance hulls [11] and to simulate both upwind and downwind sail configurations [12,13,14,15,16,17]. Usually, the flow in upwind sailing is mostly attached and the turbulence effect is quite limited; the flow regime around a downwind sail can be highly turbulent [7] and, consequently, simulating real sailing conditions can be a very difficult task. For this reason, in the past, many studies have concentrated on upwind sails [6], and in recent years, many researchers have focused on simulating the complex downwind sailing conditions [6,18,19]. The first numerical studies on downwind sails frequently treated the topic in a simplified numerical way, considering only the flow effects around a rigid structure [17]. However, the physics of downwind sails is by far more complex than the upwind ones, not only because of the highly turbulent flow, but also because these kinds of sails have an inherent unsteadiness, even in quite stable conditions [7]. This happens mainly because downwind sails are made of very light and flexible cloth, and are attached to the yacht’s structure at three points, namely, the head, the clew, and the tack [20]. The shape of a downwind sail is formed by self-generated aerodynamic forces that are strongly affected by the sail shape itself [21]. Moreover, the shape can change remarkably when sailing, depending on the trim settings and wind conditions [7]. All of these aspects must be taken into account when simulating the real sailing conditions of these kinds of sails, and the sail shape must be accurately measured in the flying condition [10,21]. For this purpose, nowadays, several authors focus on the issue of predicting the flying shape a sail develops under the impact of flow forces. In this context, it is well known that the fluid–structure interaction (FSI) is a key feature to study these kinds of problems [11,16,22,23]. Today, sail designers use specific software to define the constructed (molded) sail shape based on their experience, but they do not have any certainty about the real flying shape [6]. The possibility of predicting the real flying shape of a downwind sail could allow designers to improve the sailing performances of yachts [5,7,24,25], for instance, by maximizing the driving forces. However, knowing the real flying shape could also allow for optimizing the mechanical characteristics or the manufacturing process of the sails, in order to improve their stability, to reduce their weight, and to optimize their mechanical behavior [18,22,26,27]. Many papers have focused on evaluating the real shape of sails only to predict sailing performances [20,22,25,28]. However, to the best of our knowledge, no relevant papers have tried to measure the flying shape and the loads of a sail, taking into account both the point of view of sailors, whose primary interest is the sail performance, and of the sailmakers, who are also responsible for suitably dimensioning and manufacturing the sail.
In this work, a new approach based on the perspective of both sailors and sailmakers is proposed. In particular, the authors propose calculating the flying shape and the loads on downwind sails, not only to estimate propulsive forces [2,10], but also to help sailmakers in dimensioning and defining the best arrangement of the sail reinforcements. For this purpose, a weakly coupled FSI procedure [16] is used to find the real flying shape and to estimate the loads on the head, the clew, and the tack of the sail when the sailing conditions change. On the basis of the numerical results obtained, an analytical formulation for the calculation of the forces on the sail is also proposed.

2. Materials and Methods

In this work, to determine the flying shape and to evaluate the loads on downwind sails, a numerical FSI approach was used. Starting from a previous work of the authors [16], a new weakly coupled method was developed and tested. The new procedure is parametric and fully automated, and it only needs some input data to calculate the flying shape and the loads on the sail. As it is known, to perform a numerical FSI analysis, two embedded problems, aerodynamic and structural, need to be solved together [29]. In this study, to solve the aerodynamic and structural problems, the commercial CFD solver Ansys CFX and the Ansys Static Structural solver were used, respectively. JavaScript and a python-based programming language were used to develop the algorithm that managed the whole procedure, including the data exchange between the workbench and spreadsheet, which was developed ad hoc.

2.1. FSI Analysis: Setup of CFD Environment

The proposed procedure was developed to solve the FSI problem of a downwind sail but, with the aim to simulate the flow as accurately as possible, and so a mainsail was also considered in the CFD simulations. As mentioned, the procedure is parametric, in order to study different sailing configurations in a very simple and fast way. The main CFD parameters that must be first set to run the FSI analysis are reported in Table 1.
The downwind sail trim angle, Δαg, is a parameter that has been introduced to simulate different sail trims. In particular, the sheet angle of the downwind sail, that is the angle between the tack–clew line and the centerline of the boat (αg in Figure 1), is defined as αg = αg_design + Δαg; where αg_design is the design sheet angle of the downwind sail. In this way, once the specific sail type and its design sheet angle have been chosen, different sail trims can be simulated by simply changing the value of parameter Δαg.
The domain used for the CFD simulation was made parametric by setting the positions of the main boundary surfaces (inlet, outlet, and top) as a function of the maximum sail chord, Cm, as defined in Table 2.
Following a largely used approach, the boundary conditions at the inlet surfaces (Figure 2) were set following the formulation proposed by Richards and Hoxey [30]. Concerning the other boundary surfaces, the outlet was set as “opening condition” [16], the bottom surface was set as a rough wall, and the top surface was set as the inlet. The flow velocity at the top surface was calculated using the same formulation of the velocity profile proposed in [30]. The k-ω SST formulation [31] was chosen as the turbulence model.
Regarding the shapes of the sails, the geometry of the mainsail was fixed while the downwind sail was updated at each simulation according to the deformed shape calculated using the FEM analysis. To optimize the computational costs, the domain was divided into an inner region (grey in Figure 2), meshed using a mixed hexa-tetrahedral mesh with prismatic layers on the surfaces of the sails, and an outer region, where hexahedral elements were used. Mesh refinements were applied to the cells adjacent and close to the sail surfaces. The height of the first prismatic element attached to the surfaces of the sails was set to achieve the condition of y + ≈ 1.
A convergence study was performed to quantify the influence of the grid resolution on the CFD results. The ratio between the distances of the grid nodes is equal to 2 3 . Four similar grids with about 1.5 M (base), 0.75 M, 0.38 M, and 0.18 M cells were investigated. The relative step ratio of the i-th grid is h i = N b a s e N i 3 , with N being the number of cells, resulting in h = 1.00 (base), 1.26, 1.59, and 2.00, respectively. The convergence trend is approximated with the following equation [32]:
S h i = S 0 + C × h i p
The values of coefficients S0, C, and p are computed with the least-squares method. Figure 3 and Table 3 show the drive force computed for the different analyzed grids and the results of the uncertainty estimation. A monotonic convergence can be observed in Figure 3. The value of the estimated uncertainty, Ug, for h = 1 (corresponding to 1.5 M cells grid) is 2.1%; this value can be considered acceptable.
Basing on the grid verification results and the good performance hardware available, we used the 1.5 M cells grid. As a result of the unavailability of experimental data, it is not possible, to date, to present a complete validation of the numerical model based on a comparison of specific CFD results with corresponding experimental data. However, in Section 4, the results obtained with the proposed method are compared with others in the literature. A good level of agreement between the numerical and the experimental results has been noted; this demonstrates that the proposed method is able to give sufficiently reliable results. Of course, a complete validation will be carried out in future work.
The meshed domain and the details of mesh around the sail are shown in Figure 4.
The CFD analyses have been set with these stop conditions:
  • values of the residuals of the continuity and momentum equations lower than 1.0 × 10−5;
  • maximum number of iterations equal to 200.
For all of the analyzed configurations, it was observed that the residuals rapidly decreased under the threshold value (1.0 × 10−5), and the aerodynamic quantities (forces, pressure, etc.) converged towards stable values after a reduced number of iterations. As an example, for AWA = 90° and Δαg = 0°, after 45 iterations, the values of the driving force remained stable around a value of about 47 N; moreover, the standard deviation of the values of the driving force calculated in the last 30 iterations was lower than 0.5%. Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 show some results obtained for AWA = 90° and AWA = 115°. In particular, Figure 5 and Figure 6 show the velocity streamlines that develop from the inlet to the outflow boundary surface. For AWA = 90°, the flow was quite clean around the body of the gennaker, even if some divergences occurred near the head and the leech, indicating flow separation in these areas. For AWA = 115°, the flow was slightly less regular than AWA = 90°, and a larger recirculation region could be evidenced at the downstream.
This phenomenon is also supported by the horizontal velocity contour plots (Figure 7 and Figure 8) at the half downwind sail height. In fact, for AWA = 115° (Figure 8), a larger area with a low flow velocity at the downstream was observed compared with AWA = 90° (Figure 7).
To better visualize the flow structures, the second invariant of the velocity gradient [4], Q, was used. Figure 9 and Figure 10 show the iso-surfaces of Q = 500 s−2. In both cases, it can a strongly turbulent fluid at the downstream can be observed; for AWA = 115° (Figure 10), a large flow vorticity can also be seen at the upper part of the sail.

2.2. FSI Analysis: Setup of FEM Environment

As a result of the assumption of a rigid mainsail, FEM analysis was performed only on the downwind sail. The FEM model is parametric; the main parameters are the sail geometry, the Young’s modulus, the Poisson ratio, and the sail thickness. The geometry of the sail was automatically meshed using eight-node shell elements (Shell 281). This kind of shell element has six degrees of freedom at each node, and it is characterized by a mathematical formulation suitable for simulating thin structures (like sails). Suitable mesh refinements were applied on the head, the clew, and the tack.
A convergence study was performed with five structural meshes, with the following average sizes for the shell elements: 200, 100, 5, 30, and 20 mm. The different analyzed meshes had about 1540, 5700, 9200, 12,870, and 18,574 elements, respectively. The result of the convergence study is shown in Figure 11, where the value of the maximum displacement is plotted as a function of the average element size.
From Figure 11, it can be observed that convergence is obtained (relative differences between two consecutive meshes lower than 5%) for values of the element size lower than 100 mm. Given the result of this convergence study, the mesh with an average size of 50 mm was chosen. This represents a good compromise between the quality of the results and the computing time.
The chosen mesh, composed of about 9200 elements, is shown in Figure 12.
To calculate the loads on the rig, the fluid-dynamic pressure distribution calculated at each step by the CFD module was applied over both sail surfaces. The following boundary conditions were applied on the head and the tack (corners “B” and “C”, respectively, in Figure 12): the displacements along the x, y, and z directions were fully constrained and only the rotations along the x, y, and z axes were allowed. The clew (corner “A” in Figure 12) was joined by a spherical joint to a link element, represented by a black dashed line in Figure 12, which simulates a rope. The end of the rope connected to the sail can move freely along the x, y, and z directions; the other end, instead, can only rotate, but no displacement is allowed. A multipoint constraint approach [16] was used to distribute the constraints over many nodes of the head, tack, and clew regions so to avoid peaks of stress due to the application of a load in a single node.

2.3. FSI Procedure

The flow chart of the developed procedure is shown in Figure 13. To launch the procedure, in addition to the CFD and FEM parameters, the following input parameters must be preliminarily defined:
-
the sail design shape, Sd;
-
the number of steps, Ns, that the FSI analysis is subdivided into;
-
the displacement threshold of the current deformed shape, Td.
Sd represents the shape of the sail as defined by the designer. The geometry (shape) file can be loaded in different file formats (iges, sat, dwg, etc.). As “classic” weak coupled procedures could be very unstable because of the large deformations of the sail under the wind load, in the new proposed procedure, the wind velocity was gradually increased in different subsequent steps, until the design value was reached. In this way, the sail shape gradually changed as the wind speed increased and, consequently, no convergence problems of the structural simulations arose. The parameter Ns is used to define the preliminary number of steps and influences the wind velocity, Vs, imposed at each step. The value of the velocity at the i-th step, Vs, is calculated as follows:
Vs = AWS {1 − 0.5[(Nsi)/(Ns − 1)]2}.
Moreover, to better simulate the progressive deformation of the sail, thus avoiding any abrupt change in geometry between two subsequent steps, a check on the maximum value of the displacement was also introduced at the end of each step. For this purpose, a threshold parameter, Td, has been defined. Td represents the admissible maximum value of the displacement of the deformed shape measured at each step.
A simple user interface developed in MS Excel (Figure 14) allows for setting the numerical parameters and launching the FSI analysis.
As soon as it is started, the implemented procedure works iteratively in the following way: the sail geometry, which during the first iteration is the designed one (Sd), subjected to Vs < AWS, is analyzed through the CFD module. The calculated fluid-dynamic pressure distribution is used as the boundary condition for the subsequent structural simulation performed by the FEM module. When the FEM simulation is completed, the maximum value of the total displacement of the sail, Dmax, is calculated and compared with the threshold value, Td. If Dmax is higher than Td, it means that the sail is deformed excessively, so the current velocity is decreased and the procedure restarts from the CFD simulation. If Dmax is lower than Td, the deformed shape of the sail is extracted from the FEM solver and is used as the updated input geometry for a new CFD analysis. Before starting the new CFD simulation, the current wind velocity, Vs, is increased by an amount that depends on the value of Ns. When Vs = AWS, the procedure stops and the flying shape of the sail is obtained.

2.4. Case Study

A common all-purpose gennaker typically used to equip small size sailing boats, like dinghies, was studied. The sail area was 14.2 m2, the luff edge length was 5.95 m, and the maximum chord (Cm) length was 3.20 m. A nylon typically used for this kind of applications, characterized by a Young modulus of 5500 MPa, was chosen for this sail. As mentioned in Section 1, the aim of this paper is not only to present a new FSI approach to measure the flying shape of a downwind sail, but also to test the proposed method to evaluate the structural loads. In particular, it was investigated how the loads on the clew, the head, and the tack vary depending of different sailing conditions and, on the basis of the numerical results obtained, an attempt was made to find a possible analytical formulation that would allow for estimating loads as the wind speed varied. Based on the typical sailing conditions of these types of sails, we chose to study the gennaker at different values of AWA ranging between 90° and 130°, with 5° increments. Moreover, for each analyzed AWA value, we decided to investigate different configurations of the gennaker by varying its trim angle (Δαg) between 0° and 20° (with increments of 5°), in order to find the best trim depending on the wind angle. This resulted in about 60 different configurations to simulate. For this reason, to limit the computational time, a reduced wind velocity was chosen (AWS = 3 m/s). This choice allowed for investigating all of the sailing configurations in a reduced time using a Dual Intel Xeon E5-2670 processor (12 cores) workstation with 128 GB RAM. The mainsail boom angle (αm) was set to be constant and equal to 35°. The values of all of the parameters are presented in Table 4.

3. Results

3.1. Flying Shapes

Figure 15 and Figure 16 show the flying shapes of the gennaker obtained for two different sailing configurations. In particular, the presented results are the ones obtained for AWA = 90° (Δαg = 0°) and AWA = 110° (Δαg = 15°), respectively. The design shapes are colored in green and the flying ones are blue. The maps of the 3D deviations of the flying shapes from the design ones are also presented in Figure 15 and Figure 16. Different point of views, two from the leech and the luff sides and one frontal, are presented in order to better detect the differences between the real and the flying shape.
It can be noticed that the largest positive displacements occur for AWA = 110° and mainly occur along the leech edge. Moreover, it can be observed that for AWA = 90°, the largest positive displacements mainly involve the upper part, while for AWA = 110°, the lower part of the gennaker is the most deformed.

3.2. Loads on the Corners of the Sail: An Analytical Formulation

For all of the analyzed configurations, the driving force and the loads on the corners of the sail, which are the head, the tack, and the clew, were calculated during all of the steps of the FSI simulation. In this way, it was possible to know how these forces varied as the wind speed increased until the chosen value of AWS. To better investigate the structural loads on the head, the tack, and the clew of the sail, for each of the analyzed configurations, the polynomial trendlines of the corner loads and their analytical formulations were extracted by a specific tool of MS Excel. For example, Figure 17 shows the plots of the corner loads over the wind speed for AWA = 120° and Δαg = 20°. The polynomial trendlines and their analytical formulations are also reported.
As soon as the analytical formulations of the trendlines were calculated for all of the configurations, a careful analysis of the values of the first and second order coefficients was performed. From this analysis, it emerged that, in all of the analyzed configurations, the first order coefficient was much smaller than the second order one. In particular, it was calculated that the first order coefficient was, on average, 3.7% of the second order coefficient. By taking this into account and in order to find a simplified analytical formulation for the calculation of the loads on the corners, we neglected the first order coefficient. For this reason, the following equation has been preliminarily proposed for the calculation of the load F on the corners:
F = CC∙AWS2
where CC [kg/m] is the load coefficient that has to be defined.
To determine the most suitable value of CC to make the Equation (3) usable in different sailing conditions, the second order coefficients were analyzed for all of the configurations. Figure 18 shows the values of the second order coefficients calculated for the clew, the tack, and the head in all of the analyzed configurations.
It can be seen that for each corner (clew, tack, and head), the coefficient values change as AWA and Δαg vary. For example, it can be noted that for the clew corner at AWA = 100° (red circle in Figure 18), the highest value of the coefficient was calculated when Δαg = 0°, while for the head corner at AWA = 120° (blue circle in Figure 18), the highest value was obtained when Δαg = 20°. To simplify the data analysis and to determine the most suitable values of CC to be used in Equation (2), we considered only the maximum values of the second order coefficients at different values of AWA. Figure 19 shows the plots of the maximum values of the second order coefficients for the clew, the head, and the tack at different values of AWA.
From Figure 19, it can be observed that the trends of the maximum values of the second order coefficients are substantially constant up to the value of AWA≈110°; beyond this value, a constant decrease can be. Based on this data, adopting a conservative approach, it was decided to use the average value of the second order coefficients calculated between AWA = 90° and AWA = 110° as the load coefficient CC for Equation (3). Table 5 reports the values of the second order coefficients and their average values calculated for the clew, the head, and the tack.
Following the proposed approach, the loads on the head, the tack, and the clew could be estimated using the following analytical formulations:
Fclew = CC_clew∙AWS2 = 2.64∙AWS2
Ftack = CC_tack∙AWS2 = 4.51∙AWS2
Fhead = CC_head∙AWS2 = 5.51∙AWS2
The values of the loads evaluated through CFD simulations and by Equations (4)–(6) are plotted in Figure 20.
It can be seen that for AWA < 110°, the values of the loads evaluated through the proposed analytical formulations are quite similar to the corresponding ones calculated using the CFD simulations. For values of AWA > 110°, as expected, the proposed equations returned slightly overestimated values for the loads, but, considering that a correct design approach should take into account the highest values of the loads on a structure, this could be acceptable.

4. Discussion and Conclusions

The study of the fluid–structure interaction is still one of the most interesting topics in order to predict the performances of sails. A sail, in fact, is usually made of thin plastic material or fabric-based composite, and it forms a three-dimensional camber shape with wind pressure. The shape of the sail camber can remarkably change as the wind conditions vary, and this could result in a varied lift and drag performance by the sail [26]. These phenomena are relevant above all in the case of downwind sails, which are only attached to the yacht’s structure at three points—the head, the clew, and the tack.
In this work, the issue of the FSI problem of downwind sails has been addressed by developing and testing a new weakly coupled procedure. Regarding the presented procedure, one of the main innovative features is related to the possibility of simulating the progressive deformation of the sail, through the gradual application of aerodynamic loads. This feature allows for reducing convergence problems and to avoid solution instabilities, in order to overcome the typical drawbacks of classic weak coupled procedures [19,33]. Another key feature of the new procedure is that it is parametric and fully automated. In fact, both CFD and FEM modules have been setup in such a way as to automatically update the models and the boundary conditions when simulation parameters change. Of course, the procedure is parametric within certain limits; if the type of sail to be analyzed changes remarkably, it is necessary to manually create the preliminary setup of the numerical models. Concerning the automation of the procedure, thanks to a specifically developed algorithm that fully manages the process and the data exchange between the CFD and the FEM modules, different sailing conditions can be simulated by simply introducing some parameters in a MS Excel user-interface.
The new procedure has been used to study a gennaker of a small sailing boat to evaluate its flying shape and to calculate the forces on the corners of the sail. The FSI study has been completed with no interruption or convergence problems within the considered sailing conditions. This demonstrates that, for the specific analyzed test case, the new procedure is effective and stable.
An extensive literature search was carried out to find experimental data to validate the numerical proposed methodology. Unfortunately, to the best of our knowledge, in the literature, there are no experimental studies on sails identical or very similar to the one studied here. However, some experimental studies performed on downwind sails show that the results obtained with the proposed methodology can be considered sufficiently reliable. As evidence of this, the values of the drive force coefficient, CFx, have been calculated at different AWAs using the numerical method proposed here (Figure 21) and compared with the values measured experimentally by Motta et al. [34]. Some similarities can be evidenced.
The trends of the coefficient CFX as the AWA varies are very similar; in fact, they are almost constant with a slight increase around 100°–110°. Moreover, the order of magnitude of the coefficients is comparable, with 0.73 being the average value calculated with the numerical methodology proposed here and about 0.58 being the average value measured experimentally. Of course, as the compared sails are different, it is not possible to have identical or very similar values of CFX, but considering that the order of magnitude is identical, it could be deduced that the numerical results are sufficiently reliable. With reference to the loads calculated numerically with the proposed methodology, good levels of agreement were found with other case studies in the literature. In particular, it can be observed that the ratios of the loads on the sail corners and the trends of the forces, when AWA varies, are very similar to those found in other studies [2,20,34]. For example, the force ratios Fclew/Fhead and Ftack/Fhead calculated at AWA = 90° using the new procedure are about 48% and 81%, respectively; the same ratios have been experimentally measured by Deparday et al. [2] for a similar case study and are equal to about 51% and 85%, respectively. Regarding the trends of the loads, it can be observed that, similarly to what was found with the numerical procedure proposed here, as well as in other experimental studies on similar kinds of sails [2,34], a decrease in the forces was detected when AWA > 110°. This good level of agreement between the numerical data presented and other experimental data, even if it cannot be considered a direct validation of the numerical results, demonstrates the proposed method could give enough reliable results and can be used as a useful tool for sail design.
With regard to the analytical formulation for calculating the sail corner loads discussed in Section 3.2, the proposed Equation (3) of course cannot be used to calculate the corner loads for all types of sails. It applies only for the specific kind of gennaker; however, in our opinion, the proposed approach could be interesting because it could be effectively repeated to find other formulations for different sail types. Still, nowadays, in fact, most sailmakers design the reinforcements to apply on the head, the tack, and the clew, based on their experience. Analytical formulations, grouped by similar types of sails, could be an inexpensive and easy to use method to suitably dimension the reinforcements on the sail corners. Moreover, knowledge of the loads of the sail corners could allow for a better comprehension on wind/rig/sail interaction [2,10]. Such an approach could allow for remarkably improving the design and manufacturing process of sails, not just for competitive sailing yachts, but also for cruising ones.
For future work, the authors believe that further developments could be addressed to the setup and validation of a database of analytical formulations in order to evaluate the loads on the sail corners and to the improvements of the numerical models. A database of equations for an initial evaluation of the sail loads could be developed by analyzing numerous sails different in shape and size; in this way, sailmakers could use a simple tool to design the most suitable reinforcements for the sails, based not only on their experience, but also on numerical data. CFD simulations could also be improved by investigating unsteady sailing conditions, while FEM models could be enhanced by introducing additional information on different types of sail materials regarding structure non linearities, and by evaluating the rig deformation. In this way, very accurate and complete information could be obtained and the sail performance could be optimized. In fact, as the developed procedure is based on a weak coupled FSI method, the FEM module could be effectively interfaced with an optimization module to find the best structural parameters for the sails (e.g., optimal thickness of different regions of the sail and most suitable layout of composite materials).
In conclusion, in our opinion, the novel proposed methodology is surely promising, and further enhancements could allow for getting more information on the dynamics of downwind sails and could lead to the widespread use of this procedure in the design and manufacturing of yacht sails.

Author Contributions

Conceptualization, A.C., T.I., A.M., V.N. and D.T.; data curation, A.C.; methodology, A.C., T.I., A.M., V.N. and D.T.; software, A.C. and T.I.; validation, A.C.; writing—original draft, A.C. and T.I.; writing—review and editing, A.C., T.I., A.M., V.N. and D.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors gratefully acknowledge ANSYS Inc. for the support and the issue of academic licenses.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CFDComputational fluid dynamics
FSIFluid–structure interaction
AoAAngle of attack
AWAApparent wind angle
AWSApparent wind speed
TWATrue wind angle
TWSTrue wind speed
CmMaximum sail chord
αmMainsail boom angle
ΔαgDownwind sail trim angle
αg_designDesign sheet angle of the downwind sail
CFXDrive force coefficient

References

  1. Deparday, J.; Bot, P.; Hauville, F.; Augier, B.; Rabaud, M. Full-scale flying shape measurement of offwind yacht sails with photogrammetry. Ocean Eng. 2016, 127, 135–143. [Google Scholar] [CrossRef]
  2. Deparday, J.; Bot, P.; Hauville, F.; Motta, D.; Le Pelley, D.J.; Flay, R.G. Dynamic measurements of pressures, sail shape and forces on a full-scale spinnaker. In Proceedings of the 23rd HISWA Symposium on Yacht Design and Yacht Construction, Amsterdam, The Netherlands, 17–18 November 2014; pp. 61–73. [Google Scholar]
  3. Campbell, I.M.C. A comparison of downwind sail coefficients from tests in different wind tunnels. Ocean Eng. 2014, 90, 62–71. [Google Scholar] [CrossRef]
  4. Viola, I.M.; Bartesaghi, S.; Van-Renterghem, T.; Ponzini, R. Detached Eddy Simulation of a sailing yacht. Ocean Eng. 2014, 90, 93–103. [Google Scholar] [CrossRef] [Green Version]
  5. Ranzenbach, R.; Armitage, D.; Carrau, A. Mainsail Planform Optimization for IRC 52 Using Fluid Structure Interaction. In Proceedings of the 21st Chesapeake Sailing Yacht Symposium, Annapolis, MD, USA, 15–16 March 2013. [Google Scholar]
  6. Durand, M.; Leroyer, A.; Lothodé, C.; Hauville, F.; Visonneau, M.; Floch, R.; Guillaume, L. FSI investigation on stability of downwind sails with an automatic dynamic trimming. Ocean Eng. 2014, 90, 129–139. [Google Scholar] [CrossRef] [Green Version]
  7. Ramolini, A. Implementation of a Fluid-Structure Interaction Solver for a Spinnaker Sail. J. Sail. Technol. 2019, 4, 1–16. [Google Scholar] [CrossRef]
  8. Le Pelley, D.; Morris, D.; Richards, P.; Motta, D. Aerodynamic force deduction on yacht sails using pressure and shape measurements in real time. Int. J. Small Craft Technol. 2016, 158, 123–124. [Google Scholar] [CrossRef]
  9. Renzsch, H.; Graf, K. An experimental validation case for fluid-structure interaction simulations of downwind sails. In Proceedings of the 21th Chesapeake Sailing Yacht Symposium, Annapolis, MD, USA, 15–16 March 2013. [Google Scholar]
  10. Augier, B.; Bot, P.; Hauville, F.; Durand, M. Experimental validation of unsteady models for fluid structure interaction: Application to yacht sails and rigs. J. Wind Eng. Ind. Aerodyn. 2012, 101, 53–66. [Google Scholar] [CrossRef] [Green Version]
  11. Miyata, H.; Akimoto, H.; Hiroshima, F. CFD performance prediction simulation for hull-form design of sailing boats. J. Mar. Sci. Technol. 1997, 2, 257–267. [Google Scholar] [CrossRef]
  12. Augier, B.; Bot, P.; Hauville, F.; Durand, M. Dynamic behaviour of a flexible yacht sail plan. Ocean Eng. 2013, 66, 32–43. [Google Scholar] [CrossRef] [Green Version]
  13. Viola, I.M.; Bot, P.; Riotte, M. Upwind sail aerodynamics: A RANS numerical investigation validated with wind tunnel pressure measurements. Int. J. Heat Fluid Flow 2013, 39, 90–101. [Google Scholar] [CrossRef] [Green Version]
  14. Ciortan, C.; Soares, C.G. Computational study of sail performance in upwind condition. Ocean Eng. 2007, 34, 2198–2206. [Google Scholar] [CrossRef]
  15. Richter, H.J.; Horrigan, K.C.; Braun, J.B. Computational Fluid Dynamics for Downwind Sails. In Proceedings of the SNAME 16th Chesapeake Sailing Yacht Symposium, Annapolis, MD, USA, 1 March 2003. [Google Scholar]
  16. Cirello, A.; Cucinotta, F.; Ingrassia, T.; Nigrelli, V.; Sfravara, F. Fluid–structure interaction of downwind sails: A new computational method. J. Mar. Sci. Technol. 2019, 24, 86–97. [Google Scholar] [CrossRef]
  17. Viola, I.M. Downwind sail aerodynamics: A CFD investigation with high grid resolution. Ocean Eng. 2009, 36, 974–984. [Google Scholar] [CrossRef] [Green Version]
  18. Graf, K.; Renzsch, H. Ranse investigations of downwind sails and integration into sailing yacht design processes. In Proceedings of the 2nd High Performance Yacht Design Conference, Auckland, New Zealand, 14–16 February 2006; pp. 1–12. [Google Scholar]
  19. Lombardi, M.; Cremonesi, M.; Giampieri, A.; Parolini, N.; Quarteroni, A. A strongly coupled fluid-structure interaction model for wind-sail simulation. In Proceedings of the 4th High Performance Yacht Design Conference, Auckland, New Zealand, 12–14 March 2012; HPYD 2012. pp. 212–221. [Google Scholar]
  20. Le Pelley, D.J.L.; Richards, P.J.; Berthier, A. Development of a directional load cell to measure flying sail aerodynamic loads. In Proceedings of the 5th High Performance Yacht Design Conference, Auckland, New Zealand, 8–12 March 2015; HPYD 2015. pp. 66–75. [Google Scholar]
  21. Tahara, Y.; Masuyama, Y.; Fukasawa, T.; Katori, M. Sail Performance Analysis of Sailing Yachts by Numerical Calculations and Experiments. Fluid Dyn. Comput. Model. Appl. 2012, 91–118. [Google Scholar] [CrossRef] [Green Version]
  22. Deparday, J.; Augier, B.; Bot, P. Experimental analysis of a strong fluid–structure interaction on a soft membrane—Application to the flapping of a yacht downwind sail. J. Fluids Struct. 2018, 81, 547–564. [Google Scholar] [CrossRef] [Green Version]
  23. Renzsch, H.; Müller, O.; Graf, K. Flexsail—A fluid structure interaction program for the investigation of spinnakers. In Proceedings of the RINA—International Conference—Innovation in High Performance Sailing Yachts, Southampton, UK, 28–29 October 2015; pp. 65–78. [Google Scholar]
  24. Bergsma, F.M.J.; Moerke, N.; Zaaijer, K.S.; Hoeijmakers, H.W.M. Development of Computational Fluid-Structure Interaction Method for Yacht Sails. In Proceedings of the Third International Conference on Innovation in High Performance Sailing Yachts (INNOVSAIL), Lorient, France, 26–28 June 2013. [Google Scholar]
  25. Abel, A.G.; Viola, I.M. Force generation mechanisms of downwind sails. In Proceedings of the 7th High Performance Yacht Design Conference, Auckland, New Zealand, 11–12 March 2021. [Google Scholar]
  26. Bak, S.; Yoo, J. FSI analysis on the sail performance of a yacht with rig deformation. Int. J. Nav. Arch. Ocean Eng. 2019, 11, 648–661. [Google Scholar] [CrossRef]
  27. Lombardi, M. Numerical Simulation of a Sailing Boat: Free Surface, Fluid Structure Interaction and Shape Optimization. Ph.D. Thesis, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 2012. [Google Scholar]
  28. Motta, D.; Flay, R.; Richards, P.; Pelley, D.L.; Bot, P.; Deparday, J. An investigation of the dynamic behaviour of asymmetric spinnakers at full-scale. In Proceedings of the 5th High Performance Yacht Design Conference, Auckland, New Zealand, 8–12 March 2015; HPYD 2015. pp. 76–85. [Google Scholar]
  29. Peri, D.; Parolini, N.; Fossati, F. Multidisciplinary design optimization of a sailplan. In Proceedings of the MARINE 2015—Computational Methods in Marine Engineering VI, Rome, Italy, 15–17 June 2015; pp. 177–187. [Google Scholar]
  30. Richards, P.; Hoxey, R. Appropriate boundary conditions for computational wind engineering models using the k-ϵ turbulence model. J. Wind Eng. Ind. Aerodyn. 1996, 46–47, 145–153. [Google Scholar]
  31. Collie, S.J.; Gerritsen, M.; Jackson, P. A Review of Turbulence Modelling for Use in Sail Flow Analysis; Report No. 603; Department of Engineering Science, University of Auckland: Auckland, New Zealand, 2001; Available online: https://www.library.auckland.ac.nz/ (accessed on 26 May 2021).
  32. Viola, I.M.; Bot, P.; Riotte, M. On the uncertainty of CFD in sail aerodynamics. Int. J. Numer. Methods Fluids 2013, 72, 1146–1164. [Google Scholar] [CrossRef] [Green Version]
  33. Jackins, C.L.; Tanimoto, S.L. Oct-trees and their use in representing three-dimensional objects. Comput. Graph. Image Process. 1980, 14, 249–270. [Google Scholar] [CrossRef]
  34. Motta, D.; Flay, R.; Richards, P.; Le Pelley, D.; Deparday, J.; Bot, P. Experimental investigation of asymmetric spinnaker aerodynamics using pressure and sail shape measurements. Ocean Eng. 2014, 90, 104–118. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Main CFD parameters and sailing geometric variables.
Figure 1. Main CFD parameters and sailing geometric variables.
Jmse 09 00624 g001
Figure 2. CFD domain: boundary surfaces (top is fully transparent, bottom is white), and inner (grey) and outer regions.
Figure 2. CFD domain: boundary surfaces (top is fully transparent, bottom is white), and inner (grey) and outer regions.
Jmse 09 00624 g002
Figure 3. Drive force vs. relative grid size.
Figure 3. Drive force vs. relative grid size.
Jmse 09 00624 g003
Figure 4. CFD meshed domain and enlarged views of the mesh around the sail.
Figure 4. CFD meshed domain and enlarged views of the mesh around the sail.
Jmse 09 00624 g004
Figure 5. Velocity streamlines for AWA = 90°.
Figure 5. Velocity streamlines for AWA = 90°.
Jmse 09 00624 g005
Figure 6. Velocity streamlines for AWA = 115°.
Figure 6. Velocity streamlines for AWA = 115°.
Jmse 09 00624 g006
Figure 7. Velocity contour plots at the mid-height section for AWA = 90°.
Figure 7. Velocity contour plots at the mid-height section for AWA = 90°.
Jmse 09 00624 g007
Figure 8. Velocity contour plots at the mid-height section for AWA = 115°.
Figure 8. Velocity contour plots at the mid-height section for AWA = 115°.
Jmse 09 00624 g008
Figure 9. Iso-surfaces of Q-criterion for AWA = 90°.
Figure 9. Iso-surfaces of Q-criterion for AWA = 90°.
Jmse 09 00624 g009
Figure 10. Iso-surfaces of Q-criterion for AWA = 115°.
Figure 10. Iso-surfaces of Q-criterion for AWA = 115°.
Jmse 09 00624 g010
Figure 11. Plot of the maximum displacement vs. average element size.
Figure 11. Plot of the maximum displacement vs. average element size.
Jmse 09 00624 g011
Figure 12. FEM mesh of the sail.
Figure 12. FEM mesh of the sail.
Jmse 09 00624 g012
Figure 13. Flow chart of the developed fluid–structure interaction (FSI) procedure.
Figure 13. Flow chart of the developed fluid–structure interaction (FSI) procedure.
Jmse 09 00624 g013
Figure 14. MS Excel user interface.
Figure 14. MS Excel user interface.
Jmse 09 00624 g014
Figure 15. Design shape (green), flying shape (blue), and map of the 3D deviations [m] for AWA = 90° and Δαg = 0°.
Figure 15. Design shape (green), flying shape (blue), and map of the 3D deviations [m] for AWA = 90° and Δαg = 0°.
Jmse 09 00624 g015
Figure 16. Design shape (green), flying shape (blue), and map of the 3D deviations [m] for AWA = 110° and Δαg = 15°.
Figure 16. Design shape (green), flying shape (blue), and map of the 3D deviations [m] for AWA = 110° and Δαg = 15°.
Jmse 09 00624 g016
Figure 17. Plots of the corners loads vs. AWS for AWA = 120° and Δαg = 20°.
Figure 17. Plots of the corners loads vs. AWS for AWA = 120° and Δαg = 20°.
Jmse 09 00624 g017
Figure 18. Plots of second order coefficients vs. AWA at different values of Δαg.
Figure 18. Plots of second order coefficients vs. AWA at different values of Δαg.
Jmse 09 00624 g018
Figure 19. Plots of maximum value of second order coefficients vs. AWA.
Figure 19. Plots of maximum value of second order coefficients vs. AWA.
Jmse 09 00624 g019
Figure 20. Plots of the loads on the sail corners calculated by the CFD simulations and by Equations (4)–(6) for AWS = 3 m/s and different values of AWA.
Figure 20. Plots of the loads on the sail corners calculated by the CFD simulations and by Equations (4)–(6) for AWS = 3 m/s and different values of AWA.
Jmse 09 00624 g020
Figure 21. Drive force coefficients calculated with the numerical proposed method.
Figure 21. Drive force coefficients calculated with the numerical proposed method.
Jmse 09 00624 g021
Table 1. Computational fluid dynamics (CFD) analysis parameters.
Table 1. Computational fluid dynamics (CFD) analysis parameters.
ParameterUnit
Apparent wind speed (AWS)m/s
Apparent wind angle (AWA)deg
Mainsail boom angle (αm)deg
Downwind sail trim angle (Δαg)deg
Table 2. Distance of the boundary surfaces in function of Cm [m].
Table 2. Distance of the boundary surfaces in function of Cm [m].
InletOutletTop
5 × Cm10 × Cm1 × Cm
Table 3. Estimated uncertainty for different grids.
Table 3. Estimated uncertainty for different grids.
Number of CellsDrive Force [N]hUG [N]UG%
1,500,00047.051.000.992.1%
750,00047.281.261.443.1%
375,00048.211.592.545.5%
187,50048.872.003.257.0%
Table 4. Values of the parameters used for the case study.
Table 4. Values of the parameters used for the case study.
ParameterValue
Displacement threshold value (Td)0.035 m
Number of steps (Ns)5
Sail thickness0.7 mm
Young’s modulus5500 MPa
Poisson ratio0.25
Apparent wind speed (AWS)3 m/s
Apparent wind angle (AWA)90° ÷ 120°
Mainsail boom angle (αm)35°
Downwind sail trim angle (Δαg)0° ÷ 20°
Table 5. Second order coefficients: single and average values.
Table 5. Second order coefficients: single and average values.
AWA [°]ClewTackHead
902.694.535.53
952.674.515.54
1002.664.565.59
1052.624.505.48
1102.584.475.40
Average values (CC)2.644.515.51
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cirello, A.; Ingrassia, T.; Mancuso, A.; Nigrelli, V.; Tumino, D. Improving the Downwind Sail Design Process by Means of a Novel FSI Approach. J. Mar. Sci. Eng. 2021, 9, 624. https://doi.org/10.3390/jmse9060624

AMA Style

Cirello A, Ingrassia T, Mancuso A, Nigrelli V, Tumino D. Improving the Downwind Sail Design Process by Means of a Novel FSI Approach. Journal of Marine Science and Engineering. 2021; 9(6):624. https://doi.org/10.3390/jmse9060624

Chicago/Turabian Style

Cirello, Antonino, Tommaso Ingrassia, Antonio Mancuso, Vincenzo Nigrelli, and Davide Tumino. 2021. "Improving the Downwind Sail Design Process by Means of a Novel FSI Approach" Journal of Marine Science and Engineering 9, no. 6: 624. https://doi.org/10.3390/jmse9060624

APA Style

Cirello, A., Ingrassia, T., Mancuso, A., Nigrelli, V., & Tumino, D. (2021). Improving the Downwind Sail Design Process by Means of a Novel FSI Approach. Journal of Marine Science and Engineering, 9(6), 624. https://doi.org/10.3390/jmse9060624

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop