Method of designing blade of axial flow fluid machine and blade

Abstract

To provide a method of designing a blade of an axial flow fluid machine that has a blade surface whose radius of curvature is continuous at a leading edge thereof and has a high aerodynamic performance. The method includes a step of determining a pressure surface curve and a suction surface curve as curves capable of first to third order differentiations at respective connection points to a leading edge curve, that is, a pressure surface connection point and a suction surface connection point, and a step of forming the leading edge curve as a fifth order Bezier curve that is defined by a first control point, a second control point, a third control point, a fourth control point, a fifth control point and a sixth control point. The first control point is the suction surface connection point. The sixth control point is the pressure surface connection point. Provided that an intersection of a tangent to the pressure surface curve at the pressure surface connection point and a tangent to the suction surface curve at the suction surface connection point is referred to as a point ahead of the blade, the second and fifth control points are points that internally divide line segments connecting the point ahead of the blade to the suction surface connection point and the pressure surface connection point, respectively. The third and fourth control points are points having coordinates that are determined as solutions of simultaneous equations obtained by applying continuity conditions for first to third differential coefficients at the first and sixth control points to a fifth order Bezier function.

Claims

1. A blade having an airfoil that is formed by a leading edge curve, a trailing edge curve, a pressure surface curve and a suction surface curve, wherein the pressure surface curve is a curve that is connected to the leading edge curve at a pressure surface connection point and whose first to third order differential coefficients at the pressure surface connection point are known, the suction surface curve is a curve that is connected to the leading edge curve at a suction surface connection point and whose first to third order differential coefficients at the suction surface connection point are known; the leading edge curve is formed as a fifth order Bezier curve, the fifth order Bezier curve is defined by a first control point, a second control point, a third control point, a fourth control point, a fifth control point and a sixth control point, the first control point coincides with the suction surface connection point, the sixth control point coincides with the pressure surface connection point, provided that an intersection of a tangent to the pressure surface curve at the pressure surface connection point and a tangent to the suction surface curve at the suction surface connection point is referred to as a point ahead of the blade, the second control point is a point that internally divides a line segment connecting the point ahead of the blade and the suction surface connection point with a ratio of p to (1−p) (where 0<p<1), the fifth control point is a point that internally divides a line segment connecting the point ahead of the blade and the pressure surface connection point with a ratio of q to (1−q) (where 0<q<1), and the third control point and the fourth control point are points having coordinates that are determined as solutions of simultaneous equations obtained by applying continuity conditions for first to third differential coefficients at the first control point and the sixth control point to a fifth order Bezier function that defines the fifth order Bezier curve.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic diagram for illustrating an airfoil of a blade according to an embodiment of the present disclosure;

(2) FIG. 2 is a diagram for illustrating boundary conditions in a case where a third order Bezier curve is applied to a leading edge of the airfoil shown in FIG. 1;

(3) FIG. 3 is a diagram for illustrating boundary conditions in a case where a fifth order Bezier curve is applied to a leading edge of the airfoil shown in FIG. 1;

(4) FIG. 4A is a graph showing a shape of a part of a conventional airfoil including a leading edge and a distribution of the radius of curvature of a blade surface;

(5) FIG. 4B is a graph showing a shape of a part of the airfoil of the blade including the leading edge and a distribution of the radius of curvature of the blade surface of the present disclosure;

(6) FIG. 5A is a graph showing a distribution of the blade surface Mach number at a design point in the entire region in the chord direction;

(7) FIG. 5B is a graph showing a distribution of the blade surface Mach number at a design point in a region around the leading edge (the region surrounded by the dashed line in FIG. 5A);

(8) FIG. 6 is a graph showing a loss coefficient; and

(9) FIG. 7 is a diagram for illustrating a method of drawing a fifth order Bezier curve.

MODE FOR CARRYING OUT THE DISCLOSURE

(10) In the following, an embodiment of the present disclosure will be described in detail with reference to the drawings.

(11) FIG. 1 is a schematic diagram for illustrating an airfoil of a blade according to an embodiment of the present disclosure.

(12) An airfoil A of a blade according to an embodiment of the present disclosure is formed by a leading edge LE, a trailing edge TE, and a concave pressure surface PS and a convex suction surface SS that extend between the edges.

(13) In the following, curves forming the leading edge LE, the trailing edge TE, the pressure surface PS and the suction surface SS will be referred to as a leading edge curve CL, a trailing edge curve CT, a pressure surface curve CP, and a suction surface curve CS, respectively, and a closed curve formed by these four curves or, in other words, a curve forming the airfoil A will be referred to as an airfoil curve CA.

(14) In the airfoil A of the blade according to the embodiment of the present disclosure, the pressure surface curve CP and the suction surface curve CS are curves capable of first to third order differentiations at a pressure surface connection point TP and a suction surface connection point TS both described later, respectively, and have shapes independently determined in a process using CFD, for example. The trailing edge curve CT can be a circular or elliptic arc. The leading edge curve CL is formed as a Bezier curve.

(15) In the following, a process of determining the shape of the Bezier curve adopted as the leading edge curve CL will be described.

(16) First, a Bezier function that defines a Bezier curve is expressed by the formula (1).

(17) $\begin{array}{cc}P\left(t\right)=\underset{i=0}{\overset{N}{.Math.}}{X}_{j+1}{J}_{n,i}\left(t\right)0=t=1& \left(1\right)\end{array}$

(18) In this formula, N denotes the order of the Bezier function, X denotes the coordinates of a control point, and J denotes a Bernstein basis function defined by the formula (2).

(19) $\begin{array}{cc}{J}_{n,i}\left(t\right)=\left(\begin{array}{c}n\\ i\end{array}\right){t^i\left(1-t\right)}^{n-i}\left(\begin{array}{c}n\\ i\end{array}\right)=\frac{n!}{i!\left(n-i\right)!}& \left(2\right)\end{array}$

(20) The Bezier curve (leading edge curve CL) to be determined is connected to the pressure surface curve CP and the suction surface curve CS at the pressure surface connection point TP and the suction surface connection point TS, respectively. In order for the radius of curvature of the airfoil curve CA at each of these connection points to be continuous, the following boundary conditions need to be satisfied.

(21) (1) A continuity condition for the first order differential coefficient (first order differential coefficients (dy/dx) on opposite sides of each connection point are equal to each other).

(22) (2) A continuity condition for the second order differential coefficient (second order differential coefficients (d.sup.2y/dx.sup.2) on opposite sides of each connection point are equal to each other).

(23) The discontinuity in radius of curvature of the conventional airfoil whose leading edge curve is a simple curve such as a circular or elliptic arc that occurs at each of the connection point between the leading edge curve and the pressure surface curve and the connection point between the leading edge curve and the suction surface curve is resolved if the boundary conditions (1) and (2) described above are satisfied.

(24) Furthermore, according to the boundary layer equation, it is known to be desirable from the viewpoint of aerodynamic performance of the blade that

(25) (3) A continuity condition for the third order differential coefficient (third order differential coefficients (d.sup.3y/dx.sup.3) on opposite sides of each connection point are equal to each other) is satisfied on the blade surface.

(26) Thus, in determining the shape of the Bezier curve for the airfoil A of the blade according to the embodiment of the present disclosure, all the boundary conditions (1) to (3) described above are imposed.

(27) First, to obtain the Bezier curve that satisfies the boundary conditions described above, what order the Bezier curve should be of will be discussed.

(28) Here, a third order Bezier curve has four control points, and a fifth order Bezier curve has six control points. Two of the control points are end points of the Bezier curve, that is, the pressure surface connection point TP and the suction surface connection point TS, and the coordinates of the two control points are known. Therefore, the third order Bezier curve has four unknown quantities (the x coordinate and the y coordinate of the two control points other than the end points), and the fifth order Bezier curve has eight unknown quantities (the x coordinate and the y coordinate of the four control points other than the end points).

(29) As can be seen from the above description, when a high order Bezier curve is used, a large number of unknown quantities need to be handled, and the computational complexity is high.

(30) In view of this, first, whether the third order Bezier curve can be applied is examined. For the third order Bezier curve, as shown in FIG. 2, four control points X.sub.1, X.sub.2, X.sub.3 and X.sub.4 are used. Of these control points, the control points X.sub.1 and X.sub.4 are the suction surface connection point TS and the pressure surface connection point TP in FIG. 1, respectively.

(31) The formula (1) can be rewritten for the third order Bezier function as follows.

(32) $\begin{array}{cc}P\left(t\right)=\left(\begin{array}{c}3\\ 0\end{array}\right)X_1{t^0\left(1-t\right)}^3+\left(\begin{array}{c}3\\ 1\end{array}\right)X_2{t^1\left(1-t\right)}^2+\left(\begin{array}{c}3\\ 2\end{array}\right)X_3{t^2\left(1-t\right)}^1+\left(\begin{array}{c}3\\ 3\end{array}\right)X_4{t^3\left(1-t\right)}^0& \left(3\right)\\ P\left(t\right)=X_1\left(1-3t+3t^2-t^3\right)+3X_2\left(t-2t^2+t^3\right)+3X_3\left(t^2-t^3\right)+X_4{t}^3=X_1+\left(-3X_1+3X_2\right)t+\left(3X_1-6X_2+3X_3\right)t^2+\left(-X_1+3X_2-3X_3+X_4\right)t^3& \left(4\right)\end{array}$

(33) The function P(t) is differentiated with respect to a parameter t as follows.
{dot over (P)}(t)=(−3X.sub.1+3X.sub.2)+2(3X.sub.1−6X.sub.2+3X.sub.3)t+3(−X.sub.1+3X.sub.2−3X.sub.3+X.sub.4)t.sup.2  (5)

(34) This can be rewritten as follows.

(35) $\begin{array}{cc}\frac{\mathrm{dy}}{\mathrm{dx}}=\frac{\stackrel{.}{y}}{\stackrel{.}{x}}=\frac{\begin{array}{c}\left(-3y_1+3y_2\right)+2\left(3y_1-6y_2+3y_3\right)t+\\ 3\left(-y_1+3y_2-3y_3+y_4\right)t^2\end{array}}{\begin{array}{c}\left(-3x_1+3x_2\right)+2\left(3x_1-6x_2+3x_3\right)t+\\ 3\left(-x_1+3x_2-3x_3+x_4\right)t^2\end{array}}& \left(6\right)\end{array}$

(36) At the end point X.sub.1 (t=0) and the end point X.sub.4 (t=1), the following relations hold.

(37) $\begin{array}{cc}\frac{\mathrm{dy}}{\mathrm{dx}}{|}_{x=x_1}=\frac{-y_1+y_2}{-x_1+x_2}& \left(7\right)\\ \frac{\mathrm{dy}}{\mathrm{dx}}{|}_{x=x_{}}=\frac{-y_3+y_4}{-x_3+x_4}& \left(8\right)\end{array}$

(38) The formulas (7) and (8) show that the continuity condition for the first order differential coefficient (boundary condition (1)) is satisfied if the control point X.sub.2 is set on a tangent at the end point X.sub.1, that is, a suction surface tangent LT.sub.1, and the control point X.sub.3 is set on a tangent at the end point X.sub.4, that is, a pressure surface tangent LT.sub.4.

(39) By setting the control points X.sub.2 and X.sub.3 as described above, the number of unknown quantities decreases from four to two. However, in order to satisfy the continuity condition for the second order differential coefficient (boundary condition (2)) and the continuity condition for the third order differential coefficient (boundary condition (3)), at least four unknown quantities are required at this stage. This means that the order of the third order Bezier function is insufficient, and the third order Bezier curve cannot be applied.

(40) Thus, the fifth order Bezier curve is applied. For the fifth order Bezier curve, as shown in FIG. 3, six control points X.sub.1, X.sub.2, X.sub.3, X.sub.4, X.sub.5 and X.sub.6 are used. Of these control points, the control points X.sub.1 and X.sub.6 are the suction surface connection point TS and the pressure surface connection point TP in FIG. 1, respectively.

(41) The formula (1) can be rewritten for the fifth order Bezier function as follows.

(42) $\begin{array}{cc}P\left(t\right)=\left(\begin{array}{c}5\\ 0\end{array}\right)X_1{t^0\left(1-t\right)}^5+\left(\begin{array}{c}5\\ 1\end{array}\right)X_2{t^1\left(1-t\right)}^4+\left(\begin{array}{c}5\\ 2\end{array}\right)X_3{t^2\left(1-t\right)}^3+\left(\begin{array}{c}5\\ 3\end{array}\right)X_4{t^3\left(1-t\right)}^2+\left(\begin{array}{c}5\\ 4\end{array}\right)X_5{t^4\left(1-t\right)}^1+\left(\begin{array}{c}5\\ 5\end{array}\right)X_6{t^5\left(1-t\right)}^0& \left(9\right)\\ \begin{array}{c}P\left(t\right)=X_1\left(1-5t+10t^2-10t^3+5t^4-t^5\right)+5X_{2}t(1-4t+6t^2-\\ 4t^3+t^4)+10X_3{t}^2\left(1-3t+3t^2-t^3\right)+10X_4{t}^3(1-3t+\\ 3t^2-t^3)+5X_5{t}^4\left(1-t\right)+X_6{t}^5\\ =X_1+\left(-5X_1+5X_2\right)t+\left(10X_1-20X_2+10X_3\right)t^2+\\ \left(-10X_1+30X_2-30X_3+10X_4\right)t^3+(5X_1-20X_2+30X_3-\\ 20X_4+5X_5)t^4+\left(-X_1+5X_2-10X_3+10X_4-5X_5+X_6\right)t^5\end{array}& \left(10\right)\end{array}$

(43) The formulas (9) and (10) are functions that express the pressure surface curve CP in FIG. 3. This is because the sign of the radius of curvature depends on how the parameter of a parametric function is taken. The suction surface curve CS can be expressed by the following two formulas.

(44) $\begin{array}{cc}P\left(t\right)=\left(\begin{array}{c}5\\ 0\end{array}\right)X_6{t^0\left(1-t\right)}^5+\left(\begin{array}{c}5\\ 1\end{array}\right)X_5{t^1\left(1-t\right)}^4+\left(\begin{array}{c}5\\ 2\end{array}\right)X_4{t^2\left(1-t\right)}^3+\left(\begin{array}{c}5\\ 3\end{array}\right)X_3{t^3\left(1-t\right)}^2+\left(\begin{array}{c}5\\ 4\end{array}\right)X_2{t^4\left(1-t\right)}^1+\left(\begin{array}{c}5\\ 5\end{array}\right)X_1{t^5\left(1-t\right)}^0& \left(11\right)\\ \begin{array}{c}P\left(t\right)=X_6\left(1-5t+10t^2-10t^3+5t^4-t^5\right)+5X_{5}t(1-4t+6t^2-\\ 4t^3+t^4)+10X_4{t}^2\left(1-3t+3t^2-t^3\right)+10X_3{t}^3(1-3t+\\ 3t^2-t^3)+5X_2{t}^4\left(1-t\right)+X_1{t}^5\\ =X_6+\left(-5X_6+5X_5\right)t+\left(10X_6-20X_5+10X_4\right)t^2+\\ \left(-10X_6+30X_5-30X_4+10X_3\right)t^3+(5X_6-20X_5+30X_4-\\ 20X_3+5X_2)t^4+\left(-X_6+5X_5-10X_4+10X_3-5X_2+X_1\right)t^5\end{array}& \left(12\right)\end{array}$

(45) Here, at the end point X.sub.6, the following relation is obtained by substituting 1 to t in the formula (10).

(46) $\begin{array}{cc}\frac{\mathrm{dy}}{\mathrm{dx}}{|}_{x=}=\frac{\mathrm{dy}/\mathrm{dt}}{\mathrm{dx}/\mathrm{dt}}=\frac{-5y_6+5y_5}{-5x_6+5x_5}=\frac{-y_6+y_5}{-x_6+x_5}& \left(13\right)\end{array}$

(47) At the end point X.sub.1, the following relation is obtained by substituting 1 to t in the formula (12).

(48) $\begin{array}{cc}\frac{\mathrm{dy}}{\mathrm{dx}}{|}_{x=x_1}=\frac{\mathrm{dy}/\mathrm{dt}}{\mathrm{dx}/\mathrm{dt}}=\frac{-5y_2+5y_1}{-5x_2+5x_1}=\frac{-y_2+y_1}{-x_2+x_1}& \left(14\right)\end{array}$

(49) As in the case of the third order Bezier function, the formulas (13) and (14) show that the continuity condition for the first order differential coefficient (boundary condition (1)) is satisfied if the control point X.sub.2 is set on a tangent at the end point X.sub.1, that is, a suction surface tangent LT.sub.1, and the control point X.sub.5 is set on a tangent at the end point X.sub.6, that is, a pressure surface tangent LT.sub.6.

(50) Next, the continuity condition for the second order differential coefficient (boundary condition (2)) will be examined.

(51) Here, it is assumed that the suction surface tangent LT.sub.1 and the pressure surface tangent LT.sub.6 intersect with each other at a point X.sub.0 ahead of the blade, a point that internally divides the line segment X.sub.0X.sub.1 with a ratio of (1−p) to p is designated as the control point X.sub.2, and a point that internally divides the line segment X.sub.0X.sub.6 with a ratio of (1−q) to q is designated as the control point X.sub.5 (where 0<p<1, 0<q<1).

(52) The second order differential coefficient at the end point X.sub.1 is expressed by the following formula.

(53) 0$\begin{array}{cc}\frac{d^{2}y}{{\mathrm{dx}}^2}{|}_{x=x_1}=\frac{\stackrel{\mathrm{.Math.}}{y}\stackrel{.}{x}-\stackrel{.}{y}\stackrel{\mathrm{.Math.}}{x}}{{\stackrel{.}{x}}^3}{|}_{x=x_1}& \left(15\right)\end{array}$

(54) In this formula, the following relations hold.
ü|.sub.x=x.sub.1=20x.sub.3−40x.sub.2+20x.sub.1  (16)
ÿ|.sub.x=x.sub.1=20y.sub.3−40y.sub.2+20y.sub.1  (17)
{dot over (x)}|.sub.x=x.sub.1=−5x.sub.2+5x.sub.1  (18)
{dot over (y)}|.sub.x=x.sub.1=−5y.sub.2+5y.sub.1  (19)

(55) However, the coordinates of the control point X.sub.2 depends on the parameter p and therefore can be regarded as constants.

(56) Here, substituting the formulas (16) to (19) into the formula (15) results in the following formulas.

(57) $\begin{array}{cc}\frac{d^{2}y}{{\mathrm{dx}}^2}{|}_{x=x_1}{\stackrel{.}{x}}^3{|}_{x=x_1}=\left(20y_3-40y_2+20y_1\right)\stackrel{.}{x}{|}_{x=x_1}-\left(20x_3-40x_2+20x_1\right)\stackrel{.}{y}{|}_{x=x_1}& \left(20\right)\\ y_3=\frac{\stackrel{.}{y}{|}_{x=x_1}}{\stackrel{.}{x}{|}_{x=x_1}}{x}_3+\frac{1}{20\stackrel{.}{x}{|}_{x=x_2}}\left\{\frac{d^{2}y}{{\mathrm{dx}}^2}{|}_{x=x_1}{\stackrel{.}{x}}^3{|}_{x=x_1}-20\left(\stackrel{.}{x}{|}_{x=x_1}{y}_1-\stackrel{.}{y}{|}_{x=x_1}{x}_1\right)+40\left(\stackrel{.}{x}{|}_{x=x_1}{y}_2-\stackrel{.}{y}{|}_{x=x_2}{x}_2\right)\right\}& \left(21\right)\end{array}$

(58) The formula (21) shows that, in FIG. 3, the control point X.sub.3 must lie on a straight line LP.sub.1 that is parallel to the suction surface tangent LT.sub.1.

(59) Calculating the second order differential coefficient at the end point X.sub.6 in the same manner results the following.

(60) $\begin{array}{cc}\frac{d^{2}y}{{\mathrm{dx}}^2}{|}_{x=x_6}{\stackrel{.}{x}}^3{|}_{x=x_6}=\left(20y_4-40y_5+20y_6\right)\stackrel{.}{x}{|}_{x=x_6}-\left(20x_4-40x_5+20x_6\right)\stackrel{.}{y}{|}_{x=x_6}& \left(22\right)\\ y_4=\frac{\stackrel{.}{y}{|}_{x=x_6}}{\stackrel{.}{x}{|}_{x=x_6}}{x}_4+\frac{1}{20\stackrel{.}{x}{|}_{x=x_6}}\left\{\frac{d^{2}y}{{\mathrm{dx}}^2}{|}_{x=x_6}{\stackrel{.}{x}}^3{|}_{x=x_6}-20\left(\stackrel{.}{x}{|}_{x=x_6}{y}_6-\stackrel{.}{y}{|}_{x=x_6}{x}_6\right)+40\left(\stackrel{.}{x}{|}_{x=x_6}{y}_5-\stackrel{.}{y}{|}_{x=x_6}{x}_5\right)\right\}& \left(23\right)\end{array}$

(61) The formula (23) shows that, in FIG. 3, the control point X.sub.4 must lie on a straight line LP.sub.6 that is parallel to the pressure surface tangent LT.sub.6.

(62) Finally, the continuity condition for the third order differential coefficient (boundary condition (3)) will be examined.

(63) The third order differential coefficient is expressed by the following formula.

(64) $\begin{array}{cc}\frac{d^{2}y}{{\mathrm{dx}}^2}=\frac{\stackrel{.}{x}\stackrel{\mathrm{.Math.}}{y}-\stackrel{\mathrm{.Math.}}{x}\stackrel{.}{y}}{{\stackrel{.}{x}}^3}& \\ \begin{array}{c}\frac{d^{3}y}{{\mathrm{dx}}^3}=\frac{d}{\mathrm{dt}}\left(\frac{d^{2}y}{{\mathrm{dx}}^2}\right)\frac{\mathrm{dt}}{\mathrm{dx}}=\frac{d}{\mathrm{dt}}\left(\frac{\stackrel{.}{x}\stackrel{\mathrm{.Math.}}{y}-\stackrel{\mathrm{.Math.}}{x}\stackrel{.}{y}}{{\stackrel{.}{x}}^3}\right)\stackrel{.}{x}\\ =\frac{\left\{{\stackrel{.}{x}}^2\stackrel{...}{y}-3\stackrel{\mathrm{.Math.}}{x}\stackrel{.}{x}\stackrel{\mathrm{.Math.}}{y}+\left(3{\stackrel{\mathrm{.Math.}}{x}}^2-\stackrel{\mathrm{.Math.}}{x}\stackrel{...}{x}\right)\stackrel{.}{y}\right\}}{{\stackrel{.}{x}}^5}\end{array}& \left(24\right)\end{array}$

(65) Equations that have to be satisfied at the end point X.sub.1 on the suction surface curve CS are as follows.

(66) $\begin{array}{cc}f\left(x_3,x_4\right)=\frac{d^{3}y}{{\mathrm{dx}}^3}{|}_{x=x_1}-\frac{\left\{x^2\stackrel{...}{y}-3\stackrel{\mathrm{.Math.}}{x}\stackrel{.}{x}\stackrel{\mathrm{.Math.}}{y}+\left(3{\stackrel{\mathrm{.Math.}}{x}}^2-\stackrel{...}{x}\stackrel{\mathrm{.Math.}}{x}\right)\stackrel{.}{y}\right\}}{{\stackrel{.}{x}}^5}{|}_{x=x_1}& \left(25\right)\\ \stackrel{.}{x}{|}_{x=x_1}=-5x_2+5x_1& \left(26\right)\\ \stackrel{.}{y}{|}_{x=x_1}=-5y_2+5y_1& \left(27\right)\\ \stackrel{\mathrm{.Math.}}{x}{|}_{x=x_1}=20x_3-40x_2+20x_1& \left(28\right)\\ \stackrel{\mathrm{.Math.}}{y}{|}_{x=x_1}=20y_3-40y_2+20y_1& \left(29\right)\\ \stackrel{...}{x}{|}_{x=x_1}=-60x_4+180x_3-180x_2+60x_1& \left(30\right)\\ \stackrel{...}{y}{|}_{x=x_1}=-60y_4+180y_3-180y_2+60y_1& \left(31\right)\end{array}$

(67) On the other hand, equations that have to be satisfied at the end point X.sub.6 on the pressure surface curve CP are as follows.

(68) $\begin{array}{cc}g\left(x_3,x_4\right)=\frac{d^{3}y}{{\mathrm{dx}}^3}{|}_{x=x_6}-\frac{\left\{x^2\stackrel{...}{y}-3\stackrel{\mathrm{.Math.}}{x}\stackrel{.}{x}\stackrel{\mathrm{.Math.}}{y}+\left(3{\stackrel{\mathrm{.Math.}}{x}}^2-\stackrel{...}{x}\stackrel{\mathrm{.Math.}}{x}\right)\stackrel{.}{y}\right\}}{{\stackrel{.}{x}}^5}{|}_{x=x_6}& \left(32\right)\\ \stackrel{.}{x}{|}_{x=x_6}=-5x_5+5x_6& \left(33\right)\\ \stackrel{.}{y}{|}_{x=x_6}=-5y_5+5y_6& \left(34\right)\\ \stackrel{\mathrm{.Math.}}{x}{|}_{x=x_6}=20x_4-40x_5+20x_6& \left(35\right)\\ \stackrel{\mathrm{.Math.}}{y}{|}_{x=x_6}=20y_4-40y_5+20y_6& \left(36\right)\\ \stackrel{...}{x}{|}_{x=x_6}=-60x_3+180x_4-180x_5+60x_6& \left(37\right)\\ \stackrel{...}{y}{|}_{x=x_6}=-60y_3+180y_4-180y_5+60y_6& \left(38\right)\end{array}$

(69) If the coordinates of the control point X.sub.2 and X.sub.5 are regarded as constants since the coordinates depend on the parameters p and q, the remaining unknown quantities are only the X coordinates of the control points X.sub.3 and X.sub.4, so that the simultaneous equations including the formulas (25) and (32) can be solved.

(70) Once all the control points X.sub.1, X.sub.2, X.sub.3, X.sub.4, X.sub.5 and X.sub.6 are determined as described above, one Bezier curve, that is, the leading edge curve CL can be drawn in the method described earlier with reference to FIG. 7.

(71) However, the parameters p and q are not determined yet. Therefore, Bezier curves of various shapes, that is, leading edge curves CL of various shapes can be drawn by changing the parameters p and q between 0 and 1. In designing the airfoil, the airfoil can be optimized by individually evaluating the aerodynamic performance of the blades having the leading edge curves CL of the various shapes determined as described above.

(72) Here, summarizing the method of designing the airfoil A of the blade according to the embodiment of the present disclosure described above, the method includes the following steps.

(73) (A) A step of determining a pressure surface curve CP as a curve capable of first to third order differentiations at a pressure surface connection point TP.

(74) (B) A step of determining a suction surface curve CS as a curve capable of first to third order differentiations at a suction surface connection point TS.

(75) (C) A step of selecting a first control point X.sub.1, a second control point X.sub.2, a third control point X.sub.3, a fourth control point X.sub.4, a fifth control point X.sub.5 and a sixth control point X.sub.6 that define a fifth order Bezier curve forming a leading edge curve CL. The step (C) includes the following sub-steps.

(76) (C-1) A sub-step of designating the suction surface connection point TS as the first control point X.sub.1.

(77) (C-2) A sub-step of designating the pressure surface connection point TP as the sixth control point X.sub.6.

(78) Provided that an intersection of a tangent to the pressure surface curve CP at the pressure surface connection point TP and a tangent to the suction surface curve CS at the suction surface connection point TS is referred to as a point X.sub.0 ahead of a blade, the step (C) further includes the following sub-steps.

(79) (C-3) A sub-step of designating, as the second control point X.sub.2, a point that internally divides a line segment connecting the point X.sub.0 ahead of the blade and the suction surface connection point TS with a ratio of p to (1−p) (where 0<p<1).

(80) (C-4) A sub-step of designating, as the fifth control point X.sub.5, a point that internally divides a line segment connecting the point X.sub.0 ahead of the blade and the pressure surface connection point TP with a ratio of q to (1-q) (where 0<q<1).

(81) (C-5) A sub-step of designating, as the third control point X.sub.3 and the fourth control point X.sub.4, points having coordinates that are determined as solutions of simultaneous equations obtained by applying continuity conditions for first to third differential coefficients at the first control point X.sub.1 and the sixth control point X.sub.6 to a fifth order Bezier function that defines the fifth order Bezier curve.

(82) The shape of the airfoil A of the blade according to the embodiment of the present disclosure described above can be summarized as follows.

(83) A pressure surface curve CP is a curve whose first to third order differential coefficients at a pressure surface connection point TP are known.

(84) A suction surface curve CS is a curve whose first to third order differential coefficients at a suction surface connection point TS are known.

(85) A leading edge curve CL is formed as a fifth order Bezier curve that is defined by a first control point X.sub.1, a second control point X.sub.2, a third control point X.sub.3, a fourth control point X.sub.4, a fifth control point X.sub.5 and a sixth control point X.sub.6.

(86) The first control point X.sub.1 coincides with the suction surface connection point TS.

(87) The sixth control point X.sub.6 coincides with the pressure surface connection point TP.

(88) Provided that an intersection of a tangent to the suction surface curve CS at the first control point X.sub.1 and a tangent to the pressure surface curve CP at the sixth control point X.sub.6 is referred to as a point X.sub.0 ahead of a blade,

(89) the second control point X.sub.2 is a point that internally divides a line segment X.sub.0X.sub.1 with a ratio of (1−p) to p (where 0<p<1).

(90) The fifth control point X.sub.5 is a point that internally divides a line segment X.sub.0X.sub.6 with a ratio of (1−q) to q (where 0<q<1).

(91) The third control point X.sub.3 and the fourth control point X.sub.4 are points having coordinates that are determined as solutions of simultaneous equations (formulas (25) and (32)) obtained by applying continuity conditions for first to third differential coefficients at the first control point X.sub.1 and the sixth control point X.sub.6 to a fifth order Bezier function that defines the fifth order Bezier curve.

(92) As described above, with the airfoil A of the blade according to the embodiment of the present disclosure, since the Bezier curve adopted as the leading edge curve CL is defined so that the first to third order differential coefficients thereof are continuous at the pressure surface connection point TP and the suction surface connection point TS, the radius of curvature of the blade surface is continuous at the connection points, and the aerodynamic performance of the blade is improved.

(93) With the airfoil A of the blade according to the embodiment of the present disclosure, the order of the Bezier curve adopted as the leading edge curve CL is the fifth order. However, by applying the continuity conditions for the first to third differential coefficients described above, the number of unknown quantities is significantly reduced, and the control points of the Bezier curve can be determined using only two variables (p and q). Therefore, compared with the case where control points selected by trial and error are used, that is, where eight unknown quantities (the x coordinate and the y coordinate of the four control points other than the control points at the opposite ends) are handled, a Bezier curve having a desired shape can be quickly determined, and the design productivity can be improved. In addition, since the two variables (p and q) remain as unknown quantities, the degree of freedom of the design of the airfoil is ensured.

(94) Next, effects and advantages of the blade according to the embodiment of the present disclosure will be described below.

(95) FIGS. 4A and 4B are graphs showing shapes of a part of the airfoil including the leading edge and distributions of the radius of curvature of the blade surface. FIG. 4A shows a conventional airfoil, and FIG. 4B shows the airfoil A of the blade according to the embodiment of the present disclosure. In each graph, the horizontal axis indicates the coordinate x in the chord direction, the vertical axis (left) indicates the blade thickness y, and the vertical axis (right) indicates the radius of curvature R. FIG. 4A shows only the suction surface of the conventional airfoil, and FIG. 4B shows both the pressure surface and the suction surface of the airfoil A.

(96) As shown in FIG. 4A, with the conventional airfoil whose leading edge curve CL is elliptic, the radius of curvature of the blade surface is discontinuous at a location (x=0) in the chord direction that corresponds to the connection point (suction surface connection point TS) between the leading edge curve CL and the suction surface curve CS.

(97) To the contrary, with the airfoil A of the blade according to the embodiment of the present disclosure, as shown in FIG. 4B, it is found that the radius of curvature of the blade surface is continuous at the locations (x=0) in the chord direction that correspond to the suction surface connection point TS and the pressure surface connection point TP.

(98) FIGS. 5 and 6 show results of analysis of both the airfoils using CFD and comparison of the airfoils in terms of aerodynamic performance of the blade.

(99) FIGS. 5A and 5B are graphs showing distributions of the blade surface Mach number at a design point. FIG. 5A shows a distribution over the entire region in the chord direction, and FIG. 5B shows a distribution over a region around the leading edge (the region surrounded by the dashed line in FIG. 5A). In each graph, the horizontal axis indicates the distance in the chord direction (shown as non-dimensional numbers with respect to the chord length), the vertical axis indicates the blade surface Mach number, upper plots show a distribution on the suction surface, and lower plots show a distribution on the pressure surface.

(100) As shown in FIG. 5A, the distributions of the blade surface Mach number for the conventional airfoil and the airfoil A of the blade according to the embodiment of the present disclosure exhibit no marked difference in the macroscopic view over the entire region in the chord direction. However, the specific view of the region around the leading edge (the region surrounded by the dashed line in FIG. 5A) in FIG. 5B shows that, in the distribution of the blade surface Mach number for the airfoil A of the blade according to the embodiment of the present disclosure, the irregularities in the distribution of the blade surface Mach number for the conventional airfoil are eliminated, and the blade surface Mach number changes smoothly.

(101) That is, with the airfoil A of the blade according to the embodiment of the present disclosure, acceleration and deceleration of the flow over the blade surface are reduced compared with the conventional airfoil, and the fan or compressor incorporating the blades with the airfoil A is expected to have a higher efficiency and a wider stable operating range.

(102) FIG. 6 is a graph showing a loss coefficient. The vertical axis indicates the loss coefficient (a non-dimensional index of a total pressure drop at the time of passing over the blade with respect to a dynamic pressure), and the horizontal axis indicates the flowrate (the amount of intake air of the fan or compressor incorporating the blades). The flowrate corresponds to the incidence (inflow angle) at the inlet of the blade. The incidence is low when the flowrate is high, and is high when the flowrate is low.

(103) As shown in the drawing, it is found that when the flowrate is high (the incidence is low), no marked difference is observed. However, when the flowrate is low (the incidence is high), the airfoil A of the blade according to the embodiment of the present disclosure provides smaller loss coefficients than the conventional airfoil.

(104) With the conventional airfoil, the loss coefficient is high when the flowrate is low (the incidence is high), and in order to prevent occurrence of a stall, the operating range needs to be limited to a range of higher flowrates (lower incidences).

(105) To the contrary, with the airfoil A of the blade according to the embodiment of the present disclosure, the loss coefficient is low even when the flowrate is low (the incidence is high), so that the operating range does not need to be limited to a range of higher flowrates (lower incidences). That is, the fan or compressor incorporating the blades according to the embodiment of the present disclosure can have a stable operating range expanded to lower flowrates compared with the fan or compressor incorporating the blades with the conventional airfoil.

(106) Although the airfoil only the leading edge curve of which is formed as a Bezier curve has been described above, not only the leading edge curve but also the trailing edge curve can be formed as a Bezier curve.

(107) Although the blade according to the present disclosure is particularly advantageous when the blade is used as a rotor blade or stator vane of a fan or compressor of a turbo-fan engine of an aircraft which has a high blade surface Mach number and experiences a significant change in incidence (inflow angle) during operation, the application of the blade is not limited thereto. The blade according to the present disclosure can be widely used as a rotor blade or stator vane of an axial flow fluid machine, such as a compressor of a gas turbine other than the turbo-fan engine for an aircraft, or a fan, a compressor or a turbine as a stand-alone device.

EXPLANATION OF REFERENCE SIGNS

(108) A airfoil LE leading edge CL leading edge curve TE trailing edge CT trailing edge curve PS pressure surface CP pressure surface curve SS suction surface CS suction surface curve TP pressure surface connection point TS suction surface connection point X0 point ahead of blade X1 first control point X2 second control point X3 third control point X4 fourth control point X5 fifth control point X6 sixth control point