1 式(5.2.38)
5. CFD を用いたプロペラ性能評価
68
69
Fig.5.1.3 Grid distribution on propeller surface.
Fig.5.1.4 Velocity distribution
keyblade
0.50
0.50 0.50
0.50
0.50
0.50
0.50 0.55 0.55
0.55
0.55
0.55
0.55
0.55 0.60
0.60
0.60 0.60
0.60
0.60 0.60 0.60
0.60
0.60
0.60
0.60
Y
Z
-0.1 -0.05 0 0.05 0.1
-0.1 -0.05 0 0.05 0.1
0.50
0.50 0.50
0.50
0.50
0.50
0.50 0.55 0.55
0.55
0.55
0.55
0.55 0.60
0.60 0.60
0.60
0.60
0.60
0.60 0.60 0.60
0.60
0.60
0.60
0.60
Y
Z
-0.1 -0.05 0 0.05 0.1
-0.1 -0.05 0 0.05 0.1
70
5. CFD を用いたプロペラ性能評価
5.2.4 圧力-速度連成
Fig.5.2.1 Flowchart of pressure based solver
Fig.5.2.2 Relationship of coordinate system.
Updating physical property.
Sequentially solving for each property .
(u, v, w, k…)
Solving equation of pressure correction.
Updating mass flux, pressure and velocity
Solving scalar equation.
(Turbulance, energy…)
Converge?
Finish calculation.
No
Yes
X Y
Z
O
x y
z
o
P
r0
r
Absolute coordinate system
Moving coordinate system
71
6. 最適化システムの有効性確認
6.2 最適化計算結果
6.2.1 プロペラ性能比較
Fig.6.2.1 Optimization history of SQP steps.
6.2.2 プロペラ形状比較
Fig.6.2.2 Comparison of pitch distribution between Case0 and Case1.
0.00%
0.50%
1.00%
1.50%
2.00%
2.50%
0 1 2 3 4 5 6 7 8 9 10 improvement of propeller open efficiecy(%from original prop.)
SQP Steps
0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.500 0.600 0.700 0.800
r/R
H/Dp
Case0 Case1
72
Fig.6.2.3 Comparison of camber distribution between Case0 and Case1.
Fig.6.2.4 Comparison of cord length distribution between Case0 and Case1.
Fig. 6.2.5 Comparison of blade section at 0.70R between Case0 and Case1.
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.000 0.020 0.040 0.060 0.080 0.100 0.120
r/R
f/C
Case0 Case1
0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.000 0.050 0.100 0.150 0.200 0.250 0.300
r/R
C/Dp
Case0 Case1
0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00
0.0 10.0 20.0 30.0 40.0 50.0 60.0
Y
X
Case0 Case0
73
(a) Case0 (b) Case1
Fig.6.2.6 Comparison of blade shapes between Case0 and Case1.
6.2.3 プロペラ表面圧力比較
(a) Case0 (b) Case1
Fig.6.2.7 Comparison of pressure distribution on back side between Case0 and Case1.
74
(a) Case0 (b) Case1
Fig.6.2.8 Comparison of pressure distribution on face side between Case0 and Case1.
6.3 水槽試験による効果の確認
6.3.1 水槽試験条件
(a) Case0 (b) Case1
Fig.6.3.1 Photos of model propeller.
75
6.3.2 水槽試験結果
Fig.6.3.2 Comparison of thrust coefficient between Case0 and Case1(Tank test result).
Fig.6.3.3 Comparison of torque coefficient between Case0 and Case1(Tank test result).
0.1500 0.1600 0.1700 0.1800 0.1900 0.2000 0.2100 0.2200 0.2300 0.2400 0.2500
0.200 0.250 0.300 0.350 0.400
KT
J
Case0 Case1
0.1900 0.2000 0.2100 0.2200 0.2300 0.2400 0.2500 0.2600 0.2700
0.200 0.250 0.300 0.350 0.400
10KQ
J
Case0 Case1
76
Fig.6.3.4 Comparison of propeller open efficiency between Case0 and Case1 (Tank test result).
0.3000 0.3500 0.4000 0.4500 0.5000 0.5500
0.200 0.250 0.300 0.350 0.400
o
J
Case0 Case1
77
7. 伴流中最適化
7.1 供試船型
Fig.7.1.1 Body plan of Ship A
Fig.7.1.2 Wake distribution at propeller position of ship A.
C.L.
W.L.
0.7
0.6
0.6
0.6 0.6 0.7
0.7
0.7
0.5 0.5
0.4 0.4
0.3 0.3
Y
Z
-100 -50 0 50 100
-100 -50 0 50 100
78
Fig. 7.1.3 Averaged nominal wake distribution of ship A.
7.2 母型プロペラ
Fig.7.2.1 MAU design result for ship A at each design point.
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.000 0.200 0.400 0.600 0.800 1.000
r/R
1-Wn
0.1780 0.1790 0.1800 0.1810 0.1820 0.1830 0.1840 0.1850 0.1860 0.1870 0.1880
350.0 400.0 450.0 500.0 550.0 600.0
BHP[PS]
Rotation Speed at Design Point [RPM]
Case0
Best Design Point
79
7.3 伴流中プロペラ性能計算
Fig.7.3.1 Estimation result of thrust characteristic of case0.
Fig.7.3.1 Estimation result of torque characteristic of case0.
0.2000 0.2050 0.2100 0.2150 0.2200 0.2250 0.2300 0.2350 0.2400 0.2450 0.2500
0.150 0.200 0.250 0.300 0.350
KT
J
Case0 (Uniform flow) Case0 (Wake flow)
0.2250 0.2300 0.2350 0.2400 0.2450 0.2500 0.2550 0.2600 0.2650 0.2700 0.2750
0.150 0.200 0.250 0.300 0.350
10KQ
J
Case0 (Uniform flow) Case0 (Wake flow)
80
Fig.7.3.1 Estimation result of propeller open efficiency characteristic of case0.
7.5 最適化計算結果
Fig.7.5.1 Optimization history of SQP steps. (From Case0 to Case3)
7.5.1 プロペラ形状比較
0.2500 0.2700 0.2900 0.3100 0.3300 0.3500 0.3700 0.3900 0.4100 0.4300 0.4500
0.150 0.200 0.250 0.300 0.350
O
J Case0 (Uniform flow) Case0 (Wake flow)
0.0%
0.5%
1.0%
1.5%
2.0%
2.5%
3.0%
0 2 4 6 8 10 12 14
Improvement of propeller performance from case0
SQP Step
Case0→Case2 Case2→Case3
81
Fig.7.5.2 Comparison of pitch distribution between Case0, case2 and Case3.
Fig. 7.5.3 Comparison of camber distribution between Case0, case2 and Case3.
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.550 0.600 0.650 0.700 0.750 0.800
r/R
H/D
Case0 Case2 Case3
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
0.000 0.020 0.040 0.060 0.080 0.100 0.120 0.140
r/R
f/C
Case0 Case2 Case3
82
Fig. 7.5.4 Comparison of thickness reduction ratio distribution between Case0, case2 and Case3.
Fig. 7.5.5 Comparison of blade section at 0.70R between Case0, Case2 and Case3.
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
-100.0% -50.0% 0.0% 50.0% 100.0%
r/R
Modificate Ratio of Trailingside Thickness
Case0 Case2 Case3
0.00.5 1.01.5 2.02.5 3.0 3.54.0 4.5 5.0
0.0 10.0 20.0 30.0 40.0 50.0 60.0
Y
X
Case0 Case2 Case3
83
7.5.2 プロペラ性能比較
Fig.7.5.6 Comparison of thrust coefficient between Case0, case2 and Case3.
Fig.7.5.7 Comparison of torque coefficient between Case0, case2 and Case3.
0.2000 0.2050 0.2100 0.2150 0.2200 0.2250 0.2300 0.2350 0.2400 0.2450 0.2500
0.150 0.200 0.250 0.300 0.350
KT
J
Case0 Case2 Case3
0.2250 0.2300 0.2350 0.2400 0.2450 0.2500 0.2550 0.2600 0.2650 0.2700 0.2750
0.150 0.200 0.250 0.300 0.350
10KQ
J
Case0 Case2 Case3
84
Fig.7.5.8 Comparison of propeller open efficiency between Case0, case2 and Case3.
(a) Case0 (b) Case2 (c) Case3
Fig.7.5.9 Comparison of pressure distribution on back side.
0.2500 0.2750 0.3000 0.3250 0.3500 0.3750 0.4000 0.4250 0.4500
0.150 0.200 0.250 0.300 0.350
ηo
J
Case0 Case2 Case3
85
(a) Case0 (b) Case2 (c) Case3
Fig.7.5.10 Comparison of pressure distribution on face side.
7.6 水槽試験結果
(a) Case0 (b) Case3
Fig.7.6.1 Photo of model propellers.
86
7.6.1 伴流計測結果
Fig.7.6.1 Wake measurement results of ship A.
7.6.2 プロペラ単独試験結果
Fig.7.6.2 Comparison of thrust coefficient between Case0 and Case3. (Tank test)
0.1 0.2
0.2
0.3
0.3 0.4
0.4 0.5
0.5 0.6
0.6 0.7
0.6
0.6 0.7
0.6
0.7 0.7
0.1 0.2
0.2 0.3
0.3
0.4
0.4
0.5
0.5 0.6
0.6 RESULTS OF WAKE SURVEY
: MSNO CONDITION
VM
D.FULL (with STUD) Fn
MEASURING POSITION Wn
1637
0.1559 ( Vs= 17.13 knots ) 1.379 m/s
151.6 mm fore from A.P :
: : : :
FIG.3.2.4 WAKE CONTOURS
1 STEP = 0.10 20.0 mm
PORT STBD
Boss
Prop. Disk
CL
0.570
0.1
0.2
0.2 0.3
0.3 0.4
0.4 0.5
0.5
0.5
0.6
0.6
0.5 0.6 0.6
0.7 0.8
0.5 0.3
0.1 0.2
0.2 0.3
0.3
0.4
0.4
0.5
0.5 0.6
0.6 RESULTS OF WAKE SURVEY
: MSNO CONDITION
VM
D.FULL Fn
MEASURING POSITION Wn
1637
0.1760 ( Vs= 19.34 knots ) 1.557 m/s
151.6 mm fore from A.P :
: : : :
FIG.3.4.1 WAKE CONTOURS
1 STEP = 0.10 20.0 mm
PORT STBD
Boss
Prop. Disk
CL
-0.2000 0.2100 0.2200 0.2300 0.2400 0.2500
0.150 0.200 0.250 0.300 0.350
KT
J
Case0 Case3
87
Fig.7.6.3 Comparison of torque coefficient between Case0 and Case3. (Tank test)
Fig.7.6.4 Comparison of propeller open efficiency between Case0 and Case3. (Tank test) 0.2200
0.2300 0.2400 0.2500 0.2600 0.2700
0.150 0.200 0.250 0.300 0.350
10KQ
J
Case0 Case3
0.2000 0.2500 0.3000 0.3500 0.4000 0.4500 0.5000
0.150 0.200 0.250 0.300 0.350
ηo
J
Case0 Case3
88
7.6.3 自航試験結果
Fig.7.6.5 Total resistance curve of ship A
Fig.7.6.6 Self-propulsion test results. (O) 3.000
3.500 4.000 4.500 5.000 5.500 6.000
0.080 0.100 0.120 0.140 0.160 0.180 0.200 CTM
Fn
0.3000 0.3250 0.3500 0.3750 0.4000 0.4250 0.4500 0.4750 0.5000
0.130 0.140 0.150 0.160 0.170 0.180
O
Fn
Case0 Case3
89
Fig.7.6.7 Self-propulsion test results. (R)
Fig.7.6.8 Self-propulsion test results. (1-WQ)
Fig.7.6.9 Self-propulsion test results. (1-t) 0.9000
0.9250 0.9500 0.9750 1.0000 1.0250 1.0500 1.0750 1.1000
0.130 0.140 0.150 0.160 0.170 0.180
R
Fn
Case0 Case3
0.4000 0.4250 0.4500 0.4750 0.5000 0.5250 0.5500 0.5750 0.6000
0.130 0.140 0.150 0.160 0.170 0.180 1-WQ
Fn
Case0 Case3
0.7000 0.7250 0.7500 0.7750 0.8000 0.8250 0.8500 0.8750 0.9000
0.130 0.140 0.150 0.160 0.170 0.180
1-t
Fn
Case0 Case3
90
8. キャビテーション性能を考慮した最適化
8.1 Case2 のキャビテーション性能
Fig.8.1.1 Pressure time series at reference point on r/R=0.900.
Fig.8.1.2 Pressure time series at reference point on r/R=0.800.
-2.2000 -2.0000 -1.8000 -1.6000 -1.4000 -1.2000
0 30 60 90 120 150 180 210 240 270 300 330 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.900
Case0 Case2
-2.5000 -2.3000 -2.1000 -1.9000 -1.7000 -1.5000 -1.3000 -1.1000
0 30 60 90 120 150 180 210 240 270 300 330 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.800
Case0 Case2
91
Fig.8.1.3 Pressure time series at reference point on r/R=0.700.
8.2 圧力最小値の制約条件
Fig.8.2.1 Pressure time series at reference point on r/R=0.900 with n. -2.2000
-2.0000 -1.8000 -1.6000 -1.4000 -1.2000
0 30 60 90 120 150 180 210 240 270 300 330 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.700
Case0 Case2
-2.2000 -2.0000 -1.8000 -1.6000 -1.4000 -1.2000
0 45 90 135 180 225 270 315 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.900
Case0 Case2 σ
92
Fig.8.2.2 Pressure time series at reference point on r/R=0.800 with n.
Fig.8.2.3 Pressure time series at reference point on r/R=0.800 with n.
8.3 圧力制約を付加した最適化計算結果
-2.2000 -2.0000 -1.8000 -1.6000 -1.4000 -1.2000
0 45 90 135 180 225 270 315 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.800
Case0 Case2 σ
-2.2000 -2.0000 -1.8000 -1.6000 -1.4000 -1.2000
0 45 90 135 180 225 270 315 360
Pressure Coeficient
Rotation Angle[deg]
r/R=0.700
Case0 Case2 σ
93
Fig.8.3.1 SQP history. ( with pressure constraint)
8.3.1 翼表面圧力
Fig.8.3.2 Pressure time series of at reference point on r/R=0.900 with n.