동요하는 날개를 이용한 파랑 중 선박의 추력발생

Propulsion by Oscillating Foil Attached to Ship in Waves

Article information

J. Ocean Eng. Technol. 2013;27(1):31-36
최윤락
Corresponding author Yoon-Rak Choi: +82-52-259-2158, yrchoi@ulsan.ac.kr
Received 2012 November 02; Accepted 2013 February 14; Revised 2012 December 18.

Trans Abstract

In this paper, the effects of a foil attached to a ship on the ship motion, added resistance, and thrust generation in waves are analyzed. The unsteady theory for a two-dimensional foil is introduced to determine the coupled motion responses of the ship and foil. The thrust caused by the oscillating foil is evaluated and compared to the added resistance of the ship, so that a positive net thrust can be possible in waves. A comparison of the results of unsteady, quasi-steady, and experimental analyses is performed.

1. 서 론

선박에 날개를 부착하여 파랑 중 선체운동을 감소시키거나 날개로부터 추력을 얻고자하는 연구가 수행되어 왔다. 이를 위해서는 동요하는 날개에 의한 양력 및 추력발생에 대한 연구와 날개에 의한 힘을 선체운동과 연성시켜 해석하는 연구가 필요하다. Bessho and Kyozuka(1984)는 준정적 2차원 날개 해석법을 사용하여 날개에 의한 선박의 동요저감과 부가저항의 감소에 대해 연구하였다. Nagahama et al.(1986)도 준정적 해석법을 사용하여 선박의 동요저감 및 날개에 의한 추력발생을 해석하였고 모형시험을 수행한 바 있다. 최근에 Choi(2012)는 비정상 2차원 날개이론에 근거하여 전파하는 수직변동유장 중 동요하며 전진하는 날개에 의한 양력과 추력을 해석하였다.

본 연구에서는 Choi(2012)의 비정상 날개해석법을 사용하여 날개운동과 선박운동의 연성방정식을 도출하고, 날개에 의한 선박운동의 영향을 고찰한다. 날개영향을 고려한 선박과 날개의 운동응답을 이용하여 선박에 작용하는 파랑 부가저항과 날개에서 발생하는 추력을 해석하여 파랑에너지로부터 부가저항보다 큰 추력을 획득할 수 있는지를 검토한다. 이와 더불어 준정적 해석법과 비정상 해석법의 결과를 비교검토한다.

2. 선박과 날개의 연성 운동 방정식

파랑 중 선박은 파랑의 기진력으로 인해 6자유도 운동을 하게 된다. 본 연구에서는 Fig. 1과 같은 병진운동(η1, η2, η3)과 회전 운동(η4, η5, η6)을 정의하여 해석을 수행한다. 그림에서 선박 의 평균위치에 고정된 X축과 Y축이 이루는 평면은 정수면과 일치하고 양의 Z축은 선박의 무게중심을 지나는 수직상방의 축이다. 선박은 음의 X축방향으로 U의 속도로 전진하며 원주파수 ν의 입사파는 음의 X축과 χ의 각도를 가지고 전파한다. 이때 파장 λ는 무한수심에서의 분산관계식으로 결정된다. 이러한 입사파와 선속에 의해 선박은 다음과 같은 조우원주파수(Encounter circular frequency)로 조화운동을 하게 된다.

Fig. 1

Coordinate system for incident waves and ship motion

여기서 K는 입사파의 파수이며 ξj는 각 운동모드별 복소진폭이다.

이러한 선박에 Fig. 2 와 같이 선박 좌표계 원점으로부터 선수방향으로 거리 p, 수심 d의 위치에 날개가 부착되어 있다. 날개의 길이(Chord length)는 l이며, 부착점의 위치는 날개 중앙점(Midchord)으로부터 앞쪽으로 a의 거리에 위치하고 있다. 그리고 그림에서 xG는 날개의 무게중심의 위치이다. 날개는 선박의 상하동요, 종동요 그리고 파랑에 의해 날개 부착점에서 상하운동(q)과 회전운동(α)을 한다. 이때 회전운동은 회전 스프링(kp)에 의해 복원력을 가진다.

Fig. 2

Position of foil and coordinate system for foil

날개에 작용하는 유체력 및 날개운동을 해석하기 위해 다음과 같은 가정을 도입한다.

날개는 조우원주파수로 조화운동을 하며 운동의 진폭은 작다.

날개의 두께는 길이에 비해 매우 작으며 상하 대칭인 형상이다.

파랑에 의한 유체속도는 선속, 즉 날개의 파랑 중 전진속도에 비해 크기가 매우 작다.

날개와 자유표면과의 상호작용은 무시한다. 파랑에 의한 수직 변동류의 크기는 날개가 잠긴 수심에서의 수직방향 유체속도를 이용한다.

날개 단면에 대해 2차원적 해석을 수행하며 날개폭(Span, s)의 효과는 날개폭에 수정계수를 곱하여 보정한다. 이때 날개단면에 대한 2차원적 해석은 선체중심선(Center line, Z=0)에 대하여 수행한다.

이상의 가정에 따라 날개의 평균 캠버선(Mean camber line)에 대하여 선형 포텐셜이론을 적용할 수 있다(Choi, 2012). 이때 평균 캠버선은 상하 대칭인 날개이므로 평판날개에 해당한다. 그리고 날개의 조화운동은 다음과 같이 표현된다.

여기서 q0α0는 각각 날개의 상하운동과 회전운동의 복소진폭이다.

Choi(2012)는 전파하는 수직변동 유장 중 전진하며 동요하는 날개에 대한 양력 및 추력을 해석하였는데 선형중첩원리에 의해 정수 중 전진하며 동요하는 평판날개문제와 전파하는 변동 유장 중 전진하는 고정된 평판날개문제로 분리하여 선형문제를 해석하였다.

여기서 LM은 각각 양력과 회전점에 대한 모멘트이고, 밑첨자 mg는 각각 정수 중 동요하는 날개와 전파하는 수직변동 유장 중 고정된 날개에 의한 성분들을 의미한다.

먼저, 정수 중 동요하는 평판날개에 의한 양력과 모멘트는 식 (3)의 조화운동에 따라 다음과 같이 표현된다(Choi, 2012).

여기서 *는 공액복소수를 의미하고, ρ는 유체의 밀도이며, C(k)는 날개후류와도(Wake vortex)에 의한 시간기억효과를 나타내는 Theodorsen 함수이다(Theodorsen, 1935). k는 Reduced frequency로 다음과 같이 정의된다.

또한 식 (5)와 (7)에서 변위 위의 두 점 및 한 점은 각각 변위 방향의 가속도와 속도를 의미한다. 그리고 aqaα는 각각 평판 날개의 상하운동에 의한 부가질량과 회전운동에 의한 부가관성 질량모멘트이다.

이들 식에서 사용된 는 2차원 물리량들을 3차원의 물리량들로 확장하기 위해 곱해진 유효 날개폭으로 다음과 같이 가정하였다(Nagahama et al., 1986).

여기서 AR은 날개의 종횡비(Aspect ratio)로 AR = s/l이다.

전파하는 수직변동유장 중 전진하는 2차원 고정날개문제에서 입사각 χ로 전파하는 파랑은 날개길이방향(Chordwise direction)의 파장성분(λx)을 가진 파랑으로 해석할 수 있다.

파랑의 날개길이방향 전파속도와 날개의 전진속도간의 상대속도를 고려하여 Choi(2012)의 해석결과에 적용하면 LgMg는 다음과 같이 표현되는데, 본 연구에서는 선수파 및 선수사파인 cosχ< 0인 경우에만 해석하였다.

여기서 J0J1는 제 1종의 Bessel 함수들이며 μμ=πl/λx이다. 그리고 W는 날개위치에서의 파랑에 의한 수직 변동유속의 복소크기이다.

여기서 A는 입사파의 복소 파진폭이다.

이상의 날개에 의한 양력과 모멘트를 선박의 운동방정식과 연성시켜서 운동응답을 구할 수 있는데 날개의 상하운동(q)은 선박의 상하동요(η3)와 종동요(η5)의해 결정된다.

본 연구에서는 상하동요와 종동요 그리고 날개의 회전운동을 서로 연성시켜 해석하였다.

여기서 mijbij는 각각 선박의 부가질량과 파랑감쇠계수이고, Fexi는 선박에 작용하는 파강제력이다. 그리고 mf는 날개의 질량이며 If는 날개 회전점에 대한 날개의 관성질량모멘트이다. 식 (17)-(19)에서 LmcMmc는 식 (6)과 (8)에서 보듯이 선박의 상하 동요, 종동요 그리고 날개의 회전운동에 대한 선형함수이다. 따라서 식 (17)-(19)은 선형연립방정식이며 이 연립방정식의 해로부터 선박의 상하동요 및 종동요 그리고 날개의 회전운동을 구할 수있다.

3. 날개에 의한 추력

날개에 의한 시간평균 추력()은 2차항의 힘으로 Choi(2012)의 연구결과를 사용하여 추정할 수 있다. 추력은 6개의 성분으로 표현되며 날개 회전운동, 날개 상하운동 그리고 수직 변동유속의 크기를 사용하여 2차전달함수로 표현된다(Choi, 2012). cosχ<0인 경우 추력은 다음과 같다.

여기서 FG는 각각 Theodorsen 함수의 실수부 F(k)와 허수부G(k)를 의미한다. 그리고 θq , θα 그리고 θW는 다음과 같은 위상각이다.

식 (17)-(19)의 해와 식 (16)의 관계로부터 α0q0를 구하고 식 (15)의 W를 사용하면 식 (20)의 추력을 구할 수 있다.

4. 파랑에 의한 부가저항

선박에 날개를 부착함으로써 파랑 중에 추가적인 추력을 얻기 위해서는 파랑에 의해 선박에 작용하는 부가저항을 고려하여야 한다. 이는 날개에 의한 추력으로 선박의 파랑 중 부가저항을 일정부분 상쇠하는 효과도 있지만, 파랑에너지를 이용한 추진시스템으로 이라는 능동적 관점에서는 부가저항 보다 큰 날개추력이 요구된다.

본 연구에서는 Gerritsma and Beukelman(1972)의 방법에 따라 선박에 작용하는 파랑부가저항(RAW)을 추정하였다.

여기서 적분은 선미로부터 선수방향으로의 적분이고 VHVS는 각각 선체 단면의 파랑에 대한 수직방향 및 좌우방향 상대속도이다.

여기서 ζω는 입사파의 자유표면이고 VωY는 입사파에 의한 Y방향유속이다.

5. 해석 예 및 토의

본 연구결과를 적용하기 위해 내항성 분야의 공시선형으로 알려진 ITTC S-175 컨테이너선을 채택하여 해석을 수행하였다. 이 선박의 수선간 길이(Lpp)는 175m이며 폭(B)는 25.4m이다. S-175선박의 자세한 제원은 Fujii and Takahashi(1975)에 나와 있다. 그리고 이 선박에 부착한 날개의 제원은 Table 1과 같다.

Table 1

Principal particulars of foil

Table 1의 날개는 Nagahama et al.(1986)이 S-175선박에 부착하여 실험한 날개와 동일하며 본 연구의 해석결과와 비교하기 위해 채택되었다. 날개가 잠긴 수심은 17.24m이며, a값이 날개 길이의 절반보다 크므로 회전점이 날개의 앞전 전방에 위치한다(Fig. 2 참조). 본 연구결과에서 비교하는 실험결과는 모두 Nagahama et al.(1986)의 결과이다.

Fig. 3Fig. 4에 각각 입사파의 파진폭(A)으로 무차원화된 상하동요와 종동요 운동응답을 도시하였는데 이때 후루드수는 0.15이고 파향은 선수파이다. 본 연구에서는 스트립법을 사용하여 부가질량, 파랑감쇠계수 그리고 파강제력을 계산하였다. 그림들을 보면, 날개에 의해 계산치 및 실험치 모두 상하동요와 선수동요의 크기가 감소한다. 따라서 선박에 부착된 날개가 선박의 동요를 저감시킬 수 있다. 이러한 운동저감효과는 Bessho and Kyozuka(1984)에서도 찾아볼 수 있다. 계산치에서 Q.S.로 표기된 준정적 해석값은 Theodorsen 함수의 실수부를 F(k) = 1, 허수부를 G(k) = 0으로 해석한 결과로 날개후류와도에 의한 시간기억효과를 무시한 경우이다(Newman, 1977). 선형응답인 상하동요 및 종동요에서는 준정적 해석과 비정상 해석과의 차이가 매우 적다.

Fig. 3

Heave responses of the ship in head waves

Fig. 4

Pitch responses of the ship in head waves

Fig. 5에 선수파 중 날개의 회전운동응답을 도시하였다. 실험 결과와 계산결과사이에 상당한 차이를 보이나 정성적인 경향은 유사하다. 선형물리량인 날개 회전운동 역시 준정적 해석과 비정상 해석간의 차이는 매우 적다.

Fig. 5

Responses of rotational motion of foil in head waves

본 연구에서 도입한 Gerritsma and Beukelman(1972)의 부가저항 해석법의 유용성을 살펴보기 위해 날개가 없을 때의 부가저항 해석결과를 실험치와 비교하였다 (Fig. 6). 부가저항이 2차항의 힘임에도 불구하고 계산결과가 실험결과와 비교적 잘 일치하여 본 연구에서 사용한 부가저항해석법이 유용함을 알 수 있다.

Fig. 6

Added resistances of the ship (w/ foil)

Fig. 7에 날개를 부착하였을 때의 선수파 중 부가저항을 도시하였다. 날개부착에 의해 부가저항이 상당히 감소함을 알 수 있다. 이는 Fig. 34에서 보듯이 날개에 의해 선박의 상하동요와 종동요가 감소하여 부가저항에 가장 큰 영향을 미치는 수직상대운동이 감소하였기 때문이다. 선박의 운동은 준정적 해석과 비정상 해석간의 차이가 거의 없으므로 부가저항 역시 차이가 매우 적다.

Fig. 7

Added resistances of the ship w/o and w/ foil in head waves

후루드수가 0.15이고 선수파 중인 경우에서 날개에 의한 추력(T), 선박에 작용하는 부가저항(RAW) 그리고 순저항(RAW - T)을 Fig. 8에 도시하였다. 그림에서 순저항이 음의 값을 가지면 날개에 의한 추력이 선박에 작용하는 부가저항을 이기고 선박에 부가적인 추력을 제공함을 의미하는데, λ/Lpp가 약 0.8이상의 파랑범위에서는 이러한 상태가 됨을 알 수 있다. λ/Lpp이 약 0.8이상에서는 앞의 결과들에서 알 수 있듯이 선박의 운동과 날개의 회전운동이 큰 범위에 해당하며, 이로부터 파랑과 날개의 운동에 의해 상당량의 추력이 발생함을 알 수 있다. 준정적 해석과 비정상 해석의 결과를 비교하면 준정적 해석은 추력을 과도하게 추정하여 실험치와의 정량적 일치도가 비정상 해석결과보다 떨어진다. 선박의 동요 및 날개의 회전운동 등 선형물리량에서는 두 방법의 차이가 거의 없음을 고려하면 2차항 힘인 추력에서는 후류와도에 의한 시간기억효과의 영향이 큼을 알 수 있다.

Fig. 8

Added resistance, thrust, and net resistance in head waves (Fn=0.15)

후루드수가 0.2이고 선수파 및 선수사파 중인 경우, 순저항값의 비정상 해석 결과를 Fig. 9에 도시하였다. 선수파에서 횡파 방향으로 파랑 입사각도가 변하면 추력획득효과가 감소함을 알 수 있다.

Fig. 9

Net resistances in various heading angles (Fn=0.2)

6. 결 론

본 연구에서는 파랑 중 선박에 부착된 동요하는 날개가 선박의 운동에 미치는 영향을 해석하였고, 날개에 의한 추력을 추정하였다. Choi(2012)의 연구결과를 토대로 파랑 중 동요하는 날개에 의한 양력 산출식을 선박의 운동방정식과 연성시켜 선박의 운동 응답과 날개의 운동을 해석하였다. 이를 토대로 날개에서 발생하는 추력을 추정하였고, 파랑 중 선박에 작용하는 파랑부가저항과 비교하여 파랑에너지로부터 순추력의 획득가능성을 검토하였다.

S-175선박에 특정 날개를 적용하여 해석해 본 결과 날개 부착에 의해 선박의 상하동요와 종동요가 감소하였고, 운동의 감소에 따라 부가저항도 감소하였다. 날개에 의한 추력과 선박에 작용하는 파랑부가저항의 비교를 통해 운동응답이 큰 파랑영역에서 순추력이 발생함을 확인하였다. 이러한 결과들은 실험치에서도 확인할 수 있었다. 준정적 해석법과 비정상 해석법 결과들을 비교해보면 선형운동응답 및 선박부가저항에서는 거의 동일한 결과를 보이는데 반하여 날개의 추력에서는 상당한 차이를 보이고 있으며 비정상 해석 결과가 실험치에 더 근접한 결과를 보였다. 따라서 추력해석을 위해서는 후류와도에 의한 시간기억효과를 엄밀히 고려하는 비정상 해석이 요구된다.

본 연구에서는 2차원 해석을 통한 날개효과를 3차원으로 확장하는데 있어 날개 종횡비를 이용한 약산식을 사용하였으나, 좀 더 엄밀한 해석을 위해서는 날개의 3차원효과를 고려한 양력 및 추력해석법이 필요하리라 생각되며, 날개와 자유표면의 상호작용도 연구할 필요가 있다. 또한 날개의 회전운동이 실험치와 큰 차이를 보이는 원인에 대해 심도있는 검토가 필요하리라 사료된다.

Notes

It is noted that this paper is revised edition based on proceeding of Annual Autumn Meeting of KSOE 2012 in Busan

Acknowledgements

본 연구는 '울산대학교 조선해양공학부 세계일류화 사업(SOTOP)'의 지원을 받아 수행되었음.

References

Bessho, M., Kyozuka, Y., 1984. On the Ship Motion Reduction by Anti-Pitching Fins in Head Seas. Proceeding of 15th Symposium on Naval Hydrodynamics, Hamburg W. Germany, 109-118.

Bessho M., Kyozuka Y.. On the Ship Motion Reduction by Anti-Pitching Fins in Head Seas In : Proceeding of 15th Symposium on Naval Hydrodynamics. Hamburg W. Germany; 1984. p. 109–118.

Choi, Y.-R., 2012. Thrust Caused by Oscillating Two-Dimensional Hydrofoil Moving in Propagating Unsteady Flow Field. Journal of Ocean Engineering and Technology, 26 (5), 40-46.

Choi Y.-R.. Thrust Caused by Oscillating Two-Dimensional Hydrofoil Moving in Propagating Unsteady Flow Field. Journal of Ocean Engineering and Technology 2012;26(5):40–46. 10.5574/KSOE.2012.26.5.040.

Fujii, H., Takahashi, T., 1975. Experimental Study on the Resistance Increase of a Ship in Regular Waves. Proceeding of 14th ITTC, 4, 351-360.

Fujii H., Takahashi T.. Experimental Study on the Resistance Increase of a Ship in Regular Waves In : Proceeding of 14th ITTC; 1975. 4p. 351–360.

Gerritsma, J., Beukelman, W., 1972. Analysis of the Resistance Increase in Waves of a Fast Cargo Ship. International Shipbuilding Progress, 19(217), 285-293.

Gerritsma J., Beukelman W.. Analysis of the Resistance Increase in Waves of a Fast Cargo Ship. International Shipbuilding Progress 1972;19(217):285–293.

Nagahama, M., Murakami, M., Isshiki, H., 1986. Effects of a Foil Attached to a Ship in Waves. Hitachi Zosen Technical Review, 47(3-4), 38-43.

Nagahama M., Murakami M., Isshiki H.. Effects of a Foil Attached to a Ship in Waves. Hitachi Zosen Technical Review 1986;47(3-4):38–43.

Newman, J.N., 1977. Marine Hydrodynamics. The MIT Press, Cambridge, MA.

Newman J.N.. Marine Hydrodynamics The MIT Press. Cambridge, MA: 1977.

Theodorsen, T., 1935. General Theory of Aerodynamic Instability and the Mechanism of Flutter. N.A.C.A. Report, No 496.

Theodorsen T.. General Theory of Aerodynamic Instability and the Mechanism of Flutter. N.A.C.A. 1935.

Article information Continued

Fig. 1

Coordinate system for incident waves and ship motion

Fig. 2

Position of foil and coordinate system for foil

Table 1

Principal particulars of foil

Table 1

Fig. 3

Heave responses of the ship in head waves

Fig. 4

Pitch responses of the ship in head waves

Fig. 5

Responses of rotational motion of foil in head waves

Fig. 6

Added resistances of the ship (w/ foil)

Fig. 7

Added resistances of the ship w/o and w/ foil in head waves

Fig. 8

Added resistance, thrust, and net resistance in head waves (Fn=0.15)

Fig. 9

Net resistances in various heading angles (Fn=0.2)