CENTRIFUGAL OR MIXED FLOW TURBOMACHINERY

Abstract: An impeller (6) in a turbomachinery has blades (3) designed such that reduced static pressure difference DELTA Cp between the hub (2) and the shroud (4) on the suction surface (3c) of the blade (3) shows a remarkably decreasing tendency in the vicinity of the impeller exit (6b) as it approaches the impeller exit (6b) between the impeller inlet (6a) and the impeller exit (6b).

Description: CENTRIFUGAL OR MIXED FLOW TURBOMACHINERY The present invention relates to a turbomachinery including a centrifugal pump or a mixed flow pump for pumping liquid, a blower or a compressor for compression of gas, and more particularly to a turbomachinery having an impeller which has a fluid dynamically improved blade profile for suppressing meridional component of secondary flow. Conventionally, in flow passages of an impeller in a centrifugal or a mixed flow turbomachinery, main flow flowing along the flow passages are affected by secondary flow generated by movement of low energy fluid in boundary layers on wall surfaces due to static pressure gradients in the flow passages. This phenomenon leads to the formation of streamwise vortices or flow having non-uniform velocity in the flow passage, which in turn results in a substantial fluid energy loss not only in the impeller but also in the diffuser or guide vanes downstream of the impeller. The secondary flow is defined as flow which has a velocity component perpendicular to the main flow. The total energy loss caused by the secondary flows is referred to as the secondary flow loss. The low energy fluid accumulated at a certain region in the flow passage may cause a flow separation in a large scale, thus producing positively sloped characteristic curve and hence preventing the stable operation of the turbo achine. There are two known approaches for suppressing the secondary flows in a turbomachine, one of which is to make the impeller have a specific flow passage geometry, the other of which is to supply energy from the outside. As an example of the former approach using a specific flow passage geometry, there is a known method in which blades of the impeller in an axial turbomachine are leaned towards the circumferential direction thereof or the direction of suction or discharge side (L.H.Smith and H.Yeh, "Sweep and Dihedral Effects in Axial Flow Turbomachinery", Trans ASME, Journal of Basic Engineering, Vol.85, No.3, 1963, pp.401- 416), a method in which a radial rotor has a blade curvature in the spanwise direction with a convex blade pressure surface and/or a concave blade suction surface(GB2224083A), or a method in which blades in a turbine cascade are leaned or curved toward a circumferential direction thereof (W.Zhongqi, et al., "An Experimental Investigation into the Reasons of Reducing Secondary Flow Losses by Using Leaned Blades in Rectangular Turbine Cascades with Incidence Angle", ASME Paper 88-GT-4). These methods are known to have a favorable influence upon the secondary flows in the cascade if applied appropriately. However, since the influence of the profile of a blade camber line or a blade cross-section upon the secondary flow has not been essentially known, the effect of blade lean or spanwise blade curvature is utilized under a certain limitation without changing the blade camber line or the blade cross-section substantially. Further, the Japanese laid-open Patent Publication No. 63-10281 discloses a structure in which a projecting portion is provided at the corner of a hub surface and a blade surface in a turbomachine to reduce the secondary flow loss. Since such flow passage profile is a specific blade profile having a nonaxisymmetric hub surface, it is difficult to manufacture the impeller. In all cases of the above prior arts, the method of achieving the effect universally has not been sufficiently studied. Therefore, the universal methods of suppressing the secondary flows under different design conditions and for the different types of turbomachine have not been established. Under these circumstances, there are many cases that the above effect is reduced, or to make matters worse, undesirable effects are obtained. As a result, as of now, there is no standard design criterion for reducing secondary flow by using specific flow passage geometry. Thus the three-dimensional geometry of the impeller has been designed by trial and error to find the optimum profile of the impeller for suppressing the secondary flow. As an example of the latter approach, in which energy is supplied from the outside for suppressing secondary flow, in order to avoid instability of a turbomachine, there has been proposed a method in which secondary flow in a turbine cascade is controlled by blowing fluid into an inlet of the cascade (T.E.Biesinger and D.G.Gregory-Smith, "Reduction in Secondary Flows and Losses in a Turbine Cascade by Upstream Boundary Layer Blowing" ASME Paper 93-GT-114), or a method in which secondary flow in an impeller is controlled by blowing jet into an inlet of the impeller (PCT/JP92/01280). However, the proposed methods are disadvantageous in that an energy source for blowing fluid or jet and incidental devices thereof are required. Further, the methods are also disadvantageous in that regular energy consumption is required to suppress secondary flows. The present invention relates to the former approach, that is the method using a specific flow passage profile. It is apparent from the secondary flow theory that the secondary flow in the impeller results from the action of Corioli's force caused by the rotation of the impeller and the effects of the streamline curvature. The secondary flow in the impeller is divided broadly into two categories, one of which is blade-to-blade secondary flow generated along a shroud surface or a hub surface, the other of which is the meridional component of secondary flow generated along the pressure surface or the suction surface of a blade. It is known that the blade-to-blade secondary flow can be minimized by making the blade profile to be backswept. Regarding the other type of secondary flow, that is, the meridional component of secondary flow, it is difficult to weaken or eliminate it easily. If we wish to weaken or eliminate the meridional component of secondary flow, it is necessary to optimize the three-dimensional geometry of the flow passage very carefully. The purpose of the present invention is to suppress the meridional component of secondary flow in a centrifugal or a mixed flow turbomachine. As an example of a typical impeller in the turbomachinery to which the present invention is applied, the three-dimensional geometry of a closed type impeller is schematically shown in Figs. 1(A) and 1(B) in such a state that most part of shroud surface is removed. Fig. 1(A) is a perspective view partly in section, and Fig. 1(B) is a cross-sectional view taken along a line A-A' which is a meridional cross-sectional view. In Figs. 1(A) and 1(B), a hub surface 2 extends radially outwardly from a rotating shaft 1 so that it has a curved surface similar to a corn surface. A plurality of blades 3 are provided on the hub surface 2 so that they extend radially outward from the rotating shaft 1 and are disposed at equal intervals in the circumferential direction. The blade tip 3a of the blades 3 are covered with a shroud surface 4 as shown in Fig. 1(B). A flow passage is defined by two blades 3 in confrontation with each other, the hub surface 2 and the shroud surface 4 so that fluid flows from an impeller inlet 6a toward an impeller exit 6b. When the impeller 6 is rotated about an axis of the rotating shaft 1 at an angular velocity ω, fluid flowing into the flow passage form the impeller inlet 6a is delivered toward the impeller exit 6b of the impeller 6. In this case, the surface facing the rotational direction is the pressure surface 3b, and the opposite side of the pressure surface 3b is the suction surface 3c. In the case of open type impeller, there is no independent part for forming the shroud surface 4, but a casing (not shown in the drawing) for enclosing the impeller 6 serves as the shroud surface 4. Therefore, there is no basic fluid dynamical difference between the open type impeller and the closed type impeller in terms of the generation and the suppression of the meridional component of secondary flows, thus only the closed type impeller will be described below. The impeller 6 having a plurality of blades 3 is incorporated as a main component, the rotating shaft 1 is coupled to a driving source, thereby jointly constituting a turbomachine. Fluid is introduced into the impeller inlet 6a through a suction pipe, pumped by the impeller 6 and discharged from the impeller exit 6b, and then delivered through a discharge pipe to the outside of the turbomachine. The unsolved serious problem in connection with the impeller of turbomachine is the suppression of the meridional component of secondary flow. The mechanism of generation of meridional component of secondary flow, whose suppression is the purpose of this invention, is explained as follows: As shown in Fig. 1(B), with regard to the relative flow, the reduced static pressure distribution, defined as p * =p-0.5pu 2 , is formed by the action of a centrifugal force W 2 /R due to streamline curvature of the main flow and the action of Corioli's force 2ωW, due to the rotation of the impeller, where W is the relative velocity of flow, R is the radius of streamline curvature, ω is the angular velocity of the impeller, W β is the component in the circumferential direction of W relative to the rotating shaft 1, p * is reduced static pressure, p is static pressure, p is density of fluid, u is peripheral velocity at a certain radius r from the rotating shaft 1. The reduced static pressure p * has such a distribution in which the pressure is high at the hub side and low at the shroud side, so that the pressure gradient balances the centrifugal force W 2 /R and the Corioli's force 2ωw, directed toward the hub side. In the boundary layer along the blade surface, since the relative velocity W is reduced by the influence of the wall surface, the centrifugal force W 2 /R and the Corioli's force 2coW, acting on the fluid in the boundary layer become small. As a result, they cannot balance the reduced static pressure gradient of the main flow, and low energy fluid in the boundary layer flows towards an area of low reduced static pressure p * , thus generating the meridional component of secondary flow. That is, as shown in broken lines on the pressure surface 3b and in solid lines on the suction surface 3c in Fig. 1(A), fluid moves along the blade surface from the hub side towards the shroud side on the pressure surface 3b and the suction surface 3c forming meridional component of secondary flow. The meridional component of secondary flow is generated on both surfaces of the suction surface 3c and the pressure surface 3b. In general, since the boundary layer on the suction surface 3c is thicker than that on the pressure surface 3b, the secondary flow on the suction surface 3c has a greater influence on performance characteristics of turbomachinery. The purpose of the present invention is to suppress the meridional component of secondary flow on the suction surface of the blade. When low energy fluid in the boundary layer moves from the hub side to the shroud side, fluid flow is formed from the shroud side to the hub side at around the midpoint location to compensate for fluid flow rate which has moved. As a result, as shown schematically in Fig. 2(B) which is a cross-sectional view taken along a line B-B' in Fig. 2(A), a pair of vortices which have a different swirl direction from each other are formed in the flow passage between two blades as the flow goes towards exit. These vortices are referred to as secondary vortices. Low energy fluid in the flow passage is accumulated due to these vortices at a certain location of the impeller towards the exit where the reduced static pressure p * is lowest, and this low energy fluid is mixed with fluid which flows steadily in the flow passage, resulting in generation of a great flow loss. Furthermore, when the non-uniform flow generated by insufficient mixing of a low relative velocity (high loss) fluid and a high relative velocity (high loss) fluid is discharged to the downstream flow passage of the blades, a great flow loss is generated when both fluids are mixed. Such a non-uniform flow leaving the impeller makes the velocity triangle unfavorable at the inlet of the diffuser and causes flow separation on diffuser vanes or a reverse flow within a vaneless diffuser, resulting in a substantial decrease of the overall performance of the turbomachine. Furthermore, in the area of high loss fluid accumulated at a certain location in the flow passage, a large scale reverse flow is liable to occur, thus producing a positively sloped characteristics curve. As a result, surging, vibration, noise and the like are generated, and the turbomachinery cannot be stably operated. Therefore, in order to improve the performance of centrifugal or mixed flow turbomachinery and realize stable operation of turbomachinery, it is necessary to design the three- dimensional geometry of the flow passage for suppressing the secondary flow as much as possible. Thereby the formation of secondary vortices, the resulting non-uniform flow, and large scale flow separation or the like will be prevented. It is therefore an object of the present invention to overcome the drawbacks of increase of loss and unstable operation of turbomachinery caused by insufficient suppression of meridional component of secondary flow in the impeller, and to provide the following four design aspects which can reduce the above loss and improve stability of operation of turbomachinery. (1) According to the first aspect of the present invention, there is provided a centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that the reduced static pressure difference ΔCp or the relative Mach number difference ΔM between the hub and the shroud on the suction surface of a blade shows a remarkable decreasing tendency between the location of non-dimensional meridional distance 0 (impeller inlet) and the location of non- dimensional meridional distance 1.0 (impeller exit) as said non- dimensional meridional distance increases. Also, between the impeller inlet and the impeller exit, the reduced static pressure difference ΔCp (or Mach number difference ΔM) becomes nearly 0 or changes sign from positive to negative as it approaches the impeller exit. By designing the distribution of the reduced static pressure difference ΔCp as in the above manner, the meridional component of secondary flow is remarkably suppressed at a location where the reduced static pressure difference ΔCp or relative Mach number difference ΔM shows a remarkably decreasing tendency. As a result, the meridional component of secondary flow over the entire impeller can be effectively suppressed. In this case, if the degree of the remarkably decreasing tendency of the reduced static pressure difference ΔCp and the location where it is generated are optimum, the suppression of meridional component of secondary flow is maximized. The degree of the remarkably decreasing tendency is optimum when it is arranged such that the difference between a minimum value of the reduced static pressure difference ΔCp (or the relative Mach number difference ΔM) and the value of reduced static pressure difference ΔCp (or relative Mach number difference ΔM) at the location, corresponding to non-dimensional meridional distance obtained by subtracting non-dimensional meridional distance 0.4 from non-dimensional meridional distance representing said minimum value, is not less than 0.20 (or 0.15 for relative Much number). The location, where the remarkably decreasing tendency emerges, is optimum when it is placed at non- dimensional meridional distance of not less than m=0.4. Here, non-dimensional meridional distance is defined on the meridional plane of the impeller as shown in Fig. 1(B). At the shroud, the non-dimensional meridional distance m is defined as m=l s /l τ s , which represents the ratio of meridional distance l s , measured from the blade inlet 6a along the shroud, to the meridional distance l τ a , between the impeller inlet 6a and the impeller exit 6b measured along the shroud. Similarly, at the hub, the non- dimensional meridional distance m is defined as m=l H /l τ H , which represents the ratio of meridional distance 1 H , measured from the blade inlet 6a along the hub, to the meridional distance l τ H , between the impeller inlet 6a and the impeller exit 6b measured along the hub. So, m»0 corresponds to the impeller inlet 6a, and m=1.0 the impeller exit 6b. The meridional, component of secondary flow is suppressed when nearly zero value of ΔCp (or ΔM) or the change of sign from positive to negative of ΔCp (or ΔM) occurs at the location of not less than m=0.6. The suppression of meridional component of secondary flow is maximized when it occurs at the location in the range of non-dimensional meridional distance m=0.65-0.9. Incidentally, the reduced static pressure difference ΔCp is used with respect to incompressible fluid as in a liquid pump, while the relative Mach number difference ΔM is used with respect to compressible fluid as in a compressor. The influence of them upon generation of secondary flow is fluid dynamically equivalent to each other, therefore the following description is made using only the reduced static pressure difference ΔCp. (2) According to the second aspect of the present invention, there is provided a centrifugal or mixed flow turbomachinery having an impeller, characterized by the distribution of the meridional derivative of r ~ ϊ , i.e. d ( τ~V t )/Bm along the meridional distance. Here, r is the radial co-ordinate of a cylindrical polar co¬ ordinate system as shown in Fig. 1(B), ~ V t the tangentially averaged value of the circumferential velocity component V, of the absolute velocity V of the fluid. Impellers are designed so that the maximum value of B ( τ~V β ) /Bm at the shroud appears at the fore part of the blade, and the maximum value of d ( τ ~ V t )/θm at the hub appears at the aft part of the blade. Also, impellers are designed so that d ( rV t )/θm is larger at the shroud than that at the hub at the location in the vicinity of non-dimensional meridional distance 0 (impeller inlet side), and the said meridional derivative of r?„ i.e. 3(rV,)/ό>m is smaller at the shroud than that at the hub at the location in the vicinity of non-dimensional meridional distance 1.0 (impeller exit), and the distribution of said meridional derivative of rV, at the hub crosses over that at the shroud between the said location of non-dimensional meridional distance 0 (impeller inlet) and the said location of non-dimensional meridional distance 1.0 (impeller exit). By designing the impeller so as to realize the above distribution of d ( τ ' V t ) /dm, the blade loading can be large at the shroud and small at the hub at the impeller inlet side, and the blade loading can be small at the shroud and large at the hub at the impeller exit side. Consequently, the above is reflected on the first aspect of the present invention, as a relatively large reduced static pressure difference ΔCp is realized at the impeller inlet and a relatively small reduced static pressure difference ΔCp is realized at the impeller exit. As a result, a remarkable decrease of the reduced static pressure difference ΔCp can be realized, thus suppressing the meridional component of secondary flow effectively. In this case, if the location of the said crossing of the distributions of meridional derivative of r ' V, between the shroud and the hub is optimum, the suppression of meridional component of secondary flow is maximized. The favorable location is not less than non-dimensional meridional distance m - 0.5, and the optimum location is in the range of non-dimensional meridional distance m - 0.55-0.85. The location, where B ( τV. ) /3m at the shroud has a maximum value in the fore part of the impeller, is preferable when it is in the range of non-dimensional meridional distance m = 0-0.3, and the location, where 3 ( τV t ) /Bm at the hub has a maximum value in the aft part of the impeller, is preferable when it is in the range of non-dimensional meridional distance m = 0.6-1.0. (3) According to the third aspect of the present invention, the blade is leaned towards a circumferential direction in the whole or a part of the impeller between the inlet and the exit so that the blade at the hub precedes the blade at the shroud in the direction of impeller rotation, and also the blade lean angle γ, does not remain constant and is designed to have a decreasing tendency towards the impeller exit. Here, the blade lean angle γ is defined as an angle between the shroud surface and the blade center line on the cross-section of the flow passage, and Fig.11(C) shows the blade lean angle γ on the cross-section at the impeller exit viewed from E-E' of Fig. 10. In addition to this, on a plan view of the impeller, viewed from the inlet side of the impeller, the value of the blade exit angle o_- E (Fig. 11(C), defined as an angle between the blade camber line at the exit and the tangential direction, decreases from that at the hub α. E H to that at the shroud __- E>s , and also the blade angle __„ (Fig. 11(C), defined as an angle between the blade camber line and the tangential direction, at the hub on the plan view of the impeller is designed to have the maximum or the local maximum value at the exit side of the impeller. Namely, the blade geometry on the hub is modified so as to reduce the degree of backsweep in the latter part of the impeller and increases the blade angle α„ at the hub relative to the blade angle α. st the shroud towards the impeller exit, thus designing the flow passage of the impeller such that the hub side precedes the shroud side relative to the direction of the impeller revolution. By designing the blades having blade lean in the circumferential direction in this way, the reduced static pressure p * at the shroud side increases, whereby the pressure coefficient Cp decreases thereat. On the other hand, the reduced static pressure p * at the hub side decreases, whereby the pressure coefficient Cp increases thereat. The above is reflected on the first aspect of the present invention, and a decreasing tendency of the reduced static pressure difference ΔCp which emerges between the impeller inlet and impeller exit as it approaches the impeller exit is emphasized. In this case, the effects of secondary flow suppression is maximized when the geometry of the blade lean in the circumferential direction is appropriate. Namely, on the plan view of the impeller viewed from the inlet side of the impeller, the location of the reflection point, where the blade profile at the hub has maximum blade angle α-, is not less than the radius radio of r * =(r-r 1H )/(r 2H -r 1H )=0.6 or not less than the non- dimensional meridional distance of m=0.7. The degree of blade backsweep towards the impeller exit is higher at the shroud than at the hub (Fig. 11(C)), and the following relation is satisfied; {d(rθ)/dr}s > {d(rθ)/dr} H . Here, r and θ are the radial and angular co-ordinates of a cylindrical polar co-ordinate system respectively. The angular co-ordinate θ is defined as positive when it is measured in the counter-rotating direction of the impeller. A blade lean ratio S is defined by S=(Θ TE S -Θ TE H )/(Θ TE M - Θ E . M ) as a parameter which shows overall degree of the blade lean, and the optimum blade lean ratio is not less than S=0.14. Here, Θ TE s represents the angular co-ordinate of the impeller exit at the shroud, Θ TE H that of the impeller exit at the hub, Θ TE „ that of the impeller exit at the midspan location, i.e. the middle point between the shroud and the hub, Θ LE M that of the impeller inlet at the said midspan location (Fig. 11(C)). When the secondary flow to be suppressed are strong, and only the third aspect of the present invention is used, in order to create a sufficiently decreasing tendency of the reduced static pressure difference ΔCp between the shroud and the hub to suppress the meridional component of secondary flow, the degree of blade lean in the circumferential direction is intensified to such an extent which may lead to difficulties in manufacturing. In this case, the third aspect of the present invention has an effect for suppressing meridional component of secondary flow effectively in cooperation with the second aspect without leading to difficulties in manufacturing. (4) According to the fourth aspect of the present invention, there is provided a centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that blade angle distribution measured from the circumferential direction is arranged such that blade angle difference β H - β 3 , between blade angle β H at the hub and blade angle β 3 at the shroud, shows a remarkably increasing tendency from the location of non-dimensional meridional distance 0 (impeller inlet) to the location of non-dimensional meridional distance 1.0 (impeller exit) as said non-dimensional meridional distance increases. Here, the blade angle is defined on a blade-to-blade stream surface, viewed normal to the curved hub surface or on the shroud surface, and measured form the tangential direction. Namely, the blade angle β H represents the angle between the blade camber line at the hub and the tangential direction at radius r, and the blade angle β s the angle between the blade camber line at the shroud and the tangential direction at radius r (Figs. 12(B) and 12(C)). By designing the distribution of the blade angle difference P H ~ β s as n "the above manner, the blade angle β H becomes larger at the hub as it approaches the impeller exit, and the blade loading can be larger at the hub towards the impeller exit. On the other hand, the blade angle β s becomes smaller at the shroud as it approaches the impeller exit, and the blade loading can be smaller at the shroud towards the impeller exit. Consequently, the above is reflected in the first aspect of the present invention, and a relatively small reduced static pressure difference ΔCp is realized at the impeller exit. That is, a remarkable decrease of the reduced static pressure difference ΔCp between the impeller inlet and the impeller exit enables the meridional component of secondary flow to be suppressed effectively. In this case, if the maximum or the local maximum value of the blade angle difference β H - β s and the location where its maximum or a local maximum value emerges are optimum, the effects of the suppression of meridional component of secondary flow is maximized. The maximum value or the local maximum value of not less than 20 degrees is optimum, and the difference between the maximum or the local maximum blade angle difference and the mean blade angle difference, averaged in the range of the meridional distance m=0-0.2, is preferably not less than 10 degrees. The location where the maximum value emerges is favorable at the location not less than the non-dimensional meridional distance of m=0.5 and is optimum in the range of non-dimensional meridional distance m=0.7-1.0. There is a case where the distribution of blade angle difference β H - β s shows a decreasing tendency in a part of the non-dimensional distance between m=0-1.0 followed by a remarkably increasing tendency as the non-dimensional meridional distance increases. However, the effects of such remarkable increase of the blade angle difference on the suppression of secondary flow are basically the same. In the accompanying drawings: Figs. 1 and 2 are views for explaining the background art; Figs. 1(A) and 1(B) are views for explaining meridional component of secondary flow in three-dimensional geometry of a closed type impeller, Fig. 1(A) is a perspective view partly in section, and Fig. 1(B) is a meridional cross-sectional view taken along a line A-A' of Fig. 1(A); Figs. 2(A) and 2(B) are" views for explaining secondary vortices caused by meridional component of secondary flow in the closed type impeller, Fig. 2(A) is a perspective view partly in section, and Fig. 2(B) is a cross-sectional view taken along a line B-B' of Fig. 2(A); Figs. 3 through 14 are views showing best mode for carrying out the present invention; Figs. 3 through 6 are graphs showing the relationship between the non-dimensional meridional distance m and the pressure coefficient Cp; Figs. 7(A), 7(B) and 8 are graphs showing the relationship between the non-dimensional meridional distance m and 9(r " V β )/r3m, θ(rT β )/θm being meridional derivative of τ " V t (angular momentum) which is product of radius r in the flow passage and tangentially averaged value of the circumferential velocity component V β of absolute velocity V at the radius r; Figs. 9 through 11 are views for explaining blade lean in a circumferential direction; Fig. 9(A) is a perspective view of an impeller, partly in section, and Fig. 9(B) is a cross-sectional view taken along a line C-C of Fig. 9(A); Fig. 10 is a meridional cross-sectional view of the impeller; Lower and upper figures of Figs. 11(A), 11(B) and 11(C) are views as viewed from D-D' and E-E' of Fig. 10 respectively, Fig. 11(A) is a view showing the impeller designed by a conventional design method, Fig.11(B) is a view showing the impeller having lean blade designed by a conventional design method and Fig. 11(C) is a view showing the impeller having lean blade designed by the present invention; Figs. 12(A), 12(B) and 12(C) are views for explaining blade angle defined on the blade-to-blade stream surfaces, Fig. 12(A) is a perspective view of the impeller, Fig. 12(B) is a view of blade-to-blade stream surface along the hub 2 taken from F direction of Fig. 12(A), and Fig. 12(C) is a view of blade-to- blade stream surface along the shroud 4 taken from G direction of Fig. 12(A); Figs. 13(A) and 13(B) are graphs showing the relationship between the non-dimensional meridional distance m and blade angle, Fig. 13(A) shows blade angle β H curves at the hub, and Fig. 13(B) shows blade angle β s curves at the shroud; Fig. 14 is a graph showing curves of blade angle difference P H ~ β s between blade angle β H at the hub and blade angle β s at the shroud; Figs. 15 through 74 show verification data in the best mode for carrying out the present invention; Figs. 15 through 26 show those for impellers of a centrifugal pump with a low specific speed, Figs. 15, 18, 21 and 24 show the pressure coefficient Cp curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 16, 19, 22 and 25 show blade loading θ(rV,)/θm curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 17, 20, 23 and 26 show secondary flow vector diagrams, Figs. 15 through 17 show the case designed by the conventional design method, Figs. 18 through 20 show the case designed by using only blade lean in a circumferential direction according to the present invention, Figs. 21 through 23 show the case designed by combination of blade lean and blade loading d ( τ " v i t ) /dm according to the present invention, and Figs. 24 through 26 show the case designed by using only blade loading d ( τV t ) /dm according to the present invention; Figs. 27 through 38 show verification data for impellers of a mixed flow pump with a medium specific speed, Figs. 27, 30, 33 and 36 show the pressure coefficient Cp curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 28, 31, 34 and 37 show blade loading d(rT )/dm curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 29, 32, 35 and 38 show secondary flow velocity vector diagrams, Figs. 27 through 29 show the case designed by the conventional design method, Figs. 30 through 32 show the case designed by using only blade lean in a circumferential direction according to the present invention, Figs. 33 through 35 show the case designed by combination of blade lean and blade loading d ( τ " V. ) /dm according to the present invention, and Figs. 36 through 38 show the case designed by using only blade loading d ( τ ~ V t ) /dm according to the present invention; Figs. 39 through 50 are verification data for impellers of a centrifugal compressor, Figs. 39, 42, 45 and 48 show the relative Mach number M curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 40, 43, 46 and 49 show blade loading θ(r " V,)/θm curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 41, 44, 47 and 50 show secondary flow vector diagrams, Figs. 39 through 41 show the case designed by the conventional design method, Figs. 42 through 44 show the case designed by using only blade lean in a circumferential direction according to the present invention, Figs. 45 through 47 show the case designed by combination of blade lean and blade loading θ(r " V β )/dm according to the present invention, and Figs. 48 through 50 show the case designed by using only blade loading θ(r " V.)/c»π. according to the present invention; Figs. 51 through 62 are verification data for impellers of a mixed flow compressor, Figs. 51, 54, 57 and 60 show the relative Mach number M curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 52, 55, 58 and 61 show blade loading θ(r " V β )/r3m curves (solid lines represent curves at the hub and broken lines represent curves at the shroud), Figs. 53, 56, 59 and 62 show secondary flow vector diagrams, Figs. 51 through 53 show the case designed by the conventional design method, Figs. 54 through 56 show the case designed by using only blade lean in a circumferential direction according to the present invention, Figs. 57 through 59 show the case designed by combination of blade lean and blade loading _>(r " V β )/θm according to the present invention, and Figs. 60 through 62 show the case designed by using only blade loading θ(r " V β )/r3m according to the present invention; Figs. 63 through 74 are graphs showing the relationship between meridional distance and blade angle, comparing the case designed by the conventional design method and the case designed by using only blade lean in a circumferential direction according to the present invention, or combination of blade lean in a circumferential direction and blade loading d ( τV t ) /dm according to the present invention, or only blade loading d ( τ ~ V t )/dm according to the present invention; Figs. 63 through 65 are those for impellers of a centrifugal pump with a low specific speed, Fig. 63 shows blade angle β H curves at the hub, Fig. 64 shows blade angle β s curves at the shroud, and Fig. 65 shows blade angle difference β H - β 3 curves; Figs. 66 through 68 are those for impellers of a mixed flow pump with a medium specific speed, Fig. 66 shows blade angle β H curves at the hub, Fig. 67 shows blade angle β s curves at the shroud, Fig. 68 shows blade angle difference β H - β g curves; Figs. 69 through 71 are those for impellers of a centrifugal compressor, Fig. 69 shows blade angle β H curves at the hub, Fig. 70 shows blade angle β s curves at the shroud, and Fig. 71 shows blade angle difference β H - β s curves; Figs. 72 through 74 are those for impellers of a mixed flow compressor, Fig. 72 shows blade angle β H curves at the hub, Fig. 73 shows blade angle β 3 curves at the shroud, and Fig. 74 shows blade angle difference β H - β 3 curves; Figs. 75(A) and 75(B) are perspective views showing the difference of the geometry of blade lean in a mixed flow pump impeller, Fig. 75(A) shows the case designed by using only blade lean in a circumferential direction according to the present invention, and Fig. 75(B) shows combination of blade lean and blade loading d ( τV. ) /dm according to the present invention; Fig. 76 is a flow chart showing the flow of numerical calculation in the three-dimensional inverse design method used to acquire the verification data for the turbomachinery impellers of the present invention. An embodiment according to the first aspect of the present invention will be described below. The influence of viscosity can be neglected for main flow of the relative flow in the flow passages of an impeller, therefore the following formula is approximately satisfied in incompressible flow as in a liquid pump. p τrei = P * + °« 5 P W 2 ~ constant where P Trβl is relative stagnation pressure at the upstream of the impeller. Next, as non-dimensional quantity of reduced static pressure p * on the blade surface, pressure coefficient Cp is defined by the following equation: Cp = (P Tr _ι " P * ) / (0.5 p U 2 ) = (W / U) 2 where U represents the peripheral speed at the impeller exit. Fig. 3 shows the relationship between the non-dimensional meridional distance m and the pressure coefficient Cp, at the hub and the shroud of the blade. As is apparent from the above equation, the pressure coefficient Cp is large at the shroud where reduced static pressure p * is low, and is small at the hub where reduced static pressure p * is high. As mentioned above, since the meridional component of secondary flow on the blade suction surface is directed to the shroud side having low reduced static pressure p * from the hub side having high reduced static pressure p * , suppression of meridional component of secondary flow can be expected by reducing pressure difference ΔCp between them. Incidentally, in case of incompressible fluid, the pressure coefficient Cp is equal to (W / U) 2 , where W is relative velocity. In compressible fluid as in a compressor, the physical variable being related to the behavior of secondary flow is relative Mach number. In order to simplify the description, only the distribution of the pressure coefficient Cp will be described below. The influence of distribution of the pressure coefficient Cp in incompressible flow upon meridional component of secondary flow is equivalent to that of the relative Mach number M in compressible flow. Since the boundary layers on the blade surfaces which develop along the wall of the flow passage in the impeller increase their thickness cumulatively from the impeller inlet toward the impeller exit, the present invention proposes the structure for suppressing meridional component of secondary flow on the suction surface of the blade, considering distribution of the pressure coefficient Cp in the latter half of the impeller. That is, the blade profile is designed so as to have the pressure distribution as shown in Fig. 4 in order that the pressure difference ΔCp between the hub side and the shroud side on the suction surface shows a remarkably decreasing tendency toward the impeller exit. The pressure distribution on the blade in which the pressure difference ΔCp decreases remarkably toward the impeller exit is achieved by the following: (a) increasing the blade loading at the hub side as shown in Fig. 5, i.e. the pressure difference between the pressure surface and the suction surface of the blade at the hub side is remarkably increased towards the impeller exit. (b) decreasing the blade loading remarkably at the shroud side towards the impeller exit as shown in Fig. 6. (c) combining above described (a) and (b). In any case, it is essential to have the pressure distribution on blades such that the pressure difference ΔCp decreases remarkably towards the impeller exit on the suction surface of the blade as shown in Fig. 4. According to a number of verification data, the degree of above mentioned remarkably decreasing tendency is optimum when it is arranged such that the difference between a minimum value of reduced static pressure difference ΔCp (or relative Mach number difference ΔM) and a value of reduced static pressure difference ΔCp (or relative Mach number difference ΔM) at the location, corresponding to non-dimensional meridional distance obtained by subtracting non-dimensional meridional distance 0.4 from that representing the minimum value, is not less than 0.2 (0.15 in the case of relative Mach number difference ΔM). Also the verification data show that the optimum location where remarkably decreasing tendency emerges is in the range of non- dimensional meridional distance m=0.4-1.0. In a centrifugal impeller or a mixed flow impeller having a low specific speed, width of the impeller exit is narrow and the pressure difference between the hub surface and the shroud surface is small. Therefore the pressure difference ΔCp is usually small at the impeller exit as in Fig. 4, i.e. at the location of the non-dimensional meridional distance m-1.0. On the other hand, at the impeller inlet, the pressure coefficient Cp is larger at the shroud than that at the hub. Accordingly, the pressure difference ΔCp on the suction surface of the blade in the centrifugal impeller has a tendency to decrease towards the impeller exit in the similar manner as in Fig. 4. However, in order to suppress secondary flow in such impellers, it is necessary that the pressure difference ΔCp becomes almost zero, or the pressure difference ΔCp becomes negative, that is, the decreasing tendency of the pressure difference ΔCp is enhanced so that the pressure coefficient Cp at the hub becomes larger than that at the shroud. According to a number of verification data, the location where ΔCp becomes nearly zero or changes its sign from positive to negative is favorable when it is in the location not less than non-dimensional meridional distance of m=0.6, and is optimum when the location is in the range of non- dimensional meridional distance m=0.65-0.9. In the conventional design, since such specific consideration is not made, the meridional component of secondary flow cannot be suppressed effectively. Next, an embodiment according to the second aspect of the present invention will be described below. In order to realize the above distribution of the reduced static pressure p * , the present invention proposes a structure in which the blade profile is determined based on the meridional derivative in the direction of non-dimensional meridional distance m of the circumferentially averaged angular momentum r " V, of fluid given by the work of the impeller. In general, the blade loading of the impeller in turbomachinery, i.e. the pressure difference p * (+) - p * (-) between the pressure surface and the suction surface is expressed by the following equation. P * (+)/p - P * (-)/P - 2π(W bl VrV%)/B where W bl is relative velocity at the blade surface location, B is the number of blades, V is a differential multiplier which gives a primary derivative in space co-ordinate, and V, represents the tangentially averaged value of circumferential velocity component V β of the absolute velocity V as shown for example in the blade-to-blade plane at the hub in Fig. 12(c). From the above equation, it is acknowledged that the pressure difference between the pressure surface and the suction surface is closely related to the meridional derivative of τ ~ V t (angular momentum) which is product of radius r in the flow passage and the tangentially averaged value of circumferential velocity component V β of the absolute velocity at the radius. Therefore, as shown in Figs. 5 and 6, in order to increase or decrease the pressure coefficient Cp, the meridional derivative of τ ~ V. should be set to a proper value. In the present invention, the meridional derivative of r ~ V. is given as follows: The distribution of d ( τV t ) /dm at the hub is given such that d ( τ " V t ) /dm has a small value near the impeller inlet and a maximum value near the impeller exit as shown in Fig. 7(A), whereby the blade loading p * (+) - p * (-) on the hub surface is increased in the aft part of the impeller and is reduced in the fore part of the blade. Next, distribution on the shroud surface is given such that d ( r ~ V t )/θm has a maximum value near the impeller inlet and a low value near the impeller exit as shown in Fig. 7(B), whereby the blade loading p * (+) - p * (-) on the shroud surface is increased in the fore part of the impeller and decreased in the aft part of the blade. In Fig. 8, by combining d ( τ ~ V. ) /dm distributions on the hub surface and the shroud surface, the meridional derivative of r " V, in the direction of the meridional distance can be higher at the shroud than that at the hub near the impeller inlet, and it can be higher at the hub than that at the shroud near the impeller exit. In this case, the distribution of θ(r " V,/_»m) at the hub crosses over that at the shroud at the location between the non-dimensional meridional distance of m=0-1.0. By using this type of r3(r " V,)/θm distribution, it is possible to enhance substantially the decreasing tendency of the reduced static pressure difference ΔCp towards the impeller exit in the meridional direction, because this type of distribution results in a relatively large value of ΔCp near the inlet of the impeller and relatively small value of ΔCp near the exit of the impeller. Thereby the meridional component of secondary flow on the suction surface of blade can be effectively suppressed. Many verification data prove that the location at which both curves cross is preferably not less than non-dimensional meridional distance of m=0.5 and is optimum at in the range of non-dimensional meridional distance m=0.55-0.85. The location, where d(r )/θm at the shroud has a maximum value in the fore part of the impeller, is preferable when it is in the range of non-dimensional meridional distance m - 0-0.3, and the location, where θ(r " V„)/3m at the hub has a maximum value in the aft part of the impeller, is preferable when it is in the range of non-dimensional meridional distance m ■ 0.6-1.0. If a proper distribution of loading distribution 3(rV,)/θm is determined, the three-dimensional geometry of the impeller which realize such distribution can be determined by trial and error modification of the blade profile using a flow analysis method. However, more efficiently, such three-dimensional geometry can also be determined by using three-dimensional inverse design approach, which is widely known as a method to calculate numerically the three-dimensional blade profile which will realize specified blade loading distribution. Different theories of three-dimensional inverse design method have been proposed; for example, Zangeneh, M. , 1991, "A Compressible Three Dimensional Blade Design Method for Radial and Mixed Flow Turbomachinery Blades", International Journal of Numerical Methods in Fluids, Vol.13, pp.599-624., Borges, J.E., 1990, "A Three-Dimensional Inverse Method for Turbomachinery: Part I - Theory", Transaction of the ASME, Journal of Turbomachinery, Vol.112, pp.346-354, Yang, Y.L., Tan, C.S. and Hawthorne, W.R., 1992, "Aerodynamic Design of Turbomachinery Blading in Three- Dimensional Flow: An Application to Radial Inflow Turbines", ASME Paper 92-GT-74, Dang, T.Q., 1993, "A Fully Three-Dimensional Inverse Method for Turbomachinery Blading in Transonic Flows", Transactions of the ASME, Journal of Turbomachinery, Vol.115, pp.354-361. The method, proposed by M. Zangeneh (1991) and used in the present invention, is briefly explained as follows. Blades are represented by sheets of vorticity which is distributed along the blade camber line, and the strength of the vorticity is determined by a specified distribution of τ ~ V t which is the product of tangentially averaged value of circumferential velocity V β and radius r. V t can be easily determined from the blade loading distribution θ(r " V,)/θm and the meridional geometry by integration of θ(r " V β )/θm along the non-dimensional distance m. The blade profile for realizing the blade loading distribution d ( r ~ V t ) /dm specified in the present invention can be determined by using the condition that velocity induced by sheets of vorticity in the flow field must be aligned to the blade surfaces. The flow of numerical calculation in the three- dimensional inverse method is shown in Fig. 76 as a flow chart and can be briefly explained as follows. (Step 1) Calculate τV t distribution on the meridional plane by integrating the specified loading distribution d ( τ~V t ) /dm along the non-dimensional meridional distance m. (Step 2) Estimate meridional component of the velocity Vm based on the specified design mass flow rate of the impeller and the cross-sectional area of the specified meridional geometry. (Step 3) Estimate tangential component of the velocity V, by dividing τ ~ V t by the radius r at each location within the impeller. At this stage, the velocity field obtained is an initial approximation of the actual velocity field as the effects of blade shape are not included accurately. (Step 4) Calculate the blade shape at each location by using the velocity field and the inviscid slip condition which implies that the relative flow must be aligned to the blade surfaces. The blade shape is calculated by solving a first order hyperbolic partial differential equation, and so it is possible to give blade lean at the impeller exit by specifying the θ distribution along the blade exit as initial value for integration. (Step 5) Calculate the velocity distribution throughout the blade region based on the calculated τ ~ V. distribution obtained at Step 1 and the blade shape at this stage obtained at Step 4. (Step 6) Calculate the velocity field by solving the governing partial differential equations using the velocity distribution obtained at Step 5 and the calculated blade shape obtained at Step 4. (Step 7) Go back to Step 4 with revised velocity field obtained at Step 6, and repeat the process from Step 4 to Step 7 until changes in blade shapes between iterations fall below a certain given tolerance. Next, an embodiment according to the third aspect of the present invention will be described below. According to the present invention, in order to optimize the distribution of the reduced static pressure p * in the impeller, the blades are leaned in the circumferential direction in such a way that the blade at the hub precedes that at the shroud in the rotating direction of the impeller. Fig. 9(A) shows schematically the perspective view of the three-dimensional geometry of a closed type impeller partly in section, similar to Fig. 1(A). Fig. 9(B) shows a cross-sectional view taken along a line C-C in Fig. 9(A). Here, the symbols in Figs. 9(A) and 9(B) represent the same parts shown in Fig. 1(A) having the same symbols. As shown schematically in Fig. 9(B), due to the pressure difference between the blade pressure surface 3b and the blade suction surface 3c, a fluid force acts on the blade nearly normal to the blade surface, and a reaction force acts on the fluid in the opposite direction towards the blade suction surface 3c nearly normal to the blade 3. When there is no blade lean in the circumferential direction as shown with broken lines in Fig. 9(B), the force acts on the fluid in the circumferential direction and has no force component normal to the hub or the shroud. On the other hand, if the blade 3 is inclined such that the blade at the hub 3CH precedes that at the shroud 3CS in the rotating direction of the impeller, a force having a component normal to the shroud surface acts on the fluid as reaction force directed from the pressure side to the suction side of the flow passage as shown in Fig. 9(B). As a result, a reduced static pressure field is produced in the flow passage to balance the above described blade force acting on fluid, and the reduced static pressure p * has a higher value at the shroud surface and a lower value at the hub surface as compared to the case with no blade lean. Therefore, the pressure difference ΔCp between the hub surface and the shroud surface of pressure coefficient Cp distribution in Fig. 3 decreases, thus the pressure field suitable for secondary flow suppression can be formed as shown in Fig. 4. The influence of blade lean in the circumferential direction upon the meridional component of secondary flow has been studied qualitatively by conventional reseaches on axial flow turbomachines. However, since the influence of the blade geometry on the secondary flows has not been essentially and quantitatively understood, effects of blade lean have been utilized under a certain limitation in the conventional blade lean by simply shifting the shroud side of the blade in the circumferential direction keeping the blade lean angle γ the same between the blade inlet and the blade exit, or by stacking conventional blade profiles along a curved spanwise stacking line. Fig. 10 schematically shows the meridional cross-section of a closed impeller, and the plan view of the impeller viewed in the D-D' direction of Fig. 10 from the impeller inlet 6a and the plan view of the impeller outlet viewed in the E-E' direction of Fig. 10 from the impeller exit 6b are presented together in the lower and upper figures in Figs. 11(A), 11(B) and 11(C) respectively. When a conventional blade lean, disclosed for example in the Japanese laid-open Patent Publication No.55- 134798, is applied to a conventional impeller shown in Fig. 11(A), the impeller having the plan view of Fig. 11(B) is obtained. In this conventional blade lean, the blade in the shroud side is simply shifted in the circumferential direction without modifying blade geometry at the shroud. In this case, the location of reflection point q on the blade camber line at the hub 2, where the blade curvature changes its sign, is the same with that of the conventional impeller shown in Fig. 11(A). Also the upper figure of Fig. 11(B) shows that the blade lean angle γ, defined as the angle between the center line of the blade 3 and the shroud surface 4, remains nearly constant near the impeller exit. On the other hand, in the third aspect of the present invention, the blade geometry at the hub is modified more actively by shifting the location of the reflection point q radially outward and thus reducing the degree of backsweep at the latter part of the blade. Namely, the blade angle d H at the hub is increased toward the impeller exit, and the blade geometry of the flow passage is designed so that the hub side precede the shroud side in the direction of the impeller revolution. In this case, as shown in Fig. 11(C), the reflection point q of the blade camber line at the hub is moved toward the exit side of the impeller compared to the case of conventional blade lean shown in Fig. 11(B). Also, impellers are designed such that the degree of the backsweep, i.e. {d(rθ)/dr}, towards impeller exit becomes higher at the shroud side, giving {d(rθ)/dr} 3 > {d(rθ)/dr} H . Where, suffix "S" represents a value at the shroud and "H" at the hub. Also, it is favorable to design the blade lean angle γ such that it shows decreasing tendency towards impeller exit as shown in Fig. 11(C). Here, it is necessary not to keep the value of the blade lean angle γ constant. In this case, the pressure difference ΔCp is kept relatively high near the inlet because of the small effects of the blade lean having relatively large value of the blade lean angle γ, and the pressure difference ΔCp is reduced towards exit because of the large effects of the blade lean having small value of the blade lean angle γ towards exit of the impeller, thus remarkably reducing tendency of ΔCp is obtained. In addition to these, as shown in Fig. 11(C), it is preferable to design the impeller such that the blade exit angle α. E , defined on the plan view of the impeller viewed from the impeller inlet side, has decreasing tendency from the angle of α- E H at the hub 2 to the angle of a.. 3 at the shroud 4, and the blade angle α„ at the hub has a maximum or a local maximum value at the exit side of the impeller. The effects of blade lean become even more effective by giving such special consideration on the blade geometry in the hub side in addition to the conventional blade lean created by moving the shroud side in the opposite direction of the impeller revolution. Many verification data prove that the favorable radius position of reflection point, where the blade angle a^ becomes maximum and the blade profile at the hub changes the sign of its curvature, is located at the location not less than the radius ratio of r * =(r-r 1H )/(r 2H -r 1H )=0.6 or after the non- dimensional meridional distance of m=0.7. Where, r 1H represents the impeller inlet radius at the hub and r 2H the impeller exit radius at the hub as shown in Fig. 10. Also it was proved that the effects of blade lean become most effective when the blade lean ratio S is more than S*(Θ TE 3 -θ TE H )/(Θ TE M -Θ LE „) -0.14. Where, θ represents the circumferential co-ordinate of a cylindrical polar co-ordinate system as shown in Fig. 11(C), suffix "TE" the value at the blade exit (blade trailing edge), "LE" the value at the blade inlet (blade leading edge) , "S" the value at the shroud 4, "H" the value on the hub 2, and "M" the value at the midspan location. However, in order to completely suppress meridional component of secondary flow in a centrifugal or mixed flow turbomachinery, which shows much stronger secondary flow than the case of an axial turbomachinery, the required blade lean in a circumferential direction becomes excessive, and the use of such lean becomes unpractical because of difficulties in manufacturing, and structural and strength problems. On the other hand, according to the present invention, a method is proposed to suppress the meridional component of secondary flow by using the second aspect (optimization of the blade loading distribution θ(rl β )/dm) together with the blade lean without manufacturing difficulties, structural and strength problems. Next, an embodiment according to the fourth aspect of the present invention will be described below. With regard to optimization of d ( τ ~ V t )/dm distribution, or blade lean in a circumferential direction in which the blade is inclined towards a circumferential direction, or combination of both of them according to the present invention, the three- dimensional blade profile is calculated by the three-dimensional inverse design method following the flow chart in Fig. 76. And it is proved that the impellers in which the meridional component of secondary flow is effectively suppressed have characteristic blade angle distributions, i.e. the fourth aspect of the present invention is to design the impeller which has the following blade angle distribution. (a) Blade angle distribution with respect to non-dimensional meridional distance at the hub side increases remarkably towards the impeller exit as shown in Fig. 13(A). (b) Blade angle distribution with respect to non-dimensional meridional distance at the shroud surface side decreases remarkably towards the impeller exit as shown in Fig. 13(B). (c) Blade angle difference between the hub and the shroud in the blade angle distribution with respect to non-dimensional meridional distance increases remarkably towards the impeller exit as shown in Fig. 14. In this case, when the location of maximum value of blade angle difference is not less than the non-dimensional meridional distance of m=0.5 or most effectively in the range of non- dimensional meridional distance 0.7-1.0 and the maximum blade angle difference is not less than 20 degrees, it is ascertained in many verification data that the effect of suppression of meridional component of secondary flow is remarkable. It is also ascertained that the maximum value of blade angle difference is larger than the average blade angle difference in the range of non-dimensional distance m=0.0-0.2 by more than 10 degrees. The definition of above described blade angle is presented in Figs. 12(A) and 12(B). Fig. 12(A) shows a perspective view of a closed type impeller partly in section, which corresponds to that in Fig. 1(A). Fig. 12(B), viewed in the F direction in Fig. 12(A), shows the blade profile on hub surface 2 viewed normal to the hub 2 along the curved hub surface 2. The blade angle β H at the hub 2 is defined as an angle between the camber line of the blade 3 and the tangential direction at a radial distance of r from the axis 1 of rotation. Fig. 12(C), viewed in the G direction in Fig. 12(A), shows the blade profile on shroud surface 4 viewed normal to the shroud 4 along the curved shroud surface 4. The blade angle β 3 at the shroud 4 is defined as an angle between the camber line of the blade 3 and the tangential direction at a radial distance r from the axis 1 of the rotation. The blade angle difference β H -β 3 is defined as the angle difference between the blade angle β H at the hub 2 and the blade angle β 3 at the shroud 4. The physical meaning of the characteristic distribution in blade angles can be understood considering the basic structure for suppressing meridional component of secondary flow as shown in Fig. 4. As the blade angle (defined from a circumferential direction) becomes large, fluid flow discharged from the impeller has large swirl velocity, thus giving a large impeller work on the fluid flow. At this time, the pressure difference between the suction surface and the pressure surface, generating the impeller work, also increases, and hence the blade loading increases. That is, as the blade angle on the hub surface increases towards the impeller exit as shown in Fig. 13(A), the corresponding blade loading (the area surrounded by a solid and a broken lines in Fig. 5) on the hub surface increases towards the impeller exit as shown in Fig. 5. On the other hand, as the blade angle on the shroud surface decreases toward the impeller exit as shown in Fig. 13(B), the corresponding blade loading on the shroud surface decreases toward the impeller exit as shown in Fig. 6. As a result, it can be understood that by designing the impeller to have such a blade angle distribution, the pressure distribution on blades for effectively suppressing meridional component of secondary flow can be realized. As is apparent from the description of Fig. 4, since the effect of suppression of meridional component of secondary flow is determined by the difference of reduced static pressure distribution between the hub and the shroud on the blade suction surface, the blade angle difference between the hub and the shroud is physically important. So, the feature of the blade profile of impeller which is advantageous to suppress meridional component of secondary flow is reflected to the distribution of the blade angle difference shown in Fig. 14. There is a case where the distribution of blade angle difference β„-β s shows a decreasing tendency in a part of the non- dimensional distance between m=0-1.0 followed by a remarkably increasing tendency as the non-dimensional meridional distance increases. However, the effects of such remarkable increase of the blade angle difference on the suppression of secondary flow are basically the same. The above characteristic feature of the blade angle distribution is reflected on the third aspect of the present invention,i.e. a characteristic feature of the blade lean in the circumferential direction. The location of the reflection point, where the hub blade profile changes the sign of its curvature as shown in Fig. 11(c), gives characteristic feature of the third aspect of the present invention, and is given as the location of minimum (rθ) change in the radial direction, i.e. the location where d(rθ)/dr becomes minimum. On the other hand, the blade angle β H defined in Fig. 12(B) is given by tanβ H ={dm/d(rθ)} H . Here, tanα H ={dr/d(rθ)} H ={dra/d(rθ)} H (dr/dm) H =tanβ H (dr/dm) H and the value of (dr/dm) H is calculated once the meridional geometry is specified. Consequently, the location of the minimum {d(rθ)/dr} H is closely correlated with the location of the maximum β H . Because of this, it is easily acknowledged that the most effective location of such point of reflection on the hub, which is the characteristic feature of the third aspect of the present invention, is not less than the non-dimensional radius ratio of r * =0.6 or not less than the non-dimensional distance of m=0.7. Next, verification data for ascertaining the above embodiments will be described below. Figs. 15 through 26 and Figs. 27 through 38 show verification data for pump impellers, Figs. 15 through 26 show those for a centrifugal impeller with a low specific speed and Figs. 27 through 38 are those for a mixed flow impeller with a medium specific speed. Figs. 39 through 50 and Figs. 51 through 62 show verification data for compressor impellers, Figs. 39 through 50 are those for a centrifugal impeller and Figs. 51 through 62 are those for a mixed flow impeller. In order to ascertain the effect of suppression of meridional component of secondary flow, the result of numerical analysis by a well established three-dimensional viscous calculation method (Dawes, W.N.,1988,"Development of a 3D Navier-Stores Solver for Application to all Types of Turbomachinery" ASME Paper No. 88-GT- 70.) was used. Figs. 17, 20, 23, 26, 29, 32, 35, 38, 41, 44, 47, 50, 53, 56, 59, 62 show flow velocity vector in the vicinity of the suction surface of the blade predicted by numerical analysis and the flow pattern of secondary flow in the boundary layer along the suction surface. The invention will now be more particularly described based on the verification data of Figs. 15 through 26 for centrifugal pump impellers. In the design of conventional impeller, there is no consideration on the suppression of secondary flows in the impeller, and the meridional distribution of the pressure coefficients Cp towards such impeller exit shows larger pressure coefficient difference of ΔCp between the hub and the shroud on the blade suction surface, showing insufficient decrease of ΔCp to suppress secondary flows as shown in Fig. 15. The distribution of 6(rVΘ)/dm in Fig. 16 shows relatively flat distribution both at the shroud and the hub and is different from the distribution described in the second aspect of the present invention which is favorable to suppress secondary flows. So, it is confirmed that there is no consideration in the loading distribution of the conventional design to suppress secondary flows. As a result, the results of numerical analysis predicted strong secondary flows on the blade suction surface from the hub to the shroud as presented in Fig. 17 showing unfavorable flow pattern. On the other hand, in the case of the verification data of Figs. 18 through 20, where the third aspect of the present invention was applied, the first aspect of the present invention is realized due to the effects of the blade lean in the circumferential direction. Fig. 18 shows a remarkable decreasing tendency of the pressure coefficient difference ΔCp on the suction surface between the hub and the shroud towards impeller exit, and ΔCp changes its sign at the non-dimensional meridional distance of about 0.7, giving sufficient tendency to suppress secondary flows. In this case, although the maximum loading at the shroud occurs in the fore part of the blade, the maximum loading at the hub does not occur at the latter part of the blade and the second aspect of the present invention is not sufficiently satisfied. However, the third aspect of the present invention is fully utilized by using the blade lean ratio of S=0.142, and so the favorable pressure distribution is obtained as shown in Fig. 18. As a result of this, the secondary flows from the hub to the shroud on the suction surface are completely suppressed and a favorable flow pattern is achieved as shown in the result of numerical analysis of Fig. 20. Next, Figs. 21 though 23 show the case which was obtained by the combination of the third and the second aspects of the present invention. Because of the synergism of the effects of the third aspect of blade lean in the circumferential direction and the second aspect of the blade loading, a remarkable decreasing tendency of the pressure coefficient difference ΔCp on the suction surface between the hub and the shroud towards impeller exit is realized, and ΔCp changes it sign at the non- dimensional meridional distance of about 0.75, giving sufficient tendency to suppress secondary flows. As a result of this, the secondary flows from the hub to the shroud on the suction surface are suppressed and a favorable flow pattern is achieved as shown in the result of numerical analysis of Fig. 23. In this case, the blade lean ratio is S=0.078, which is about the half of that obtained for the case of Figs. 18 through 20 designed using only the third aspect. As a result of this, a three dimensional flow passage geometry, which is more favorable for manufacturing, has been successfully obtained. Figs. 24 though 26 show the case where only the second aspect of the present invention is applied. Because of the favorable distribution of the blade loading d ( rW. ) /dm of Fig. 25, the pressure coefficient difference ΔCp on the suction surface between the hub and the shroud shows a remarkably decreasing tendency towards impeller exit compared to that for conventional design, having a favorable tendency to suppress secondary flows. As a result of this, the secondary flows from the hub to the shroud on the suction surface are fully suppressed as shown in the result of numerical analysis of Fig. 26 compared to that of conventional design case of Fig. 17. However, as is apparent from the comparison between Figs. 18 and 24, the decreasing tendency, in Fig. 24, of the pressure coefficient on the suction surface between the hub and the shroud is not as remarkable as that of Fig. 18. As a result of this, and as is apparent from the comparison between Figs. 20 and 26, the secondary flow is suppressed less effectively than the previous design using the third aspect of the present invention. Figs. 27 through 38 show the verification data for mixed flow pump impellers, and the effects of secondary flow suppression similar to those for centrifugal pump impellers described above (Figs. 15 through 26), and the applicability of the present invention can be confirmed also for mixed flow pump impellers. Next, the verification data of Figs. 39 through 50 for centrifugal compressor impellers will be described in detail. As is described in the previous section, in compressible flows, the physical variable, which is closely related to the behavior of secondary flows and is used in the first aspect of the present invention, is the relative Mach number distribution. In the conventionally designed impellers, there is no consideration on the suppression of secondary flows, and the meridional distribution of the relative Mach number M towards impeller exit does not show sufficient decreasing tendency of relative Mach number M on the suction surface between the hub and the shroud as shown in Fig. 39. Furthermore, the blade loading in Fig. 40 shows relatively flat distribution both at the hub and the shroud and is different from the favorable distribution to suppress secondary flows. So, it is confirmed that there is no consideration in the loading distribution of the conventional design to suppress secondary flows. Thus, the results of the numerical analysis 1 presented in Fig. 41, predicts strong secondary flows on the blade suction surface from the hub to the shroud towards impeller exit. On the other hand, in the case of the verification data of Figs. 42 through 44, where the third aspect of the present invention was applied, the first aspect of the present invention is realized due to the effects of the blade lean in the circumferential direction. Fig.42 shows a remarkable decreasing tendency of the relative Mach number difference ΔM on the suction surface between the hub and the shroud towards impeller exit, and ΔM changes it sign at the non-dimensional meridional distance of about 0.85, giving sufficient tendency to suppress secondary flows. In this case, the maximum loading at the shroud occurs in the latter part of the blade and the second aspect of the present invention is not sufficiently satisfied. However, the third aspect of the present invention is fully utilized, and so the favorable relative Mach number distribution is obtained as shown in Fig. 42. As a result of this, the secondary flows from the hub to the shroud on the suction surface are completely suppressed and a favorable flow pattern is achieved as shown in the result of numerical analysis of Fig. 44. Next, Figs. 45 though 47 show the case which was obtained by the combination of the third and the second aspects of the present invention. Because of the synergism of the effects of the third aspect of blade lean in the circumferential direction and the second aspect of the blade loading, a remarkable decreasing tendency of the relative Mach number difference ΔM on the suction surface between the hub and the shroud towards impeller exit, can be seen in Fig. 45, which shows that ΔM changes its sign at the non-dimensional meridional distance of about 0.75, giving sufficient tendency to suppress secondary flows. As a result of this, the secondary flows from the hub to the shroud on the suction surface are suppressed and a favorable flow pattern is achieved as shown in the result of numerical analysis of Fig. 47. In this case, the blade lean ratio is similar to that applied for the case of Figs. 42 through 44 designed using only the third aspect. However, as is apparent from the comparison between Figs. 42 and 45, the remarkable decreasing tendency of the relative Mach number is achieved by using the second aspect of the present invention in combination. Figs. 48 though 50 show the case where only the second aspect of the present invention is applied. Because of favorable distribution of the blade loading ->(rVΘ)/r3m presented in Fig. 49, the relative Mach number difference ΔM on the suction surface between the hub and the shroud shows a decreasing tendency towards impeller exit, and ΔM becomes very close to zero at the non-dimensional meridional distance of about 0.90, giving favorable distribution to suppress secondary flows. As a result of this, the secondary flows from the hub to the shroud on the suction surface are completely suppressed as shown in the result of numerical analysis of Fig. 50 and a favorable flow pattern is achieved. Figs. 51 through 62 show the verification data for mixed flow compressor impellers, and the effects of secondary flow suppression similar to those for centrifugal compressor impellers described above (Figs. 39 through 50), and the applicability of the present invention can be confirmed also for mixed flow compressor impellers. Figs. 63 through 74 show a characteristic blade angle distribution according to the forth aspect of the present invention. Through the comparison with the predicted results of meridional component of secondary flow by numerical analysis shown in Figs. 15 through 38, it is ascertained that the effect of suppression of the meridional component of secondary flow increases as the characteristic feature of blade angle distribution becomes remarkable. Namely, the results of the numerical analysis for the centrifugal pump impellers show the increase of the effects of secondary flow suppression in the order of Figs. 20, 23 and 26, and the corresponding distribution of the blade angle difference in Figs. 65 shows the tendency of the increasing blade angle difference towards the impeller exit more remarkably in the same order. So, the effectiveness of the fourth aspect of the present invention can be confirmed very clearly. Also, the blade angle distributions on the hub and the shroud show characteristic feature of the fourth aspect of the present invention described in relation to Figs. 13 (A) and (B). In the verification data for mixed flow pump impellers, the characteristic feature of the blade angles β H and β 3 can be observed in Figs. 66 through 68. It can be confirmed that, similar to those of centrifugal pump impeller cases, the fourth aspect of the present invention of Fig. 68 is well correlated with the degree of secondary flow suppression presented in Figs. 32, 35 and 38. In the above drawings, the blade angle distribution of the impeller designed by a conventional design method are shown by broken lines. The difference between the conventional method and the present invention can be clearly acknowledged. According to the conventional design method of a centrifugal or mixed flow pump impeller, as described in A. J.Stepanoff, "Centrifugal and Axial Flow Pumps", 2nd ed., John Wiley & Sons, New York, 1957, pp. 95-104, or J.L.Dianas, "Vertical Turbine, Mixed Flow and Propeller Pumps", MacGraw-Hill, New York, 1989, pp.305-311, once the blade angles at the impeller inlet and the impeller exit are determined from design specification, the blade angle distribution between the impeller inlet and the impeller exit is generally determined by connecting both of them by a curved line which changes gradually and smoothly. So, the consideration for suppressing secondary flow is not made at all in the conventional design method. The broken lines in Figs. 63 through 68 show the blade angle distributions of the impellers designed by such a conventional method. Because of this, the secondary flows can not be suppressed in the conventional impellers, and this is confirmed by the results of numerical analysis of Figs. 17 and 29, which correspond to the conventional impellers presented with broken lines in Figs. 65 and 68. Figs. 69 through 74 show the verification data for compressor impellers having the characteristic feature of the blade angle distribution of the fourth aspect of the present invention. Through the comparison with the predicted meridional secondary flows by the numerical analysis in Figs. 44, 47, 50, 56, 59, 69, it can be confirmed clearly that the effectiveness of the secondary flow suppression increases as the characteristic feature of the blade angle distribution of the fourth aspect becomes increasingly remarkable. Here, broken lines in Figs. 69 through 74 show the blade angle distributions of the conventional impellers of a typical compressor. In the compressor, since deceleration of the flow along the shroud surface is completed in the fore part of the blade in order to reduce friction loss, the blade angle distribution is different from that of the pump impeller. However, since the consideration for secondary flow suppression is not made at all, the meridional component of secondary flow cannot be suppressed. This can be confirmed by the corresponding results of numerical analysis as shown in Figs. 41 and 53. Fig. 75 shows the effectiveness to combine the third aspect of the blade lean and the second aspect of the blade loading in the case of mixed flow pump impeller. The geometry of the impeller, designed by using only the third aspect of the present invention, has the blade lean ratio of S-0.217 showing substantial blade lean in the circumferential direction with very small blade lean angle γ as shown in Fig. 75(A), and the impeller has the geometry which is difficult to manufacture. On the other hand, the impeller, designed by using the second aspect of the present invention in combination with the third aspect, has the halved blade lean ratio of S=0.10, relatively large blade lean angle γ as shown in Fig. 75(B), thus solving the difficulty in manufacturing, and the effectiveness of the present invention is ascertained. Industrial Applicability The present invention is characterized by the following aspects in designing an impeller in a turbomachinery. (1) The reduced static pressure difference ΔCp or the relative Mach number difference ΔM shows a remarkably decreasing tendency with respect to non-dimensional meridional distance. (2) The meridional derivative of V. t i.e. d ( τ ~ V t )/d which is derivative of product ( τ ~ V t ) of radius r in the impeller and the tangentially averaged value of the circumferential velocity component V β of absolute velocity V of fluid at the radius r, has such distribution along the meridional distance that the peak value of d(r " V„)/θm at the shroud emerges on the fore part of the blade and/or that at the hub on the latter part of the blade. (3) The blade is leaned toward a circumferential direction with a non constant blade lean angle γ so that the blade at the hub precedes the blade at the shroud in a rotational direction of the impeller. (4) The blade angle difference β H - β 3 between blade angle β H at the hub and blade angle β 3 at the shroud shows a remarkably increasing tendency with respect to non-dimensional meridional distance. According to the above aspects, since the meridional component of secondary flow can be effectively suppressed, a loss which occurs in the turbomachinery or the downstream flow passage can be reduced, emergence of positively sloped characteristic curve can be avoided, and stability of operation can be improved. Therefore, the present invention has a great utility value in industry. With regard to the above four aspects, those propriety is substantiated by verification data in accordance with the inverse design method and the flow analysis techniques, which are now widely available, therefore the present invention can be utilized effectively in industry.

Claims: CLAIMS 1. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that reduced static pressure difference ΔCp or relative Mach number difference ΔM between the hub and the shroud on the suction surface of a blade shows a remarkably decreasing tendency between the location of non-dimensional meridional distance m*0 (impeller inlet) and the location of non-dimensional meridional distance m=1.0 (impeller exit) as said non-dimensional meridional distance increases. 2. The centrifugal or mixed flow turbomachinery as recited in claim 1, wherein said remarkably decreasing tendency is arranged such that the difference between a minimum value of reduced static pressure difference ΔCp and a value of reduced static pressure difference ΔCp at the location corresponding to non-dimensional meridional distance obtained by subtracting non- dimensional meridional distance 0.4 from non-dimensional meridional distance representing said minimum value is not less than 0.2, or the difference between a minimum value of relative Mach number difference ΔM and a value of relative Mach number difference ΔM at the location corresponding to non-dimensional meridional distance obtained by subtracting non-dimensional meridional distance 0.4 from non-dimensional meridional distance representing said minimum value is not less than 0.15. 3. The centrifugal or mixed flow turbomachinery as recited in claim 1 or 2, wherein said remarkably decreasing tendency emerges at the location in the range of non-dimensional meridional distance m»0.4-1.0. 4. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that reduced static pressure difference ΔCp or relative Mach number difference ΔM between the hub and the shroud on the suction surface of a blade becomes nearly zero or changes its sign from positive to negative between the location of non-dimensional meridional distance 0 (impeller inlet) and the location of non- dimensional meridional distance 1.0 (Impeller exit). 5. The centrifugal or mixed flow turbomachinery as recited in claim 4, wherein said nearly zero value of ΔCp (or ΔM) or said change sign from positive to negative emerges at the location of non-dimensional meridional distance of not less than 0.6 (impeller exit side) . 6. The centrifugal or mixed flow turbomachinery as recited in claim 4, wherein said nearly zero value of ΔCp(or ΔM) or said change sign from positive to negative emerges at the location in the range of non-dimensional meridional distance m=0.65-0.9. 7. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed such that meridional derivative of rT , i.e. d(rV„)/dm which is the derivative of (r " V,), product of radius r in the impeller and the tangentially averaged value of the circumferential velocity component V, of the absolute velocity V of the fluid in a flow passage in the said impeller at the said radius r, has a maximum value at the shroud at the location in the fore part of the impeller, and/or said meridional derivative of r?„i.e. d ( τ ~ V. ) /dm has a maximum value at the hub at the location in the aft part of the impeller. 8. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed such that meridional derivative of τV t , i.e. d ( τV t ) /dm which is derivative of (r?,) ( product of radius r in the impeller and the tangentially averaged value of the circumferential velocity component V, of the absolute velocity V of the fluid in a flow passage in the said impeller at the said radius r, is larger at the shroud of the blade than that at the hub of the blade at the location in the vicinity of non-dimensional meridional distance m=0 (impeller inlet side), and the said meridional derivative of r?„ i.e. d ( r ~ V t ) /dm, is smaller at the shroud than that at the hub at the location in the vicinity of non-dimensional meridional distance m=1.0 (impeller exit side), and the distribution of said meridional derivative of τV t at the hub crosses over that at the shroud between said location of non-dimensional meridional distance 0 (impeller inlet) and said location of non-dimensional meridional distance m=1.0 (impeller exit). 9. The centrifugal or mixed flow turbomachinery as recited in claim 8, wherein said crossing emerges at the location in the range of non-dimensional meridional distance of not less than 0.50 (impeller exit side). 10. The centrifugal or mixed flow turbomachinery as recited in claim 8, wherein said crossing emerges at the location in the range of non-dimensional meridional distance 0.55-0.85 (impeller exit side). 11. The centrifugal or mixed flow turbomachinery having a impeller, wherein said blade is leaned toward a circumferential direction, between the location of non-dimensional meridional distance m=0 (impeller inlet) and the location of non-dimensional distance m=1.0 (impeller exit), so that the blade at the hub side precedes the blade at the shroud side in a rotational direction of said impeller and the blade lean angle γ, defined as an angle between the shroud surface and the blade center line on the cross-sectional view of the flow passage, shows a decreasing tendency towards the exit of the said impeller. 12. The centrifugal or mixed flow turbomachinery as recited in claim 11, wherein the blade angle α_,- defined as an angle between the blade camber line at the exit and the tangential direction on the plan view of the said impeller, viewed from the impeller inlet side, decreases from the value of . κ „ at the hub towards the value of a~ z 3 at the shroud. 13. The centrifugal or mixed flow turbomachinery as recited in claim 12, wherein the blade angle a^, defined as an angle between the blade camber line and the tangential direction at the hub on the plan view of said impeller, viewed from the impeller inlet side, has a maximum or a local maximum value at the exit side of the impeller. 14. The centrifugal or mixed flow turbomachinery as recited in claim 13, wherein said maximum or local maximum value of the blade angle __„ emerges at the location of non-dimensional meridional distance of not less than m=0.7 or at the location of non-dimensional radius ratio of not less than r*=0.6 (impeller exit side). 15. The centrifugal or mixed flow turbomachinery having a impeller, wherein said blade is leaned toward circumferential direction between the location of non-dimensional meridional distance m=0 (impeller inlet) and the location of non-dimensional distance m=1.0 (impeller exit) so that the blade at the hub side precedes the blade at the shroud side in the rotational direction of the said impeller and the backsweep of the blade becomes significant at the shroud towards impeller exit giving {d(rθ)/dr} s >{d(rθ)/dr}„. 16. The centrifugal or mixed flow turbomachinery as recited in any one of claims 11 to 15, wherein said blade lean in the circumferential direction is designed such that the blade lean ratio S, which represents overall degree of the blade lean and defined as S= ( Θ TE . S ~ΘT_. H )/(ΘT E . M ~€»L E .M) : ,where Θ TE g represents the angular co-ordinate of the impeller exit at the shroud, Θ TE H that of the impeller exit at the hub, Θ TE M that of the impeller exit at the middle point between the shroud and the hub, Θ LE M that of the impeller inlet at the middle point between the shroud and the hub, is not less than S=0.14. 17. The centrifugal or mixed flow turbomachinery as recited in any one of claims 1 to 10, wherein the said blade is leaned towards the circumferential direction so that the blade at the hub precedes the blade at the shroud in the rotational direction of the impeller. 18. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that the distribution of blade angle measured form the circumferential direction is arranged such that blade angle difference β H -β s , between the blade angle β H at the hub and the blade angle β s at the shroud, shows a remarkably increasing tendency from the location of non-dimensional meridional distance 0 (impeller inlet) to the location of non-dimensional meridional distance 1.0 (impeller exit) as the said non-dimensional meridional distance increases. 19. A centrifugal or mixed flow turbomachinery having an impeller, characterized in that said impeller is designed so that the distribution of blade angle measured form the circumferential direction is arranged such that blade angle difference β H -β s , between blade angle β H at the hub and blade angle β s at the shroud shows a decreasing tendency in a part of the non- dimensional distance between the location of non-dimensional meridional distance 0 (impeller inlet) and the location of non- dimensional meridional distance 1.0 (impeller exit) followed by a remarkably increasing tendency as said non-dimensional meridional distance increases. 20. The centrifugal or mixed flow turbomachinery as recited in claim 18 or 19, wherein said blade angle difference β H -β 3 has a maximum value or a local maximum value which is not less than 20 degrees. 21. The centrifugal or mixed flow turbomachinery as recited in claim 18 or 19, wherein said blade angle difference β H -β s has a maximum value or a local maximum value which is larger than the averaged value of blade angle difference between the location of non-dimensional meridional distance of m«0-0.2 (impeller inlet side) by not less than 10 degrees. 22. The centrifugal or mixed flow turbomachinery as recited in claim 20 or 21, wherein said maximum value of said blade angle difference β H -β s emerges at the location in the range of non- dimensional meridional distance not less than m-0.5 (in the impeller exit side). 23. The centrifugal or mixed flow turbomachinery as recited in claim 20 or 21, wherein said maximum value of said blade angle difference β H -β s emerges at the location in the range of non- dimensional meridional distance m=0.7-1.0 (near the impeller exit).

Inventors: GOTO AKIRA, HARADA HIDEOMI, ZANGENEH MEHRDAD, GOTO, AKIRA, HARADA, HIDEOMI, ZANGENEH, MEHRDAD
Applicants: EBARA CORP, EBARA RES CO LTD, GOTO AKIRA, HARADA HIDEOMI, UNIV LONDON, ZANGENEH MEHRDAD, EBARA CORPORATION, EBARA RESEARCH CO., LTD., GOTO, AKIRA, HARADA, HIDEOMI, UNIVERSITY COLLEGE LONDON, ZANGENEH, MEHRDAD
Application year: 1994
Main classification: F01D5
IPC codes: F01D5, F04D29