Research Article  Open Access
Tomomi Uchiyama, Haruki Fukuhara, Shouichiro Iio, Toshihiko Ikeda, "Numerical Simulation of Water Flow through a NanoHydraulic Turbine of WaterfallType by Particle Method", International Journal of Rotating Machinery, vol. 2013, Article ID 473842, 8 pages, 2013. https://doi.org/10.1155/2013/473842
Numerical Simulation of Water Flow through a NanoHydraulic Turbine of WaterfallType by Particle Method
Abstract
This study simulates the flow through an impulsetype smallscale hydraulic turbine utilizing a waterfall of extralow head. The twodimensional Moving Particle Semiimplicit (MPS) method is employed for the simulation. The fluid is discretized by particles, and the flow is computed by the Lagrangian calculation for the particle motion. When the distance between the particles discretizing the waterfall of a width , , is set at , the flow can be simulated with the sufficiently high spatial resolution, and the rotor performance can also be favorably predicted. The present simulation also successfully analyzes the effect of the rotational speed of rotor on the flow and the turbine performance.
1. Introduction
Hydropower is one of the promising renewable energy resources. It is converted to electric energy through hydraulic turbines. In Japan in 2010, approximately 10% of all supplying electric energy is provided by hydraulic power generations. The percentage is expected to increase steadily on the basis of a government policy promoting the development of renewable energy. As largescale hydroelectric plants require huge dams and long conduits, the places for the construction are hardly remained. Thus, expectations for the development of a smallscale hydropower, of which output is less than 1000 kW, have been increasing. When natural disaster occurs, the largescale centralized hydraulic power plants may lose the power supplying ability due to the collapse of the power grid. Since the smallscale hydropower, existing in smallscale rivers, irrigation canals, and industrial drainages, distributes widely in Japan, it realizes the smallscale distributed power generation. The smallscale hydropower can contribute the local production for local consumption of electric power, which is more resistant to disaster. Consequently, the development is also of great worth from the viewpoint of constructing a disasterresistant society.
To exploit effectively the smallscale hydropower, various types of hydraulic turbine have been presented [1–5]. Ikeda et al. [6] developed an impulsetype hydraulic turbine utilizing waterfalls of extralow head which is 2 m or less in small rivers and agricultural canals and so forth. The nanohydraulic turbine, of which output power is less than several kW, can be easily carried to the places where they are necessary, and produces electric power easily without damaging the environment. Ikeda et al. [6] investigated the power characteristic of the turbine through a laboratory experiment and made clear the flow inside the rotor by the experimental visualization.
This study proposes the simulation method for the flow in the nanohydraulic turbine of Ikeda et al. [6], which can be favorably employed for the prediction of the turbine performance. The turbine performance is generally affected by the geometric conditions, such as the blade shape, the blade installation angle, and the number of blades. It is also influenced by the position of turbine relative to the water flow. The numerical simulation can examine the effect of each condition individually and promises to yield design guidelines for highperformance turbine. The nanohydraulic turbine is driven by a waterfall. The falling water collides with the rotating blades, and it is scattered around the rotor. Therefore, the flow through the turbine includes free surfaces, and it is very complicated. The simulation of this study is based on the Moving Particle Semiimplicit (MPS) method [7–9], which is one of the particle methods for freesurface flows. The fluid is discretized with particles, and the flow is simulated by the Lagrangian computation of the particle motion. There are few numerical simulations of the flow through smallscale hydraulic turbines utilizing smallscale hydropower, except for a MPS simulation of an impulsive turbine driven by a nozzle jet [5]. But the parameters for the MPS method, such as the distance between the particles, were not examined, and accordingly the sufficient knowledge on the simulation of smallscale hydraulic turbines was not obtained. First, this study elucidates the appropriate distance between the particles for the flow and performance simulations of the hydraulic turbine. Secondly, it is demonstrated that the present simulation can analyze the effect of the rotational speed of rotor on the turbine performance.
2. Basic Equations and Simulation Method
2.1. Conservation Equations for Flow
If the flow through the hydraulic turbine is incompressible, it is governed by the mass and momentum conservation equations: where is the density, is the time, is the velocity, is the pressure, is the kinematic viscosity, and is the external forces such as the gravitational force and the surface tension.
This study analyzes (1) and (2) by the Moving Particle Semiimplicit (MPS) method [7–9], which is one of the particle methods. In the MPS method, the fluid is discretized by particles, and the particle motion is computed by the Lagrangian method. Equations (1) and (2) are discretized through the interactions between the particles.
2.2. Particle Interaction Model
The interactions between the particles are modeled with a weight function defined by the following equation: where is the distance between two particles and stands for a kernel size.
The particle number density at the position of the th particle, , is defined as where is the position vector of the th particle.
As the fluid density remains unaltered in the incompressible flow, the particle number density is required to be constant. The incompressible flow condition in the MPS method is satisfied by maintaining at a constant value .
The Laplacian operator, expressing the viscous term on the righthand side of (2), is modeled with the weight function. The Laplacian operator at the position of the th particle is given as where is a physical quantity. A parameter is introduced so that the variance increase is equal to the analytical solution:
The gradient operator, expressing the pressure gradient term on the righthand side of (2), is also modeled. The gradient operator at the position of the th particle is modeled by setting the interparticle force at the repulsion to ensure the numerical stability [7–9]. Thus, the interacting pressure forces between two particles are not antisymmetric, and the momentum is not always conserved. To resolve this problem, this simulation employs the following model presented by Khayyer and Gotoh [10]: where is the number of space dimensions and is defined as
2.3. Simulation Method
Equations (1) and (2) are solved by a semiimplicit method, which is used in the SMAC method [11]. If the particle velocity and position at time are known, the flow at time is simulated by the following two steps.
In the first step, the temporal velocity and position for the particle, and , respectively, are calculated from (2) without considering the pressure gradient term. Then, the temporal particle number density is computed by using (4).
In the second step, the following Poisson equation is solved for the pressure so that the mass conservation is satisfied or is made to coincide with :
Then, the temporal velocity and position for the particle are corrected by the obtained pressure gradient: where
The viscous term in (2) and the lefthand side of (9) are computed by the Laplacian operator (5). The righthand side of (12) is calculated by the gradient operator (7).
3. Simulation Conditions
The flow through an impulsetype nanohydraulic turbine developed by Ikeda et al. [6] is simulated. The power characteristic of the turbine was investigated by a laboratory experiment, and the flow field was made clear by the experimental visualization [6]. Figure 1 shows the rotor and blade specifications. The rotor diameter is 200 mm, the axis diameter is 15 mm, and the number of blades is 12. The blade has a circulararc shape. The curvature radius is 23 mm, the thickness is 3 mm, and the chord length is 40 mm. To visualize the flow inside the rotor, the blades are sandwiched between two circular plates made of transparent vinyl chloride.
In the experiment, the rotor is placed at the bottom of a waterfall as shown in Figure 2. The origin of coordinates is set at the initial point of the waterfall. The axis is horizontal, and the axis is vertical. The flow through the rotor is simulated by a twodimensional MPS method. The water flow rate is 0.0035 m^{3}/s, and the head of waterfall is 570 mm. The nondimensional horizontal distance between the initial point of the waterfall and the colliding point of the water with the blade , , is 1.33.
The square region of 500 mm 500 mm around the rotor axis is chosen as the computational domain as shown in Figure 3, where the coordinates of the axis are supposed to be . The simulation is performed at the conditions of , where is defined by the rotor tip speed and the impact velocity of waterfall with the blade .
The waterfall flows into the computational domain from the location A on the upper boundary, as indicated in Figure 3. The particles, discretizing the falling water, are released from the location A into the domain. The releasing velocity and angle are determined from the condition of the free fall. It is required that the released particles are arranged at a uniform interval in the horizontal and vertical directions to ensure the computational accuracy. The particles are arranged at an interval in the horizontal direction and released into the domain at a constant time interval , so that the width of the waterfall at the location A corresponds to the measured value 8.4 mm. The value is prescribed so that the vertical distance between the particles coincides with . The distance corresponds to the grid width for gridbased simulations such as a finite difference method. Therefore, the simulated results are considered to depend on . This study investigates the effect of on the flow and the turbine performance to search for the appropriate value of .
The rotor axis and the blades, which are the solid walls, are discretized by the particles having the same angular velocity as the rotor. The distance between these particles is also set at . The Neumann boundary condition for the pressure gradient is imposed on the particles contacting with the fluid.
The pressure is set at zero on the free surface. The position of the free surface is detected according to the value of the particle number density. When the particle number density obtained by the firststep calculation of each time step satisfies the following relation, the th particle is decided to be on the free surface: where is a parameter of .
The time increment is determined from the maximum particle velocity at each computational time step. In this simulation, the initial value is set at 0.001 s, resulting in the maximum value of 0.001 s. The value of in (3) is generally chosen at [7, 8]. The value is for the particle number density and the gradient operator, while it is for the Laplacian operator [8]. The parameter in (13) is set at 0.97. It is reported that the simulation of a fragmentation of fluid scarcely depends on the value in the case of [7].
The simulation conditions are listed in Table 1.

4. Results and Discussion
4.1. Results at Tip Speed Ratio of
Figure 4 shows the time variation of the particle distribution for the fully developed flow, where the tip speed ratio is 0.56. The distributions at five time points during a period for the interaction between the waterfall and a blade are presented. The distance between the particles discretizing waterfall is set at . When the time is as in Figure 4(a), the waterfall collides directly with the tip of the concave surface for the blade A. At the subsequent time points of Figures 4(b), 4(c) and 4(d), the collision point moves toward the center of the concave surface for the blade A as the rotor rotates. When the time is as in Figure 4(e), a part of the waterfall gets in contact with the tip of the subsequent blade B. But the collision between the waterfall and the concave surface for the blade A is still maintained. The collision point is closest to the rotor axis during the one period for the interaction between the waterfall and a blade. On the concave surface of the blade A, the number of particles or the mass of water increases with the passage of time during the one period. When the time is as in Figure 4(d), the water on the concave surface of the blade A flows toward the subsequent blade B, and it collides with the convex surface of the blade B at the time of Figure 4(e). It should be noted that the blade C corresponds to a blade, which has just finished colliding with the waterfall. The water on the concave surface flows toward the inside and outside of the rotor. The water flowing toward the inside collides with the concave surface of the subsequent blade A, and then it enters into the rotor. The water directing toward the outside disperses markedly in the radial direction as the rotor rotates. The flow rate toward the outside is much larger. The particles always exist on the convex surface of the blades, which are on the opposite side of the waterfall. This demonstrates the stagnation of water inside the rotor. The blades give the angular momentum to the stagnant water and flick the water away from the rotor, causing the deterioration of the rotor performance.
(a)
(b)
(c)
(d)
(e)
Figure 5 shows the superposition for the particle distributions at every time interval 3 during 30 time period, where and . One can grasp the water dispersion around the rotor and the water stagnation in the rotor. The water scatters mainly toward the lower left direction and the waterfall direction.
The timeaveraged water velocity for is shown in Figure 6. The water dispersion toward the outside of the rotor and the flow into the rotor are reconfirmed.
The flow pattern inside and around the rotor, experimentally visualized by Ikeda et al. [6] at , is presented in Figure 7. It was acquired by using a CCD camera and a strobe light sheet shaped through a 2 mm wide slit. The image visualizes vividly the water flow along the concave surface of the blades, the water dispersion and flick toward the outside of the rotor, and the water stagnation inside the rotor. The simulated flow, shown in Figure 4, agrees well with the experimental visualization, demonstrating the validity of the present simulation.
The distance between the particles discretizing the waterfall or the initial distance between the particles corresponds to the grid width in gridbased simulation methods, such as a finite difference method. The simulation with smaller has superior space resolution. The timeaveraged number of particles, , for the fully developed flow is plotted against in Figure 8. For the abovementioned simulation of , is 14437. The value increases greatly with the decrement of . The improvement of the spatial resolution increases the number of particles, and therefore it causes the increment of the computational time. In the simulation of , 348 and 299 particles are used to discretize the axis and the single blade, respectively. These particles are included in Figure 8. For the flow simulation during the one revolution of the rotor, about 6.4 hours are required on a workstation (Processor: Intel Xeon X5660 2.8 GHz × 6, Memory: 12 GB).
The superimposed particles at and 0.12 distribute as shown in Figure 9, where . In the case of , the number of particles is low, and the water dispersion outside the rotor and the stagnation inside the rotor are not fully resolved. The simulation of , composed of more particles, yields almost the same result at shown in Figure 5. It is discovered that the simulation of can ensure a sufficiently high spatial resolution.
(a) In case of
(b) In case of
Calculating the rotor power by estimating the water kinematic energy at the rotor inlet and outlet, the nondimensional value changes as the function of as plotted in Figure 10. When the spatial resolution is high enough (), the value remains almost unaltered. It is slightly larger than the experimental result of . But it is considered to be valid in due consideration of the twodimensional simulation. Ikeda et al. [6] calculated the value at by estimating the force on a blade from the visualized flow pattern. It is also plotted in Figure 10, being in good agreement with the present simulation. Consequently, it is found that the simulation of can accurately predict the value.
Ikeda et al. [6] made it clear by their experiment that the horizontal distance between the initial point of the waterfall and the blade, affects the rotor performance . The distance varies with the water flow rate . For the practical use of the hydraulic turbine, it is desirable that is always set at the optimal value irrespective of . Ikeda et al. [6] proposed a method to control the value by installing a flat plate along the waterfall. The current simulation is effectively employed to search for the applicability of the method.
4.2. Results at Tip Speed Ratios of and 0.7
The flows at and 0.7 are simulated for the condition of . Figure 11 depicts the relation between and . The value decreases with the increment of . This change agrees with the experimental result of Ikeda et al. [6], indicating that the effect of on is successfully analyzed by the simulation. The simulated value is slightly larger. This is because the current simulation does not sufficiently resolve the turbulent flow, and therefore it ignores the losses caused by the turbulence on the blade surface as well as inside the rotor. It may be also because the current simulation employs the twodimensional MPS method. Using the threedimensional MPS method, the flow and the value would be simulated more accurately. But the number of particles increases, and a longer computational time is required.
The flow fields at are shown in Figures 12 and 13. The waterfall disperses due to the collision with the rotor. When compared with the results at (Figures 5 and 6), the divergence angle of the dispersed water lessens. One can grasp the decrement of the rotor angular momentum obtained from the waterfall. The amount of water inside the rotor is less than that at .
5. Conclusions
The flow through an impulsetype smallscale hydraulic turbine utilizing a waterfall of extralow head is simulated by a twodimensional MPS method. The rotor performance is also analyzed by using the simulated flow field. The results are summarized as follows.(1)When the distance between the particles discretizing the waterfall of a width , , is set at , the flow simulated at the tip speed ratio is confirmed to agree well with the experimentally visualized one. Thus, the simulation of has a sufficiently high spatial resolution.(2)For the simulation of at , it is also confirmed that the simulated power coefficient agrees nearly with the experiment. Thus, the simulation of can also favorably predict the rotor performance.(3)The values simulated at and 0.7 agree almost with the experimental results. Therefore, the present simulation can successfully analyze the effect of on the flow and the rotor performance.
Nomenclatures
:  Thickness of waterfall 
:  Power coefficient 
:  Number of space dimensions 
:  Diameter of rotor 
:  External force 
:  Head of waterfall 
:  Distance between particles discretizing waterfall 
:  Horizontal distance between blade and initial point of waterfall 
:  Particle number density 
:  Timeaveraged number of particles 
:  Pressure 
:  Power output from rotor 
:  Water flow rate 
:  Position vector of particle 
:  Time 
:  Velocity 
:  Impactvelocity of waterfall with blade 
:  Rotor tip speed 
:  Weight function 
:  Spatial coordinates 
:  Time increment 
:  Tip speed ratio 
:  Kinematic viscosity 
:  Density 
:  Angular velocity of rotor. 
References
 Y. Takamatsu, A. Furukawa, K. Okuma, and K. Takenouchi, “Experimental studies on a preferable blade profile for high efficiency and the blade characteristics of Darrieustype crossflow water turbines,” JSME International Journal, vol. 34, no. 2, pp. 149–156, 1991. View at: Google Scholar
 M. Nakajima, S. Iio, and T. Ikeda, “Performance of doublestep Savonius rotor for environmentally friendly hydraulic turbine,” Journal of Fluid Science and Technology, vol. 3, no. 3, pp. 410–419, 2008. View at: Google Scholar
 M. Nakajima, S. Iio, and T. Ikeda, “Performance of Savonius rotor for environ¬mentally friendly hydraulic turbine,” Journal of Fluid Science and Technology, vol. 3, pp. 420–429, 2008. View at: Google Scholar
 S. Derakhshan and A. Nourbakhsh, “Experimental study of characteristic curves of centrifugal pumps working as turbines in different specific speeds,” Experimental Thermal and Fluid Science, vol. 32, no. 3, pp. 800–807, 2008. View at: Publisher Site  Google Scholar
 Y. Nakanishi, S. Iio, Y. Takahashi, A. Kato, and T. Ikeda, “Development of a simple impulse turbine for nano hydropower,” Journal of Fluid Science and Technology, vol. 4, pp. 567–577, 2009. View at: Google Scholar
 T. Ikeda, S. Iio, and K. Tatsuno, “Performance of nanohydraulic turbine utilizing waterfalls,” Renewable Energy, vol. 35, no. 1, pp. 293–300, 2010. View at: Publisher Site  Google Scholar
 S. Koshizuka and Y. Oka, “Movingparticle semiimplicit method for fragmentation of incompressible fluid,” Nuclear Science and Engineering, vol. 123, no. 3, pp. 421–434, 1996. View at: Google Scholar
 S. Koshizuka, A. Nobe, and Y. Oka, “Numerical analysis of breaking waves using the moving particle semiimplicit method,” International Journal for Numerical Methods in Fluids, vol. 26, no. 7, pp. 751–769, 1998. View at: Google Scholar
 K. Shibata, S. Koshizuka, M. Sakai, and K. Tanizawa, “Lagrangian simulations of shipwave interactions in rough seas,” Ocean Engineering, vol. 42, pp. 13–25, 2012. View at: Publisher Site  Google Scholar
 A. Khayyer and H. Gotoh, “Development of CMPS method for accurate watersurface tracking in breaking waves,” Coastal Engineering Journal, vol. 50, no. 2, pp. 179–207, 2008. View at: Publisher Site  Google Scholar
 A. A. Amsden and F. H. Harlow, “The SMAC method: a numerical technique for calculating incompressible fluid flows,” Los Alamos Scientific Laboratory Report LA4370, 1970. View at: Google Scholar
Copyright
Copyright © 2013 Tomomi Uchiyama et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.