DrawnOnward
Junior Member
- Joined
- Jul 23, 2015
- Messages
- 18
- Reaction score
- 0
- Location
- Brisbane
I have designed the S410 foil in XFLR5 with the intention of using it as a lifting T foil in the Re 2m to 5m range on the 16.8m proa I'm building. The modelling suggests that this foil exhibits properties that might make it a good choice at Re as low as 0.5m. Of course this is only a computer model. The grey areas surrounding the real world application of ncrit values notwithstanding, I hope to establish, as at best a lay hydrodynamicist, that the design path I have followed has merit. That said, I'll run with the data. Critiques appreciated.
The computer likes the aft loading but the real world may not. I can run an FEA on the foil to determine that a particular layup will work for a given aspect ratio, but haven't the tools or expertise to run CFD. Of course, some other more conservative foil shape may not perform as well in 2d, but be stiffer and consequently more suited to higher aspect ratio loading and thus more efficient.
Figure 1. 3 foils.
Alpha plots require a little too much effort to compare positive lift angles, explaining why this Cp min analysis was conducted for positive values of Cl, rather than a fixed alpha range. It'd be nice to be able to plot foil polars relative their zero lift angle. From it we get an idea of the range of angles of attack the foils can sustain at a given Cp (and thus a given speed). From the second graph, the lift coefficient range at various Cp levels can be found. The S410 was designed to provide the widest possible angle and lift range in the -1 < Cpmin < -0.6 range. This is sufficient for my proa to hit speeds of 25 to 35 knots, explaining why I haven't targeted a lower Cp. I have no desire for more and, and with hindsight, most likely a desire for less speed.
Figure 2. Cpmin values plotted against alpha and Cl.
Cp min graphs seem pretty much universal across various Re and ncrit. One variation found is in the form of a blip in the Cp curve appearing on the plot (between 0 and -1 degrees) due to the laminar flow transition. This feature is evident on the Eppler and the S410, but not the H105. I think I can attribute this to a larger separation bubble. Tom Speer has observed in the excellent 'Foiler Design' thread (http://www.boatdesign.net/forums/sailboats/foiler-design-2447.html) that one of the goals for foil design in the laminar flow Re range was a more uniform rate of transition of the bubble across the surface of the foil, as the flatter transition (refer to the xtr1 graph in Fig. 3) seen in the Eppler and S410 will result in a longer separation bubble and increased turbulence - a problem at low Re.
Figure 3 compares the 3 foils at Re of 3m. This isn't the ideal comparison point for the H105, but I'm interested in thoughts regarding sacrificing L/D for slower bubble transition and at what Re the benefit diminishes. I presume that the Cl/Cd graph for the relevant Re contains that information. Refer to figure 5 which compares the S410 and H105 at lower Re. The H105 looks superior through the majority of the useful lift range from something like 0.5m Re down. I'm not concerned, given that 0.5m Re is roughly 3 knots for my craft and that foiling won't have higher L/D than displacement for a 16.8m craft until something like 15 knots of boat speed (the windward foil may be deployed sooner as this hull is only 11.2m long). This is however relevant for craft with smaller chorded foils and lower takeoff speeds - and perhaps unless the foil is also retractable, which is something I'm considering.
Figure 3. 3 foils at Re 3m
Figure 4 compares the L/D of the S410 with the Eppler 817. A small uniform advantage at 2m Re is lost at higher Re in the lift range targeted at higher speeds. This a direct tradeoff on my part for a wider Cp min bucket.
Figure 4. Eppler 817 and S410. Re 2 to 5m.
Figure 5. H105 and S410. Re 0.4m to 1.2m.
The configuration I'm considering is rotating T foils close to each leeward bow and a single surface piercing foil to windward. As speed increases (along with the wind speed increasing) the apparent wind moves forward and the heeling force increases. This increased heeling force: increases the load on the leeward T foils, thereby reducing the required Cl reduction; and decreases the load on the windward foil, making it relatively more suited to a tapered surface piercing foil. The loads upon these T foils and the supporting structure will not be trivial. So all this is still subject to an engineering analysis. I'm mentioning this as the reader might be curious concerning the details of the application, but it is presently only a subject for speculation. Once I've selected a suitable foil section and span, I'll run an FEA on the foil and supporting structure.
Whereas the surface piercer would not require a flap (but perhaps benefit from some pitch adjustment), the T foils will.
Figure 6. The range of flap angles under consideration. Flap hinge at 80%.
This is a small flap. I trialled larger flaps up to 33% of chord in size. They tended to be more draggy and have higher Cp peaks for positive angle. This reduced the benefit of somewhat less negative angle required for a given lift reduction.
It is interesting to note (Fig 7) that positive flap angles feature a very low minimum drag. The sharpening drag increase with progressively more positive flap I imagine is attributable to the rapid transit of the correspondingly larger separation bubble (refer to xtr1 in Fig 7). The biggest problem with positive flap however is the significant increase in Cp min (Fig 7). Consequently, I consider it undesirable to use flap angles greater than +4° with this foil. I can't comment on whether a bow down induced high flap angle on a moth at high speed will cause cavitation without analysing some other foil shapes, but the sharp junction at +8° doesn't look good.
In contrast, the negative flap angles are progressively more draggy, but maintain low Cp min and feature a more benign separation bubble. What is particularly of interest is the Cp min bucket shifting with the zero lift angle - effectively widening the Cp min bucket by an amount equivalent to almost a 6° pitch adjustment.
Figure 7. 2d analysis of flapped foil.
The next step in this process involved 3d foil analysis. From Fig. 12 in the NACA theoretical and experimental analysis of lift and drag http://naca.central.cranfield.ac.uk/reports/1955/naca-report-1232.pdf, which Tom Speer has referred to, it can be taken that a biplane modelled in XFLR5 can account for the free surface effects experienced by a submerged hydrofoil - where the biplane foil separation is twice that of the foil and surface separation. At least the lift can. The NACA report develops a specific set of formulae to account for drag.
First, lift. In figure 8, I have modelled 5 different chord depths below the free surface. The smallest (0.0283c) is the smallest separation that produced useful results. The single foiler is the base line, i.e. infinite depth of immersion. In the Cl vs alpha graph the foil closest to the surface generates a tad more than half the lift of the single foil - excepting the cross over at low positive values of Cl! I'm curious as to why this is happening. Perhaps the result would be more consistent with the theory if the foil was more symmetrical? Perhaps if the foil wasn't tapered? Note that the chord depth is measured relative to the chord at the centre of the foil, which tapers from 420 to 280 mm. I'm aware that the NACA report 'neglects' geometric camber in its calculation of the position of the image foil and an "infinite array of images is, of course, required to give an exact value". Does XFLR5 handle biplane theory in the same way?
Drag. I can't really be sure that the drag results accord with the empirical data from the NACA report, although they have similarites in magnitude (allowing for large Re differences). The comparison with the NACA study could obviously be made clearer by using the same foil sections and chords at the same Re as found in the study. If necessary, I can do this but I'd rather not bother if I'm on the right track. If the Cl vs Cd graph in Fig 8 is a reasonable representation of the effect that proximity of a free surface has upon the foil then none of these questions are crucial. My goal is merely to design a foil that will operate successfully at the targeted Re allowing for free surface effects. If the free surface adjustments merely entail lift and drag penalties of the order of 10 to 20% below 0.5 chords depth, then this adjustment is probably less sensitive than the selection of ncrit and accuracy of manufacture.
Figure 8. Hydrofoil depth analysis
The Cl/Cd graph is barely influence by changes in speed in the target velocity range - the only change is due to changes in viscous drag.
Fig 8a. Viscous drag at 8 and 18 m/sec and various depths
The final figure displays the results of a fixed lift analysis of the 3d aspect ratio 9 flapped biplane foil, without strut, at a depth of 1 metre (2.38 x max chord).
I'll just make a couple of points at this stage. Along with Cp min from the 2d analysis, the Cl/Cd max for various flap angles tends to track with angle of attack. That is, for a given angle of attack, the flap angle required to generate the fixed quantity of lift, will be operating close to its Cl/Cd max and Cp min. The polars in the Vx vs alpha plot track for a given flap angle, the speed required to generate a fixed amount of lift, across the range of angles of attack. Points in the graph that lie on the concave side a polar, will generate lift in excess of that required for that flap angle in static equilibrium.
Figure 9 3d flapped foil analysis at a depth below the free surface of 2.38 x max chord
The computer likes the aft loading but the real world may not. I can run an FEA on the foil to determine that a particular layup will work for a given aspect ratio, but haven't the tools or expertise to run CFD. Of course, some other more conservative foil shape may not perform as well in 2d, but be stiffer and consequently more suited to higher aspect ratio loading and thus more efficient.
Figure 1. 3 foils.
Alpha plots require a little too much effort to compare positive lift angles, explaining why this Cp min analysis was conducted for positive values of Cl, rather than a fixed alpha range. It'd be nice to be able to plot foil polars relative their zero lift angle. From it we get an idea of the range of angles of attack the foils can sustain at a given Cp (and thus a given speed). From the second graph, the lift coefficient range at various Cp levels can be found. The S410 was designed to provide the widest possible angle and lift range in the -1 < Cpmin < -0.6 range. This is sufficient for my proa to hit speeds of 25 to 35 knots, explaining why I haven't targeted a lower Cp. I have no desire for more and, and with hindsight, most likely a desire for less speed.
Figure 2. Cpmin values plotted against alpha and Cl.
Cp min graphs seem pretty much universal across various Re and ncrit. One variation found is in the form of a blip in the Cp curve appearing on the plot (between 0 and -1 degrees) due to the laminar flow transition. This feature is evident on the Eppler and the S410, but not the H105. I think I can attribute this to a larger separation bubble. Tom Speer has observed in the excellent 'Foiler Design' thread (http://www.boatdesign.net/forums/sailboats/foiler-design-2447.html) that one of the goals for foil design in the laminar flow Re range was a more uniform rate of transition of the bubble across the surface of the foil, as the flatter transition (refer to the xtr1 graph in Fig. 3) seen in the Eppler and S410 will result in a longer separation bubble and increased turbulence - a problem at low Re.
Figure 3 compares the 3 foils at Re of 3m. This isn't the ideal comparison point for the H105, but I'm interested in thoughts regarding sacrificing L/D for slower bubble transition and at what Re the benefit diminishes. I presume that the Cl/Cd graph for the relevant Re contains that information. Refer to figure 5 which compares the S410 and H105 at lower Re. The H105 looks superior through the majority of the useful lift range from something like 0.5m Re down. I'm not concerned, given that 0.5m Re is roughly 3 knots for my craft and that foiling won't have higher L/D than displacement for a 16.8m craft until something like 15 knots of boat speed (the windward foil may be deployed sooner as this hull is only 11.2m long). This is however relevant for craft with smaller chorded foils and lower takeoff speeds - and perhaps unless the foil is also retractable, which is something I'm considering.
Figure 3. 3 foils at Re 3m
Figure 4 compares the L/D of the S410 with the Eppler 817. A small uniform advantage at 2m Re is lost at higher Re in the lift range targeted at higher speeds. This a direct tradeoff on my part for a wider Cp min bucket.
Figure 4. Eppler 817 and S410. Re 2 to 5m.
Figure 5. H105 and S410. Re 0.4m to 1.2m.
The configuration I'm considering is rotating T foils close to each leeward bow and a single surface piercing foil to windward. As speed increases (along with the wind speed increasing) the apparent wind moves forward and the heeling force increases. This increased heeling force: increases the load on the leeward T foils, thereby reducing the required Cl reduction; and decreases the load on the windward foil, making it relatively more suited to a tapered surface piercing foil. The loads upon these T foils and the supporting structure will not be trivial. So all this is still subject to an engineering analysis. I'm mentioning this as the reader might be curious concerning the details of the application, but it is presently only a subject for speculation. Once I've selected a suitable foil section and span, I'll run an FEA on the foil and supporting structure.
Whereas the surface piercer would not require a flap (but perhaps benefit from some pitch adjustment), the T foils will.
Figure 6. The range of flap angles under consideration. Flap hinge at 80%.
This is a small flap. I trialled larger flaps up to 33% of chord in size. They tended to be more draggy and have higher Cp peaks for positive angle. This reduced the benefit of somewhat less negative angle required for a given lift reduction.
It is interesting to note (Fig 7) that positive flap angles feature a very low minimum drag. The sharpening drag increase with progressively more positive flap I imagine is attributable to the rapid transit of the correspondingly larger separation bubble (refer to xtr1 in Fig 7). The biggest problem with positive flap however is the significant increase in Cp min (Fig 7). Consequently, I consider it undesirable to use flap angles greater than +4° with this foil. I can't comment on whether a bow down induced high flap angle on a moth at high speed will cause cavitation without analysing some other foil shapes, but the sharp junction at +8° doesn't look good.
In contrast, the negative flap angles are progressively more draggy, but maintain low Cp min and feature a more benign separation bubble. What is particularly of interest is the Cp min bucket shifting with the zero lift angle - effectively widening the Cp min bucket by an amount equivalent to almost a 6° pitch adjustment.
Figure 7. 2d analysis of flapped foil.
The next step in this process involved 3d foil analysis. From Fig. 12 in the NACA theoretical and experimental analysis of lift and drag http://naca.central.cranfield.ac.uk/reports/1955/naca-report-1232.pdf, which Tom Speer has referred to, it can be taken that a biplane modelled in XFLR5 can account for the free surface effects experienced by a submerged hydrofoil - where the biplane foil separation is twice that of the foil and surface separation. At least the lift can. The NACA report develops a specific set of formulae to account for drag.
First, lift. In figure 8, I have modelled 5 different chord depths below the free surface. The smallest (0.0283c) is the smallest separation that produced useful results. The single foiler is the base line, i.e. infinite depth of immersion. In the Cl vs alpha graph the foil closest to the surface generates a tad more than half the lift of the single foil - excepting the cross over at low positive values of Cl! I'm curious as to why this is happening. Perhaps the result would be more consistent with the theory if the foil was more symmetrical? Perhaps if the foil wasn't tapered? Note that the chord depth is measured relative to the chord at the centre of the foil, which tapers from 420 to 280 mm. I'm aware that the NACA report 'neglects' geometric camber in its calculation of the position of the image foil and an "infinite array of images is, of course, required to give an exact value". Does XFLR5 handle biplane theory in the same way?
Drag. I can't really be sure that the drag results accord with the empirical data from the NACA report, although they have similarites in magnitude (allowing for large Re differences). The comparison with the NACA study could obviously be made clearer by using the same foil sections and chords at the same Re as found in the study. If necessary, I can do this but I'd rather not bother if I'm on the right track. If the Cl vs Cd graph in Fig 8 is a reasonable representation of the effect that proximity of a free surface has upon the foil then none of these questions are crucial. My goal is merely to design a foil that will operate successfully at the targeted Re allowing for free surface effects. If the free surface adjustments merely entail lift and drag penalties of the order of 10 to 20% below 0.5 chords depth, then this adjustment is probably less sensitive than the selection of ncrit and accuracy of manufacture.
Figure 8. Hydrofoil depth analysis
The Cl/Cd graph is barely influence by changes in speed in the target velocity range - the only change is due to changes in viscous drag.
Fig 8a. Viscous drag at 8 and 18 m/sec and various depths
The final figure displays the results of a fixed lift analysis of the 3d aspect ratio 9 flapped biplane foil, without strut, at a depth of 1 metre (2.38 x max chord).
I'll just make a couple of points at this stage. Along with Cp min from the 2d analysis, the Cl/Cd max for various flap angles tends to track with angle of attack. That is, for a given angle of attack, the flap angle required to generate the fixed quantity of lift, will be operating close to its Cl/Cd max and Cp min. The polars in the Vx vs alpha plot track for a given flap angle, the speed required to generate a fixed amount of lift, across the range of angles of attack. Points in the graph that lie on the concave side a polar, will generate lift in excess of that required for that flap angle in static equilibrium.
Figure 9 3d flapped foil analysis at a depth below the free surface of 2.38 x max chord
Attachments
-
3 foils.jpg50 KB · Views: 622 -
3 foils y.jpg78.9 KB · Views: 710 -
3 foils analysis.jpg148.1 KB · Views: 643 -
3 foils analysis 2.jpg170.3 KB · Views: 709 -
3 foils analysis 3.jpg128.7 KB · Views: 703 -
3 foils analysis 4.jpg195.6 KB · Views: 814 -
80pc flap foil shapes.jpg51 KB · Views: 689 -
80pc flap 2d foil charts.jpg262.6 KB · Views: 829 -
Foil depth analysis - fixed speed.jpg184.7 KB · Views: 784 -
Foil depth analysis -2 speeds viscous drag.jpg116.7 KB · Views: 858 -
Foil lift analysis - fixed lift.jpg199.3 KB · Views: 884
Last edited: