# Numerical Simulation of Water Flow through a Nano-Hydraulic Turbine of Waterfall-Type by Particle Method

Authors: Tomomi Uchiyama; Haruki Fukuhara; Shouichiro Iio; Toshihiko Ikeda
Journal: International Journal of Rotating Machinery
 (2013)

---

## Abstract

This study simulates the flow through an impulse-type small-scale hydraulic turbine utilizing a waterfall of extra-low head. The two-dimensional Moving Particle Semi-implicit (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 widthB, l0, is set at l0/B≤0.084, 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.


---

## Body

## 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 large-scale hydroelectric plants require huge dams and long conduits, the places for the construction are hardly remained. Thus, expectations for the development of a small-scale hydropower, of which output is less than 1000 kW, have been increasing. When natural disaster occurs, the large-scale centralized hydraulic power plants may lose the power supplying ability due to the collapse of the power grid. Since the small-scale hydropower, existing in small-scale rivers, irrigation canals, and industrial drainages, distributes widely in Japan, it realizes the small-scale distributed power generation. The small-scale 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 disaster-resistant society.To exploit effectively the small-scale hydropower, various types of hydraulic turbine have been presented [1–5]. Ikeda et al. [6] developed an impulse-type hydraulic turbine utilizing waterfalls of extra-low head which is 2 m or less in small rivers and agricultural canals and so forth. The nano-hydraulic 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 nano-hydraulic 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 high-performance turbine. The nano-hydraulic 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 Semi-implicit (MPS) method [7–9], which is one of the particle methods for free-surface 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 small-scale hydraulic turbines utilizing small-scale 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 small-scale 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:(1)DρDt=0,(2)DuDt=-1ρ∇p+ν∇2u+F,
where ρ is the density, t is the time, u is the velocity, p is the pressure, ν is the kinematic viscosity, and F is the external forces such as the gravitational force and the surface tension.This study analyzes (1) and (2) by the Moving Particle Semi-implicit (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 functionw defined by the following equation:
(3)w(r)={rer-1r≤re,0r>re,
where r is the distance between two particles and re stands for a kernel size.The particle number density at the position of theith particle, 〈n〉i, is defined as
(4)〈n〉i=∑j≠i‍w(|rj-ri|),
where ri is the position vector of the ith 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〈n〉i at a constant value n0.The Laplacian operator, expressing the viscous term on the right-hand side of (2), is modeled with the weight function. The Laplacian operator at the position of the ith particle is given as
(5)〈∇2ϕ〉i=2dn0Λ∑j≠i‍[(ϕj-ϕi)w(|rj-ri|)],
where ϕ is a physical quantity. A parameter Λ is introduced so that the variance increase is equal to the analytical solution:
(6)Λ=∫V‍w(r)r2dv∫V‍w(r)dv.The gradient operator, expressing the pressure gradient term on the right-hand side of (2), is also modeled. The gradient operator at the position of the ith 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]:
(7)〈∇ϕ〉i=dn0∑j≠i‍[(ϕi+ϕj)-(ϕ^i+ϕ^j)|rj-ri|2(rj-ri)w(|rj-ri|)],
where d is the number of space dimensions and ϕ^i is defined as
(8)ϕi^=minj∈J(ϕi,ϕj),J={j:w(|rj-ri|)≠0}.


### 2.3. Simulation Method

Equations (1) and (2) are solved by a semi-implicit method, which is used in the SMAC method [11]. If the particle velocity uik and position rik at time t=kΔt are known, the flow at time t=(k+1)Δt is simulated by the following two steps.In the first step, the temporal velocity and position for the particle,ui* and ri*, respectively, are calculated from (2) without considering the pressure gradient term. Then, the temporal particle number density 〈n*〉i is computed by using (4).In the second step, the following Poisson equation is solved for the pressurepk+1 so that the mass conservation is satisfied or 〈n*〉i is made to coincide with n0:
(9)〈∇2pk+1〉i=-ρΔt2〈n*〉i-n0n0.Then, the temporal velocity and position for the particle are corrected by the obtained pressure gradient:(10)uik+1=ui*+ui′(11)rik+1=ri*+ui′Δt,
where
(12)ui′=-Δtρ〈∇pk+1〉i.The viscous term in (2) and the left-hand side of (9) are computed by the Laplacian operator (5). The right-hand side of (12) is calculated by the gradient operator (7).


## 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:(1)DρDt=0,(2)DuDt=-1ρ∇p+ν∇2u+F,
where ρ is the density, t is the time, u is the velocity, p is the pressure, ν is the kinematic viscosity, and F is the external forces such as the gravitational force and the surface tension.This study analyzes (1) and (2) by the Moving Particle Semi-implicit (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 functionw defined by the following equation:
(3)w(r)={rer-1r≤re,0r>re,
where r is the distance between two particles and re stands for a kernel size.The particle number density at the position of theith particle, 〈n〉i, is defined as
(4)〈n〉i=∑j≠i‍w(|rj-ri|),
where ri is the position vector of the ith 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〈n〉i at a constant value n0.The Laplacian operator, expressing the viscous term on the right-hand side of (2), is modeled with the weight function. The Laplacian operator at the position of the ith particle is given as
(5)〈∇2ϕ〉i=2dn0Λ∑j≠i‍[(ϕj-ϕi)w(|rj-ri|)],
where ϕ is a physical quantity. A parameter Λ is introduced so that the variance increase is equal to the analytical solution:
(6)Λ=∫V‍w(r)r2dv∫V‍w(r)dv.The gradient operator, expressing the pressure gradient term on the right-hand side of (2), is also modeled. The gradient operator at the position of the ith 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]:
(7)〈∇ϕ〉i=dn0∑j≠i‍[(ϕi+ϕj)-(ϕ^i+ϕ^j)|rj-ri|2(rj-ri)w(|rj-ri|)],
where d is the number of space dimensions and ϕ^i is defined as
(8)ϕi^=minj∈J(ϕi,ϕj),J={j:w(|rj-ri|)≠0}.


## 2.3. Simulation Method

Equations (1) and (2) are solved by a semi-implicit method, which is used in the SMAC method [11]. If the particle velocity uik and position rik at time t=kΔt are known, the flow at time t=(k+1)Δt is simulated by the following two steps.In the first step, the temporal velocity and position for the particle,ui* and ri*, respectively, are calculated from (2) without considering the pressure gradient term. Then, the temporal particle number density 〈n*〉i is computed by using (4).In the second step, the following Poisson equation is solved for the pressurepk+1 so that the mass conservation is satisfied or 〈n*〉i is made to coincide with n0:
(9)〈∇2pk+1〉i=-ρΔt2〈n*〉i-n0n0.Then, the temporal velocity and position for the particle are corrected by the obtained pressure gradient:(10)uik+1=ui*+ui′(11)rik+1=ri*+ui′Δt,
where
(12)ui′=-Δtρ〈∇pk+1〉i.The viscous term in (2) and the left-hand side of (9) are computed by the Laplacian operator (5). The right-hand side of (12) is calculated by the gradient operator (7).


## 3. Simulation Conditions

The flow through an impulse-type nano-hydraulic 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 D is 200 mm, the axis diameter is 15 mm, and the number of blades is 12. The blade has a circular-arc 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.Figure 1
Rotor and blade specifications.In the experiment, the rotor is placed at the bottom of a waterfall as shown in Figure2. The origin of coordinates is set at the initial point of the waterfall. The x-axis is horizontal, and the z-axis is vertical. The flow through the rotor is simulated by a two-dimensional MPS method. The water flow rate Q is 0.0035 m3/s, and the head of waterfall HF is 570 mm. The nondimensional horizontal distance between the initial point of the waterfall and the colliding point of the water with the blade LF, LF/D, is 1.33.Figure 2
Positional relation between waterfall and rotor.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 (xa,za). The simulation is performed at the conditions of 0.56≤λ≤0.7, where λ(=Vt/UF) is defined by the rotor tip speed Vt and the impact velocity of waterfall with the blade UF.Figure 3
Computational domain.The waterfall flows into the computational domain from the location A on the upper boundary, as indicated in Figure3. 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 l0 in the horizontal direction and released into the domain at a constant time interval Δtp, so that the width of the waterfall B at the location A corresponds to the measured value 8.4 mm. The Δtp value is prescribed so that the vertical distance between the particles coincides with l0. The distance l0 corresponds to the grid width for grid-based simulations such as a finite difference method. Therefore, the simulated results are considered to depend on l0. This study investigates the effect of l0 on the flow and the turbine performance to search for the appropriate value of l0.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 atl0. 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〈n*〉i obtained by the first-step calculation of each time step satisfies the following relation, the ith particle is decided to be on the free surface:
(13)〈n*〉i<βn0,
where β is a parameter of β<1.The time incrementΔt 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 re in (3) is generally chosen at 2≤re/l0≤4 [7, 8]. The value is 2.1l0 for the particle number density and the gradient operator, while it is 4l0 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 0.8≤β≤0.99 [7].The simulation conditions are listed in Table1.Table 1
Simulation conditions.
Rotor diameterD
0.2 m
Number of blades
12
Water flow rateQ
0.0035m3/s
Vertical velocity of waterfall colliding with bladeUF
3.34 m/s
Head of waterfallHF
0.57 m
Horizontal distance between blade and initial point of waterfall
1.33D
Distance between particles discretizing waterfalll0
0.069B–0.12B
Maximum of time incrementΔt
0.001 s


## 4. Results and Discussion

### 4.1. Results at Tip Speed Ratio ofλ=0.56

Figure4 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 l0 is set at l0/B=0.084. 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.Time variation for particle distribution whenλ=0.56.
(a)
t
=
t
1
[
s
]
(b)
t
=
t
1
+
0.005
(c)
t
=
t
1
+
0.01
(d)
t
=
t
1
+
0.015
(e)
t
=
t
1
+
0.02Figure5 shows the superposition for the particle distributions at every time interval 3Δt during 30Δt time period, where λ=0.56 and l0/B=0.084. 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.Figure 5
Distribution of superimposed particles whenλ=0.56.The time-averaged water velocity forλ=0.56 is shown in Figure 6. The water dispersion toward the outside of the rotor and the flow into the rotor are reconfirmed.Figure 6
Time-averaged velocity distribution whenλ=0.56.The flow pattern inside and around the rotor, experimentally visualized by Ikeda et al. [6] at λ=0.56, 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.Figure 7
Experimentally visualized flow pattern whenλ=0.56.The distance between the particles discretizing the waterfalll0 or the initial distance between the particles corresponds to the grid width in grid-based simulation methods, such as a finite difference method. The simulation with smaller l0 has superior space resolution. The time-averaged number of particles, N, for the fully developed flow is plotted against l0 in Figure 8. For the abovementioned simulation of l0/B=0.084, N is 14437. The N value increases greatly with the decrement of l0. 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 l0/B=0.084, 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).Figure 8
Change in number of particlesN due to distance between particles l0 when λ=0.56.The superimposed particles atl0/B=0.069 and 0.12 distribute as shown in Figure 9, where λ=0.56. In the case of l0/B=0.12, 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 l0/B=0.069, composed of more particles, yields almost the same result at l0/B=0.084 shown in Figure 5. It is discovered that the simulation of l0/B≤0.084 can ensure a sufficiently high spatial resolution.Effect of distance between particlesl0 on particle distribution when λ=0.56.
(a)
In case ofl0/B=0.12
(b)
In case ofl0/B=0.069Calculating the rotor powerP by estimating the water kinematic energy at the rotor inlet and outlet, the nondimensional value Cp changes as the function of l0/B as plotted in Figure 10. When the spatial resolution is high enough (l0/B≤0.084), the Cp value remains almost unaltered. It is slightly larger than the experimental result of Cp=0.66. But it is considered to be valid in due consideration of the two-dimensional simulation. Ikeda et al. [6] calculated the Cp value at λ=0.56 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 l0/B≤0.084 can accurately predict the Cp value.Figure 10
Change in power coefficientCp due to distance between particles l0 when λ=0.56.Ikeda et al. [6] made it clear by their experiment that the horizontal distance between the initial point of the waterfall and the blade, LF affects the rotor performance Cp. The distance LF varies with the water flow rate Q. For the practical use of the hydraulic turbine, it is desirable that LF is always set at the optimal value irrespective of Q. Ikeda et al. [6] proposed a method to control the LF 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λ=0.6 and 0.7

The flows atλ=0.6 and 0.7 are simulated for the condition of l0/B=0.084. Figure 11 depicts the relation between λ and Cp. The Cp value decreases with the increment of λ. This change agrees with the experimental result of Ikeda et al. [6], indicating that the effect of λ on Cp is successfully analyzed by the simulation. The simulated Cp 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 two-dimensional MPS method. Using the three-dimensional MPS method, the flow and the Cp value would be simulated more accurately. But the number of particles increases, and a longer computational time is required.Figure 11
Relation between power coefficientCp and tip speed ratio λ.The flow fields atλ=0.7 are shown in Figures 12 and 13. The waterfall disperses due to the collision with the rotor. When compared with the results at λ=0.56 (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 λ=0.56.Figure 12
Distribution of superimposed particles whenλ=0.7.Figure 13
Time-averaged velocity distribution whenλ=0.7.


## 4.1. Results at Tip Speed Ratio ofλ=0.56

Figure4 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 l0 is set at l0/B=0.084. 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.Time variation for particle distribution whenλ=0.56.
(a)
t
=
t
1
[
s
]
(b)
t
=
t
1
+
0.005
(c)
t
=
t
1
+
0.01
(d)
t
=
t
1
+
0.015
(e)
t
=
t
1
+
0.02Figure5 shows the superposition for the particle distributions at every time interval 3Δt during 30Δt time period, where λ=0.56 and l0/B=0.084. 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.Figure 5
Distribution of superimposed particles whenλ=0.56.The time-averaged water velocity forλ=0.56 is shown in Figure 6. The water dispersion toward the outside of the rotor and the flow into the rotor are reconfirmed.Figure 6
Time-averaged velocity distribution whenλ=0.56.The flow pattern inside and around the rotor, experimentally visualized by Ikeda et al. [6] at λ=0.56, 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.Figure 7
Experimentally visualized flow pattern whenλ=0.56.The distance between the particles discretizing the waterfalll0 or the initial distance between the particles corresponds to the grid width in grid-based simulation methods, such as a finite difference method. The simulation with smaller l0 has superior space resolution. The time-averaged number of particles, N, for the fully developed flow is plotted against l0 in Figure 8. For the abovementioned simulation of l0/B=0.084, N is 14437. The N value increases greatly with the decrement of l0. 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 l0/B=0.084, 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).Figure 8
Change in number of particlesN due to distance between particles l0 when λ=0.56.The superimposed particles atl0/B=0.069 and 0.12 distribute as shown in Figure 9, where λ=0.56. In the case of l0/B=0.12, 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 l0/B=0.069, composed of more particles, yields almost the same result at l0/B=0.084 shown in Figure 5. It is discovered that the simulation of l0/B≤0.084 can ensure a sufficiently high spatial resolution.Effect of distance between particlesl0 on particle distribution when λ=0.56.
(a)
In case ofl0/B=0.12
(b)
In case ofl0/B=0.069Calculating the rotor powerP by estimating the water kinematic energy at the rotor inlet and outlet, the nondimensional value Cp changes as the function of l0/B as plotted in Figure 10. When the spatial resolution is high enough (l0/B≤0.084), the Cp value remains almost unaltered. It is slightly larger than the experimental result of Cp=0.66. But it is considered to be valid in due consideration of the two-dimensional simulation. Ikeda et al. [6] calculated the Cp value at λ=0.56 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 l0/B≤0.084 can accurately predict the Cp value.Figure 10
Change in power coefficientCp due to distance between particles l0 when λ=0.56.Ikeda et al. [6] made it clear by their experiment that the horizontal distance between the initial point of the waterfall and the blade, LF affects the rotor performance Cp. The distance LF varies with the water flow rate Q. For the practical use of the hydraulic turbine, it is desirable that LF is always set at the optimal value irrespective of Q. Ikeda et al. [6] proposed a method to control the LF 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λ=0.6 and 0.7

The flows atλ=0.6 and 0.7 are simulated for the condition of l0/B=0.084. Figure 11 depicts the relation between λ and Cp. The Cp value decreases with the increment of λ. This change agrees with the experimental result of Ikeda et al. [6], indicating that the effect of λ on Cp is successfully analyzed by the simulation. The simulated Cp 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 two-dimensional MPS method. Using the three-dimensional MPS method, the flow and the Cp value would be simulated more accurately. But the number of particles increases, and a longer computational time is required.Figure 11
Relation between power coefficientCp and tip speed ratio λ.The flow fields atλ=0.7 are shown in Figures 12 and 13. The waterfall disperses due to the collision with the rotor. When compared with the results at λ=0.56 (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 λ=0.56.Figure 12
Distribution of superimposed particles whenλ=0.7.Figure 13
Time-averaged velocity distribution whenλ=0.7.


## 5. Conclusions

The flow through an impulse-type small-scale hydraulic turbine utilizing a waterfall of extra-low head is simulated by a two-dimensional 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 widthB, l0, is set at l0/B≤0.084, the flow simulated at the tip speed ratio λ=0.56 is confirmed to agree well with the experimentally visualized one. Thus, the simulation of l0/B≤0.084 has a sufficiently high spatial resolution.(2)
For the simulation ofl0/B≤0.084 at λ=0.56, it is also confirmed that the simulated power coefficient Cp agrees nearly with the experiment. Thus, the simulation of l0/B≤0.084 can also favorably predict the rotor performance.(3)
TheCp values simulated at λ=0.6 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.



---

*Source: 473842-2013-11-18.xml*