Pintakitkan mittauslaitteen virtauskentän simulointi

Työn tarkoitus on yrittää selvittää CFD-simuloinnin avulla APSO-nimisen pintakitkamittarin kokeissa ilmenneitä ristiriitaisuuksia. Neljä valittua tuulitunnelitapausta simuloitiin Reynolds-keskiarvotetuilla Navier-Stokes-yhtälöillä käyttäen Menterin k - omega BSL-turbulenssimallia. Ajasta riippumaton symmetrinen ratkaisu iteroitiin FINFLO-virtausratkaisijalla. Geometriaksi valittiin APSO:sta yksinkertaistettu sylinteri tasolevyn pinnalla. Laskentahilat sisältävät rajakerroksen sekä 0,4 m x 0,4 m alueen tasolevystä pintaesteen ympäriltä. Hiloissa on noin kaksi miljoonaa koppia kummassakin ja pienimmät kopit ovat kooltaan 4 µm x 4 µm x 123 µm. Ennalta lasketut rajakerrosvirtaukset on asetettu ylävirran puolelle sisäänvirtausreunaehdoksi, joka kiinnittää nopeusprofiilin lisäksi turbulenssiasteen ja pyörreviskositeetin. Kaksi ylimääräistä tapausta simuloitiin käyttäen mitattua nopeusprofiilia reunaehtona. Diplomityön dokumentti alkaa valikoivalla virtausmekaniikan kertauksella, johon on valittu perusvirtaustapauksia ja rajakerrosteoriaa taustaksi simulointituloksia varten. Tämän jälkeen on esitelty valikoima pintakitkan mittaustapoja, minkä jälkeen APSO ja siihen liittyvät tuulitunnelitulokset esitellään. Tämän työn laskennallinen simulointi käydään läpi ja sen tulokset esitellään. Mukana on visualisointia virtauskentästä sekä paine- ja kitkajakaumia. Simuloitujen tapausten kalibrointiyhteydet kitkan ja paineen välille esitetään yhdessä tuulitunnelitulosten kalibrointikäyrien kanssa. Simuloidut tapaukset edustavat kahta eri rajakerrosta ja kahta eri pintaesteen korkeutta. kahdessa tapauksessa esteen korkeus suhteessa rajakerrokseen on sama. Kun simuloidut paine-erot ja pintakitkat skaalataan kalibrointimuuttujiksi ja piirretään kalibrointikäyriksi, huomataan, että ne asettuvat melko hyvin samalle suoralle logaritmiasteikolla ja mitään merkittävää siirtymää tai ryhmittymistä ei havaita perustuen esteen suhteelliseen korkeuteen. Simulointitulokset asettuvat kahden tuulitunnelisarjan välille, joten ne eivät erityisesti tue kummankaan mittaussarjan oikeellisuutta. Tapauksissa, joissa käytettiin mitattua nopeusprofiilia, tulokset osuvat lähelle vastaavia puhtaasti laskennallisia tapauksia. Kaikki tapaukset edustavat puristumatonta virtausta laajalta kalibrointialueelta. Kolmiulotteinen tiheän hilan malli erikoisella sisäänvirtausreunaehtojärjestelyllä osoittautui melko raskaaksi geometria- ja virtausparametrien variointiin. Kaksiulotteinen malli olisi suositeltavampi kattavaan parametrianalyysiin

Modelling load and vibrations due to iced turbine operation

Wind energy in icing and low-temperature climate has a huge growth potential, but rotor icing effects on turbine dynamics and lifetime are not well known and simulations with iced rotor are not required in current IEC 61400-1 turbine design standard. In this article, simulations with iced rotor are compared to measured mechanical loads. The dynamic behaviour of the wind turbine was simulated with FLEX5 aeroelastic code for Senvion MM92 2 MW wind turbine. Simulations with typical iced airfoil lift and drag coefficients, aerodynamic and mass imbalances for iced rotor were performed and compared to measured iced turbine loads. Resulting iced turbine simulation parameters can be used in defining new design load cases for cold climate turbines. The most representative simulation parameter combination was achieved with a symmetric aerodynamic penalty applied on all blades and an asymmetric rotor mass imbalance of 166 kg ice load on two blades and 83 kg ice load on one blade