Generally, global performance analysis in offshore platforms is performed using potential-based numerical tools, which neglect hydrodynamic viscous effects. In comparison with the potential theory, computational fluid dynamics (CFD) methods can take into account the viscous effects by solving the Navier-Stokes equation using the finite-volume method. The open-source field operation and manipulation (OpenFOAM) C++ libraries are employed for a finite volume method (FVM) numerical analysis. In this study, in order to apply CFD to the global performance analysis of a hull-mooring coupled system, we developed a numerical wave basin to analyze the global performance problem of a floating body with a catenary mooring system under regular wave conditions. The mooring system was modeled using a catenary equation and solved in a quasi-static condition, which excluded the dynamics of the mooring lines such as the inertia and drag effects. To demonstrate the capability of the numerical basin, the global performance of a barge with four mooring lines was simulated under regular wave conditions. The simulation results were compared to the analysis results from a commercial mooring analysis program, Orcaflex. The comparison included the motion of the barge, catenary shape, and tension in the mooring lines. The study found good agreement between the results from the developed CFD-based numerical calculation and commercial software.
특정 해역에 머무르며 오랜 기간 작업을 수행하는 부유식 해양구조물의 특성상 정확한 운동응답특성 및 위치유치성능이 요구된다. 특히 위치 유지 성능은 구조물의 가동률이나 작업능률에 크게 영향을 미치기 때문에 매우 중요한 요소이다. 본 논문에서는 파랑 중 해양구조물의 위치를 유지하기 위해 계류 시스템을 사용하였으며, 그 중에서 현수선 계류(catenary mooring)를 사용하였다.
부유식 해상구조물 및 계류 시스템의 설계 및 성능 검증을 위해서는 일반적으로 모형실험 및 포텐셜(potential)기반의 수치 기법을 이용한 연구가 보편적이다. Hong et al.(1988)은 계류된 반잠수식 시추선의 규칙파 중 운동특성과 현수선 계류선의 인장력을 실험을 통해 계산하고 이를 3차원 포텐셜 이론의 수치해석을 통해 비교하였다. Sirnivas et al.(2016)은 부유체와 계류 시스템의 연성운동에 대해 포텐셜 기반인 프로그램을 개발하고 이를 상용프로그램과 실험을 통해 비교 및 검증하였다.
최근 전산유체역학의 발달로 인해 유체동역학(Hydrodynamic) 분야에서 모형실험의 한계를 보완하고 포텐셜 기반의 수치해석과 비교하여 보다 정도 높은 해석을 위한 접근이 가능해졌다. Ok et al.(2017)은 오픈소스 CFD 프로그램을 사용하여 파랑 중 단일선체 및 병렬 배치된 두 부유체의 운동에 관해 수치해석을 수행하였다. Wu et al.(2016)은 자체 계류 시스템 모듈을 개발하고 이를 상용 CFD 프로그램인 STAR CCM+과 연성을 통한 부유체와 계류 시스템의 운동특성을 파악하였다.
본 연구에서는 현수선 계류 시스템의 준정적 해석 라이브러리를 개발하고 이를 오픈소스 CFD 프로그램인 Open-source Field Operation And Manipulation(OpenFOAM)과의 연성을 통해 부유체와 계류 시스템 간의 양방향 연성 모듈을 개발하였다. 여기서 계류선의 관성과 양력 효과는 고려되지 않았다. 개발된 모듈에 대한 비교 및 검증을 위해 부유체 운동 계산을 위한 상용프로그램인 AQWA와 계류선의 장력 계산을 위한 상용프로그램인 Orcaflex를 이용하여 수치해석을 수행하였다. 두 프로그램을 선택한 이유는 AQWA는 포텐셜 기반의 경계요소법을 이용한 파랑 중 부유체의 운동계산 프로그램으로서 유사한 기능을 가진 다른 프로그램에 비해 범용성이 가장 높기 때문이고, Orcaflex는 계류선의 동적 거동을 계산하는 프로그램 중에 비선형 요소들을 추가할 수 있는 확장성이 가장 높기 때문에 추후 본 연구의 확장에 대응하기 위함이다.
본 유동장의 유체는 비압축성·점성유체이며 유동장내의 유동에 대한 속도와 압력을 계산하기 위해 연속방정식과 Navier-Stokes 방정식을 사용하였다. 각각의 방정식은 식 (1), (2)와 같다.
여기서
여기서 아래 첨자
강체로 가정한 부유체의 6자유도 운동은 다음과 같이 나타낼 수 있다.
식 (6)은
병진운동을 기준으로
Fig. 1에서는 입사파, 부유체, 계류 시스템 간의 상관관계에 따른 모식도를 나타낸다. 입사파는 선형규칙파이고, 파랑으로 야기된 부유체의 표류운동으로부터 위치제어를 위해 현수선형태의 계류 시스템을 적용하였다. 즉, 시간 영역에서 파랑 중 부유체의 거동에 대하여 페어리드(Fairlead) 좌표가 현수선 계류 시스템 라이브러리(Library)의 입력값이 되고 부유체의 표류운동 제어를 위한 계류선의 장력이 현수선 계류 시스템 라이브러리의 출력값이 되는 알고리즘이다.
Fig. 2는 적용된 현수선 방정식의 형상을 나타내며, 식(10), (11)과 같이 정의된다.
여기서
Step. 1
식 (10), (11)을 각각
Step. 2
Newton raphson & Jacobian 방법을 이용하여 매 시간 간격 마다 변화하는 페어리드 좌표로부터 계류선의 장력을 계산하기 위한 수치해석을 수행한다.
단, 매 타임스텝마다 장력계산을 위한 반복계산의 범위는 이다.
Step. 3
좌표변환 단계이다. 앞서 계산된
유체의 VOF 이송 방정식과 지배방정식은 유한 체적법(Finite Volume Method)에 의하여 이산화 된다. 지배방정식 중 시간 항은 1차 정확도의 Euler scheme, 공간 항은 2차 정확도의 Linear upwind scheme을 사용하였다. 속도와 압력의 연성을 위해서 SIMPLE 알고리즘과 PISO 알고리즘을 결합한 PIMPLE 알고리즘을 사용하였다. 수치해석 과정은 Fig. 4와 같다.
해석 도메인은 25
[Table 1.] Regular wave condition
Regular wave condition
파랑에 의한 부유체 운동응답의 해석에 앞서 규칙파에 대한 격자 민감도 테스트를 수행하였다. Fig. 5는 주기가 1.0
Fig. 6으로부터 입력한 파 진폭과 계산된 파 진폭을 비교한 결과, 격자 크기가 가장 큰 Test1에서는 약 4%의 오차가 발생하였고, 격자 크기가 가장 작은 Test4에서는 약 0.5%오차가 발생하였다. 계산 정도와 계산시간의 합리적인 절충안으로써 약 0.7%의 오차를 보이는 Test3 격자 크기를 선정하여 수치해석을 수행하였다. Fig. 7은 유동장의 격자분포를 나타낸다. 자유수면의 정도 높은 계산을 위해 파가 지나가는 영역에 가장 조밀한 격자를 생성하였으며 앞서 수행한 격자 수렴도 테스트를 바탕으로 한 파장을 기준으로
부유체의 운동은 경계요소법(Boundary Element Method, BEM)을 이용하는 ANSYS AQWA를 사용하여 계산하였다. 이상유체를 가정한 선형 포텐셜 이론을 기반으로 단위진폭을 가지는 선형파에 대한 부유체의 운동 응답 및 파랑 하중을 계산 하였다. 주파수 영역에서 부가질량계수(Added mass coefficient), 방사감쇠계수(Radiation damping coefficient), 파진기력(Wave exciting force) 등을 산출하고 이를 통해 입사파의 주파수에 대한 운동의 크기를 나타내는 진폭 응답 함수(Motion RAO)를 계산하였다.
유한요소법(Finite Element Method, FEM)을 이용하는 상용프로그램인 Orcaflex를 통하여 부유체와 계류 시스템 간의 시간 영역 연성해석을 수행하였다. Orcaflex를 통한 수치해석은 전체 계류선의 세그먼트에 따라 노드점과 스프링의 개수가 정해지고, 이를 통해, 각 노드점에 질량을 부과하는 집중질량법(Lumped mass)을 기반으로 한다. Fig. 8과 같이 총 4개의 계류선이 선수 선미로부터 좌현, 우현으로 향하도록 계류선을 배치하였고 페어리드와 앵커링 지점의 수평 방향 거리와 수심의 비를 의미하는 계류선반경(Foot print radius)은 2배로 설계하였다.
또한, 계류선의 양력 및 항력계수는 각각 2.4, 1.5로 계산하였다. 이는 DNVGL-OS-E301(2015)에서 제공하는 값을 기준으로 산정하였다.
파랑 하중을 적용하기에 앞서 정적 평형상태(Static equilibrium state)에서 계류선의 형상 및 페어리드에 작용하는 장력을 비교하였다. Fig. 9는 페어리드에서 앵커링 지점까지 총 50개의 세그먼트로 나눈 계류선의 형상을 나타내고 있다. 개발된 OpenFOAM과 상용프로그램인 Orcaflex 간의 전체적인 형상이 잘 일치함을 확인할 수 있다.
Fig. 10은 페어리드에 작용하는 계류선의 초기장력을 비교한 그림이다. Line1, Line3은 선수, 선미의 대각선 방향에 위치한 계류선을 나타낸다. 현수선 방정식으로부터 계산된 OpenFOAM 계류 시스템의 초기장력은 Orcaflex의 결과와 잘 일치하는 것을 확인할 수 있다.
입사파의 각도가 0°(선수파)인 조건에서 입사파 주기에 따른 부유체 운동 응답 특성 및 계류선의 장력변화에 대한 비교를 수행하였다(Davidson and Ringwood, 2017).
Fig. 11은 주기가 1.6s인 조건에서 OpenFOAM, Orcaflex로부터 계산된 부유체의 운동 응답에 대한 시계열을 나타낸다. 정량적 비교를 위해 두 수치해석 간의 위상(Phase)을 맞추어 비교하였다. Heave와 Pitch의 경우 두 개의 시그널이 정량적으로 잘 일치하는 것을 확인할 수 있다. Surge 응답에서는 정량적으로 일부 상이한 결과가 나타나지만, 개발된 계류모듈 코드는 Orcaflex의 결과와 정량적으로 유사한 결과를 보여주는 것을 확인하였다. 또한, Surge 응답에서 부유체와 계류시스템 간의 연성 운동으로 부터 기인하는 하모닉(Harmonic) 현상이 관측되고 있는데, 이를 해석하고자 각각의 시간이력을 파워스펙트럼으로 Fig. 12에 나타내었다.
Fig. 12는 입사파의 주기에 따른 Surge 운동의 파워스펙트럼을 나타낸다. 스펙트럼마다 2개의 피크가 확연히 구분되어 관찰되는데 상대적으로 고주파수에 위치한 피크는 입사파의 주파수를 의미하고 저주파수인 2.2rad/s에 위치하는 피크는 계류시스템의 고유진동수와 일치하는 응답을 의미한다. OpenFOAM과 Orcaflex의 결과들 간에 진폭의 차이는 존재하지만 각각의 파 주파수와 계류시스템의 고유진동수에 의한 피크위치는 정확히 예상하는 결과와 일치하는 것을 확인 할 수 있다.
Fig. 13, 14는 계류선의 장력에 대해 평균값과 최대값을 나타내었다. 장력의 평균값은 두 방법 간에 큰 차이를 보이지 않았으나 장력의 최대값은 특정 파주기에서 약간의 차이를 보이고 있는데, 이는 개발된 계류시스템용 모듈이 현수선방정식을 사용함으로써 평면내의 운동만 표현이 가능하고, 동역학적인 계류선의 거동을 표현하는데 한계가 있어서 생기는 것으로 판단된다. 이는 추후 모형실험을 통해 추가적인 비교를 수행하거나, 연성모듈의 현수선방정식을 집중질량법으로 개량하는 방법을 통해 검증이 이루어질 예정이다.
본 연구에서는 오픈소스 CFD 프로그램인 OpenFOAM를 이용하여 부유체와 계류 시스템간의 양방향 연성 모듈을 개발하고 이를 상용프로그램의 수치해석결과와 비교하여 검증하였다. 정적평형 상태에서 계류선의 형상 및 초기장력을 비교했으며, 선수파 조건의 규칙파 조건에서 부유체의 거동 및 페어리드 지점의 장력을 주파수 영역에서의 스펙트럼 및 시간영역에서의 최대값과 평균값을 통해 정량적으로 평가하였다. 본 연구로부터 다음과 같은 결론을 얻을 수 있다.
(1) 개발된 OpenFOAM 코드의 Heave, Pitch 운동은 상용프로그램인 Orcaflex간에 정량적으로 잘 일치하는 것을 확인하였다.
(2) 개발된 OpenFOAM 코드의 Surge 운동 및 계류선의 장력은 상용프로그램인 Orcaflex와 정량적으로 허용 범위 내의 차이를 보이며 같은 경향을 보임을 확인하였다.
(3) 본 논문에서는 점성효과를 고려하는 CFD를 이용하여 부유체와 계류 시스템간의 연성 운동을 구현하였다. 추후 모형실험 또는 연성모듈의 개선을 통해 추가적인 비교 및 검증이 필요하며, 더 나아가 이를 바탕으로 극한 환경 조건과 계류선의 동적 효과를 고려한 연구가 요구된다.
