In the present investigation, the hydrodynamic characteristics of a vertically floating hollow cylinder in regular waves have been studied. The potential theory for solving the diffraction and radiation problem was employed by assuming that the heave response motion was linear. By using the matched eigenfunction expansion method, the characteristics of the exciting forces, hydrodynamic coefficients, and heave motion responses were investigated with various system parameters such as the radius and draft of a hollow cylinder. In the present analytical model, two resonances are identified: the system resonance of a hollow cylinder and the piston-mode resonance in the confined inner fluid region. The piston resonance mode is especially important in the motion response of a hollow circular cylinder. In many cases, the heave response at the piston resonance mode is large, and its resonant frequency can be predicted using the empirical formula of Fukuda (1977). The present design tool can be applied to analyze the motion response of a spar offshore structure with a moon pool.
내부가 빈 부유식 원기둥 구조물은 해양공학의 다양한 분야에 활용이 가능한 구조물로 최근에 많은 관심을 받고 있다. 대표적인 활용 분야로 석유와 가스를 채굴하는 문풀(Moon pool)을 갖는 스파형(Spar) 구조물, 파도로부터 계류시설과 연안 시설물을 보호하는 부유식 방파제, 내부 유체의 공진을 이용하여 파랑 에너지를 증폭시키는 파랑 에너지 변환장치, 태양 에너지의 흡수와 축적을 동시에 실행하는 태양 연못(Solar pond) 등이 있다.
특히 탐사나 시추 작업을 위하여 스파형 구조물 내부에 조성된 문풀은 갑판부터 바닥까지 관통된 개구부로 주로 원형이나 사각형 형상을 갖는다. 문풀은 유체와 공기가 만나는 자유 수면이 존재하므로 입사파의 주파수와 문풀 내부의 고유주파수가 일치하면 공진 현상이 발생하여 문풀 내의 수면은 크게 증폭한다. 문풀 내부 유체 영역에서 크게 증폭된 수면 변동은 스파형 구조물의 운동에 영향을 주어 구조물의 작업성 및 안정성을 저하시키거나 계류시스템에 치명적인 손상을 가져온다. 문풀 내부 유체 영역에서 발생하는 공진모드는 전후방향 거동인 슬로싱(Sloshing) 모드와 상하방향 거동인 피스톤(Piston) 모드로 구분된다. 특히 문풀 내부에서 형태의 변화 없이 상하로 진동하는 피스톤 모드의 공진 주파수는 해양에서 빈번히 발생하는 해양파의 주파수 범위내에 위치하므로 문풀 내의 피스톤 공진모드는 스파형 해양구조물을 설계할 때 반드시 고려되어야 한다.
Faltinsen(1990)는 문풀을 감쇠가 없는 질량-스프링 시스템으로 가정하여 피스톤 모드의 공진 주파수를 구하는 간단한 식을 제안하였다. Fukuda(1977)는 공진시 상하로 진동하는 내부 유체의 부가질량 효과를 고려하여 피스톤 모드의 공진 주파수에 대한 공식을 제안하였다. Molin(2001)은 선형 포텐셜 이론을 이용하여 사각형 문풀에 대한 슬로싱과 피스톤 모드의 공진 주파수를 추정하기 위한 간략식을 제안하였다.
피스톤 공진 모드에서의 부가질량(Added mass)과 감쇠계수(Damping coefficient)는 흥미로운 거동특성을 보여준다. 즉, 부가질량은 공진 주파수에서 갑자기 증가하다가 급격히 작아져 음의 값을 갖는다. 이때 방사감쇠계수는 피크값을 보인다. 해양유체역학 분야에서 유체동역학적 공진으로 공진 주파수에서 부가질량이음의 값을 갖는 예들은 많다. 안벽에 계류된 선박이나 병렬로 계류된 선박에서 두 물체 사이에는 공진을 유발하는 2차원 내부 유체 영역이 존재한다. Cho and Pyo(2009)는 부분 반사를 허용하는 안벽 앞에 계류된 부유체의 운동응답 특성을 고유함수전개법을 사용하여 해석하였다. Hong et al.(2005)은 병렬로 계류된 선체에 미치는 파력과 선체 사이의 내부 유체 영역내의 파형을 경계요소법(Boundary element method)을 사용하여 계산하였다. Faltinsen et al.(2007)은 두 선체간의 병렬계류와 쌍동선 등에서 발생하는 피스톤 모드의 공진 현상을 실험과 함께 선형 포텐셜 이론을 적용하여 연구하였다. 또한 Kristiansen and Faltinsen(2008)은 보텍스(Vortex) 추적법을 통하여 부유체 사이의 내부 유체 영역에서 발생하는 물입자의 와류의 형태 및 진행경로를 수치 계산하였다. Koo and Lee(2011)은 비선형 수치 파동 수조를 이용하여 2차원 주상체의 수직운동에 따른 내부 유체 영역내의 파형과 수평/수직 방향 동유체력을 시간영역에서 해석하였다.
축대칭 구조를 갖는 내부가 빈 원기둥 형상에 대한 연구로 Cho(2002)은 진동수주형(Oscillating water column) 파력발전장치에 대한 회절문제와 방사문제를 풀어 파랑에너지 흡수효율을 계산하였다. Mavrakos(1985, 1988)는 내부가 빈 원기둥에 대한 회절문제와 방사문제를 고유함수전개법을 이용하여 해석하였고, Sioris and Memos(1999)는 Green 함수법을 사용하였는데 Mavrakos(1985, 1988)의 결과와는 차이를 보여주었다. Chau and Yeung(2010)는 파랑에너지 변환장치로 내부가 빈 원기둥 모델을 채택하여 수직 원기둥 운동에 의한 방사문제를 고유함수전개법으로 풀어 부가질량과 방사감쇠계수를 구하였다.
본 연구에서는 내부가 빈 원기둥의 수직 운동특성을 주파수 영역에서 해석하였다. 선형포텐셜 이론에 근거한 경계치 문제를 회절문제와 방사문제로 각각 나누어 고유함수전개법을 사용하여 파기진력과 동유체력을 구하고 이 결과를 운동방정식에 대입하여 수직운동 변위와 내부 유체의 파형과 유량을 살펴보았다. 특히 피스톤 공진 모드에서의 내부가 빈 원기둥의 수직운동 특성을 살펴보았고, 원기둥의 형상(반경, 흘수)이 부가질량과 방사감쇠계수, 그리고 수직운동 변위에 미치는 영향을 조사하였다.
본 연구에서는 외경
여기서
고유함수전개법을 적용하기 위하여 유체영역을 Fig. 1과 같이 영역 (I), (II), (III)으로 나눈다. 영역 (I)은
Laplace방정식과 함께 자유표면 및 해저면 경계조건 그리고 방사조건을 만족하는 영역 (I)의 회절포텐셜(
여기서
영역 (I)에서의 고유값(
식 (3)에 정의된 고유함수
여기서
영역 (II)는 Laplace방정식과 원기둥 바닥과 해저면 바닥에서의 경계조건식을 만족한다. 주어진 경계조건식들을 만족하는 영역 (II)에서의 속도포텐셜은 다음과 같다.
여기서
여기서
영역 (III)에서의 회절포텐셜은 다음과 같다.
식 (2), (5), (7)의 미지수
먼저
여기서
같은 방법으로
여기서
여기서
식 (10)의 양변에
같은 방법으로
여기서
식 (8)과 (9)를 식 (11)과 (12)에 대입하여 정리하면 미지수
여기서
식 (13)에 주어진 대수방정식을 풀어 미지수
내부가 빈 원기둥에 작용하는 수직방향의 파기진력 (
원기둥의 수직운동에 의한 방사문제는 회절문제와 같은 방법으로 유체영역을 영역 (I), (II), (III)으로 나누어 각 영역에서의 해당 경계조건식들을 만족하는 속도포텐셜을 구한다.
식 (15)에서 특별해는 아래와 같이 주어진다.
앞에서 살펴본 회절문제와 같은 방법으로
여기서
영역 (II)의 미지수 , 는 아래 식을 이용하여 구할 수 있다.
원기둥의 수직운동에 의한 동유체력은 영역 (II)의 동압력을 원기둥 바닥면에 대하여 적분하여 다음과 같이 구할 수 있다.
여기서
원기둥 내부의 자유 수면의 변동유량(
여기서
파기진력과 동유체력을 아래 주어진 운동방정식에 대입하여 원기둥의 수직운동 변위를 구한다.
여기서
Fig. 3은 입사파의 진폭으로 무차원화시킨 내부가 빈 원기둥 수직운동의 진폭 비인 RAO(|
고유함수전개법을 사용하여 얻은 해석해의 수렴도를 살펴보기 위하여 고유함수의 개수(
Convergence test of (Cf, ) with the number of eigenfunctions (N) for M = 10, h/a = 5.0, d/a = 0.5, b/a = 0.5
본 해석결과의 타당성을 검증하기 위하여 Chau and Yeung(2010)의 계산결과와 본 계산결과를 Fig. 4에서 비교하였다. 내부가 빈 원기둥의 형상은
Fig. 5는 Fig. 4와 같은 계산 조건하에서 파기진력을 무차원화된 파장(
Fig. 6은
원기둥의 흘수는 내부 유체 영역이 외부영역과 차단되는 정도를 나타내는 변수이다. 원기둥의 흘수가 커질수록 내부 유체 영역으로 들어온 파가 외부로 빠져나지 못하여 방사 감쇠력은 줄어든다. 이러한 현상은 내부 유체가 피스톤과 같이 상하로 진동하는 피스톤 공진 모드에서 더욱 뚜렷이 나타난다. Fig. 7은 원기둥 흘수변화(
원기둥 흘수변화(
Fig. 9는 내부 유체 영역내 수면의 변동유량을
내부가 빈 원기둥의 동유체력(부가질량, 감쇠계수)과 파기진력, 그리고 수직운동 변위를 원기둥의 반경과 흘수를 변화시키면서 살펴본 결과 다음과 같은 결론을 얻었다.
(1) 자유감쇠실험을 통하여 구한 점성 감쇠계수를 포함한 수직운동 RAO의 해석결과는 실험결과와 잘 일치하였다.
(2) 계산 주파수 영역내에는 2개의 공진이 존재한다. 즉, 입사파의 주파수와 원기둥의 고유 주파수와 일치될 때 나타나는 시스템 공진과 원기둥 내부 유체의 공진으로 생긴 피스톤 공진이다. 피스톤 공진 모드에서 수직운동 변위는 시스템 공진보다 상대적으로 큰 값을 보인다.
(3) Fukuda(1997)가 제안한 피스톤 모드에서의 공진 주파수를 구하는 간략식은 현재 해석결과와 모형실험에서 나타난 공진주파수를 잘 따라감을 확인하였다.
(4) 피스톤 공진 모드에서 부가질량은 음의 값을 가지며 그 때 감쇠계수는 피크값을 갖는다. 특히, 원기둥의 흘수가 증가할수록 원기둥 내부로 들어온 파는 외부로 빠져나가지 못하여 방사감쇠계수는 줄어들며 이에 따라 수직운동의 RAO가 줄어드는 것을 볼 수 있다.