Αντικείμενο της παρούσας διδακτορικής διατριβής είναι η πειραματική και υπολογιστική διερεύνηση αεροδυναμικής συμπεριφοράς πτερύγων σε διφασική ροή αέρα–νερού και η εφαρμογή σε πτερύγια ανεμοκινητήρων. Αρχικά, γίνεται πειραματική και υπολογιστική μελέτη μονοφασικής ροής αέρα γύρω από αεροτομές, πτέρυγες και πτερύγιο ανεμοκινητήρα και στη συνέχεια μελέτη διφασικής ροής αέρα-νερού γύρω από τα ίδια σώματα. Η σύγκριση μεταξύ των αποτελεσμάτων της μονοφασικής ροής με τα αντίστοιχα της διφασικής ροής αέρα-νερού είναι αναγκαία ώστε να μελετηθούν οι επιπτώσεις της διφασικής ροής αέρα–νερού στην αεροδυναμική απόδοση. Η πειραματική ανάλυση αφορά τη διεξαγωγή πειραμάτων για τη μελέτη της αεροδυναμικής συμπεριφοράς αεροτομών και πτερύγων σε συνθήκες μονοφασικής και διφασικής ροής. Για την προσομοίωση συνθηκών διφασικής ροής αέρα-νερού τροποποιείται η αεροσήραγγα που διαθέτει ήδη το Εργαστήριο με την προσαρμογή ειδικών ακροφυσίων ψεκασμού νερού (συνθήκες βροχής). Για τις ανάγκες των πειραμάτων χρησιμοποιούνται τα μοντέλα αεροτομών και πτερύγων NACA 0012 που συνοδεύουν την αεροσήραγγα και κατασκευάζονται αεροτομή και πτέρυγες S809. Τα πειράματα μονοφασικής και διφασικής ροής γίνονται για την ίδια ταχύτητα αέρα. Για τη διφασική ροή αέρα-νερού εξετάστηκαν τέσσερις διαφορετικές πυκνότητες περιεχόμενης βροχής. Η υπολογιστική ανάλυση γίνεται με το υπολογιστικό πακέτο ANSYS CFD-Fluent. Αρχικά, γίνονται προσομοιώσεις για μονοφασική ροή αέρα γύρω από την αεροτομή NACA 0012, για την οποία υπάρχει πλήθος δημοσιευμένων αποτελεσμάτων, με τρία διαφορετικά μοντέλα τύρβης ώστε να βρεθεί το καταλληλότερο. Ο συντελεστής άνωσης υπολογίζεται με μεγάλη ακρίβεια, σε αντίθεση με το συντελεστή αντίστασης. Το πρόβλημα αυτό οφείλεται στην αδυναμία του Fluent να υπολογίσει το σημείο μετάβασης του οριακού στρώματος από στρωτό σε τυρβώδες. Κρίνεται επομένως αναγκαίο να γίνει σύγκριση του συντελεστή αντίστασης με πειραματικά δεδομένα για πλήρως τυρβώδες οριακό στρώμα. Για ακόμα πιο ακριβή αποτελέσματα αναπτύσσεται αλγόριθμος για τον υπολογισμό του σημείου μετάβασης από στρωτό σε τυρβώδες οριακό στρώμα και γίνονται προσομοιώσεις ορίζοντας την περιοχή αριστερά από το σημείο μετάβασης ως στρωτή και δεξιά από αυτό ως τυρβώδη. Υπολογίζονται οι κατανομές πίεσης και ταχύτητας γύρω από την αεροτομή, καθώς επίσης και τα σημεία ανακοπής, μέγιστης ταχύτητας, αποκόλλησης και επανακόλλησης του οριακού στρώματος. Παρουσιάζονται επίσης οι ροϊκές γραμμές και τα διανύσματα της ταχύτητας γύρω από την αεροτομή. Αντίστοιχες προσομοιώσεις γίνονται και για την αεροτομή S809. Για τη μελέτη του τρισδιάστατου χαρακτήρα της ροής, γίνονται προσομοιώσεις γύρω από πτέρυγα S809. Υπολογίζονται οι συντελεστές άνωσης και αντίστασης, τα σημεία ανακοπής, μέγιστης ταχύτητας, αποκόλλησης και επανακόλλησης του οριακού στρώματος. Επίσης παρουσιάζονται κατανομές της έντασης της τύρβης στην άνω επιφάνεια της πτέρυγας και της συνισταμένης ταχύτητας, της ταχύτητας στη z-διεύθυνση, της έντασης της τύρβης και της επιτάχυνσης της ροής πίσω από την πτέρυγα. Για τη μελέτη της ροής γύρω από περιστρεφόμενο πτερύγιο γίνονται προσομοιώσεις γύρω από το πτερύγιο Phase IV της NREL. Γίνεται μελέτη της κατανομής της αξονικής ταχύτητας πίσω από το δρομέα, της κατανομής της στατικής πίεσης και της έντασης της τύρβης πάνω στην επιφάνεια του πτερυγίου και της κατανομής της στατικής πίεσης σε διάφορα σημεία πάνω στο πτερύγιο. Η υπολογιστική μελέτη της διφασικής ροής αέρα-νερού γίνεται αρχικά για την αεροτομή NACA 0012 με πυκνότητα περιεχόμενης βροχής LWC=30 g/m³, επειδή υπάρχουν αντίστοιχα έγκυρα πειραματικά αποτελέσματα ώστε να γίνει σύγκριση για την εγκυρότητα της διαδικασίας της προσομοίωσης. Στη συνέχεια γίνονται προσομοιώσεις για διφασική ροή αέρα-νερού γύρω από την αεροτομή S809, την πτέρυγα S809 και το περιστρεφόμενο πτερύγιο Phase IV της NREL. Προσομοιώσεις γίνονται επίσης για διαφορετικές πυκνότητες περιεχόμενης βροχής για τη ροή γύρω από τις αεροτομές σε χαμηλό αριθμό Reynolds. Τα αποτελέσματα της διφασικής ροής αέρα-νερού συγκρίνονται με τα αντίστοιχα της μονοφασικής ροής ώστε να προκύψουν συμπεράσματα για τις επιπτώσεις της βροχής στην αεροδυναμική απόδοση. Γίνεται επίσης υπολογισμός του συντελεστή ισχύος του ανεμοκινητήρα σε συνθήκες μονοφασικής ροής αέρα και διφασικής ροής αέρα-νερού. Σε συνθήκες διφασικής ροής αέρα-νερού παρατηρείται υποβάθμιση της αεροδυναμικής απόδοσης, συγκεκριμένα μείωση της άνωσης με παράλληλη αύξηση της αντίστασης. Δυο είναι οι βασικοί μηχανισμοί που επικρατούν και έχουν ως αποτέλεσμα την υποβάθμιση αυτή. Στην επιφάνεια της αεροτομής δημιουργείται ανομοιόμορφο φιλμ νερού που αυξάνει την τραχύτητα και το πάχος της αεροτομής. Τα σταγονίδια καθώς προσκρούουν πάνω στο φιλμ νερού δημιουργούν «κρατήρες» αυξάνοντας την τραχύτητα της αεροτομής. Επίσης, τα σωματίδια νερού διασπώνται κατά την πρόσκρουσή τους πάνω στην αεροτομή σε άλλα σταγονίδια μικρότερης διαμέτρου και μειωμένης ταχύτητας. Αυτό έχει ως αποτέλεσμα τα σταγονίδια αυτά, επαναεπιταχυνόμενα από τη ροή του αέρα να αποσπούν ποσό ενέργειας από το οριακό στρώμα καθιστώντας το πιο ευάλωτο σε αποκόλληση. Στόχος της μελέτης της αεροδυναμικής συμπεριφοράς των πτερυγίων σε διφασική ροή αέρα-νερού είναι η κατασκευή ανεμοκινητήρων υψηλού βαθμού απόδοσης και η παραγωγή φθηνής ενέργειας από την όσο το δυνατόν καλύτερη αξιοποίηση της αιολικής ενέργειας. / The aim of the present doctoral thesis is the experimental and computational study of the aerodynamic behavior of wings in two-phase flow and the application on wind turbine blades. First of all, experimental and computational study of one-phase flow over airfoils, wings and wind turbine blade and afterwards study of two-phase flow over the same bodies is conducted. The comparison of the results between dry and wet conditions is necessary in order to show the effects of two-phase flow at the aerodynamic performance. Wind tunnel tests were conducted to show the aerodynamic behavior of airfoils and wings in one-phase and two-phase flows. To simulate two-phase flow, the wind tunnel of the Fluid Mechanics Laboratory has to be configured with adding commercial rain simulated nozzles. For the experiments NACA 0012 airfoils and wings which come along the wind turbine are utilized and airfoil and wings S809 are constructed. The experiments of one-phase flow and two-phase flow are conducted for the same air velocity. For the two-phase flow four different Liquid Water Contents are examined. For the computational analysis the commercial CFD code ANSYS Fluent is used. In first place, simulations of one-phase flow over the NACA 0012 airfoil are done with three different turbulence models. The NACA0012 airfoil is chosen because it has been studied in depth and has a precise data base to compare the results of the simulation with. The lift coefficients are computed with accuracy in contrast to the drag coefficient. The overprediction of drag is expected since the actual airfoil has laminar flow over the forward half. The turbulence models cannot calculate the transition point from laminar to turbulent and consider that the boundary layer is turbulent throughout its length. Therefore, it is necessary to compare the computational results with experimental data of a fully turbulent boundary layer. In order to get more accurate results, the computational domain could be split into two different domains to run mixed laminar and turbulent flow. The contours of pressure and velocity over the airfoil are presented, as well as stagnation, maximum velocity, detachment and reattachment points of the boundary layer are computed. Streamlines and velocity vectors over the airfoil are also presented. Similar simulations are conducted for the S809 airfoil. In order to study the tree-dimensional effects of the flow, simulations over the S809 wing are made. Lift and drag coefficients, stagnation, maximum velocity, detachment and reattachment points of the boundary layer are computed. Moreover, contours of turbulent intensity on the upper surface of the wing and velocity, z-velocity, turbulence intensity and helicity behind the wing are presented. Simulations over the Phase IV blade of NREL are also conducted. The axial velocity behind the rotor, the static pressure and the turbulence intensity contribution on the blade’s surface and the static pressure contours at several blade cross-sections are studied. First of all, the computational study of the two-phase flow over a NACA 0012 airfoil and Liquid Water Content LWC=30 g/m3 is conducted, because there are published experimental data for comparison, in order to validate the CFD developed model. After that, simulations of two-phase flow over the S809 airfoil, S809 wing and Phase IV blade are made. In addition, computational study of the effects of different Liquid Water Content on the aerodynamic performance of NACA 0012 and S809 airfoil at low Reynolds number is made. The results from two-phase flow are compared with the corresponding results from one-phase flow in order to show the effects of two-phase flow at the aerodynamic performance. The influence of two-phase flow on the power coefficient of a wind turbine is also investigated. The results show that the aerodynamic performance degrades when encountering rain, especially lift is degreased and drag is increased. The aerodynamic degradation is caused by the water film formation on the airfoil’s surface and the cratering effects from the raindrops impact. The presence of uneven water film on the airfoil surface roughens the airfoil surface and increases the airfoil thickness. The cratering effects from the water droplets impact on the water film layer increase also the airfoil thickness. Moreover, the droplets splash-back when they impact the airfoil and as a result droplets with smaller diameter and velocity are formed. The acceleration of the splashed-back droplets by the air flowfield acts as a momentum sink, deenergizing the boundary layer and leaving it more susceptible to separation. The aim of the study of the aerodynamic behavior of blades in two-phase flow is the construction of wind turbines with greater efficiency and the production of energy from wind with low cost.

Με τη διατριβή έγινε ανασκόπηση των πιο γνωστών θεωριών που διέπουν τις πολυφασικές ροές. Ακολούθησε συλλογή και κριτική αξιολόγηση των θεωρητικών μοντέλων. Εκπονήθηκαν δύο προγράμματα σε Η/Υ, τόσο για την ομογενή ροή όσο και τη χωριστή ροή. Έγινε σύγκριση των θεωρητικών με υπάρχοντα πειραματικά αποτελέσματα. Σχεδιάστηκε και κατασκευάστηκε πειραματική διάταξη. Λήφθηκαν πειραματικές μετρήσεις και έγινε αξιολόγησή τους και σύγκριση με θεωρητικά αποτελέσματα. Προτάθηκε μέθοδος διόρθωσης της πρόβλεψης λειτουργίας μιας υδροπνευματικής αντλίας με τη βοήθεια του μοντέλου ‘drift-flux’ (μέθοδος CoSM). Προτάθηκαν νέοι ροϊκοί χάρτες κατάλληλοι για την παρουσίαση της λειτουργίας των υδροπνευματικών αντλιών, τον υπολογισμό του κλάσματος κενού και την μετάβαση των ροϊκών καταστάσεων. Προτάθηκε τρόπος προσδιορισμού της μετάβασης ‘slug-churn’ με απλή φωτογραφική μέθοδο (camera) και χρήση του μοντέλου ‘drift-flux’. Τα αποτελέσματα της διατριβής (πειράματα, μέθοδος CoSM, ροϊκοί χάρτες και μετάβαση slug-churn) παρουσιάστηκαν σε συνέδρια και επιστημονικά περιοδικά. / This PhD thesis deals with multiphase flows and air-lift pumps. All well-known theories concerning these two scientific fields are presented and analyzed in detail. Two computational codes were developed for homogeneous two-phase and separated three-phase flow. A comparison between theoretical results and experimental data followed. An experimental investigation was performed in a lab scale air-lift pump installation at Fluid Mechanics Laboratory, University of Patras. A new method for the precise prediction of the performance of a two-phase air-lift pump with the aid of drift-flux model was presented (CoSM method). New regime maps were introduced suitable for air-lift pump presentations. That means the direct view of the flow behaviour inside the air-lift pump, void fraction calculation and the regime transitions. An experimental method was presented for the prediction of slug-churn transition in two-phase flow, using a camera and ‘drift-flux’ model. The results of this work were presented in International Conferences and Journals.