| Home | E-Submission | Sitemap | Contact Us |  
Cho: Wave-Energy Extraction by a Compact Circular Array of Buoys

Abstract

The wave power, extracted from a circular array of small power buoys, is investigated under the potential theory. It is assumed that the buoy's radius, the draft, and the separation distance are much smaller than the water depth, the wave length, and the radius of a circular deployment area. The boundary value problem involving the macro-scale boundary condition on the mean surface covered by buoys is solved using the eigenfunction expansion method. The capture width, which is defined as the ratio of the extracted power to the wave power per unit length of the incident wave crest, is assessed for various combinations of packing ratio, radius of a circular array, and PTO damping coefficient. It is found that the circular array deployment is more effective in the viewpoint of efficiency than the single large buoy of the same total displaced volume.

Abstract

크기가 작은 다수의 부이들이 원형으로 배열되어 있을 때 부이를 통한 추출 파워를 포텐셜이론 통하여 살펴보았다. 부이의 크기와 배열 간격이 수심과 파장 그리고 설치면적의 반경에 비하여 작다는 가정 아래에서 다수부이로 채워진 수면에서의 경계조건식을 유도하고 이를 포함한 경계치문제를 고유함수전개법을 이용하여 해석하였다. 추출 파워와 입사파 에너지의 비인 취득 폭을 packing ratio, 배열 반경, 그리고 PTO 감쇠계수를 바꿔가면서 살펴보았다. 또한 다수의 소형 부이를 원형 배열하였을 때 추출된 파력이 같은 잠긴 체적을 갖는 대형 단일 부이보다 효율적임을 계산을 통하여 밝혔다.

1. 서 론

해양 선진국인 노르웨이와 영국에서는 반잠수식 해양구조물(semi-submersible rig) 아래에 매달인 여러개의 소형 부이들로 구성된 FO3(http://www.abb.com)와 맨체스터 Bobber(http://www.manchesterbobber.com)라 불리는 파력발전장치를 개발하였다. 진행하는 파는 수면 위에 떠 있는 부이의 수직운동을 일으키고 이러한 운동에너지는 유압방식으로 작동하는 발전기를 통하여 전기에너지로 변환된다. 현재 축척비 1/3인 연구모델(12 m × 12 m ×8 m)을 제작하여 실해역에서 실험중에 있으며 실물로 제작되었을 때 유의파고 6 m, 피크주기(peak period) 9초인 해역에서 2.52 MW 전기를 생산할 것으로 예측하고 있다. 최근 국내에서도 한국해양과학기술원을 중심으로 반잠수식 해양구조물 각 모서리에 2 MW 발전량을 가진 풍력터빈 4기를 설치하여 총 8 MW의 전기를 풍력발전으로 생산하고 구조물 아래 수면에 여러개의 파력발전장치를 설치하여 총 10 MW급 부유식 파력-해상풍력 연계형 발전시스템 설계기술 개발을 진행하고 있다. 파랑에너지로부터 1차 변환과정을 통하여 얻은 부이의 운동에너지에서 전기에너지를 추출하기 위한 2차 변환방식으로 유압방식이 아닌 영구자석과 코일로 구성된 선형발전기 방식을 채택하였다. 노르웨이와 영국에서 연구중인 파력발전장치내의 1차 변환장치는 다수의 소형 부이를 사용한다. 따라서 부이는 입사파의 피크 주파수와 공진을 일으키지 않고 파를 따라 수직운동을 하는 반면에 국내에서 채택한 1차 변환장치는 입사파 스펙트럼의 피크 주파수와 공진을 일으키도록 부이를 설계하므로 부이의 크기는 대형이며 개수가 상대적으로 줄어든다.
다양한 형태의 단면을 갖는 부이에 대한 회절과 방사 문제는 선형포텐셜이론을 이용하여 Black et al.(1971), Garrett(1971), Yeung(1982) 등에 의해 다루어졌다. 원형 부이에 의한 파랑 에너지 흡수는 Budal and Falnes(1975)에 의해 연구되었다. 다수의 부이들이 배열된 경우에 부이 간의 상호간섭이 약하다는 가정 아래에서 다수 부이들에 의해 흡수된 파력을 구하는 근사 이론들이 개발되었다(Falnes, 1980, 1984; Falnes and Budal, 1982; Falcäo, 2002). 여러개의 부유체간의 상호작용을 정확히 고려한 회절과 방사 문제는 Linton and Evans(1990), Linton and McIver(1996), Chamberlain(2007)등에 의해 연구되었는데 연구결과를 파랑에너지 추출에 적용하지는 않았다. Kweon and Lee(2012)는 파랑에너지 추출을 위하여 배열된 부이들에 의한 산란문제를 천수역 방정식인 완경사 방정식(mild-slope equation)을 수치적으로 풀어 해석하였다. Garnaud and Mei(2009a,b)는 다수의 소형 부이들이 원형과 사각형으로 배열된 경우에 다척도(multi-scale) 근사기법을 사용하여 입사파와 배열된 부이간의 상호작용 문제를 풀었다.
본 연구에서는 다수의 소형 부이들이 일정한 간격을 갖고 원형으로 배열되어 있을 때 Garnaud and Mei(2009a,b)가 제안한 다수 부이들로 채워진 수면에서 만족하는 경계조건식을 유도하고 이를 포함한 경계치문제를 고유함수전개법(eigenfunction expansion method)으로 풀어 추출 파워를 구하였다. 다수의 부이가 채워진 수면에서의 경계조건식은 부이의 바닥면에서의 운동학적 경계조건과 동역학적 경계조건을 결합한 근사식으로 수심, 파장, 배열 폭이 부이의 반경, 흘수, 간격에 비하여 아주 길다는 가정에 기반을 두고 있다. 각 부이에는 PTO(Power take-off) 감쇠력으로 모델링된 선형발전기가 연결되어 파랑에너지를 흡수한다. 따라서 다수 부이가 배열된 수면에서 만족하는 경계조건식내에 PTO 감쇠력이 포함되어 있어 고유치 문제를 풀어 구하는 고유값은 복소수 값을 갖는다. 고유값을 구하는 수치해석 방법으로 Steffensen의 축차법(iteration method)을 사용하였다(Cho, 2013). 부이가 채워진 수면 전체 면적과 부이의 면적의 비인 packing ratio과 PTO 감쇠력, 그리고 배열 반경을 바꿔가면서 원형으로 배열된 다수 부이들에 의한 에너지 플럭스(energy flux)를 구하고, 이 값을 단위 폭당 입사파의 에너지로 나눈 취득 폭(capture width)을 여러 가지 설계변수들을 바꿔가면서 살펴보았다. 다수의 소형 부이를 원형 배열하였을 때 얻은 파력 추출은 같은 잠긴 체적을 갖는 단일 대형 부이보다 효율적임을 밝혔다.
Fig. 1.

Definition sketch of a circular array of small buoys.

JKSCOE_2014_v26n2_103_f001.jpg

2. 문제의 정식화

2.1 단일 부이

반경(as)과 흘수(ds)가 같은 원형 단면을 갖는 부이가 수심(h)인 수면 위에 떠 있다고 가정하자. 부이의 운동에너지로부터 전기에너지를 추출하기 위해서는 선형발전기를 연결하는데 이때 부이와 함께 움직이는 영구자석의 운동방향과 반대방향으로 작용하는 PTO 감쇠력이 존재한다. 본 연구에서는 편의상 PTO 감쇠력이 영구자석의 운동속도에 선형적으로 비례 한다고 가정한다. 선형발전기와 연성된 부이의 수직운동 변위(ξ)를 입사파의 진폭(A)으로 나눈 RAO(response amplitude operator)는 아래식과 같다.
JKSCOE_2014_v26n2_103_s001.jpg
여기서 m은 부이의 질량이며, a33, b33, X3 는 각각 부가질량, 방사감쇠계수, 파기진력으로 부이의 형상과 주파수(ω)의 함수이다. PTO 감쇠계수인 cPTO는 선형발전기의 고정자(코일)와 가동자(영구자석)의 형상과 코일의 감긴 모양, 구동 주파수와 밀접한 관련이 있다.
시간평균 파력은 아래와 같다.
JKSCOE_2014_v26n2_103_s002.jpg
최대 파워는 ① 공진조건 pgπas2 - ω2(m+a33) = 0 에서와 ② PTO 감쇠계수가 방사감쇠계수와 같을 때(cPTO = b33 )발생한다. 식 (1)과 (2)로부터 시간평균 최대 파워는 다음과 같다.
JKSCOE_2014_v26n2_103_s003.jpg
한편, 시간평균 파워 JKSCOE_2014_v26n2_103_ms001.jpg 와 단위 폭을 갖는 입사파의 평균파워의 비를 취득 폭이라 하며 기호로 Cw로 쓴다.
JKSCOE_2014_v26n2_103_s004.jpg
여기서 JKSCOE_2014_v26n2_103_ms002.jpg 는 군속도이다. ηmech(=JKSCOE_2014_v26n2_103_ms001.jpg/JKSCOE_2014_v26n2_103_ms001.jpgmax)는 추출된 파워와 최대 파워의 비로 기계적인 효율을 나타내며 최적 조건인 ① 과 ②를 만족하면 ηmech = 1 이 된다. 길이 차원을 갖는 취득 폭에 입사파의 파수를 곱하여 무차원화시키면 kCw이다.

2.2 원형배열

부이(반지름 a, 흘수 d )가 일정한 간격 s을 가지고 반경 R인 설치 영역내에 배열되어 있다고 생각하자(Fig. 1 참조). 배열된 다수 부이들에 의해 추출되는 파랑에너지를 구하기 위하여 선형포텐셜이론을 가정하여 속도포텐셜을 도입한다. 파동장이 주파수 ω을 갖고 조화운동을 한다고 가정하면 속도포텐셜을 시간과 공간의 함수로 분리하여 Φ(r,θ, z, t) = Re {φ(r,θ, z)eiωt} 로 쓸 수 있다. 원통좌표계에서 살펴본 속도포텐셜 φ(r, θ, z) 로부터 변수분리법을 이용하여 θ의 고유함수를 분리하면 아래와 같이 표현된다.
JKSCOE_2014_v26n2_103_s005.jpg
고유함수전개법을 사용하기 위하여 유체영역을 Fig. 1과 같이 영역 (I)과 영역 (II)로 나눈다. 영역 (I)은 r ≥ R, 영역(II)는 r ≤ R으로 정의된다. 영역 (I)의 속도포텐셜 φl(1) (r, z)은 아래와 같이 입사포텐셜과 산란포텐셜의 합으로 표현된다.
JKSCOE_2014_v26n2_103_s006.jpg
여기서 I0 = N10 ⁄ cosh kh 이며, l = 0 일 때 βl = 1 이며, l ≥ 1 일 때 βl 2(i)l 이다. n = 0 는 진행파(propagating waves)성분을 나타내며, n ≥ 1 은 배열된 부이 주위에만 존재하며 r이 커짐에 따라 급격히 사라지는 비진행파(evanescent waves) 성분을 나타낸다. JlKl는 각각 제 1종 Bessel함수, Kelvin함수이다. 영역 (I)에서의 고유값(k0 = –ikkn,n = 1,2,…)은 잘 알려진 선형분산식(kntanknh = –ω2g)을 만족하며, 고유함수 fn(z)은 다음과 같이 정의된다.
JKSCOE_2014_v26n2_103_s007.jpg
영역 (I)에서의 고유함수 fn(z)는 아래와 같은 직교성(orthogonality)을 만족한다.
JKSCOE_2014_v26n2_103_s008.jpg
여기서 δmnm = n 일 때 δmn = 1, mn일 때 δmn = 0 으로 정의된 Kronecker delta이다.
다수의 부이들이 원형으로 배열되어 있을 때 부이간의 상호작용을 정확히 고려하기 위해서는 복잡한 수식전개과정과 많은 계산시간이 필요하나 본 연구에서는 근사해법을 이용하여 다수 부이들에 의하여 추출되는 파워를 간략하게 구하였다. 이를 위하여 거시적 규모(macro-scale)의 길이 척도(length-scale)를 갖는 입사파의 파장(l), 수심(h), 배열 반경(R)이 미시적 규모(micro-scale)의 길이 척도를 갖는 반지름(a), 흘수(d), 배열 간격(s)에 비하여 훨씬 크다고 가정한다. 이를 수식으로 정리하면 아래 식과 같다.
JKSCOE_2014_v26n2_103_s009.jpg
이와 같이 해석모델 내에 크기가 전혀 다른 2개의 길이가 공존할 때 다척도 기법을 적용하여 근사해를 구한다. 거시적 규모의 길이 척도를 가지고 문제를 재구성하였을 때 부이가 채워진 수면에서의 경계조건식을 유도할 수 있고, 이를 포함한 경계치문제를 풀어 원형 영역내에 배열된 부이들을 통하여 흡수되는 파워를 구하였다.
다수의 부이들이 일정한 간격으로 배열된 수면에서 만족하는 경계조건식은 다음과 같다(Garnaud and Mei, 2009a). 부록 A에 자세한 유도과정을 소개하였다.
JKSCOE_2014_v26n2_103_s010.jpg
여기서 cPTO = cPTO/ ρgπa2 는 부이의 복원력 계수로 나눈 PTO 감쇠계수이며, p(= πa2/s2) 는 다수의 부이가 놓인 수면의 면적 중 부이가 차지하는 면적의 비로 packing ratio라 부른다. 반경의 범위는 0 < a < s ⁄ 2 이므로, p의 범위는 JKSCOE_2014_v26n2_103_ms003.jpg 이다. 또한 p→0는 수면 위에 부이가 없는 경우를 의미하며, 이 경우 식 (10)은 잘 알려진 자유표면 경계조건식 JKSCOE_2014_v26n2_103_ms004.jpg으로 바뀐다.
Laplace 방정식과 수면에서의 경계조건식(10) 그리고 해저면 경계조건식(∂φ∂z = 0, at z = –h)을 만족하는 영역 (II)에서의 속도포텐셜은 다음과 같다.
JKSCOE_2014_v26n2_103_s011.jpg
영역 (II)의 고유함수는 아래 주어진 Stüm-Liouville형 고유치 문제를 만족한다.
JKSCOE_2014_v26n2_103_s012.jpg
여기서 JKSCOE_2014_v26n2_103_ms005.jpg 이다. 식 (12a)와 (12c)를 만족하는 Fn(z) = Cncoshλn(z+h) 이며, 이 함수를 식 (12b)에 대입하면 고유값 λn, (n=1,2,...)은 아래 주어진 비선형 분산식을 만족한다.
JKSCOE_2014_v26n2_103_s013.jpg
식 (13)에 주어진 비선형 방정식의 복소해를 얻기 위하여 Steffensen법을 사용하였다(Mathews and Fink, 2004). Steffensen법은 3개의 초기값 (λ0 W (λ0)), (λ1 W (λ1)), (λ2 W (λ2)) 을 가지고 출발하는 축차법의 일종인 Muller법으로 비선형방정식의 복소수 해를 구하는데 많이 이용된다. 첫 번째 초기값을 가정하고 나머지 초기값을 Newton-Rapson 방법을 이용한다(Cho, 2013).
고유함수의 직교성JKSCOE_2014_v26n2_103_ms006.jpg 을 만족하도록 정규화된 고유함수를 구하면 다음과 같다.
JKSCOE_2014_v26n2_103_s014.jpg
여기서 JKSCOE_2014_v26n2_103_ms007.jpg 이다.
영역 (I),(II)의 미지수 Aln, Bln, (l,n=0,1,2,…)는 두 영역이 만나는 r = R 에서 속도포텐셜 (φl) 과 법선방향 속도(∂φl/∂r)가 서로 같다는 정합조건식으로부터 구해진다.
r = R 에서 속도포텐셜이 서로 같다는 정합조건식과 함께 고유함수 fn(z)의 직교성을 이용하면 아래와 같은 관계식을 얻을 수 있다.
JKSCOE_2014_v26n2_103_s015.jpg
여기서 JKSCOE_2014_v26n2_103_ms008.jpg 이다.
r = R에서 영역 (I)과 영역(II)의 법선방향 속도가 서로 같다는 정합조건식을 적용한 뒤 양변에 Fm(z)를 곱하고 −h에서 0까지 z에 대해 적분하면 다음식을 얻는다.
JKSCOE_2014_v26n2_103_s016.jpg
식 (16)을 식 (15)에 대입하여 정리하면 아래와 같은 대수방정식을 유도할 수 있다.
JKSCOE_2014_v26n2_103_s017.jpg
식 (17)에 주어진 대수방정식을 수치적으로 풀기 위하여 고유함수의 개수(N)와 θ방향의 모드 개수를 유한개(L)로 잡아 주어진 l = 0, 1, 2,…, L 에 대하여 (N + 1) × (N + 1) 대수방정식을 풀어 미지수 Alm, (l = 0, 1, 2,…, L, m = 0, 1, 2,…,N)을 수치적으로 구한다. Alm이 구해지면 영역 (II)에서의 미지수 Blm는 식 (16)으로부터 구할 수 있다.
배열된 부이 주변에서의 파형식은 다음과 같다.
JKSCOE_2014_v26n2_103_s018.jpg
배열된 부이들로부터 추출되는 파워는 배열된 부이들로부터 멀리 떨어진 가상의 원기둥 표면(S)을 통한 시간평균 에너지 플럭스와 같다. S 를 통과하는 시간평균 에너지 플럭스는 다음과 같이 정의된다.
JKSCOE_2014_v26n2_103_s019.jpg
여기서 윗첨자 (*)는 공액 복소수를 나타낸다.
식 (19)내의 속도포텐셜은 입사포텐셜과 산란포텐셜의 합이다. 실제로 적분값에 기여하는 산란포텐셜의 r→∞ 한 근사식을 구하면 다음과 같다.
JKSCOE_2014_v26n2_103_s020.jpg
여기서 JKSCOE_2014_v26n2_103_ms009.jpg(θ)는 가상의 원기둥 (S)을 통하여 빠져나가는 산란파의 진폭으로 JKSCOE_2014_v26n2_103_ms010.jpg이다.
식 (20)를 식 (19)에 대입하고 stationary phase 법을 이용하면 시간평균 에너지 플럭스는 아래와 같다.
JKSCOE_2014_v26n2_103_s021.jpg
시간평균 파워 JKSCOE_2014_v26n2_103_ms011.jpg와 단위 폭을 갖는 입사파의 평균 파워의 비를 취득 폭이라 부르며 기호로 Cw로 쓴다.
JKSCOE_2014_v26n2_103_s022.jpg
길이 차원을 갖는 취득 폭을 무차원화된 값으로 만들기 위한 방법으로 입사파의 파수를 곱한 kCw를 사용하거나 배열된 면적의 직경 2R로 나눈 Cw/2R 을 사용한다. 두 값 모두 배열된 부이에 의한 에너지 흡수 효율을 나타내는 척도이다.

3. 계산결과 및 고찰

Fig. 2는 원기둥의 반경(as)과 흘수(ds)는 같고 선형발전기의 PTO 감쇠계수(cPTO)는 최대 파워가 추출되는 최적조건(cPTO= b33)일 때의 계산결과이다. 시간평균 추출 파워(JKSCOE_2014_v26n2_103_ms001.jpg)와 최대 파워(JKSCOE_2014_v26n2_103_ms001.jpgmax)의 비인 기계적 효율(ηmech)과 추출 파워 JKSCOE_2014_v26n2_103_ms001.jpg와 단위 폭을 갖는 입사파의 평균 파워의 비인 취득 폭(Cw)을 살펴보았다. 질량에 비하여 부가질량이 상대적으로 작은 경우, 공진주파수는 JKSCOE_2014_v26n2_103_ms012.jpg로 근사된다. 따라서 공진주파수는 원기둥의 흘수에 의해 결정된다. 3개의 무차원화된 흘수(ds ⁄ h = 0.2, 0.3, 0.4)에 대한 공진을 일으키는 kh는 각각 3.2, 2.2, 1.7이다. 예상대로 Fig. 2a에서 기계적 효율이 1인 주파수는 kh = 3.2, 2.2, 1.7 와 일치한다. 또한 흘수가 증가함에 따라 기계적 효율의 피크값이 저주파수 영역으로 이동하며 피크값을 나타내는 주파수의 범위가 점차 줄어드는 경향을 보여주고 있다. Fig. 2b에서 흘수가 증가함에 따라 취득폭은 커지나 기계적 효율과 마찬가지로 최대값을 주는 주파수의 범위는 줄어드는 경향을 보여주고 있다.
Fig. 2.

Mechanical efficiency (a) capture width (b) for a single buoy of different draft as a function of non-dimensional wavelength for as = ds.

JKSCOE_2014_v26n2_103_f002.jpg
Muller법과 Newton-Raphson법을 결합한 Steffensen법을 사용하여 식 (13)에 주어진 비선형 방정식을 풀어 구한 고유값 λn, (n = 1, 2,…, 10)을 cPTO = 1.0, p = 0.2 일 때 3개의 무차원화된 주파수JKSCOE_2014_v26n2_103_ms013.jpg에 대하여 Table 1에 정리하였다. 이 계산결과는 Garnaud and Mei (2009a)의 수치계산결과와 일치한다. 그들은 축차법 알고리즘을 사용하여 해를 구하였는데 축차법의 초기값을 얻기 위하여 식 (12)에 주어진 고유치 문제를 유한요소법을 이용하여 수치 계산하였다.
Table 1.

First ten eigenvalues for cPTO* = 1.0 and p = 0.2

Kn JKSCOE_2014_v26n2_103_t1-s001.jpg JKSCOE_2014_v26n2_103_t1-s002.jpg JKSCOE_2014_v26n2_103_t1-s003.jpg
1 0.5107+0.0230i 1.1165+0.0835i 3.3669+0.3159i
2 6.6600e-3+3.0634i 0.0357+2.8342i 0.0545+2.1332i
3 3.2200e-3+6.2448i 0.0163+6.1376i 0.0449+5.7538i
4 2.1300e-3+9.3993i 0.0107+9.3286i 0.0322+9.0697i
5 1.6000e-3+12.5473i 8.0100e-3+12.4945i 0.0247+12.2996i
6 1.2800e-3+15.6927i 6.3900e-3+15.6505i 0.0200+15.4944i
7 1.0600e-3+18.8368i 5.3200e-3+18.8017i 0.0168+18.6715i
8 9.1000e-4+21.9802i 4.5600e-3+21.9502i 0.0144+21.8385i
9 8.0000e-4+25.1232i 3.9900e-3+25.0969i 0.0126+24.9991i
10 7.1000e-4+28.2658i 3.5400e-3+28.2425i 0.0113+28.1556i
Fig. 3R ⁄ h = 1.0, p = 0.2일 때 균일하게 배열된 부이들에 의한 무차원화된 취득 폭 kCwCw ⁄ 2R 를 입사파의 파장(kh)과 PTO 감쇠계수(cPTO*)를 바꿔가면서 살펴보았다. 취득폭 kCwCw ⁄ 2RcPTO* = 0.5 일 때 파장에 관계없이 가장 큰 값을 보이고 있다. 단일 부이인 경우 최대 파워는 최적조건(cPTO = b33)일 때 일어나지만 다수 부이를 배열한 경우 부이들에 의해 흡수되는 최대 파워는 cPTO* = 0.5 일 때 발생한다.
Fig. 3.

Non-dimensional capture widths (a) kCw (b) Cw/2R for a circular array with various PTO damping coefficients cPTO* as a function of kh for R/h =1.0, p = 0.2.

JKSCOE_2014_v26n2_103_f003.jpg
다음 그림은 p = 0.2로 고정시키고 배열 반경(R)을 바꿔가면서 무차원화된 취득 폭을 그렸다. PTO 감쇠계수는 에너지를 최대로 흡수할 수 있는 cPTO* = 0.5 이다. 취득 폭은 예상대로 배열 반경이 커짐에 따라 일정하게 늘어나고 있음을 보여주고 있다. 특히 배열 직경과 같은 폭을 갖는 입사파의 파워에 대한 흡수된 파워의 비인 Cw/2R 는 배열 반경이 수심의 3배이고 kh가 2.5이상이면 입사파 에너지의 80%이상이 배열된 부이를 통하여 흡수된다. p = 0.2인 상태에서 배열 폭이 증가한다는 것은 많은 부이가 제작되어 투입됨을 의미하므로 경제성을 검토하여 적절한 개수의 부이를 배열하여야 한다.
Fig. 4.

Non-dimensional capture widths (a) kCw (b) Cw/2R for a circular array of buoys with various radius R as a function of for cPTO* = 0.5, p = 0.2.

JKSCOE_2014_v26n2_103_f004.jpg
Fig. 5는 수심으로 무차원화된 배열 반경(R/h)을 1.0으로 고정시키고 p를 0.1, 0.2, 0.4, 0.7로 바꿔가면서 취득 폭을 그렸다. p가 커짐은 주어진 설치면적내에 부이의 개수를 증가시켜 촘촘히 배열함을 의미한다. p가 커질수록 취득 폭은 증가하는 경향을 보여주고 있다. 그러나 kh = 4보다 큰 단파에서는 오히려 p = 0.4일 때의 취득 폭이 오히려 p = 0.7일 때보다 커짐을 볼 수 있다.
Fig. 5.

Non-dimensional capture widths kCw (b) Cw/2R for a circular array of buoys with various packing ratio p as a function of kh for cPTO* = 0.5, R/h = 1.0.

JKSCOE_2014_v26n2_103_f005.jpg
Fig. 6은 원형으로 배열된 부이 주변에서의 파형의 증폭비(|η(r, θ)/A|)를 등고선으로 나타내었다. 이때 p = 0.2, cPTO* = 0.5 이며, 무차원화된 배열 폭(R/h)은 1.0이다. 2개의 무차원화된 파장에 대하여 살펴보았다. kh가 작다는 것은 입사파의 파장에 비하여 배열 반경에 작음을 의미한다. 그림에서 상대적으로 큰 배열 반경에 대하여 배열된 부이 후면의 증폭비가 크게 줄어들었음을 알 수 있다.
Fig. 6.

Free surface elevation in the neighbourhood of a circular array of buoys for cPTO* = 0.5, R/h = 1.0, p = 0.2.

JKSCOE_2014_v26n2_103_f006.jpg
다음은 반경(as)과 흘수(ds)가 같은 대형 단일 부이의 취득폭과 소형 부이들을 원형 배열하였을 때의 취득 폭을 비교하였다. 이때 단일 부이의 잠긴 부분의 체적(πab3)과 다수 부이들의 체적(πpR2d) 은 서로 같다고 가정하였다. 여기서 d는 소형 부이의 흘수이다. 이러한 조건 아래에서 단일 부이의 반경과 흘수는 ab = (pR2d)1/3이다. 4개의 배열 반경 R/h = (0.5, 1, 2, 3)에 대하여 계산을 수행하였으며, 이때 p = 0.2이다. 소형 부이의 흘수를 d/h = 0.1로 하였을 때 대형 부이의 반경은 ab/h = (0.17, 0.27, 0.43, 0.56)이다. 계산에 사용한 PTO 감쇠계수는 단일 부이인 경우 최적조건(cPTO= b33)에서 계산하였고, 배열된 다수 부이인 경우는 최대 파워를 주는 cPTO* = 0.5로 놓고 계산하였다. 다수 부이를 원형으로 배열하였을 때의 무차원화된 취득 폭이 입사파의 파장에 관계없이 단일 부이를 통해 얻는 취득 폭에 비하여 훨씬 커짐을 알 수 있다. 그러나 다수 부이를 배열하는 경우 설치면적이 크게 증가하며, 부이 각각에 선형발전기를 한 개씩 설치해야 하므로 설치비용이 크게 증가한다.

4. 결 론

수면 위에 다수의 부이를 일정한 간격으로 원형 배열하였을 때 다수 부이를 통하여 추출되는 파워를 입사파의 파장, PTO 감쇠계수, packing ratio, 배열 반경을 변화시키면서 살펴보았다. 단일 부이를 통하여 추출 할 수 있는 최대 파워는 최적조건(cPTO= b33)일 때 일어났지만 다수 부이를 원형 배열한 경우 부이들에 의해 흡수되는 최대 파력은 파장에 관계없이 cPTO* ≅ 0.5 일 때 발생하였다.
고정된 packing ratio에 대하여 배열 반경을 증가시키면 무차원화된 취득 폭은 증가한다. 또한 배열 반경을 고정시키고 packing ratio을 증가시켜도 같은 결과를 준다. 그러나 두 경우 모두 에너지 추출은 증가하나 많은 양의 부이가 제작되어 설치되어야 하므로 경제성을 검토하여 적절한 개수의 부이를 선택하여야 한다. 한편 입사파의 파장에 비하여 배열 반경이 큰 경우 배열된 부이 후면에서의 파 증폭비가 크게 줄어들었다.
공진을 이용하는 대형 단일 부이와 공진이 발생하지 않는 소형 부이들을 원형 배열하였을 때의 취득 폭을 서로 비교하였다. 이때 잠긴 부분의 체적이 서로 같게 하였다. 다수 부이를 원형으로 배열하였을 때의 무차원화된 취득 폭이 입사파의 파장에 관계없이 단일 부이를 통한 취득 폭에 비하여 훨씬 컸다.
본 해석방법은 입사파의 파장, 수심, 배열 반경이 부이의 크기와 배열 간격에 비하여 아주 크다는 가정 아래에서 구한 근사식이다. 따라서 배열된 각 부이간의 상호간섭 효과를 고려하여 추출 파워를 정확히 산정하기에는 어려움이 따른다. 따라서 실해역에 적용하기에 앞서 배열된 부이들을 통하여 추출할 있는 전체 파워의 근사값을 추정하거나 부이 배치안의 기초자료로 활용하는데 도움을 줄 것이다. 앞으로 부이간의 상호작용을 정확히 고려할 수 있는 해석방법을 개발하여 정확한 추출 파워를 구하거나 다수의 부이를 효과적으로 배치하는데 활용할 것이다.

ACKNOWLEDGEMENTS

본 연구는 해양수산부의 해양청정에너지개발사업 “10MW급 부유식 파력-해상풍력 연계형 발전시스템 개발 및 인프라구축”의 지원으로 수행되었으며, 연구비 지원에 감사드립니다.

References

1. Black, J.L, Mei, C.C, Bray, M.C.G. (Radiation and scattering of water waves by rigid bodies, Journal of Fluid Mechanics, 1971). 46(1):151-164.
crossref
2. Budal, K, Falnes, J. (A resonant point absorber of ocean-wave power, Nature, 1975). 256, 478-479.
crossref
3. Chamberlain, P.G. (Water wave scattering by finite arrays of circular structures, IMA Journal of Applied Mathematics, 2007). 72(1):52-56.
crossref
4. Cho, I.H. (Reflection and transmission coefficients by a surface- mounted horizontal porous plate, Journal of Korean Society of Coastal and Ocean Engineers, 2013). 25(5):327-334. 10.9765/KSCOE.2013.25.5.327.
crossref
5. Falco, A.F de O. (Wave-power absorption by a periodic linear array of oscillating water columns, Ocean Engineering, 2002). 29(4):1164-186.

6. Falnes, J. (Radiation impedance matrix and optimum power absorption for interacting oscillators in surface waves, Applied Ocean Research, 1980). 2, 75-80.
crossref
7. Falnes, J, Budal, K. (Wave-power absorption by parallel rows of interacting oscillating bodies, Applied Ocean Research, 1982). 4(4):194-207.
crossref
8. Falnes, J. (Wave-power absorption by an array of attenuators oscillating with unconstrained amplitudes, Applied Ocean Research, 1984). 6(1):16-22.
crossref
9. Garnaud, X, Mei, C.C. (Wave power extraction by a compact array of buoys, Journal of Fluid Mechanics, 2009a). 635, 389-413.
crossref
10. Garnaud, X, Mei, C.C. (Comparison of wave power extraction by a compact array of small buoys and by a large buoy, Proceedings of the 8th European Wave and Tidal Energy Conference, Uppsala; Sweden: 2009b.
crossref
11. Garrett, C.J.R. (Wave forces on a circular dock, Journal of Fluid Mechanics, 1971). 46(1):129-139.
crossref
12. Kweon, H.-M, Lee, J.-L. (An array effect of wave energy farm buoys. Inter, Journal of Naval Architecture and Ocean Engineering, 2012). 4, 437-446.
crossref
13. Linton, C. M, Evans, D. V. (The interaction of waves with arrays of vertical cylinders, Journal of Fluid Mechanics, 1990). 215, 549-67.
crossref
14. Linton, C.M, Mclver, R. (The scattering of water waves by an array of circular cylinders in a channel, Journal of Engineering Mathematics, 1996). 30, 661-82.
crossref
15. Mathews, J.H, Fink, K.K. (Numerical Methods Using Matlab, 4th Edition, Prentice-Hall Inc: 2004.

16. Yeung, R.W. (Added mass and damping of a vertical cylinder in finite-depth waters, Applied Ocean Research, 1981). 3, 119-113.
crossref

APPENDICES

부록 A

원형 단면을 갖는 부이의 바닥면(z = −d)에서의 수직운동 속도는 물입자의 속도와 같다는 운동학적 경계조건식과 부이 바닥면에서의 전체 힘은 동유체력, 복원력, 그리고 PTO감쇠력의 합으로 표현되며 Newton의 제2법칙으로부터 유도한 동역학적 경계조건식은 다음과 같다.
JKSCOE_2014_v26n2_103_s-a001.jpg
부이 바닥면에서 위의 경계조건식들을 결합한 후 부이 바닥면에서의 속도포텐셜은 일정하다고 가정하면 아래와 같은 식을 얻을 수 있다.
JKSCOE_2014_v26n2_103_s-a002.jpg
부이가 주파수 ω을 갖고 조화운동을 한다고 가정하면 식(A2)는 아래와 같이 바뀐다.
JKSCOE_2014_v26n2_103_s-a003.jpg
거시적인 규모의 관점에서 문제를 바라보았을 때, 위의 경계조건식은 z = −d가 아닌 z = 0 에 적용된다. 부이의 중량은 부력과 같으므로 m = ρπa2d 이다. 따라서 2/ρgπa2 = ω2d/gkd = O(μ) 이므로 식 (A3)는 아래식으로 근사시킬 수 있다.
JKSCOE_2014_v26n2_103_s-a004.jpg
여기서 cPTO* = cPTO/ρgπa2이다. 다수 부이들로 채워지지 않은 수면에서의 경계조건식은 잘 알려진 자유표면경계조건식 JKSCOE_2014_v26n2_103_ms014.jpg이다. 전체 면적에서 부이가 차지하는 면적의 비를 p라 할 때 수면에서의 평균적인 경계조건식은 아래와 같다.
JKSCOE_2014_v26n2_103_s-a005.jpg
Editorial Office
Korean Society of Coastal and Ocean Engineers,
#1132, LG EClat, 71 Banpo-daero 14-gil, Seocho, Seoul, Korea
Tel: +82-2-3474-1934,   Fax: +82-2-3473-1934   E-mail : cocean@kscoe.or.kr
Copyright© Korean Society of Coastal and Ocean Engineers. All rights reserved.       powerd by M2community
About |  Browse Articles |  Current Issue |  For Authors and Reviewers