Performance Of Wells Turbine Engineering Essay

Performance Of swell Turbine Engineering EssayABSTRACTA rise up turbine has inherent disadvantages in comparison with conventional turbines relative low faculty and poor starting characteristics. In this case, the performance of wells turbine is studied on computational analysis by changing aerofoils and providing antithetic be given of relative incidence for the improvement of the turbines performance.Study is establish on analysing the lean of air on turbine victimization computational analysis at steady condition.1 INTRODUCTIONThe Ozone depletion and global warming have altered the international community and urged the need for more focus on alternative unripe showtimes of energy. Ocean wave energy is one of the renewable forms of energy which back end be utilized in response to the disturbing prospect of an exhaustible source of energy. Several wave energy devices being studied under many wave energy programs make employ of the principle of the oscillating water mains tay (OWC).Potentially the most successful device used in harnessing on wave energy has been the OWC wave energy converter. The OWC chamber, either floating or bottom standing, with the immersed devastation opened to the action of the sea. A reciprocating air precipitate is created by the action of the free surface of the water within the chamber. The conversion of this airflow into mechanical energy may be achieved by a egress of devices like-A. TAPCHAN The TAPCHAN comprises a gradually narrowing channel with wall heights typically 3 to 5 m above mean water take. The waves enter the wide end of the channel and, as they propagate down the narrowing channel, the wave height is amplified until the wave crests expire over the walls to a reservoir which provides a s shelve water supply to a conventional low head turbine. The requirements of low tidal range and come acrossting shoreline limit the world-wide installation of this device. public figure (a) TAPERED CHANNEL1 (TAPCHAN)h ttp//re.emsd.gov.hk/english/other/marine/images/marine_tech_010_2.gifB. PENDULOR The PENDULOR device consists of a rectangular box, which is open to the sea at one end. A pendulum flap is hinged over this opening, so that the action of the waves pay offs it to swing back and forth. This motion is then used to force out a hydraulic pump and generator.http//t2.gstatic.com/images?q=tbnANd9GcQ7yge9ouptnhszDgsXGA_gCvAXKqbo78BeXZHFFtPB89433p0pFig (b) PEDULOR 2C. come up TURBINE The come up turbine is one of the most suitable air turbines for energy conversion from oscillating air flow. . A schematic view of the OWC device with a Wells turbine is shown in Fig. c. The Wells turbine is an axial flow air turbine. It consists of several symmetrical aerofoil blades set around a hub. As waves Impinge on the device, they cause the water column to rise and fall in the air chamber, which alternately compresses and depressurized the trapped air. This air is allowed to flow to and from the atmos phere by a turbine which drives an electric generator.http//www.aussiestockforums.com/forums/attachment.php?attachmentid=9213stc=1d=1180172232Fig (c) WELLS TURBINE3Fig. 1.1 ceremonious of the Three Major shoreline DevicesThe Wells turbine is a self-rectifying air turbine which is expected to be widely used in wave energy devices with the OWC (Raghunathan, 1995) .It can extract power at low airflow rate, when other turbines would be inefficient. The Wells turbines for wave power conversion have less aptitude. To increase the expertness is the major quest all over the world, the flow of air through the wells turbine impeller is carried out in this reckon by using distinct size of it impeller and introducing biplane i.e. two rows of symmetrical aerofoil blades.1.1 WAVE RESOURCEShttp//www.oceanenergy.ie/images/world-map.jpgFig. 1.2 Global Distri entirelyion of Deep Water Wave power levels in kW/m crest length4 scorn the climate change phenomena, the world imagery for wave remain s very much as shown in fig. 1.2 by Dr Tom Thorpe 5. The full(prenominal)est energy waves ar strong off the western coasts in the 40o-60o latitude range north and south. The power in the wave fronts varies in these atomic number 18as between 30 and 70kW/m with peaks to 100kW/m in the Atlantic SW of Ireland, the Southern Ocean and off Cape Horn. The capability to supply electricity from this resource is such that, if harnessed appropriately, 10% of the current level of world supply could be provided 4PERIODAMPLITUDEPOWER DENSITYVELOCITY (m/s)WAVELENGTH( sulfur)(m)(kW/m)(m)Storm1414170023320Average93.56015150Calm5.50.51950Fig 1.3 Tthe nautral and technical wave energy resource for the north and west side of the UK6The techinical resource is dependent on the nautral conditions like the shape of the rock and location i.e beaches and gullies. The wave energy at calm sea is considered in this project.1.2 WELLS TURBINEThe monoplane Wells turbine i.e. the basal Wells turbine consists o f several symmetrical aerofoil blades (NACA four digit series) set around a hub at 90 details with respect to the airflow. Since its an invention by Prof. A.A. Wells in 1976, most researchers have focused on improving its efficiency and its range of efficient operation. In fact, compared to other conventional air turbines (e.g. Francis turbine) the Wells turbine has a lower efficiency and a narrow operational region. Nevertheless, it can extract power at low airflow rate, when other turbines would be inefficient.Fig. 1.4 Schematic of the Monoplane Wells Turbine7A schematic diagram of a Wells turbine is shown in Fig. 1.4. At first sight the arrangement might seem to be unlikely elbow room of energy conversion. However, once the blades have attained design speed the turbine is capable of producing a time-averaged positive power output from the cyclically reversing airflow with a fairly high efficiency. Wells turbine has low efficiency and poor starting characteristics.The Biplane We lls turbineMuhammad Mamun in the Study on the Hysteretic Characteristics of the Wells Turbine in a Deep Stall specification says the pressure drop across a mono-plane Wells turbine above is proportional to the straightforwardly of the tip speed which has to be limited if transonic effects are to avoided. For wave energy devices which call forth significantly larger pressure drops than the limit for a single plane turbine a biplane turbine can be usedFig. 1.5 Schematic of the Biplane Wells 7It has certain advantages over the conventional monoplane Wells turbine as followsI. It can operate under high loading.II. It can absorb higher wave power than the monoplane turbine if the diameter and rotational speed of the turbine are kept constant.III. The design speed is lower than that of the monoplane for the same loading.IV. It avoids the use of guide vanes and therefore the turbine would require less support and repairs. 71.3 PRINCIPLE OF OPERATIONThe principle of operation of Wells t urbine is based on the classical aerofoil theory. According to the classical aerofoil theory, an aerofoil which is set at an tap of incidence in a fluid flow generates a lift force L normal to the free stream. The aerofoil also experiences a drag force D in the statement of the free stream (relative velocity). These lift and drag forces can be resolved into tangential (in the plane of rotation) and axial (normal to the plane of rotation) components FT and FA respectively.Fig. 1.6 Notation for determining lift, drag, and axial and tangential forces onAn aerofoil7Resulting expression for axial and tangential forcesFA = Lcos + DsinFT = Lsin DcosThe axial force is absorbed but the turbine while the tangential force causes the turbine to revolve. For a symmetrical airfoil the direction of tangential force is the same for both positive and negative values of . Therefore, the direction of rotation of the rotor coil is independent of airflow direction.2 AIMS AND OBJECTIVES OF THE PROJE CTSimulation of air flow through wells turbine impeller by means of numerical rule using a CFD (Computational fluid dynamics) called unstable and check the flow process of different parameters and the factors affecting the differences. Since wells turbine is a low efficiency turbine to increase the efficiency of Wells turbine is the other aim. Path followed to meet the requirements is first calculation of efficiency theoretically. Simulating a modified design by using different angle of incidence and making biplane i.e. two rows. Comparing the results of different model and selecting the suitable design.3 LITERATURE REVIEW3.1 Types of CFD PROCESS USEDCommercial CFD code FLUENT, Star-CD, FLOW-3D, CFX/AEA, and so on interrogation CFD code Self-developedPublic domain software (PHI3D, HYDRO, and WinpipeD, etc.)Other CFD software includes the Grid extension software (e.g. Gridgen, Gambit) and flow visualization software (e.g. Tecplot, FieldView)Commercial CFD code FLUENT is used in t his project.3.2 General working on CFDTable3.1 CFD working layoutThe Processes shown in the table 3.1 is divide into pre- process and post- process viz. GAMBIET AND FLUENTGeneral grade of GAMBIT operationsInitial setupSolver selection, profit size, Defaults, etc.Geometry creation (ACIS, IGES or Mesh import)Create full geometryDecompose into take-able sectionsMeshing local anaesthetic web Edge and Boundary layersGlobal meshing Face and/or VolumeMesh examinationZone assignmentContinuum and Boundary attachmentMesh exportGeneral sequence of FLUENT operationsSelection of appropriate models.Turbulence, combustion, multiphase, etc.Define material propertiesFluidSolidMixturePrescribe operating(a) conditionsPrescribe boundary conditions at all boundaru zonesProvide and initial solutionSet up solver controlsSet up convergence monitors3.3 Grid generationGrid generation is one of the key elements in Computational Fluid Dynamics (CFC). It has now become a fairly common tool for use in the n umerical solution of partial differential equations on arbitrarily shaped regions. The numerical solution of partial differential equations requires some discretization of the field into a collection of points (nodes) or elemental volumes (cells). The differential equations are approximated by a set of algebraic equations on this collection, and this system of algebraic equations is thensolved to produce a set of discrete values which approximates the solution of the partial differential system over the field. The practice of discretizing the physical domain into a finite number of elements is called as grid generation.3.4 Grid topologiesGenerally, the governing equations may be transformed into finite element, finite difference, or finite volume equations. The cell types supported by FLUENT are followed as triangular and quadrilateral cells in 2D are accepted, and in 3D, tetrahedral, hexahedral, wedge, and pyramid cells can be usedFIG3.2 Different types of grids organise versus Un incorporate GridsThe section presents a brief description of grid generation. The grid generation techniques available at present fall into two categories, namely a) organize grid generation and b) unstructured grid generation. The structured grid generation techniques are based on the transformation of the complex physical domain into a simple computational domain, which is often chosen to be rectangular in shape (quadrilateral and hexahedron).The unstructured grid generations have been used with FEM (finite element method) procedure only, whereas structured grids have general applicability.7.3.5 Types of structured gridIn FLUENT, both single-block and multi-block structured meshes are acceptable, as well as hybrid meshes containing quadrilateral and triangular or hexahedral, tetrahedral, pyramid, and wedge cellsMultiple BlockSometimes, it is possible to combine several structured computational meshes together to fit the physical domain. Multi- locking has the advantage of the spe ed of a structured solver, without as many mapping constraints apparent in single block meshes.Single BlockIn this technique, one computational grid is mapped to fit the whole physical domain. For even moderately complex shapes, it may be practically impossible to define a transformation which will map the outer(prenominal) surface of the computational domain to the required physical shape, while ensuring that the resulting grid has desirable attributes of smoothness.73.6 Mesh qualityThe quality of mesh plays a significant section in the accuracy and stability of the numerical simulations. The attributes associated with mesh quality are density of node, cell shape, smoothness and flow-field dependency. In many cases, poor resolution in critical regions can dramatically alter the flow characteristics.3.7 The Capabilities of FLUENTThis section provides a brief introduction to FLUENT and an explanation of its capabilities 10.FLUENT used in this project is a commercial code and a stat e-of-the-art computer program for manikin single and multiphase flows, agitate and mass transfer, chemical reaction phenomena, and etc. in complex geometries. This code includes following components FLUENT, the flow solver GAMBIT, the pre-processor for geometry modelling and mesh generation pre-PD, and etc. FLUENT solver utilizes a finite-volume, pressure-based, multiphase infinite marching method (SIMPLE algorithm), for solving the governing integral equations for conservation of mass and momentum, and for energy and other scalars such as turbulence and chemical species. It has the following modelling capabilities Flows in 2D or 3D geometries using triangular/tetrahedral, quadrilateral/hexahedral, or mixed (hybrid) grids that include prisms (wedges) or pyramids In compressible or compressible flows Steady-state or transient analysis Laminar and turbulent flows Newtonian or non-Newtonian flow Convective heat transfer, including natural or forced convection Coupled conduction/con vective heat transfer Radiation heat transfer Inertial (stationary) or non-inertial (rotating) reference frame models Multiple moving reference frames, including sliding mesh interfaces and mixing planes for rotor/stator interaction modelling Chemical species mixing and reaction, including combustion sub models and surface deposition reaction models Arbitrary volumetric sources of heat, mass, momentum, turbulence, and chemical species Flow through porous media One- symmetryal fan/heat-exchanger performance models Two-phase flows, including cavitations Free-surface flows with complex surface shapesFLUENT can provides a number of boundary conditions, including Velocity or Pressure set Inlets/ subjects Stationary or Moving Walls, with or without Friction Periodic Boundary Conditions Symmetry Boundary Conditions Pressure Far-filed Boundary Conditions Outflow Boundary Conditions Inlet/Outlet Vent Boundary Conditions Intake/Exhaust Fan Boundary ConditionsAs the Well turbine has a complex geometry for modelling, a large number of modelling capabilities are required of the CFD code for the turbine. FLUENT can incorporates all of these capabilities, and is most suitable for modelling the Wells turbine.104 ANALYSIS OF TASK4.1 Theoretical calculationThe dimension used in this project is of prototype obtained from others experimental work, the model is designed and simulated by using the two different models shown in the table below.a 8 b9Table 4.1 Dimension of wells turbineThe theoretical calculation of efficiency is done using the above two different dimension, the method used to calculate the efficiency is shown below.CALCULATION FOR EFFICIENCYSimilarly,CALCULATION AT 4 peak ANGLE OF ATTACKAt = 4 degreeThe table below shows a work out efficiency at different angle of onrush calculated using the format shown above.(degree)(radians)W(relative velocity)(rads/sec)ReClCd4.000.07143.3642.451221641.450.400.0121.865.000.09114.7442.39977760.010.500.0123.816.000.1095.6742.32 815255.460.600.0125.187.000.1282.0642.24699252.150.700.0126.028.000.1471.8542.14612312.290.800.0126.449.000.1663.9242.03544748.240.900.0126.7310.000.1757.5941.91490747.451.000.0126.7811.000.1952.4141.77446611.021.100.0126.9612.000.2148.1041.62409873.061.200.0226.9713.000.2344.4541.46378826.411.300.0226.8914.000.2441.3441.29352251.701.400.0226.9015.000.2638.6441.10329254.751.300.0223.38Table 4.2 Efficiency at different angleUsing the values of efficiency and the angle of glide path from the above table (4.2) a direct relation between efficiency and the angle of attack is obtained which can be seen in the graph below (fig 4.3). Usig a Trendline option in Microsoft Excel an equation of direct relation between angle of attack and efficiency is obtained. The equation shown in the graph is a sixth order equation which is difficult to differentiate to obtained the angle at which the efficiency will be maximal,so a 2nd order equation is obtained from trendline option. Differentiating the equation gives the value of an angle at which the efficiency is max. From this procedure 12 degree is the calculated angle obtained at which the efficiency is max.Fig 4.3 Efficiency Vs. Angle of Attacky = -0.00026 + 0.00785 0.10484 + 0.70883 2.70462 + 6.5617x + 17.369when x = 12y = = 25.98 %After substituting the value on angle obtained for maximum efficiency a difference between the two values is found and it is due to the R squared value. More closer the value of R square to unity more accurate results can be obtained.Equations obtained from Microsoft Excel at different orders are shown below- distinguish 2y = -0.12842 + 1.8889x + 20.336R = 0.8797Order 3y = -0.00383 0.05482 + 1.4909x + 20.851R = 0.8848Order 4y = -0.00524 + 0.13263 1.22762 + 5.2125x + 17.578R = 0.9636Order 5y = -0.0015 + 0.02864 0.26813 + 0.86162 + 0.6489x + 20.649R = 0.9869Order 6y = -0.00026 + 0.00785 0.10484 + 0.70883 2.70462 + 6.5617x + 17.369R = 0.9945Similarly using the dimension in table 4.1 (b) the calculated efficiency is show below(degree)(radians)W(relative velocity)(rads/sec)ReClCd4.000.07143.3678.24916231.000.400.0153.725.000.09114.7478.13733319.930.500.0158.496.000.1095.6778.00611441.530.600.0161.877.000.1282.0677.85524439.060.700.0163.938.000.1471.8577.67459234.170.800.0164.969.000.1663.9277.47408561.140.900.0165.6810.000.1757.5977.24368060.551.000.0165.8011.000.1952.4176.99334958.231.100.0166.2412.000.2148.1076.72307404.761.200.0266.2813.000.2344.4576.42284119.781.300.0266.0714.000.2441.3476.10264188.751.400.0266.1015.000.2638.6475.76246941.041.300.0257.46TABLE 4.4 Efficiency at different anglesSimilarly in this case a graphical representation of Angle of Attack Vs. Efficiency is obtained which can be seen below and the equation represents a direct relation between efficiency and angle of attack.Fig 4.5 Efficiency Vs. Angle of AttackOrder 6y = -0.00066 + 0.02925 0.62024 + 6.85633 42.0712 + 139.39x 138.43R = 0.9945when x =12y = = 64.24%Similarly using the order 2 equation to find the angle at which the efficiency will be maximum. The calculate angle using the same procedure as above is 12 degree at which the efficiency is maximum.. 4.2 Gambiet (Pre Processing)-The figure below shows an impeller of wells turbine designed with blades at 0 degree angle of incidence and using the dimension from the table 4.1 (a).Fig 4.6 Impeller of wells turbineCreating a model using gambiet and then meshing the geometry for which meshing size is selected based on the Reynolds number. Since the Reynolds number lies in the transational flow at the angle in which the efficiency is maximum,using turbulence boundary layer formula =0.00269The thicknes of boundary layer is 0.003 m. The mesh size comes to be 0.001m to get three elements in one layer to get fine meshing. In case of 3-Dimensional model the mesh elemet used is Tet/Hybrid. Checking the meshing quality the Aspect Ratio lies between 1 to 4. Boundary conditions is given for impeller is moving wall and interfa ces is decided so that the fluid can be go around within this volume. The mesh is exported for post processing in Fluent4.3 Fluent (Post Processing)Steps used in fluent is as followsStep 1 source the case fileStep 2 Defining the grid interfacesStep 3 Grid checkStep 4 Defining model as viscous and using K-epsilon (2 equation )Step 5 Defining boundary conditionIn boundary condition fluid within the impeller is made to rotate at 40 rads/sec.The impeller is a moving wall rotating relative to cell zone at 0 rads/sec.Inlet velocity is 10 m/sec and the turbulence method selected is intensity and hydraulic cylinder.Step 6 Solution is converged after ilteraion5 RESULTS AND DISCUSSIONThe results shown below contains pressure contours, velocity vectors and pathlines at different cros-section of the models designed using the dimension from table 4.1 (a).Model with blades at 0 degree angle of incidence and time out flow from topModel with blades at 0 and 2 degree(+) angle of incidence a nd portal flow from topModel with blades at 0 and 2 degree(+) angle of incidence and inlet flow from bottomModel with blades at 2 degree(+) angle of incidence and inlet flow from topModel with blades at 2 degree(+) angle of incidence and inlet flow from bottomBiplane modelsModel with blades at 0 degree angle of incidece and inlet from topModel with blades at 0 degree angle of incidence and inlet flow from bottomModel with blades at 2(+)and 2(-) degree angle of incidence and inlet flow from topModel with blades at 2(+)and 2(-) degree angle of incidence and inlet flow from bottomModel with blades at 0 degree angle of incidence and inlet flow from topModel with blades at 0 and 2 degree(+) angle of incidence and inlet flow from topModel with blades at 0 and 2 degree(+) angle of incidence and inlet flow from bottomModel with blades at 2 degree(+) angle of incidence and inlet flow from topModel with blades at 2 degree(+) angle of incidence and inlet flow from bottomBiplane modelsModel wi th blades at 0 degree angle of incidece and inlet from topModel with blades at 0 degree angle of incidence and inlet flow from bottomModel with blades at 2(+)and 2(-) degree angle of incidence and inlet flow from topModel with blades at 2(+)and 2(-) degree angle of incidence and inlet flow from bottomComparing the above graphical results under a range of 0-400 for comparison except the last two model. The table below shows the value of dynamic pressure (max) in Pascals of above design.From the table it can be seen that the introduction of two rows provides a better result in terms of dynamic pressure. After giving the installatoin angle the maximun dynamic pressure obtained is 1176 pascals by which we can say that the two rows impeller with and an installaition angle is better than the single rows .Assumptions non-homogeneous assumptions made to carry out the simulation is as followsPATHLINES OF PARTICLES ON IMPELLERAOA 0 recess FROM assoilAOA 0 AND 2 DEGREE INLET AT TOP AOA 0 AND 2 DEGREE DEGREE INLET AT BOTTOMAOA 2 DEGREE INLET AT TOP AOA 2 DEGREE INLET AT BOTTOMTWO ROWSAOA 0 DEGREE INLET AT TOP AOA 0 DEGREE INLET AT BOTTOMAOA +2 -2 DEGREE INLET AT TOP AOA +2 -2 DEGREE INLET AT BOTTOMThe results shown below is of the dimension used from table 4.1 (b).modelling of the wells turbine is divided into two parts theoretical and practical