Research Article  Open Access
Takayuki Matsunuma, Takehiko Segawa, "Effects of Input Voltage on Flow Separation Control for LowPressure Turbine at Low Reynolds Number by Plasma Actuators", International Journal of Rotating Machinery, vol. 2012, Article ID 902548, 10 pages, 2012. https://doi.org/10.1155/2012/902548
Effects of Input Voltage on Flow Separation Control for LowPressure Turbine at Low Reynolds Number by Plasma Actuators
Abstract
Active flow control using dielectric barrier discharge (DBD) plasma actuators was investigated to reattach the simulated boundary layer separation on the suction surface of a turbine blade at low Reynolds number, Re = 1.7 × 10^{4}. The flow separation is induced on a curved plate installed in the test section of a lowspeed wind tunnel. Particle image velocimetry (PIV) was used to obtain instantaneous and timeaveraged twodimensional velocity measurements. The amplitude of input voltage for the plasma actuator was varied from ±2.0 kV to ±2.8 kV. The separated flow reattached on the curved wall when the input voltage was ±2.4 kV and above. The displacement thickness of the boundary layer near the trailing edge decreased by 20% at ±2.0 kV. The displacement thickness was suddenly reduced as much as 56% at ±2.2 kV, and it was reduced gradually from ±2.4 kV to ±2.8 kV (77% reduction). The total pressure loss coefficient, estimated from the boundary layer displacement thickness and momentum thickness, was 0.172 at the baseline (actuator off) condition. The total pressure loss was reduced to 0.107 (38% reduction) at ±2.2 kV and 0.078 (55% reduction) at ±2.8 kV.
1. Introduction
Blade Reynolds numbers for the lowpressure turbine of small and mediumsized gas turbines for aircraft propulsion can drop to below 2.5 × 10^{4} at high altitudes, where the air density is low [1]. At these low Reynolds numbers, the boundary layer is dominated by laminar flow and is susceptible to flow separation and strong secondary vortices that develop into increased losses, leading to reduced performance.
Newly developed highefficiency smallscale gas turbines for industrial power generation also encounter the low Reynolds number problem because of the increased viscosity caused by high turbine inlet temperatures and miniaturization of the cascade [2]. Several studies focusing on the aerodynamics of turbine cascades at low Reynolds numbers have been published [3–7].
Different passive and active flow control techniques have been developed to increase the efficiency of the turbine blade at low Reynolds numbers [8]. Passive control devices, such as dimples and fixed turbulators, are traditional approaches to trigger transition and induce reattachment, but create undesirable drag at high Reynolds number. Active control devices can be employed only when needed. Vortex generating jets (VGJs) and plasma actuators are prominent subjects of recent research in active flow control.
The application of dielectric barrier discharge (DBD) plasma actuator was demonstrated by Roth et al. [9] in 1998 and has been developed over the last decade from fundamental studies to a wide range of applications [10].
Figure 1 shows a simple schematic configuration of a DBD plasma actuator. The DBD plasma actuator consists of a dielectric layer sandwiched between top and bottom electrodes. Applying high voltages at high frequencies (approximately ) between the top and bottom electrodes forms a layer of glow discharge plasma across the dielectric surface. This induces a tangential air jet with a strong horizontal velocity component. It is known that the direction of the resultant tangential jet is constant and independent of which electrode the voltage is applied to, and independent of the polarity of the applied voltage [11].
DBD plasma actuators have also been used in a number of lowpressure turbine separation control studies [12–18]. In order to distinguish this work from these prior studies, this paper focused on the presentation of detailed PIV flow data and accompanying analysis at various input voltage of the plasma actuator.
The objective of this study is to investigate the most effective active flow control operation of DBD plasma actuators by simplified experiments. The experimental results in this study will be used as the fundamental data for considering the application of DBD plasma actuators for the annular turbine wind tunnel with a singlestage axialflow turbine [19, 20]. This paper presents the results of particle image velocimetry (PIV) measurement to understand the flow mechanisms and the effect of different amplitude of the plasma actuator input voltage.
2. Experimental Facility and Method
Figure 2 shows the measurement system. The wind tunnel is a lowspeed, opencircuit, blowertype facility with 305 mm × 85 mm × 65 mm test section.
A curved wall plate (streamwise length mm) was installed in the test section of the wind tunnel in order to simulate the separated flow on the suction surface of a turbine blade. Figure 3 shows the geometry and the designsurface velocity distribution, derived from an inviscid calculation, at the midspan of the corresponding turbine rotor blade in the annular turbine wind tunnel in AIST [19, 20]. Table 1 shows the nondimensional parameters of the corresponding turbine rotor midspan. The shape of the curved wall was designed using a simple onedimensional continuity argument to match the design surface velocity and pressure distributions of the corresponding turbine blade. A plasma actuator was mounted at the front of the adverse pressure gradient region (deceleration region) on the curved wall.

Particle image velocimetry (PIV) was employed to quantify the behavior of the flow field around the curved wall. The laser was a 25 mJ/pulse, doublepulse NdYAG laser (MiniLase II, 20 Hz, New Wave Research Co. Ltd.). Atomized Dioctyl sebacate (DOS) oil with a mean particle diameter of 1 μm was injected upstream of the test section via a pressurized oil chamber. Image pairs were taken by a camera (PIV CAM 138, TSI Inc.), which has 1,280 × 1,024 pixel resolution. TSI software calculated the velocity vectors from the peak correlation of groups of particles between frames, using conventional crosscorrelation algorithms on a 32 × 32 pixel grid. In order to calculate the averaged velocity distributions, 50 instantaneous velocity distributions were measured at each experimental condition in this study. It should be noted that the PIV data in the vicinity of the plasma actuator possesses lower reliability than that in the main flow, due to the laser reflections from the top electrode of the plasma actuator, the large velocity gradients across the induced jet, and the possibility that the seeding particles are influenced by the electric field [21, 22].
Figure 4 shows the test section geometry and PIV measurement grids. In these experiments, the rotating speed of the wind tunnel blower was kept constant. The streamwise velocity of the main flow near the trailing edge ( mm) was almost constant, m/s, for each input voltage condition, as shown later in Figure 11. Therefore, the Reynolds number based on the streamwise length of the curved wall, , and the streamwise velocity of the main flow near the trailing edge, , was Re = 1.7 × 10^{4}.
The edge of the top electrode of the plasma actuator, where the surface plasma is formed, was located at mm from the leading edge. A thin sheet of polymide (125 μm thickness, relative permittivity ) was used as the dielectric barrier of the plasma actuator. As for the electrodes, thin sheets of copper (thickness = 35 μm) were glued and pressed onto the both sides of the dielectric. A high voltage, highfrequency power supply (PG1040F, PSI Inc.), which can output bipolar sinusoidal waveforms, was used to supply input signals to the top and bottom electrodes. The bottom electrode was connected to ground, V. The amplitude of input voltage to the top electrode was varied from 2.0 kV to 2.8 kV (4.0–5.6 k). The frequency of input voltage was fixed at kHz.
3. Results and Discussion
3.1. Instantaneous Absolute Velocity Distributions
Figure 5 shows a sample of the PIV images, the instantaneous absolute velocity, and vorticity distributions for the baseline condition (plasma actuator off). Massive flow separation is visualized in the flow deceleration region in the PIV image of Figure 5(a). The large lowvelocity region, drawn by blue contours in Figure 5(b), corresponds to the flow separation. (The criteria for separation or reattachment of the boundary layer are based on the boundary layer shape factor , as shown in Figure 15 later). The separated boundary layer is stable from the streamwise location mm (starting position of the flow separation) to mm (namely 70% chord length position). The separated flow downstream of mm becomes unstable and large clockwise vortices are generated.
(a) PIV image
(b) Instantaneous absolute velocity
(c) Vorticity
Figure 6 shows a sample of the PIV images, the instantaneous absolute velocity and vorticity distributions for the flow control condition (plasma actuator on, input voltage kV, maximum input voltage in this study). The separated flow region is dramatically reduced by flow control using the plasma actuator. The flow separation from mm became unstable after and small clockwise vortices were generated. The similar phenomenon of the smaller flow structure by the effect of the plasma actuator was observed in the flow visualization by Marks et al. [18]. The flow control using plasma actuator results in attached flow before the trailing edge.
(a) PIV image
(b) Instantaneous absolute velocity
(c) Vorticity
3.2. TimeAveraged Absolute Velocity Distributions
Figure 7 shows enlarged views of the timeaveraged absolute velocity distributions in the separation region (blue dashed frame in Figure 4) at various input voltages. These distributions were calculated by averaging 50 instantaneous velocity distributions as shown in Figures 5 and 6. In these figures, streamlines were superimposed upon the velocity contours.
The large separation region observed in the baseline (actuator off) condition in Figure 7(a) is gradually reduced by the effect of active control in Figures 7(b), 7(c), 7(d), 7(e), and 7(f).
In the baseline condition in Figure 7(a), large reverse flow is observed in the separation region. In the flow control case of kV in Figure 6(b), a large recirculation exists in the separation region, and the lowvelocity region is slightly reduced near the trailing edge. In the flow control case of kV in Figure 6(c), the recirculation in the separation region become reduced in size and the position of its center moves upstream. The separated flow reattaches onto the surface of the curved wall before the trailing edge. In Figures 6(d), 6(e), and 6(f), at , as the input voltage is increased, the recirculation in the separation region become smaller and moves more upstream. The reduction of the separation region was observed in the LDV measurements by Huang et al. [13], the PIV measurements by Boxx et al. [16], total pressure measurements by Burman et al. [17], and the surface pressure measurements by Marks et al. [18].
3.3. Turbulence Intensity Distributions
Figure 8 shows turbulence intensity distributions at various input voltages. In these figures, velocity vectors were superimposed upon the turbulence intensity contours. The turbulence intensity was calculated as follows: The turbulence intensity was normalized by the averaged streamwise velocity of the main flow near the trailing edge at mm, m/s, as shown in Figure 11 later.
In the baseline condition in Figure 8(a), a high turbulence intensity region exists at the boundary line between the main flow and the separated flow. A wide high turbulence intensity region is suddenly generated after mm. This phenomenon corresponds to the fact that the separated boundary layer is stable from mm to mm and becomes unstable after mm, as shown in Figure 5. The wide high turbulence intensity region after mm is due to the unstable transitional separated boundary layer.
In the flow control case of kV in Figure 8(b), the high turbulence intensity region is generated after mm, which is further upstream than the baseline condition.
In the flow control case of kV in Figure 8(c), the high turbulence intensity region moves even further upstream and closer to the curved wall than that in kV condition, Figure 8(b).
In Figures 8(c), 8(d), 8(e), and 8(f), as the input voltage is increased, the high turbulence intensity region moves more upstream and closer to the curved wall and becomes smaller in width.
3.4. Velocity Distributions Around Plasma Actuator
Figure 9 shows the vertical distribution of the streamwise velocity near the center of the plasma actuator, mm, (red dashed line in Figure 4). The vertical position mm corresponds to the surface of the curved wall. In the baseline condition (black line in the figure), the average streamwise velocity of the main flow is 2.05 m/s. The streamwise velocity in the boundary layer at mm is close to zero due to the flow separation near the curved wall. At kV (red line), the average streamwise velocity of the main flow increases to 2.26 m/s (10% increase from the baseline condition). The streamwise velocity in the boundary layer increases to 0.73 m/s. At kV (pink), the average streamwise velocity of the main flow increases to 2.48 m/s (21% increase from the baseline) and that in the boundary layer increases to 0.89 m/s. At kV (orange), the average streamwise velocity of the main flow increases slightly to 2.51 m/s and that in the boundary layer increases to 1.27 m/s. At maximum input voltage of kV (blue), the average streamwise velocity of the main flow increases to 2.58 m/s (26% increase from the baseline) and that in the boundary layer increases to 1.66 m/s.
Figure 10 shows the vertical distribution of the vertical velocity near the plasma actuator. The positive value means upward flow in the direction and the negative value means downward flow.
In the baseline condition (black line), the vertical velocity of the main flow at mm is 0.27 m/s. As the input voltage is increased, the vertical velocity of the main flow is gradually reduced by the effect of active flow control of the plasma actuator. At maximum voltage of kV (blue), the vertical velocity of the main flow at mm becomes −0.07 m/s (0.34 m/s decrease from the baseline condition).
The vertical velocity of the boundary layer at mm is almost zero in the baseline condition. At kV (red line), the vertical velocity in the boundary layer increases to 0.3 m/s due to the induced flow by the plasma actuator. As the input voltage is increased, the vertical velocity in the boundary layer is gradually reduced.
Figure 11 shows the vertical distribution of the flow angle near the plasma actuator. The positive value means upward flow and the negative value means downward flow. As the input voltage is increased, the flow angle in both the main flow and in the boundary layer is gradually reduced by the effect of active flow control of the plasma actuator (i.e., decrease in separation angle). At kV, the flow angle is reduced 9 deg in the main flow and 23 deg in the boundary layer, respectively.
3.5. Velocity Distributions Near Trailing Edge
Figure 12 shows the vertical distribution of the streamwise velocity near the trailing edge, mm, (brown dashed line in Figure 4). The bottom of the vertical axis (vertical position mm) corresponds to the surface of the curved wall at the trailing edge. The streamwise velocity of the main flow is almost constant ( m/s) for all input voltage conditions. In the baseline condition (black line in the figure), large velocity deficit with reverse flow exists near the surface of the curved wall. The boundary layer thickness (distance between the curved wall surface and the vertical position with 99% of the main flow velocity) is approximately 21 mm (from mm to mm). At kV (red line), the reverse flow is reduced but the boundary layer thickness is the same as the baseline. At kV (pink), the velocity in the boundary layer is suddenly increased. The reverse flow region completely disappears and the boundary layer thickness is reduced to 17 mm. At kV (orange), the shape of the velocity distributions in the boundary layer becomes quite linear. At kV and kV (green and blue), the velocity distributions in the boundary layer is further improved and the shape becomes convex.
Figure 13 shows vertical distributions of the vertical velocity near the trailing edge, mm. As the input voltage is increased, the vertical velocity in both the main flow and in the boundary layer is reduced (maximum 0.42 m/s reduction at kV).
Figure 14 shows vertical distributions of the turbulence intensity near the trailing edge, mm (dashed lines in Figure 8). Although the turbulence intensity of the main flow at the baseline is 2.1% for the baseline, that of the flow control conditions increase to 2.8% at kV and to at . The slight increase in the turbulence intensity of the main flow is due to the flow fluctuation generated by the plasma actuator AC operation. In the baseline condition, the maximum turbulence intensity in the boundary layer is 17.5%. The maximum turbulence intensity increases to 20% at , and decreases to 16% and 14% at kV at kV, respectively. As the input voltage is increased, the vertical location of the maximum turbulence intensity moves closer to the wall, suggesting reattachment of the separated boundary layer. The similar phenomenon of the turbulence intensity change by the effect of the plasma actuator was observed in the PIV measurements by Boxx et al. [16].
3.6. Boundary Layer Displacement Thickness and Momentum Thickness at Trailing Edge
Figure 15 shows the displacement thickness, , of the boundary layer near the trailing edge. The displacement thickness decreased by 20% at ±2.0 kV. The displacement thickness is suddenly reduced as much as 56% at ±2.2 kV, and it is reduced gradually from ±2.4 kV to ±2.8 kV (77% reduction).
Figure 16 shows the momentum thickness, , of the boundary layer near the trailing edge. The momentum thickness is 2.45 mm at the baseline condition. The momentum thickness slightly increases to 2.61 mm (6.5% rise) at kV, where the shape of boundary layer velocity is a straight line, as shown in Figure 12. The momentum thickness then decreases to 1.95 mm (20.4% reduction) at kV.
3.7. Boundary Layer Shape Factor at Trailing Edge
Figure 17 shows the boundary layer shape factor () near the trailing edge. The shape factor is 6.40 at the baseline condition and reduces as the amplitude of input voltage increases.
An empirical singlevariable correlation using was developed for incipient and full detachment and for reattachment of turbulent boundary layers on twodimensional surfaces [23]. The criteria of the correlation are for separated region (full detachment), for intermittent detachment, and for attached boundary layer (reattachment).
From the correlation, the boundary layer at baseline and kV is separated, that at is intermittently detached, and that at , ±2.6 and ±2.8 kV is reattached.
3.8. Total Pressure Loss Estimation
Figure 18 shows total pressure loss estimation. The total pressure loss coefficient, , was estimated from the boundary layer displacement thickness, , and momentum thickness, , using the following correlation [24]: where is passage width (60 mm), is trailing edge thickness (1.9 mm), and is base pressure coefficient. Because the value is unknown in this study, the typical value of is used for this estimation. The third term including is about of ; therefore, the effect of unknown value is considered to be relatively smaller than the first and second terms.
The total pressure loss coefficient is 0.172 at the baseline condition. The total pressure loss is reduced to 0.107 (38% reduction) at ±2.2 kV, and 0.078 (55% reduction) at ±2.8 kV.
4. Conclusions
Active flow control using dielectric barrier discharge plasma actuators was investigated to reattach the simulated separation flow over the suction surface of a turbine blade at low Reynolds number, Re = 1.7 × 10^{4}. The flow separation was induced on a curved plate installed in the test section of a lowspeed wind tunnel. Particle image velocimetry (PIV) was used to obtain instantaneous and timeaveraged twodimensional velocity measurements. The amplitude of input voltage for the plasma actuator was varied from ±2.0 kV to ±2.8 kV. Due to the flow control using the plasma actuator, the separated flow was successfully reduced. The decrease in displacement thickness of the boundary layer near the trailing edge was 20% at ±2.0 kV. The displacement thickness was suddenly reduced by as much as 56% at ±2.2 kV and it was reduced gradually from ±2.4 kV to ±2.8 kV (77% reduction). The total pressure loss coefficient, estimated from the boundary layer displacement thickness and momentum thickness, was 0.172 at the baseline condition. The total pressure loss coefficient was reduced to 0.107 (38% reduction) at ±2.2 kV, and 0.078 (55% reduction) at ±2.8 kV.
Nomenclature
:  Base pressure coefficient, −0.15 (typical value) 
:  Total pressure loss coefficient 
:  Diffusion factor 
:  Frequency of input voltage 
:  Shape factor () 
:  Streamwise length of curved wall, 100 mm 
:  Total pressure at inlet 
:  Massaveraged total pressure at outlet 
:  Reynolds number 
:  Turbulence intensity 
:  Trailing edge thickness, 1.9 mm 
:  Absolute velocity 
:  Massaveraged absolute velocity at outlet 
:  Maximum absolute velocity 
:  Streamwise velocity 
:  Averaged streamwise velocity of main flow near trailing edge 
:  Random fluctuation (turbulence) component of streamwise velocity 
:  Vertical velocity 
:  Random fluctuation (turbulence) component of vertical velocity 
:  Amplitude of input voltage 
:  Flow passage width at trailing edge, 60 mm 
:  Streamwise (horizontal) distance 
:  Vertical distance 
:  Boundary layer displacement thickness 
:  Boundary layer momentum thickness 
:  Flow angle 
:  Density. 
Acknowledgments
This study was supported by GrantinAid for Young Scientists (B) of KAKENHI (no. 20760124) from The Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan. The authors would like to sincerely thank Dr. Timothy Jukes of the University of Nottingham for his helpful advice.
References
 R. Sondergaard, R. B. Rivir, and J. P. Bons, “Control of lowpressure turbine separation using vortexgenerator jets,” Journal of Propulsion and Power, vol. 18, no. 4, pp. 889–895, 2002. View at: Google Scholar
 H. Arakawa, T. Suzuki, K. Saito, S. Tamura, and S. Kishi, “Research and development of 300 kW class ceramic gas turbine project in Japan,” ASME Paper 97GT87, 1997. View at: Google Scholar
 C. G. Murawski, R. Sondergaard, R. B. Rivir, K. Vafai, T. W. Simon, and R. J. Volino, “Experimental study of the unsteady aerodynamics in a linear cascade with low reynolds number low pressure turbine blades,” ASME Paper 97GT95, 1997. View at: Google Scholar
 W. Lou and J. Hourmouziadis, “Separation bubbles under steady and periodicunsteady main flow conditions,” ASME Journal of Turbomachinery, vol. 122, no. 4, pp. 634–643, 2000. View at: Publisher Site  Google Scholar
 R. J. Howell, O. N. Ramesh, H. P. Hodson, N. W. Harvey, and V. Schulte, “High lift and aftloaded profiles for lowpressure turbines,” ASME Journal of Turbomachinery, vol. 123, no. 2, pp. 181–188, 2001. View at: Publisher Site  Google Scholar
 R. J. Volino and L. S. Hultgren, “Measurements in separated and transitional boundary layers under lowpressure turbine airfoil conditions,” ASME Journal of Turbomachinery, vol. 123, no. 2, pp. 189–197, 2001. View at: Publisher Site  Google Scholar
 K. W. Van Treuren, T. Simon, M. von Koller, A. R. Byerley, J. W. Baughn, and R. Rivir, “Measurements in a turbine cascade flow under ultra low reynolds number conditions,” ASME Journal of Turbomachinery, vol. 124, no. 1, pp. 100–106, 2002. View at: Publisher Site  Google Scholar
 R. B. Rivir, R. Sondergaard, and J. P. Bons, “Control of separation in turbine boundary layers,” AIAA Paper 20042201, 2004. View at: Google Scholar
 J. R. Roth, D. M. Sherman, and S. P. Wilkinson, “Boundary layer flow control with a one atmosphere uniform grow discharge,” AIAA Paper 980328, 1998. View at: Google Scholar
 T. C. Corke, C. L. Enloe, and S. P. Wilkinson, “Dielectric barrier discharge plasma actuators for flow control,” Annual Review of Fluid Mechanics, vol. 42, pp. 505–529, 2010. View at: Publisher Site  Google Scholar
 Y. B. Suzen and P. G. Huang, “Simulations of flow separation control using plasma actuators,” AIAA Paper 2006877, 2006. View at: Google Scholar
 J. Huang, T. C. Corke, and F. O. Thomas, “Plasma actuators for separation control of lowpressure turbine blades,” AIAA Paper 20031027, 2003. View at: Google Scholar
 J. Huang, T. C. Corke, and F. O. Thomas, “Unsteady plasma actuators for separation control of lowpressure turbine blades,” AIAA Journal, vol. 44, no. 7, pp. 1477–1487, 2006. View at: Publisher Site  Google Scholar
 D. P. Rizzetta and M. R. Visbal, “Numerical investigation of plasmabased flow control for a transitional highlyloaded lowpressure turbine,” AIAA Paper 2007938, 2007. View at: Google Scholar
 Y. B. Suzen, P. G. Huang, and D. E. Ashpis, “Numerical simulations of flow separation control in lowpressure turbines using plasma actuators,” AIAA Paper 2007937, 2007. View at: Google Scholar
 I. G. Boxx, J. M. Newcamp, M. E. Franke, N. M. Woods, and R. B. Rivir, “A PIV study of a plasma discharge flowcontrol actuator on a flat plate in an aggressive pressure induced Separation,” ASME Paper GT200691044, 2006. View at: Google Scholar
 D. Burman, T. W. Simon, U. Kortshagen, and D. Ernie, “Separation control using plasma actuators: steady flow in low pressure turbines,” ASME Paper GT201146807, 2011. View at: Google Scholar
 C. Marks, R. Sondergaard, and M. Wolff, “Experimental comparison of DBD plasma actuators for low reynolds number separation control,” ASME Paper GT201145397, 2011. View at: Google Scholar
 T. Matsunuma, “Effects of Reynolds number and freestream turbulence on turbine tip clearance flow,” ASME Journal of Turbomachinery, vol. 128, no. 1, pp. 166–177, 2006. View at: Publisher Site  Google Scholar
 T. Matsunuma, “Unsteady flow field of an axialflow turbine rotor at a low Reynolds number,” ASME Journal of Turbomachinery, vol. 129, no. 2, pp. 360–371, 2007. View at: Publisher Site  Google Scholar
 P. Magnier, D. Hong, A. LeroyChesneau, J. M. Bauchire, and J. Hureau, “Control of separated flows with the ionic wind generated by a DC corona discharge,” Experiments in Fluids, vol. 42, no. 5, pp. 815–825, 2007. View at: Publisher Site  Google Scholar
 T. N. Jukes and K.S. Choi, “Dielectricbarrierdischarge vortex generators: characterisation and optimisation for flow separation control,” Experiments in Fluids, vol. 52, no. 2, pp. 329–345, 2012. View at: Publisher Site  Google Scholar
 S. J. Kline, J. G. Bardina, and R. C. Strawn, “Correlation of the detachment of twodimensional turbulent boundary layers,” AIAA Journal, vol. 21, no. 1, pp. 68–73, 1983. View at: Google Scholar
 J. D. Denton, “Loss mechanisms in turbomachines,” ASME Journal of Turbomachinery, vol. 115, no. 4, pp. 621–656, 1993. View at: Google Scholar
Copyright
Copyright © 2012 Takayuki Matsunuma and Takehiko Segawa. 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.