About | Browse Archives | Editorial Board | For Contributors |
Sorry.
You are not permitted to access the full text of articles.
If you have any questions about permissions,
please contact the Society.
죄송합니다.
회원님은 논문 이용 권한이 없습니다.
권한 관련 문의는 학회로 부탁 드립니다.
[ Original Article ] | |
Journal of the Korean Society for Marine Environment and Energy - Vol. 21, No. 1, pp. 10-22 | |
Abbreviation: J. Korean Soc. Mar. Environ. Energy | |
ISSN: 2288-0089 (Print) 2288-081X (Online) | |
Print publication date 25 Feb 2018 | |
Received 16 Oct 2017 Revised 27 Nov 2017 Accepted 29 Nov 2017 | |
DOI: https://doi.org/10.7846/JKOSMEE.2018.21.1.10 | |
편심된 회전축을 갖는 수평 원기둥 파력발전장치의 성능해석 | |
1제주대학교 해양시스템공학과 | |
2제주대학교 풍력공학부 | |
Performance Analysis of a Horizontal Cylinder Wave Energy Converter with Off-centered Rotational Axis | |
1Department of Ocean System Engineering, Jeju National University, Jeju 63243, Korea | |
2Faculty of Wind Energy Engineering, Jeju National University, Jeju 63243, Korea | |
Correspondence to : †cho0904@jejunu.ac.kr | |
Funding Information ▼ |
횡 파로부터 에너지를 추출하는 수평 원기둥 파력발전장치의 성능은 입사파의 주파수, 원기둥의 크기와 잠긴 깊이, 회전축의 위치, 파력발전장치내의 질량분포 등에 영향을 받는다. 본 연구에서는 원기둥의 잠긴 깊이와 회전축 위치 변화에 따른 수평 원기둥 파력발전장치의 횡 운동변위와 추출파워를 주파수영역에서 살펴보았다. 수평 원기둥에 작용하는 유체력(동유체력, 파기진력)을 구하는 수치해석 방법으로 WAMIT 상용코드를 사용하였다. 이때 편심된 회전축에서의 유체력은 원기둥 중심에서 구한 유체력을 축 변환하여 얻어진다. 회전축의 위치 변화는 횡 운동모드의 고유주파수에 영향을 준다. 따라서 설치해역의 입사파의 주파수 특성에 맞게 회전축의 위치를 적절히 선택하면 공진을 일으켜 추출파워를 극대화시킬 수 있다.
The performance of a horizontal cylinder wave energy converter (WEC) for extracting wave energy in beam sea is affected by the wave frequency, the size and submerged depth of the cylinder, the position of the rotational axis, and the mass distribution in the WEC. In this study, the extracted power and roll motion of the horizontal cylinder WEC in beam sea was examined according to the change of the submergence depth and rotational axis of WEC in the frequency domain. WAMIT commercial code was used as a numerical method to obtain the fluid forces (hydrodynamic forces, wave exciting forces) acting on the horizontal cylinder WEC. The fluid forces at the off-centered rotational axis is easily obtained by using axis-transformation formulas from the fluid forces at the geometric center of the cylinder. The change of the position of the rotational axis affects the natural frequency of the rolling motion mode. Therefore, if the position of the rotational axis is appropriately selected according to the wave characteristics at the installation site, the resonance can be generated to maximize the extraction power.
Keywords: Wave Energy Converter, Roll Motion, Power Take-Off, Off-centered Rotational Axis, Extracted Power 키워드: 파력발전장치, 횡운동, 파워추출장치, 편심된 회전축, 추출 파워 |
Salter(1974)가 1974년도 Nature 학술지에 제안한 1차 변환장치는 운동 특성이 수면에 떠 있는 오리의 움직임과 같다고 하여 Salter nodding duck 또는 Edinburgh duck이라 불리면서 70년대 후반과 80년대 초반에 모형실험(Salter et al., 1975; Jeffrey et al., 1976)과 이론해석(Evans, 1976; Count, 1978, Mynett et al., 1979)등을 통하여 폭 넓게 연구되었다. Salter duck의 횡 운동에 의한 방사파(radiated waves)의 크기는 흡수된 에너지와 밀접한 관련이 있음을 밝혔다. 양쪽 방향으로 진행하는 방사파 중에서 입사파의 진행방향과 같은 방향의 방사파의 진폭(
최근에 덴마크에서 개발한 WEPTOS 파력발전시스템 내에 Salter duck 형상의 1차 변환장치를 도입한 이후 다시금 Salter duck 단면 형상이 세간의 관심을 받게 되었다. WEPTOS 파력발전시스템은 일정한 길이를 갖는 Salter duck 형상의 1차 변환장치를 V자형 부유구조물에 여러개 배열한 형태로 V자형 부유구조물을 일점 계류(turret mooring)시켜 그 점을 중심으로 회전케 하여(weather vane) 파랑하중을 줄일 수 있도록 고안되었다(Pecher and Kofoed, 2013, 2014). 또한 파랑조건에 따라 V자형 부유구조물의 열림 각(opening angle)을 바꾸면 파랑에너지의 입력값을 조절할 수 있어 입력 에너지의 변동성을 낮출 수도 있으며 극한 해상상태에서 파력발전시스템에 가해지는 충격력을 줄일 수 있는 장점을 가지고 있다.
Salter duck이 높은 추출 효율을 가지고 파랑 에너지를 흡수할 수 있다고 알려졌지만 단면 형상이 복잡하여 제작에 어려움이 따른다. 따라서 cylindrical duck이나 Bristol cylinder와 같이 형상이 단순한 수평 원기둥 1차 변환장치에 편심된 회전축을 두어 파랑에너지를 추출하는 연구가 몇몇 학자들의 관심을 받았다(Lucas et al., 2009). 이러한 원형 단면을 갖는 1차 변환장치의 성능에 중요한 설계변수로 원기둥의 크기, 잠긴 깊이, 폭, 회전축의 위치, 무게중심의 위치 등이며, 이 값들을 바꿔가면서 각 변수들이 추출파워에 미치는 영향을 분석하는 것이 중요하다.
본 연구에서는 원형 단면을 갖는 1차 변환장치의 여러 설계변수들 중에서 특히 잠긴 깊이와 회전축의 변화가 1차 변환장치의 성능에 미치는 영향에 대하여 자세히 살펴보았고, 해석결과를 불규칙파로 확장하여 파랑스펙트럼이 주어졌을 때 잠긴 깊이와 회전축의 위치에 따른 흡수 파워와 취득 폭의 대푯값을 구하였다. 수평 원기둥 1차 변환장치에 작용하는 횡 방향 유체력(동유체력, 정유체력, 파기진력)을 상용코드 WAMIT을 이용하여 구하였다. 수치계산 과정에서 임의의 회전축 위치에서 유체력을 직접 구하지 않고 원기둥 중심에서의 값들을 이용하여 회전축에서의 값들로 변환하는 방법을 사용하여 계산시간을 크게 단축시켰다(Cruz and Salter, 2006). 수평과 수직 운동이 구속된 상태에서 1자유도 횡 운동방정식을 주파수영역에서 유도하고 이를 풀어 횡 운동 RAO, 추출 파워, 취득 폭을 구하였다. 이때 시간평균 추출 파워가 최대값을 갖도록 최적의 PTO(power take-off)감쇠계수를 구하였다. 주파수영역에서의 해를 불규칙파로 확장하여 파랑스펙트럼이 주어져 있을 때 횡 운동스펙트럼과 파워스펙트럼을 구하고 각각의 스펙트럼의 면적을 적분하여 운동과 추출파워의 대푯값을 구하였다. 회전축의 위치를 바꿔가면서 수평 원기둥 1차 변환장치를 통하여 흡수되는 파워를 분석하여 최적의 회전축의 위치를 정하였다. 해석결과, 공진조건을 만족할 때 추출 파워는 최대가 되며 이때 회전축의 위치가 공진을 결정하는데 중요한 변수임을 확인하였다.
일정한 폭(W)을 갖는 반경 R인 수평 원기둥 1차 변환장치가 횡파 중에 놓여있다. 원기둥의 잠긴 깊이는 d이며, 입사파는 진폭 A, 주파수 ω을 갖고 x축의 양의 방향으로 진행한다. 좌표축의 원점은 원기둥의 중심(C)에 위치하며 z축을 연직상향으로 잡았다. 회전축의 중심은 R0이며, 원기둥 중심과 l0 떨어져 있다. 회전축 중심과 원기둥 중심간의 상대각도 α는 Fig. 1과 같으며 반시계 방향을 양의 방향으로 잡았다. 선형포텐셜 이론을 가정하여 속도포텐셜 Φ(x, y, z, t) = Re{φ(x, y, z)e−iωt}을 도입하면 전체포텐셜을 식 (1)과 같이 입사포텐셜(φ0), 산란포텐셜(φ4), 그리고 방사포텐셜(φj, j = 1,2,3)의 합으로 표현할 수 있다.
(1) |
여기서 ξj, (j = 1,2,3)는 수평(sway), 수직(heave), 횡(roll) 방향의 복소운동진폭이다.
수평 원기둥 1차 변환장치의 운동방정식은 Newton의 제 2법칙으로부터 다음과 같이 쓸 수 있다.
(2) |
여기서 ξj(t) = Re{ξje−iωt}, (j = 1,2, 3)이다.
식 (2)에서 Mji는 질량행렬의 요소를 나타내며 다음과 같다.
(3) |
여기서 m은 1차 변환장치의 전체 질량, (xg, zg)와 J는 각각 회전축(R0)에서의 무게중심의 좌표와 관성 모멘트이다. 1차 변환장치의 전체 질량(m)은 고정 질량(mH)과 가변 질량(mB)의 합과 같다.
(4) |
여기서 고정 질량은 일정한 두께를 갖는 내부가 빈 원기둥의 질량을 나타낸다. 따라서 mH의 무게중심의 좌표는 원기둥의 중심과 일치한다(xH, zH) = (0,0). 한편 가변 질량(mB)은 1차 변환장치가 일정한 잠긴 깊이를 갖도록 원기둥 내부 바닥에 설치하는 발라스트(ballast) 질량으로 무게중심의 위치는 발라스트의 위치에 따라 원기둥의 중심으로부터 이격된다. 이와 동일하게 전체 관성모멘트도 고정 질량에 의한 관성모멘트의 성분(JH)과 가변 질량에 의한 관성모멘트의 성분(JB)의 합으로 쓸 수 있다.
(5) |
평행축의 정리를 적용하여 회전축(R0)에서의 관성 모멘트(JH, JB)를 각각의 무게 중심의 위치인(xH, zH)와 (xR, zB)에서의 관성 모멘트 값(IH, IB)을 가지고 표현하면 다음과 같다.
(6) |
여기서 lH=l0과
식 (2)의 Fj(t)는 WEC에 작용하는 j(= 1,2,3)방향의 유체력으로 파기진력(
Added mass | Radiation damping | Hydrostatic force | Exciting force |
---|---|---|---|
(7) |
여기서 zB, zG는 원기둥의 중심에서의 부력 중심과 무게 중심의 z좌표이다. S00는 수선면적이며,
수평 원기둥 파력발전장치의 운동에너지로부터 전기에너지를 추출하기 위하여 2차 변환장치인 PTO장치를 설치하여야 한다. 이때 발생하는 PTO감쇠력은 운동속도와 선형적으로 비례한다고 가정하였다.
(8) |
여기서
1차 변환장치가 입사파의 주파수에 따라 조화운동을 한다고 가정한 뒤 동유체력, 정유체력, PTO감쇠력, 그리고 파기진력(
(9) |
여기서 δji는 Kronecker delta로 δji= 1, j = i, δji= 0, j≠i이다. Xj, aji, bji는 수평 원기둥 1차 변환장치의 잠긴 부분의 형상과 주파수의 함수로 회절과 방사 포텐셜로부터 아래와 같이 구할 수 있다.
(10) |
본 연구에서는 회절과 방사 포텐셜을 구하는 수치해석 기법으로 경계요소법(boundary element method)에 기반을 둔 상용코드 WAMIT을 이용하였다.
회전축에서 계산된 정유체력, 동유체력, 파기진력을 가지고 수평과 수직 운동을 구속한 상태에서 1자유도 횡 운동의 RAO(Response Amplitude Operator)를 구하면 아래 식과 같다.
(11) |
수평 원기둥 파력발전장치의 횡 운동을 통하여 얻는 시간평균 파워는 아래와 같다.
(12) |
윗식에서 시간평균 파워가 최대값을 갖도록 식(
(13) |
여기서 ωN는 횡 운동모드의 비감쇠 고유주파수(Undamped Natural Frequency)로 아래와 같다.
(14) |
식 (14)는 초기값
식 (13)을 식 (11)에 대입하면 PTO-on 상태에서의 1차 변환장치의 횡 운동 변위를 구할 수 있으며, 식 (12)에 대입하면 시간평균 추출파워의 최적 값 (
(15) |
여기서
식 (12)의 시간평균 추출파워로부터 최대 파워를 얻기 위해서는 최적의 PTO감쇠계수를 갖는 것 외에 입사파의 주파수와 횡 운동의 고유주파수가 서로 같다(ω = ωN)는 공진조건을 만족하여야 한다.
(16) |
공진조건을 만족할 때 식 (13)에서 최적의 PTO감쇠계수는
(17) |
불규칙파에 대한 결과로 확장하기 위하여 파랑스펙트럼으로 TMA 스펙트럼 을 사용하였다(Bouws et al., 1985). TMA 스펙트럼은 아래식과 같이 수심의 효과를 준다.
(18) |
여기서 Sj(ω)는 JONSWAP 스펙트럼으로 Goda (1988)에 의하면 다음과 같다.
(19) |
여기서 H1/3는 유의파고,
유한 수심의 효과를 나타내는 함수 F(ω*)는 아래와 같다.
여기서
파랑스펙트럼에 횡 운동 RAO의 제곱을 곱하여 횡 운동스펙트럼을 얻고, 그 면적을 적분하여 유의 횡 운동진폭(significant roll motion amplitude)을 아래와 같이 구한다.
(20) |
단위 폭당 불규칙 입사파가 지닌 파워는 다음과 같다.
(21) |
식 (15)로부터 구한 최적의 시간평균 파워에 입사파의 스펙트럼을 곱하면 시간평균 파워 스펙트럼을 구할 수 있고 이를 주파수에 따라 적분하면 불규칙파에서의 시간평균 파워의 대푯값 (PE)을 구할 수 있다. 시간평균 파워의 대푯값을 단위 폭당 입사파의 파워로 나누면 취득 폭의 대푯값을 구할 수 있다.
(22) |
Fig. 2는 원기둥의 중심에서 계산된 수평 원기둥 파력발전장치에 작용하는 수평과 수직 방향의 무차원화된 유체력(부가질량, 방사감쇠계수, 파기진력)을 보여주고 있다. 계산에 사용한 수평 원기둥의 폭(W)과 반경(R)은 각각 5 m, 2m이며, 잠긴 깊이(d)는 1.6 m이며 수심은 80 m이다. WAMIT 상용코드에서 사용한 총 격자수는 1,152개이다. 회전축 중심과 원기둥의 중심이 일치할 때 횡 방향의 유체력은 모두 사라진다. 그러나 회전축의 중심이 원기둥 중심으로부터 이격되면 수평, 수직 방향의 유체력 성분들이 횡 방향의 유체력 발생에 기여한다. 이때 이격 거리가 클수록 횡 방향 유체력은 증가한다.
회전축의 위치 변화에 따른 유체력을 WAMIT 상용코드를 이용하여 구할 수 있지만 매번 계산을 해야 하는 번거로움이 따른다. 계산 양을 줄이기 위하여 본 연구에서는 Cruz and Salter(2006)이 제안한 축 변환식을 이용하였다(Table 1 참조). Table 1의 축 변환식에서 윗첨자 '는 원기둥 중심에서 계산된 값들을 나타내며 윗첨자가 없는 왼쪽 항은 회전축의 중심에서의 값을 의미한다. Fig. 3은 회전축의 중심에서 WAMIT을 사용하여 직접 계산한 수치해와 원기둥 중심에서 계산된 해를 Cruz and Salter(2006)의 축 변환식을 통하여 동일한 회전축 중심에서의 값으로 바꾼 결과를 비교한 그림이다. 회전축의 위치는 원기둥 중심으로부터 l0= 0.75R, α = 300o 만큼 떨어져 있다. 두 결과는 서로 완벽히 일치하였다. 따라서 이후 모든 계산에서는 Table 1에 나타난 축 변환식을 사용하였다.
Fig. 4,5,6은 6개의 서로 다른 회전축의 중심(α = 60o, 90o, 120o, 240o, 270o, 300o, l0= 0.75R)에서 횡 방향 유체력(부가질량, 방사감쇠계수, 파기진력)과 횡 운동 RAO을 주파수에 따라 살펴보았다. 무차원화된 점성 감쇠계수는 κ = 0.01로 고정시켰다. Fig. 4는 잠긴 깊이가 가장 작은 d/R = 0.8의 결과로 이 모델에 대한 자세한 제원은 Table 2에 정리하였다. Fig. 5와 6은 각각 d/R = 1.2, 1.6에 대한 결과로 각 모델에 대한 자세한 제원은 Table 3과 4에 나타내었다. 수평 원기둥에 작용하는 부가질량과 방사감쇠계수는 대칭 구조를 갖는 원기둥 단면의 특성상, 회전축의 위치가 α = 60o, 120o, 240o, 300o일 때 모두 같은 결과를 준다. 같은 이유로 회전축의 중심이 원기둥 중심의 바로 위(α = 270o)와 아래(α = 90o)에 놓일 때도 같은 결과를 준다. 그러나 같은 파기진력을 주는 회전축의 위치는 α가 (60o, 240o), (90o, 270o), (120o, 300o)에 위치할 때이다. 1자유도 운동방정식을 풀어 회전축의 중심에서 구한 횡 운동 RAO는 Fig. 4(d), 5(d), 6(d)에 나타내었다. 회전축의 위치가 α = 60o에서 α = 300o로 이동함에 따라 횡 운동모드의 고유주파수가 ωN= 1.42, 1.17, 1.42, 1.11, 0.87, 1.11 rad/s로 변화한다. 이는 회전축의 위치가 바뀜에 따라 횡 방향 관성모멘트와 정유체력이 달라지기 때문이다. α = 60o, 120o에서 고유주파수는 1.42 rad/s이며, α = 240o, 300o에서는 고유주파수가 다소 줄어든 1.11 rad/s이다. z축에 대칭으로 위치할 때 관성모멘트, 정유체력, 부가 질량력가 서로 같기 때문이다 동일한 고유주파수를 준다. 회전축의 위치에 따른 고유주파수의 변동성을 살펴보면 회전축의 중심이 수면 아래에 또는 z축으로부터 멀리 벗어날수록 고유주파수는 커지며, 반대로 수면 위에 또는 z축 상에 가깝게 놓일수록 고유주파수는 줄어든다. 잠긴 깊이 변화에 따라 고유주파수의 정량적인 값은 달라지지만 앞서 기술한 회전축 위치에 따른 고유주파수의 특성은 잠긴 깊이에 관계없이 동일하게 적용된다.
Draft (d) | 1.6 m |
Radius of cylinder (R) | 2.0 m |
Mass of cylinder (mH) | 14405.3 kg |
Mass of ballast (mB) | 9613.3 kg |
Vertical distance from its geometric center to the center of ballast mass (ZB) | 1.84 m |
Moment of inertia (IH) at the center of mass (mH) | 56763.4 kg m2 |
Moment of inertia (IB) at the center of mass (mB) | 1551.1 kg m2 |
Heave Hydrostatic coefficient ( |
19.6 m2 |
Pitch Hydrostatic coefficient ( |
17.3 m3 |
Draft (d) | 2.4 m |
Radius of cylinder (R) | 2.0 m |
Mass of cylinder (mH) | 14405.3 kg |
Mass of ballast (mB) | 25898.7 kg |
Vertical distance from its geometric center to the center of ballast mass (ZB) | 1.73 m |
Moment of inertia (IH) at the center of mass (mH) | 56763.4 kg m2 |
Moment of inertia (IB) at the center of mass (mB) | 7998.5 kg m2 |
Heave Hydrostatic coefficient ( |
19.2 m2 |
Roll Hydrostatic coefficient ( |
42.7 m3 |
Draft (d) | 3.2 m |
Radius of cylinder (R) | 2.0 m |
Mass of cylinder (mH) | 14405.3 kg |
Mass of ballast (mB) | 40830.5 kg |
Vertical distance from its geometric center to the center of ballast mass (ZB) | 1.64 m |
Moment of inertia (IH) at the center of mass (mH) | 56763.4 kg m2 |
Moment of inertia (IB) at the center of mass (mB) | 16937.7 kg m2 |
Heave Hydrostatic coefficient ( |
20.0 m2 |
Roll Hydrostatic coefficient ( |
78.2 m3 |
Fig. 7는 시간평균 추출 파워의 최대값을 주는 식 (13)의 최적의 PTO감쇠계수(
Fig. 8은 수평 원기둥 파력발전장치로부터 추출되는 최적의 파워와 이 값을 입사파의 파워와 원기둥의 전체 폭(W)으로 나눈 무차원화된 취득 폭을 주파수에 따라 살펴보았다. 파워 곡선의 최대값을 주는 특정 주파수는 공진주파수와 일치함을 알 수 있다. 공진주파수가 저주파수 영역에 위치할 때는 추출 파워의 최대값은 다소 증가한 반면에 공진 폭은 좁아지는 특징을 가지고 있다. 그러나 공진주파수가 고주파수 영역에 위치할 때는 반대로 추출 파워의 최대값은 줄어든 반면에 공진 폭은 넓어진다. 다양한 주파수가 혼재된 실제 해역에서 많은 양의 파랑에너지를 추출하기 위해서는 공진주파수 주변에서 높은 파워를 가짐과 동시에 공진 폭도 넓어야 한다. 그러나 무엇보다도 중요한 것은 공진주파수가 파랑에너지가 밀집된 주파수 대역에 위치해야 한다. 예상대로 무차원화된 취득 폭을 나타내는 곡선은 최적의 추출 파워 곡선과 유사하며 최대값은 대략 1.4이다.
Fig. 9는 불규칙파로 확장하였을 때 최적의 회전축의 위치를 찾는 과정을 소개한 그림이다. 본 계산에서 사용한 불규칙파는 γ = 2.2인 JONSWAP 스펙트럼에 수심의 효과를 고려한 TMA스펙트럼을 사용하였다. 유의파고와 피크주기는 각각 H1/3= 2.0m, TP= 6.65s이다. 회전축의 중심과 원기둥 중심간의 거리를 l0= 0.75R로 고정시키고 상대 각도(α)를 0부터 360도까지 5도씩 증가시키면서 식 (22)에 주어진 시간평균 파워의 대푯값을 극 좌표계로 나타내었다. 단위는 kW이다. 3가지 잠긴 깊이(d/R = 0.8, 1.2, 1.6)에 대하여 살펴보았다. 먼저 잠긴 깊이에 관계없이 불규칙파중 시간평균 파워의 최대값은 α = 300o(회전축이 수면 밖에 놓일 때) 주변에서 일어나며 특히, 잠긴 깊이가 가장 낮은 d/R = 0.8에서 상대적으로 가장 큰 시간평균 파워(12.5 kW)을 주었다. d/R = 1.2, 1.6일 때도 회전축의 중심이 수면보다 위에 놓일 때가 수면 아래에 놓일 때 보다 큰 파워 값을 준다. 그러나 앞서 살펴본 d/R = 0.8보다는 전반적으로 추출 파워가 떨어진다.
Fig. 9에서 나타난 유의파고 H1/3= 2.0m와 피크주기 TP= 6.65s인 불규칙파 중에서 얻을 수 있는 시간평균 파워의 대푯값과 회전축의 중심간의 상관 관계를 명확히 설명하기 위하여 고정된 l0= 0.75R에서 상대 각도 α의 변화에 따른 횡 운동의 고유주파수의 변화를 Fig. 10a에 그렸다. 수평 실선은 피크주기 TP= 6.65s에 해당하는 피크주파수이다. 횡 운동의 고유주파수와 피크주파수가 일치할 때 공진으로 큰 운동 변위와 큰 파워 추출을 기대할 수 있다. 그림에서 무차원화된 잠긴 깊이가 d/R = 0.8에서 상대 각도 α = 250o, 290o에서 공진이 발생함을 알 수 있었다. 이는 Fig. 9에서 큰 추출 파워의 발생으로 잘 나타나 있다. 무차원화된 잠긴 깊이가 d/R = 1.2에서는 α = 270o에서 공진이 일어나지만 Fig. 9에서 α = 270o보다 다소 큰 각도에서 최대 흡수 파워를 보이고 있다. 잠긴 깊이가 가장 큰 d/R=1.6에서는 공진이 일어나지 않았다. 따라서 Fig. 9에서 보듯이 3개의 잠긴 깊이 중에서 가장 낮은 추출 파워를 보인다. Fig. 10b는 상대 각도 α = 300o에 대하여 회전축의 중심과 원기둥의 중심간의 거리를 증가시키면서 횡 운동 고유주파수를 그렸다. l0/R = 0은 회전축의 중심이 원기둥의 중심과 일치한 경우이다. 잠긴 깊이에 관계없이 거리가 증가할수록 횡 운동 고유주파수가 줄어드는 경향을 보여주고 있다. 즉, 회전축의 위치를 결정하는 2개의 중요한 설계변수(α, l0) 모두 1차 변환장치의 고유주파수에 영향을 미치며 불규칙파의 피크주파수와 횡운동의 고유주파수가 일치할 때 공진으로 인하여 추출 파워가 증가한다. 따라서 회전축의 위치를 설치해역의 파랑을 특성을 고려하여 공진이 발생하도록 적절히 잡으면 단순한 단면 형상을 갖는 1차 변환장치를 가지고도 비교적 큰 파랑에너지를 얻을 수 있다.
Table 5,6,7은 운동스펙트럼과 파워스펙트럼을 적분하여 구한 횡운동의 유의 진폭(ξ1/3)과 흡수 파워의 대푯값(PE) 그리고 취득 폭의 대푯값 (
d/R | α = 60o | α = 90o | α = 120o | α = 240o | α = 270o | α = 300o |
---|---|---|---|---|---|---|
0.8 | 0.75186 | 1.56093 | 0.91488 | 1.41569 | 1.36679 | 1.54427 |
1.2 | 0.38729 | 0.59965 | 0.47387 | 0.91088 | 1.24088 | 1.03010 |
1.6 | 0.20049 | 0.26035 | 0.23667 | 0.58471 | 0.88310 | 0.68443 |
d/R | α = 60o | α = 90o | α = 120o | α = 240o | α = 270o | α = 300o |
---|---|---|---|---|---|---|
0.8 | 6.59333 | 6.38133 | 9.33458 | 10.26533 | 3.47775 | 12.53613 |
1.2 | 4.77958 | 6.75611 | 6.94889 | 7.76469 | 7.60679 | 10.05305 |
1.6 | 3.43525 | 4.90952 | 4.75371 | 6.16085 | 8.31669 | 8.30858 |
d/R | α = 60o | α = 90o | α = 120o | α = 240o | α = 270o | α = 300o |
---|---|---|---|---|---|---|
0.8 | 0.53049 | 0.51343 | 0.75105 | 0.82593 | 0.27981 | 1.00864 |
1.2 | 0.38447 | 0.54347 | 0.55897 | 0.62460 | 0.61190 | 0.80868 |
1.6 | 0.27633 | 0.39493 | 0.38239 | 0.49558 | 0.66900 | 0.66835 |
형상이 단순한 수평 원기둥 1차 변환장치의 중요한 설계변수인 원기둥의 잠긴 깊이와 회전축의 위치 등을 바꿔가면서 이 값들이 추출 파워에 미치는 영향을 살펴보았다. 원기둥의 잠긴 깊이와 회전축의 위치 모두 관성모멘트와 정유체력을 변화시키므로 횡 운동모드의 고유주파수의 정량적인 값을 결정하는데 중요하다. 모든 파력발전장치의 1차 변환장치는 최대 추출 파워를 얻기 위해서 입사파의 주파수와 특정 운동모드의 고유주파수를 일치시켜 공진을 일으킨다. 따라서 설치해역의 파랑특성이 주어졌을 때 잠긴 깊이와 회전축의 위치를 적절히 설계하여 공진조건을 만족하도록 설계해야 한다.
주파수영역 해석법을 불규칙파(H1/3= 2.0m, TP= 6.65s)로 확장하여 1차 변환장치의 잠긴 깊이와 회전축의 위치 변화에 따른 횡 운동변위, 추출파워, 취득 폭의 대푯값을 구하였다. 3개의 잠긴 깊이(d = 1.6 m, 2.4 m, 3.2 m)에 대하여 l0= 0.75R을 고정시키고 상대 각도 α를 0부터 360o까지 변화시키면서 계산한 결과, 잠긴 깊이에 관계없이 추출 파워의 최대값은 α = 300o(회전축이 수면 밖에 놓일 때)에서 일어나며 특히, 잠긴 깊이가 가장 낮을 때 상대적으로 가장 큰 추출 파워를 주었다. 또한 상대 각도를 고정시키고 회전축의 중심과 원기둥의 중심간의 거리를 증가시키면 잠긴 깊이에 관계없이 횡 운동 고유주파수가 줄어드는 경향을 보여주었다. 즉, 회전축의 위치를 결정하는 2개의 중요한 설계변수(α, l0)모두 횡 운동의 고유주파수에 영향을 미치며, 횡 운동의 고유주파수가 불규칙파의 피크주파수와 일치할 때 공진으로 인하여 추출 파워는 크게 증가한다. 따라서 회전축의 위치를 설치해역의 파랑특성을 고려하여 공진이 발생하도록 적절히 잡으면 단순한 형상을 갖는 1차 변환장치를 통하여 비교적 큰 파랑에너지를 추출할 수 있다.
본 연구는 산업통상자원부(MOTIE)와 한국에너지기술평가원(KETEP)의 지원을 받아 수행한 연구 과제입니다(No. 20163010071690).
1. | Bouws, F., Gunther, H., Rosenthal, W., and Vincent, C.L., (1985), Similarity of the Wind Wave Spectrum in Finite Depth Water: 1. Spectral Form, J. Geophys. Res, 90(1), p975-986. |
2. | Count, B.M., (1978), On the Dynamics of Wave-power Devices, Proc. of the Royal Society of London. Series A, Mathematical and Physical Sciences, 363, p559-579. |
3. | Cruz, J., and Salter, S.H., (2006), Numerical and experimental modelling of a modified version of the Edinburgh duck wave energy device, Proc. IMechE PartM- Journal of Engineering for the Maritime Environment, 220, p129-147. |
4. | Evans, D.V., (1976), A Theory for Wave-power Absorption by Oscillating Bodies, J. Fluid Mech, 77, p1-25. |
5. | Evans, D.V., Jeffrey, D.C., Salter, S.H., and Taylor, J., (1979), Submerged Cylinder Wave Energy Device: Theory and Experiment, Appl. Ocean Res, 1(1), p3-12. |
6. | Goda, Y., (1988), Statistical Variability of Sea State Parameters as a Function of a Wave Spectrum, Coastal Eng. in Japan, 31, p39-52. |
7. | Greenhow, M., and Ahn, S.I., (1988), Added Mass and Damping of Horizontal Circular Cylinder Sections, Ocean Eng, 15(5), p495-504. |
8. | Jeffrey, D.C., Richmond, D.J.E., Salter, S.H., Taylor, J.R.M., and Young, I.A., 1976, Second year interim report on Edinburgh wave power project: study of mechanism to extract power from sea waves, Interim report, Wave-power project, University of Edinburgh, September. |
9. | Lucas, J., Salter, S.H., Cruz, J., Taylor, J., and Bryden, I., 2009, Performance optimization of a modified duck through optimal mass distribution, In Proc. of the 8th European Wave and Tidal Energy Conference, page 10, Uppsala, Sweden, September. |
10. | Mynett, A.E., Serman, D.D., and Mei, C.C., (1979), Characteristics of Salter’s Cam for Extracting Energy from Ocean Waves, Appl. Ocean Res, 1(1), p13-20. |
11. | Pizer, D.J., (1992), The numerical predictions of the performance of a solo duck, Technical report, University Of Edinburgh. |
12. | Pecher, A., and Kofoed, J.P., (2013), Experimental study on the updated PTO system of the WEPTOS wave energy converter, Aalborg University DCE Contract Report 138. |
13. | Pecher, A., and Kofoed, J.P., (2014), Experimental study on the structural and mooring loads of the WEPTOS wave energy converter, Aalborg University DCE Contract Report 142. |
14. | Salter, S.H., (1974), Wave Power, Nature, 249(249), p720-724. |
15. | Salter, S.H., Jeffrey, D.C., and Taylor, J., 1975, First year interim report on Edinburgh wave power project: study of mechanism to extract power from sea waves, Technical report, Wave-power project, University of Edinburgh, September. |
16. | Salter, S.H., Cruz, J., Lucas, J., and Pascal, R., (2007), Wave powered desalination, In International Conference on Integrated Sustainable Energy Resources in Arid Regions, Abu Dhabi. |
Copyright (c) 1997 by Korean Society for Marine Environment & Energy (KOSMEE).
Authorization to photocopy articles may be granted by the KOSMEE, provided the material is used on a personal basis only.
The Society does not consent copying for general distribution, promotion, for creating a new work or for resale. Permission
to photocopy articles must be requested to the KOSMEE main office.
(office contact : kosmee@kosmee.or.kr)