Sunday 22 July 2018

11x4.7 Aeronautic Propeller Characteristics (Using CFD) (Verified and Validated) (Update 02)

     This post presents the results from an aeronautic propeller CFD analysis.

     An 11x4.7 propeller was modelled using SolidWorks CAD package using the geometry from [1]. The simulations were run at two different rotational velocities and each rotational velocity was simulated at three advance ratios. The mesh for the 3,000 RPM rotational velocity had 206,184 total cells among which 22,103 cells were at the solid fluid boundary. While, the mesh for the 6,000 RPM rotational velocity had 357,300 total cells among which 64,012 cells were at the solid fluid boundary. A mesh control was employed to refine the mesh near the propeller geometry and at the boundary of the rotating region and the stationery domain for all of the cases simulated. This was done to ensure accuracy of the results was within an acceptable range. The results of the numerical simulations are plotted along with the experimental results [1] in Fig. 1.



Fig. 1 J= Advance Ratio, ηprop = propeller efficiency
 
     It can be seen from Fig. 1 that the trends for the propeller efficiency are in agreement with the experimental results. The fine mesh had the number of cells in each of the respective co-ordinate directions increased by a factor of 1.1. The mesh is shown in the Fig. 2.
 
Fig. 2 The computational mesh around the propeller.
 
     The computational domain size was at 2D x 2D x 2.4D, D being the propeller diameter, as shown in Fig. 3. In Fig. 3, the curved teal arrow represents the direction of rotation of the sliding mesh. The blue arrow represents the direction of free stream velocity while the brown arrow represents the force of gravity.

Fig. 3 The computational domain.

Fig. 4 The pressure distribution and the velocity vectors around the propeller.

     The CAD model and numerical simulation setup files are available here.
 
     Thank you for reading. If you'd like to collaborate on research projects, please reach out.

     [1] Brandt, J. B., & Selig, M. S., “Propeller Performance Data at Low Reynolds Numbers,” 49th AIAA Aerospace Sciences Meeting, AIAA Paper 2011-1255, Orlando, FL, 2011.
doi.org/10.2514/6.2011-1255
 

Update 01

     Mesh independent test results are now available.
 

Update 02

     CAD files for the propeller including the CFD analysis setup are now available.

Monday 9 July 2018

Desktop Computer Part List (Summer 2018-Pakistan Market)

     At the time of writing, 1 USD = Rs. 121.936. The prices mentioned in this post are based on the local market prices of computer components in Pakistan. Please remember that, prices vary from city-to-city within the country and shop-to-shop within a city. This is the reason why a price range is mentioned.

CPUs

     Intel Core i7-8700 for Rs. 41,500-Rs. 43,500. Top of the line processor from Intel. Always prefer a Core i7 or a Core i9 processor.

     Intel Core i5-8400 Processor Rs. 26,000-Rs. 28,000. Only buy this processor if there is a budget constraint.

Motherboard

     Gigabyte Z370 AORUS Gaming 3 for Rs. 21,200-Rs. 21,500. This motherboard has many USB ports and also comes equipped with a USB type-C port and M.2 slots etc. for future proofing.

Storage

     WD Blue 500GB Solid State Drive - WDS500G1B0A for Rs. 16,500-Rs. 17,500. Please do not buy a hard drive with rotating mechanism, it's 2018! When later in the year the SSD's prices go down, probably around November-December 2018, then buy another ~500 GB SSD. Do not waste money on a legacy hard drive.

Memory

     Corsair Vengeance LPX 16GB (1x16GB) DDR4 DRAM 3000MHz Rs. 25,800-Rs. 32,000. Please do not buy two sticks of 8 GB each. Save the remaining memory slots for future upgrading. Memory prices will fall significantly in October-November 2018 once the Chinese memory plants become operational. Don't fall for the shop keepers trickery. A common ploy employed by shopkeepers is that 1x16GB memory modules don't work in single channel mode for the 2400 MHz+ modules. It works perfectly well.

Casing

     Corsair Carbide Series® 100R Mid-Tower Case Rs.6,650-Rs. 7,500. This is the best option, really. Don't waste money on casing, it's just a box.

Power Supply

     Corsair VS550 - 550 Watt Power Supply Rs.5,800 -Rs. 6,000. A 550 Watt power supply for an i7 8700 CPU, 4 sticks of 1x16GB DDR4 memory modules, 1 SSD and up to a GTX 1x70 level graphic cards without any over clocking and 16 hours per day usage will be enough. May be this power supply will even be enough for GTX 1180 graphic card, as the graphic card chips are becoming more and more energy efficient. Yet again, do not fall for the shop keepers ploys.

Graphic Card

     Wait for the graphic new graphic cards from NVidia. The new cards are just around the corner. Local shops in Pakistan are selling 2 year old graphic cards, the GTX 10 series, at much inflated prices as compared to the rest of the world.

Conclusion

     This system will last at least 5 years, in terms of gaming with a GTX 1x70 level graphics card. It will even perform well for 10+ years if you keep upgrading it and take care of it cooling wise. Upgrade to a PCIe-NVMe SSD down the road when the prices drop, add more memory and update the graphic card every 5 years etc. Currently, the system will cost anywhere between Rs. 118,000 -Rs. 128,000, depending on you city and the shop, with an i7 processor.

     If you'd like to collaborate on research projects, please reach out. Thank you for reading.

Monday 11 June 2018

Marine Propeller Characteristics (Using CFD) (Verified and Validated)

     This post presents the results from a marine propeller CFD analysis. The key thing about this CFD analysis was that the propeller efficiency obtained was within 7% of the experimental results, by using only 112,081 total cells in the computational mesh, of these cells, 21,915 cells were at the solid fluid boundary. The results are, indeed, mesh independent. The software employed was Flow Simulation Premium. The results of the numerical simulations are plotted along with the experimental results in Fig. 1.

Fig. 1 KT = coefficient of thrust, 10KQ = coefficient of torque multiplied by a factor of 10, ηprop = propeller efficiency
 
     It can be seen from Fig. 1 that the trends for the thrust and the torque coefficients and the propeller efficiency are in agreement with the experimental results. The experimental data was taken from here. The flow conditions were following. The propeller diameter was at 0.254 m. Propeller rotational velocity was at 15 rev/s. The propeller inclination angle was at 12°. Fluid considered was water. To change the advance ratio, fluid flow velocity was altered. The computational domain size was at 2D x 2D x 3.2D, D being the propeller diameter, as shown in Fig. 3. The mesh is shown in the Fig. 2.

Fig. 2 The computational mesh around the propeller.
 
Fig. 3 The computational domain.
 
     In Fig. 3, the curved teal arrow represents the direction of rotation of the sliding mesh. The blue arrow represents the direction of free stream velocity while the brown arrow represents the force of gravity. The Pressure distribution on the propeller blades and the velocity streamlines are shown in Fig. 4. The streamlines were drawn using line integral convolution, relative to the rotating frame of reference.
 
Fig. 4 The pressure distribution and the velocity profile around the propeller.
 
     If you'd like to collaborate on research projects, please reach out. Thank you for reading.

Sunday 1 April 2018

Comparison of VAWT Blade Designs (Leading-Edge Tubercle, Leading and Trailing-Edge Tubercle, Unmodified) (Update 05)

Numerical simulations were run on SolidWorks Flow Simulation Premium (model files are available here) software to compare the torque characterizes of three distinct vertical-axis wind turbine blade designs shown in Fig. 1. The torque characteristics are shown in Fig. 2.

This publication was used to verify and validate the numerical methodology. The results were within 8% of the publication's results at the design point of TSR of 1.2 at 90 RPM and 7.85 m/s wind speed. The dimensions of the turbine, the  blades and the cross section used are mentioned in the publication.

Fig. 1. Top Row, L-R: VAWT with blades having tubercles at the leading edge (ten tubercles per blade span, configuration name 10T), VAWT with blades having tubercles at both the leading and the trailing edge (ten tubercles per blade span). Bottom Row, VAWT with blades having no modifications.

It is clear from the Fig. 2 that the baseline design provides the most stable torque. On average the turbine with no modifications on the blades produced 5.31 Nm torque in one complete rotation, while the turbine with tubercles at the leading edge only, produced 5.20 Nm torque. The turbine with tubercles added to both the leading and the trailing edge produced 5.09 Nm torque in one complete rotation.

The peak torque was maximum for the turbine with the leading edge tubercles, followed by the turbine with the tubercles added to both the leading and the trailing edge of the turbine blades and the turbine with no modifications on the blades at 21.59 Nm, 21.45 Nm and 20.58 Nm respectively.

Fig. 2. Top Row, L-R: Torque curves for VAWT with blades having tubercles at the leading edge, Torque curves for VAWT with blades having tubercles at both the leading and the trailing edge. Bottom Row, Torque curves for VAWT with blades having no modifications. Three colors denote each of the blades in the turbine.

CFD post processing will be added later (may be next week). The effect of leading edge tubercle geometry will be investigated next. The blade design with tubercles added to both the leading and the trailing edge will not be investigated further because it produced the lowest average torque and second highest peak torque.

Update 01:
Decreased the number of tubercles per unit length of the blade, i.e. made the wavelength of the tubercles longer, kept the sweep angle same. As a result, the average and peak torque decreased to 4.53 Nm, and 19.33 Nm, respectively. The figure is attached.


Fig. 3. T-B: Torque curves for VAWT with blades having large wavelength tubercles at the leading edge (five tubercles per blade span, configuration name 5T45). Three colors denote each of the blades in the turbine. Render of the blades.

Update 02:
Increased the number of tubercles per blade span, i.e. made the wavelength of the tubercles smaller, kept the sweep angle same. As a result, the average and peak torque increased to 5.80 Nm, and 23.36 Nm, respectively. The figure is attached.


Fig. 4. T-B: Torque curves for VAWT with blades having smaller wavelength tubercles at the leading edge (fifteen tubercles per blade span, configuration name 15T45). Three colors denote each of the blades in the turbine. Render of the blades.
Update 03:
Again, increased the number of tubercles per blade span, i.e. made the wavelength of the tubercles smaller, kept the sweep angle same. As a result, the average and peak torque increased to 6.1 Nm, and 24.12 Nm, respectively. The figure is attached.


Fig. 5. T-B: Torque curves for VAWT with blades having smaller wavelength tubercles at the leading edge (twenty tubercles per blade span, configuration name 20T45). Three colors denote each of the blades in the turbine. Render of the blades.
Update 04:
Once more, increased the number of tubercles per blade span, i.e. made the wavelength of the tubercles smaller, kept the sweep angle same. As a result, the average and peak torque increased to 6.42 Nm, and 24.63 Nm, respectively. The figure is attached.


Fig. 6. T-B: Torque curves for VAWT with blades having smaller wavelength tubercles at the leading edge (twenty-five tubercles per blade span, configuration name 25T45). Three colors denote each of the blades in the turbine. Render of the blades.
A table for the tubercle geometry is shown below.

Table 01, Tubercle Geometry
Configuration Name
Amplitude (m)
Wavelength (m)
Sweep Angle (°)
Baseline
0
0
0
5T45
0.12777778
0.25555556
45
10T45
0.06052632
0.12105263
45
15T45
0.03965517
0.07931034
45
20T45
0.02948718
0.05897436
45
25T45
0.02346939
0.04693878
45

It is evident from Table 2 that adding more tubercles to the wind turbine's blade causes an increase in both the peak and the average torque. But it is also clear from the Table 2 that the percentage difference in both the average and the peak torque from the previous configuration (less tubercles per blade span) decreases as the number of tubercles per blade span is increased. It appears to be converging to a value.
Table 02, Tubercle Efficiency
Configuration Name
Peak Torque (Nm)
Average Torque (Nm)
Percentage Difference in the Average Torque from the Previous Configuration
Percentage Difference in the Average Torque from then Baseline Configuration
Baseline
20.58
5.31
N/A
N/A
5T45
19.33
4.53
-17.22
-17.22
10T45
21.59
5.2
12.89
-2.12
15T45
23.36
5.8
10.35
8.45
20T45
24.12
6.1
4.92
12.95
25T45
24.63
6.42
4.98
17.29
I think the difference between both the peak and the average torque produced by 25T45 and 20T45 configuration is comparable, up next, a new sweep angle.

Update 05

Following are my publications relating to the subject of this post.

Butt, F.R., and Talha, T., "A Numerical Investigation of the Effect of Leading-Edge Tubercles on Propeller Performance," Journal of Aircraft. Vol. 56, No. 2 or No. 3, 2019, pp. XX. (Issue/page number(s) to assigned soon. Active DOI: https://arc.aiaa.org/doi/10.2514/1.C034845)

Butt, F.R., and Talha, T., "A Parametric Study of the Effect of the Leading-Edge Tubercles Geometry on the Performance of Aeronautic Propeller using Computational Fluid Dynamics (CFD)," Proceedings of the World Congress on Engineering, Vol. 2, Newswood Limited, Hong Kong, 2018, pp. 586-595, (active link: http://www.iaeng.org/publication/WCE2018/WCE2018_pp586-595.pdf).

Butt, F.R., and Talha, T., "Optimization of the Geometry and the Span-wise Positioning of the Leading-Edge Tubercles on a Helical Vertical-Axis Marine Turbine Blade ," AIAA Science and Technology Forum and Exposition 2019, Turbomachinery and Energy Systems, accepted for publication.

Thank you for reading.