Large-Eddy Simulation of Low-Frequency Flow Oscillations for NACA0012 and Dynarig Sail at Large Attack Angles
Abstract
:1. Introduction
2. Methods
2.1. Mathematical Modeling
2.1.1. Governing Equations
2.1.2. Lambda2 Criterion
2.2. Geometric Model
2.3. Numerical Setting
2.3.1. Solver and Mesh
2.3.2. Boundary Conditions
2.4. Fast Fourier Transform (FFT)
3. Verification and Validation
3.1. Verification
3.2. Validation
4. Results and Analysis
4.1. Force Performance of NACA0012 Sail
- When α = 15°, wiggles in CL are generally located in a small range of 1.225–1.300. Since the flow is not fully separated, the amplitude of the flow oscillations is suppressed and shows physically stable characteristics in the wake. This phenomenon also applies to angles of attack less than 15°.
- With the increase in attack angle, an increment is found in the fluctuation of the force coefficients. The proportion of oscillating amplitude for CL or CD can go up to 20%, 35%, and 50% for α = 30°, 45° and 90°.
4.2. Performance of Dynarig Sail and Comparisons with NACA0012
- When α = 15°, the flow is in a transitional state between attached and detached. The separation point for the NACA0012 sail is located towards the rear of the mid-section, whereas the flow close to the leading edge of the Dynarig sail has already fully separated. This suggests that compared to wing sails like the NACA0012, thin-arc sails like the Dynarig are more prone to flow separation as the attack angle increases. Additionally, at this angle of attack, the NACA0012 sail can accelerate the incoming flow up to 2.81-fold, which is more pronounced than the 1.52-fold acceleration caused by the Dynarig sail.
- At α = 30°, the vortex structures of the NACA0012 and Dynarig sails are similar, and their acceleration effects on the flow field are also comparable. This is different from the significant disparity observed at α = 15°, which demonstrates the impact of the sail’s geometric shape on the flow field.
- At α = 45°, the NACA0012 sail generates a more extensive low-speed wake region compared to the Dynarig sail, indicating a state of complete stall. In contrast, the low-speed region in the wake of the Dynarig sail is primarily concentrated near the surface of the sail and in the vicinity of some detached vortices. This phenomenon results in the forces on the Dynarig sail exhibiting more regular variations compared to those on the NACA0012 sail, as reflected in Figure 7c and Figure 9c.
- At α = 90°, where the wind speed is perpendicular to the chord length of the sail, the NACA0012 sail exhibits a large-scale flow-separation region at its trailing end, with no discernible regularity. In contrast, the wake region of the Dynarig sail displays a stronger integrity of vortical structures, as evident from Figure 9d, which also reveals a marked periodicity in the sail’s force fluctuations. This behavior may be attributed to the differences in geometry, specifically the fact that the Dynarig sail possesses a fore–aft symmetric shape, which is geometrically different from the NACA0012 sail.
4.3. Spectrum Analysis for NACA0012 Sail and Dynarig Sail
- For both the NACA0012 and Dynarig, with the increment in attack angle, the amplitudes of CL and CD in the frequency domain also increase. For instance, the maximum amplitude of CL of both NACA0012 and Dynarig sails at α = 90° is about 10 times than that at α = 15°.
- The dominant frequency (with the largest value of amplitude) is located in the range of 0.5–10 Hz, which clearly indicates a low-frequency oscillation.
- For the NACA0012 sail, a second local maximum-frequency region can be found at 30–60 Hz at both α = 30° and α = 45°. This phenomenon is not found for the Dynarig sail, which might be caused by the geometry difference.
- Compared to the NACA0012 sail, the Dynarig sail’s spectrum has significantly fewer local maxima. This implies that the Dynarig sail comprises fewer frequency components when subjected to low-frequency oscillations. Such a distinction can serve as a reference for frequency-based sail selection.
5. Conclusions
- For the NACA0012 sail, when α = 15°, the flow is not fully separated. The amplitude of the flow oscillations is suppressed and shows physically stable flow in the wake. With increased attack angle, an increment is also found in the amplitude of the force coefficients, where the oscillating percentage of CL or CD can rise to about 20%, 35%, and 50% for α = 30°, 45° and 90°.
- For the Dynarig sail, the amplitude of wiggles of the force coefficient increases with the rise in attack angle, which is similar to the NACA0012 sail. However, the values of CD at α = 90° are close to a state of regular change for the Dynarig sail.
- For both the NACA0012 and Dynarig, with the increment in attack angle, the amplitude of CL and CD in the frequency domain also increase. Different from the Dynarig sail, a second local maximum-frequency region can be found at 30–60 Hz.
- Compared to the NACA0012 sail, the Dynarig sail’s spectrum has significantly fewer local maxima. Such a distinction can serve as a reference for frequency-based sail selection.
- Strouhal numbers for both sails are predominantly below 0.1, signifying a relatively low level, and exhibit a decreasing trend with increasing angles of attack.
- The amplitude of force fluctuations deviates up to 30% and 17% from the time-averaged forces for Dynarig and NACA0012, respectively, which are too significant to be ignored.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
References
- Meinig, C.; Lawrence-Slavas, N.; Jenkins, R.; Tabisola, H.M. The Use of Saildrones to Examine Spring Conditions in the Bering Sea: Vehicle Specification and Mission Performance. In Proceedings of the OCEANS 2015—MTS/IEEE, Washington, DC, USA, 19–22 October 2015; pp. 1–6. [Google Scholar]
- An, Y.; Yu, J.; Zhang, J. Autonomous Sailboat Design: A Review from the Performance Perspective. Ocean Eng. 2021, 238, 109753. [Google Scholar] [CrossRef]
- Sun, Q.; Qi, W.; Liu, H.; Ji, X.; Qian, H. Toward Long-Term Sailing Robots: State of the Art From Energy Perspectives. Front. Robot. AI 2021, 8, 7253. [Google Scholar] [CrossRef] [PubMed]
- Silva, M.F.; Friebe, A.; Malheiro, B.; Guedes, P.; Ferreira, P.; Waller, M. Rigid Wing Sailboats: A State of the Art Survey. Ocean Eng. 2019, 187, 106150. [Google Scholar] [CrossRef]
- Enqvist, T.; Friebe, A.; Haug, F. Free Rotating Wingsail Arrangement for Åland Sailing Robots. In Robotic Sailing 2016; Springer: Berlin/Heidelberg, Germany, 2017; pp. 3–18. [Google Scholar] [CrossRef]
- Selig, M.S.; Donovan, J.F.; Fraser, D.B. Airfoils at Low Speeds; H.A. Stokely: Virginia Beach, VA, USA, 1989. [Google Scholar]
- Souppez, J.-B.R.G.; Arredondo-Galeana, A.; Viola, I.M. Recent Advances in Numerical and Experimental Downwind Sail Aerodynamics. J. Sail. Technol. 2019, 4, 45–65. [Google Scholar] [CrossRef]
- Almutairi, J.H.; AlQadi, I.M. Large-Eddy Simulation of Natural Low-Frequency Oscillations of Separating-Reattaching Flow near Stall Conditions. AIAA J. 2013, 51, 981–991. [Google Scholar] [CrossRef]
- Eljack, E.M.; Soria, J. Investigation of the Low-Frequency Oscillations in the Flowfield about an Airfoil. AIAA J. 2020, 58, 4271–4286. [Google Scholar] [CrossRef]
- ElJack, E. High-Fidelity Numerical Simulation of the Flow Field around a NACA-0012 Aerofoil from the Laminar Separation Bubble to a Full Stall. Int. J. Comput. Fluid Dyn. 2017, 31, 230–245. [Google Scholar] [CrossRef]
- ElAwad, Y.A.; ElJack, E.M. Numerical Investigation of the Low-Frequency Flow Oscillation over a NACA-0012 Aerofoil at the Inception of Stall. Int. J. Micro Air Veh. 2019, 11, 1756829319833687. [Google Scholar] [CrossRef]
- Arredondo-Galeana, A.; Babinsky, H.; Viola, I.M. Vortex Flow of Downwind Sails. Flow 2023, 3, E8. [Google Scholar] [CrossRef]
- William, Y.E.; Kanagalingam, S.; Mohamed, M.H. Ground Effect Investigation on the Aerodynamic Airfoil Behavior Using Large Eddy Simulation. J. Fluids Eng. 2024, 146, 031205. [Google Scholar] [CrossRef]
- Cravero, C. Aerodynamic Performance Prediction of a Profile in Ground Effect With and Without a Gurney Flap. J. Fluids Eng. 2017, 139, 031105. [Google Scholar] [CrossRef]
- Zhu, H.; Yao, H.-D.; Ringsberg, J.W. Unsteady RANS and IDDES Studies on a Telescopic Crescent-Shaped Wingsail. Ships Offshore Struct. 2024, 19, 134–147. [Google Scholar] [CrossRef]
- Dawi, A.H.; Akkermans, R.A.D. Direct Noise Computation of a Generic Vehicle Model Using a Finite Volume Method. Comput. Fluids 2019, 191, 104243. [Google Scholar] [CrossRef]
- Akkermans, R.A.D.; Ewert, R.; Moghadam, S.M.A.; Dierke, J.; Buchmann, N. Overset DNS with Application to Sound Source Prediction. In Proceedings of the Progress in Hybrid RANS-LES Modelling, College Station, TX, USA, 19–21 March 2015; Girimaji, S., Haase, W., Peng, S.-H., Schwamborn, D., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 59–68. [Google Scholar]
- Zeng, Q.; Cai, W. Simulation and Characterization of Low-Frequency Flow Oscillations Over a Symmetrical Sail at Large Attack Angles for an Unmanned Sailboat. In Proceedings of the ASME 2023 42nd International Conference on Ocean, Offshore and Arctic Engineering, Melbourne, Australia, 11–16 June 2023; American Society of Mechanical Engineers Digital Collection: New York, NY, USA, 2023. [Google Scholar]
- Jabbari, H.; Esmaeili, A.; Rabizadeh, S. Phase Portrait Analysis of Laminar Separation Bubble and Ground Clearance Interaction at Critical (Low) Reynolds Number Flow. Ocean Eng. 2021, 238, 109731. [Google Scholar] [CrossRef]
- Eljack, E.; Soria, J.; Elawad, Y.; Ohtake, T. Simulation and Characterization of the Laminar Separation Bubble over a NACA-0012 Airfoil as a Function of Angle of Attack. Phys. Rev. Fluids 2021, 6, 034701. [Google Scholar] [CrossRef]
- Moise, P.; Zauner, M.; Sandham, N.D. Connecting Transonic Buffet with Incompressible Low-Frequency Oscillations on Aerofoils. J. Fluid Mech. 2024, 981, A23. [Google Scholar] [CrossRef]
- Leonard, A. Energy Cascade in Large-Eddy Simulations of Turbulent Fluid Flows. In Advances in Geophysics; Frenkiel, F.N., Munn, R.E., Eds.; Turbulent Diffusion in Environmental Pollution; Elsevier: Amsterdam, The Netherlands, 1975; Volume 18, pp. 237–248. [Google Scholar]
- Abbott, I.H.; von Doenhoff, A.E. Theory of Wing Sections: Including a Summary of Airfoil Data; Courier Corporation: North Chelmsford, MA, USA, 2012; ISBN 978-0-486-13499-4. [Google Scholar]
- Zeng, Q.; Hekkenberg, R.; Thill, C.; Hopman, H. Scale Effects on the Wave-Making Resistance of Ships Sailing in Shallow Water. Ocean Eng. 2020, 212, 107654. [Google Scholar] [CrossRef]
- Eça, L.; Hoekstra, M. A Procedure for the Estimation of the Numerical Uncertainty of CFD Calculations Based on Grid Refinement Studies. J. Comput. Phys. 2014, 262, 104–130. [Google Scholar] [CrossRef]
- Zeng, Q.; Zhang, X.; Cai, W.; Zhou, Y. Wake Distortion Analysis of a Dynarig and Its Application in a Sail Array Design. Ocean Eng. 2023, 278, 114341. [Google Scholar] [CrossRef]
- Sheldahl, R.E.; Klimas, P.C. Aerodynamic Characteristics of Seven Symmetrical Airfoil Sections through 180-Degree Angle of Attack for Use in Aerodynamic Analysis of Vertical Axis Wind Turbines; Sandia National Labs.: Albuquerque, NM, USA, 1981. [Google Scholar]
Item | Symbol | Value | Unit | |
---|---|---|---|---|
NACA0012 | Dynarig | |||
Chord | c | 0.2 | 0.7 | m |
Camber | cam | - | 0.07 | m |
Initial velocity | U0 | 28.20 | 8.06 | m/s |
Reynolds number | Re | 3.6 × 105 | 3.6 × 105 | - |
NACA0012 | Dynarig | ||||
---|---|---|---|---|---|
Cells (Million) | CL | CD | Cells (Million) | CL | CD |
2.71 | 0.7009 | 0.7647 | 2.10 | 0.3186 | 0.2841 |
4.68 | 0.7439 | 0.8035 | 3.63 | 0.3600 | 0.3345 |
8.09 | 0.7737 | 0.8337 | 6.27 | 0.3829 | 0.3616 |
13.98 | 0.7793 | 0.8390 | 10.84 | 0.3876 | 0.3669 |
NACA0012 | Dynarig | |||
---|---|---|---|---|
Key Variable | CL | CD | CL | CD |
ϕ0 | 0.7936 | 0.8566 | 0.3956 | 0.3758 |
α | −0.0143 | −0.0176 | −0.0067 | −0.0074 |
p | 3.4187 | 3.0198 | 4.4672 | 4.6127 |
Error | 0.0429 | 0.0600 | 0.0251 | 0.0281 |
Uncertainty | 5.50% | 7.19% | 6.48% | 7.76% |
Layers | Cells (Million) | CL | CD |
---|---|---|---|
14 | 11.1 | 0.7760 | 0.8312 |
20 | 12.5 | 0.7768 | 0.8361 |
26 | 14.5 | 0.7793 | 0.8390 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Zeng, Q.; Cai, W.; Xu, J. Large-Eddy Simulation of Low-Frequency Flow Oscillations for NACA0012 and Dynarig Sail at Large Attack Angles. J. Mar. Sci. Eng. 2024, 12, 835. https://doi.org/10.3390/jmse12050835
Zeng Q, Cai W, Xu J. Large-Eddy Simulation of Low-Frequency Flow Oscillations for NACA0012 and Dynarig Sail at Large Attack Angles. Journal of Marine Science and Engineering. 2024; 12(5):835. https://doi.org/10.3390/jmse12050835
Chicago/Turabian StyleZeng, Qingsong, Wei Cai, and Junhui Xu. 2024. "Large-Eddy Simulation of Low-Frequency Flow Oscillations for NACA0012 and Dynarig Sail at Large Attack Angles" Journal of Marine Science and Engineering 12, no. 5: 835. https://doi.org/10.3390/jmse12050835
APA StyleZeng, Q., Cai, W., & Xu, J. (2024). Large-Eddy Simulation of Low-Frequency Flow Oscillations for NACA0012 and Dynarig Sail at Large Attack Angles. Journal of Marine Science and Engineering, 12(5), 835. https://doi.org/10.3390/jmse12050835