Multidisciplinär modellering av vattenkolvsoscillationer i vågenergiomvandlare : Bedömning av flödesmotstånd genom CFD-modellering med vätskestrukturinteraktioner

As the world’s need for electricity increases, so does its demand for sustainable energy with low to no greenhouse gas emissions. One of these renewable sources of energy is the Ocean which is one of the world’s largest and most predictable energy source, where the extraction of energy can be from waves or tidal current, with zero greenhouse gas emissions during production. A company which works with wave energy is Waves4Power which has developed the wave energy buoy WaveEL.WaveEL is comprised of a buoy which is eight meters in diameter with a 36 meters long vertical cylinder which goes through the buoy. In the cylinder is a piston that oscillates in pace with the waves and generates electricity. Between the piston and the cylinder wall is a gap where the water can move from one side of the piston to the other in pace with the piston’s oscillations. The gap is called the leakage clearance. The leakage clearance effect on the flow resistance is the focus of this master thesis as something which has not been studied before in scientific articles for similar wave energy buoys or in other fields.The aim of the master thesis is to improve the understanding of how the water flow, because of the leakage clearance in the WaveEL buoy, affects the force which the piston is subjected to, and in turn how much electricity can be generated. As it is a complex system the focus will be to determine the dynamic flow resistance parameter because of the leakage clearance and the changes to the dynamic flow resistance parameter as the dimensions of the piston is varied.The leakage clearance effect on the flow resistance has been studied with the help of Computational Fluid Dynamics (CFD) in the software COMSOL Multiphysics 6.0 in two different models. In the first model, model A, the piston is locked in different positions in the cylinder and the pressure at the bottom of the cylinder varies to reflect the motion of the waves. For the second model, model B, the piston is allowed to move vertically in the cylinder due to a set force which reflects the motion of the waves.The dynamic flow resistance parameter for model A is lower at higher Reynolds number and within an interval between 0.4 and 1.6 within the working area. Outside of the working area the dynamic flow resistance parameter is lower at a higher Reynolds number and higher than in the working area at an interval between 0.7 and 45.For model B, the dynamic flow resistance parameter has only been calculated for the working area and is within an interval of 0.1 and 7, the value for the dynamic flow resistance parameter is low when the Reynolds number is high. Dissimilar to model A where the piston is locked into position, the piston oscillates in model B. There is a phase shift between the velocity of the piston and the velocity of the water, which leads to the piston being subjected to a larger force than in model A at lower water velocities. This is one of the reasons why the dynamic flow resistance parameter is higher in model B than model A at low Reynolds number. As model B calculates the dynamic flow resistance parameter based on the relative velocity between the piston and the water, the dynamic flow resistance parameter becomes lower than for model A at higher Reynolds number.For the performed sensitivity analysis, the results shows that a more advantageous value on the dynamic flow resistance parameter can be achieved by altering the dimensions of the piston. A more advantageous result was achieved for example when the rounded edge on the piston became sharper or when the leakage clearance width was increased by 10%. If this master thesis work is to be extended, the studies should focus on elaborating model B either more in depth or with values derived from experiments from the WaveEL buoy for a more realistic model and thus determine a more accurate dynamic flow resistance parameter.The results from the sensitivity analysis justifies a future study where the dynamic flow resistance parameter should be investigated with greater variations of the piston diameters, as this can increase the flow resistance and thus generate more electricity. However, it should be investigated in relation to the cost of manufacture, to obtain the ultimate design which generates maximum electricity for a reasonable manufacturing cost.Allt eftersom värdens elektricitetsbehov ökar, ökar också efterfrågan på förnybar energi med låga eller inga växthusgasutsläpp. En av dessa förnybara energikällor är havet som anses vara en av världens största och mest förutsägbara energikällor, där utvinningen av energi kan ske från vågor eller tidvattenströmmar, med noll växthusgasutsläpp under produktion. Ett av företagen som arbetar med vågenergi är Waves4Power som har utvecklat vågenergibojen WaveEL.WaveEL består av en boj som är åtta meter i diameter med en 36 meter lång vertikal cylinder som går igenom bojen. I cylindern finns en kolv som oscillerar i takt med vågorna och genererar elektricitet. Mellan kolven och cylinderväggen finns en spalt där vatten kan förflytta sig från ena sidan av kolven till den andra i takt med att kolven oscillerar, denna spalt kallas läckspalten. Läckspaltens påverkan på flödesmotståndet är någonting som tidigare inte har studerats i andra vetenskapliga artiklar för likartade vågenergibojar eller inom andra branscher och är därför fokus för detta examensarbete.Syftet med examensarbetet är att förbättra förståelsen hur vattenflödet till följd av läckspalten i WaveEL-bojen påverkar kraften som kolven utsätts för, vilket i sin tur påverkar hur mycket elektricitet som genereras. Eftersom det är ett komplext system kommer focus vara att fastställa den dynamiska flödesmotstånds parametern till följd av läckspalten och förändringarna i den dynamiska flödesmotstånds parametern när dimensionerna på kolven varieras.Läckspaltens påverkan på flödesmotståndet har studerats med hjälp av Computational Fluid Dynamics (CFD) med programvaran COMSOL Multiphysics 6.0 i två olika modeller. I den första modellen, modell A, är kolven låst i olika positioner i cylindern och trycket vid botten av cylindern varierar för att efterlikna vågrörelser. För den andra modellen, modell B, tillåts kolven röra sig vertikalt i cylindern genom en bestämd kraft som efterliknar vågrörelserna.Resultatet från modell A visar att flödesmotståndet befinner sig mellan 0,4 och 1,6 inom arbetsområdet och utanför arbetsområdet är den dynamiska flödesmotståndsparametern större än i arbetsområdet, här är den dynamiska flödesmotstånds parametern mellan 0,7 och 45. Oberoende av kolvens position i bojen sjunker den dynamiska flödesmotståndsparametern i takt med att Reynolds talet ökar.För modell B har den dynamiska flödesmotstånds parametern endast beräknas för arbetsområdet och ligger inom intervallet 0,1 och 7, där ett lågt värde på den dynamiska flödesmotstånds parametern fås när Reynolds talet är högt. Tillskillnad mot modell A där kolven är låst i sin position, oscillerar kolven i modell B. Det finns en fasvridning mellan kolvens och vattnets hastighet, vilket innebär att kolven utsätts för en större kraft i modell B än i modell A vid låga vattenhastigheter. Det är en av anledningarna till att den dynamiska flödesmotstånds parametern är högre i modell B vid låga Reynolds tal. Eftersom modell B beräknar ut den dynamiska flödesmotstånds parametern med hjälp av den relativa hastigheten mellan kolven och vattnet blir den dynamiska flödesmotståndsparametern lägre än för modell A högre Reynolds tal.För den genomförda känslighetsanalysen visar resultatet att ett mer fördelaktigt värde på den dynamiska flödesmotsstånds parametern kan uppnås vid förändring av kolvens dimensioner. Exempelvis uppnåddes ett mer fördelaktigt resultat när den avrundade kanten på kolven blev skarpare eller när bredden på läckspalten ökades med 10%. Om detta examensarbete ska utökas bör studierna fokusera på att utveckla modell B antingen mer på djupet eller med värden härledda från experiment från WaveEL-bojen för en mer realistisk modell och därmed en mer exakt bestämning av den dynamiska flödesmotstånds parametern.Resultatet från känslighetsanalysen motiverar en framtida studie där den dynamiska flödesmotstånds parametern bör undersökas med kraftigare dimensions variationer på kolven, eftersom detta kan komma att öka flödesmotståndet och således generera mer elektricitet. Dock behöver de undersökas i relation med tillverkningskostnaderna, för att hitta den ultimata design som genererar maximalt med elektricitet vid en rimlig tillverkningskostnad