In this paper, the wave run.up on a vertical cylinder is presented. Various cross sections of a cylinder were simulated using the panel method for various wave periods. Two.dimensional model tests were performed in a wave flume. The simulation results are compared with the test results. The simulation is based on the linear diffraction theory.
최근 개발도상국들의 급속한 산업화로 인한 화석에너지의 수요증가와 육상 화석에너지자원 고갈의 가속화의 영향으로 국제 유가 또한 급속한 상승세를 보이고 있다. 이러한 에너지자원 문제의 해소를 위하여 해상에너지자원 개발에 많은 기술개발과 투자가 이루어지고 있고, 연근해뿐만 아니라 심해 및 극지에서의 자원개발이 점차 구체화 되고 있다. 해상에 설치 운용되는 해양구조물의 경우 육상구조물과는 달리 파랑 및 바람과 같은 환경하중의 영향을 많이 받게 된다. 이에 따라 극한의 환경하중에 대한 해양구조물의 구조적 안전성과 극한 운동성능 파악의 중요성이 대두되고 있다.
해양구조물에 작용하는 환경하중은 파랑에 의한 하중이 지배적이다. 이에 따라 해양구조물에 작용하는 파랑 하중과 Wave runup의 추정은 구조물 설계를 위한 해석기법 적용에 있어서 가장 중요한 요소 중 하나로, 이에 대한 연구가 다수 수행 되었다. Stansberg and Braaten(2002), Kristiansen et al.(2004), Sung et al. (2011)에서는 바닥 부분이 절단된 연직 원형실린더(Truncated vertical circular cylinder)에 대하여 규칙파에 의한 Run-up 현상을 연구하였다. 특히 Sung et al. (2011)은 Run-up 현상뿐만 아니라 Run-up에 의해 야기되는 파랑하중 분석에 대한 연구도 수행하였다.
본 연구에서는 해양구조물 Column의 일반적인 형태인 원형 및 사각형 단면을 가지는 실린더에 작용하는 Wave run-up 현상에 대하여 연구를 수행 하였다. 실린더 형상의 변화 및 입사파의 주기 변화에 따른 Wave run-up 값의 변화를 확인하기 위하여 Panel method을 이용한 수치 시뮬레이션을 수행 하였다. 아울러 2차원 Wave flume에서 1/100로 축소된 각각의 모형에 대하여 모형실험을 수행하였고, 모형시험 결과와 시뮬레이션 결과를 비교하였다. 특히 원형실린더에 작용하는 Wave run-up 현상에 대해서는 산란파 이론을 이용한 계산 값을 함께 비교하였다.
이상유체, 비회전성 유체 그리고 선형 파 이론을 가정하였을 때, 입사파의 속도 퍼텐셜은 식 (1)과 같이 나타낼 수 있다.
여기서
여기서 𝜖
여기서 은 오더가
입사파의 퍼텐셜과 Diffraction 퍼텐셜의 합을 Dynamic free surface boundary condition에 적용하면 실린더 주위 Elevation은 식 (4)과 같이 표현된다.
실린더 형상에 따른 Wave run-up 현상의 변화를 확인하기 위하여 Fig. 1과 같은 단면을 가지는 세 가지 실린더를 사용하여 시뮬레이션을 수행 하였다. 그림에서와 같이 Model A는 직경이 16m인 원형 단면이며, Model B는 한 변의 길이가 16m인 정사각형의 네 모서리에 2.5m의 곡면을 가지는 단면이다. 또한 Model C는 Model B가 입사파의 진행 방향에 대하여 45° 회전된 형태의 단면을 가진다.
입사파의 주기의 변화에 따른 Wave run-up 현상의 변화를 알아보기 위하여 입사파의 주기를 7sec ~ 15sec까지 변화는 6가지 주기에 대해 시뮬레이션을 수행하였다. 시뮬레이션에 사용된 입사파의 조건은 Table 1과 같다. 시뮬레이션은 실린더가 수면 아래 24m, 수면 위 20m로 고정된 상태에서 수행되었고, 이때 수심은 60m이다. 시뮬레이션에 사용된 주요 치수는 Table 2와 같다.
[Table 1] Incident wave conditions for simulation
Incident wave conditions for simulation
[Table 2] Principal characteristics of model
Principal characteristics of model
수치 시뮬레이션은 BV(Bureau Veritas)운동해석 프로그램인 Hydro-Star V6.2를 사용하였다. HydroStar는 일반적으로 사용하는 운동해석 프로그램이지만 본 연구에서는 실린더를 인위적으로 고정시켜 Incident 퍼텐셜과 Diffraction 퍼텐셜만을 계산하여 Elevation을 계산하였다. 계산에 사용된 Mesh는 Fig. 2와 같다.
4.1 Experimental Model and Facility
시뮬레이션 결과와의 비교를 위해 모형 시험을 수행하였다. 모형실험은 부산대학교에 설치된 30m×0.6m×0.95m(
[Table 3] Specifications of wave maker
Specifications of wave maker
실험모델은 Wave flume의 크기를 고려하여 실선의 1/100로 축소된 모형을 사용하였으며 모형 실린더의 형상은 Fig. 4와 같다.
Wave run-up 현상을 파악하기 위한 파고 계측을 위하여, 총 8개의 파고계를 실린더 주변에 설치하였다. 파고계는 Fig. 5에서 보이는 바와 같이 실린더 표면에서 5mm 떨어진 지점에 0°에서 180°까지 45° 간격으로 5개(WGP, Wave gauge position #1~5), 실린더 센터에서 160mm 떨어진 지점에 0°에서 90°까지 45°간격으로 3개(WGP #6~8) 설치하였다. 파고계는 Fig. 6과 같이 센서 케이블의 직경이 2mm인 용량식 상대파고계를 사용하였다. 이는 실린더의 크기(160mm)에 비해 상대적으로 작은 크기를 가지게되어 센서 설치로 인한 실험오차를 줄이도록 하였다.
입사파의 주기의 변화에 따른 Wave run-up 현상의 변화를 알아보기 위한 실험은 실선 시뮬레이션과 같이 입사파의 주기가 7sec에서 15sec 사이의 6개 주기에 대해 실험을 수행하였다. 또한 각각의 주기에서 입사파의 파형경사(Wave steepness) 변화에 따른 Wave run-up의 특징을 알아보기 위하여 파형경사(
[Table 4] Incident wave conditions and amplitude
Incident wave conditions and amplitude
입사파의 주기 변화 및 실린더 단면 형상의 변화에 따른 파고 계측은 Fig. 7에서와 같이 진행되었다. 계측된 파고 값은 Fig. 8에서와 같이 파가 완전히 안정화된 이후 10주기 이상의 통계해석을 통해 얻어진 값을 사용하였다.
Fig. 9~14까지 입사파 주기의 변화에 따른 Wave run-up 값을 각각의 실린더에 대해 반경별(WGP : #1-1 ~ #1-5)로 도시하였다. Model A의 경우 실험값과 해석 값, 시뮬레이션 값을 비교하였으며, Model B와 C는 단면 형상이 수학적으로 표현되지 않아 해석해가 존재하지 않기 때문에 실험값과 시뮬레이션 값만을 비교하였다.
Fig. 9~14의 (a)에서 보이는 것처럼 해석 값과 시뮬레이션 값은 잘 일치함을 확인 할 수 있다. 해석 값의 조건이 실린더가 바닥에 고정된 경우이고, 시뮬레이션 값은 Tuncated 조건이지만 실린더 바닥의 유체입자의 운동이 Wave elevation에는 크게 영향을 미치지 못하기 때문에 잘 일치된 결과를 보였다고 판단된다. 실험 값은 모든 경우에 대해 해석 및 시뮬레이션 값과 다소간의 차이를 보이는데 이는 수조의 폭이 실린더에 비해 상대적으로 좁고, 반사파의 영향으로 충분한 통계해석구간을 선정하지 못한 결과로 생각된다. 또한 해석 값 및 시뮬레이션 값은 선형 파 이론에 근거한 값이기 때문에 실험에서의 비선형성도 결과에 영향을 주었을 것으로 생각된다. 이에 따라 향후 비선형 및 물의 점성 등을 고려한 시뮬레이션을 수행할 예정이다. Wave run-up의 상대적인 값의 차이는 발생하였지만, 전반적인 경향은 잘 일치함을 확인할 수 있다. 또한 그림에서 확인 할 수 있는 것처럼 파형경사가 작을 경우 시뮬레이션 및 해석 값과 유사한 경향을 보이지만 파형경사가 증가하게 되면 다소간의 경향 차이가 발생하게 된다. 수치 및 해석 값은 선형이론에 기반을 둔 계산이기 때문에 파형경사가 증가하게 되면 선형이론으로 계산하기에 무리가 따르게 되는데 이러한 영향으로 인한 차이로 사료된다.
입사파의 주기변화에 따른 Wave run-up 현상의 변화를 살펴보면, 모든 경우에서 입사파의 주기가 길어짐에 따라 Run-up 현상이 감소하였다. 이는 입사파의 파장이 길어지게 되면 파장에 비해 실린더의 크기가 상대적으로 작아지게 되어 실린더로 인해 발생하는 Diffraction 퍼텐셜의 크기가 작아져서 Run-up 현상도 줄어들게 되는 것으로 생각된다.
Fig. 15에서 실린더 단면의 형상에 따른 Wave run-up 값을 비교하였다. 그림에서 보이는 것처럼 실린더의 형상에 따른 Runup은 Model A에 비해 B와 C가 다소 크게 나타나는 것으로 계산되었다. 입사파의 주기가 길어지게 되면 (a)~(c)에서 보이는 것처럼 Run-up 값의 차이가 점점 줄어들게 된다.
본 연구에서는 수직 실린더에 작용하는 Wave run-up 현상에 대한 연구를 수행하였다. 수치시뮬레이션을 통하여 실린더 형상 및 입사파의 변화에 따른 Run-up 현상의 변화를 확인하였고, 해석해 및 실험값과의 비교를 수행하였다.
결과 값들의 비교를 통하여 입사파의 주기가 짧을수록 Runup 현상이 가중됨을 확인하였다. 또한 본 연구에서 검토한 단면 형상은 Wave run-up 현상에 상대적으로 작은 영향을 미치게 되는 것을 확인하였다. 다시 말해 Wave run-up 현상은 실린더의 형상보다는 파장에 대한 실린더의 상대적인 크기에 많은 영향을 받는 것으로 판단된다.
