선택한 단어 수는 입니다.
최소 단어 이상 선택하여야 합니다.
최소 단어 이상 선택하여야 합니다.
선택한 단어 수는 30입니다.
최대 10 단어까지만 선택 가능합니다.
최대 10 단어까지만 선택 가능합니다.
문제 정의
- 본 연구에서는 현수선 계류 시스템의 준정적 해석 라이브러리를 개발하고 이를 오픈소스 CFD 프로그램인 Open-source Field Operation And Manipulation(OpenFOAM)과의 연성을 통해 부유체와 계류 시스템 간의 양방향 연성 모듈을 개발하였다. 여기서 계류선의 관성과 양력 효과는 고려되지 않았다.
제안 방법
- (3) 본 논문에서는 점성효과를 고려하는 CFD를 이용하여 부유체와 계류 시스템간의 연성 운동을 구현하였다. 추후 모형 실험 또는 연성모듈의 개선을 통해 추가적인 비교 및 검증이 필요하며, 더 나아가 이를 바탕으로 극한 환경 조건과 계류선의 동적 효과를 고려한 연구가 요구된다.
- 5%오차가 발생하였다. 계산 정도와 계산시간의 합리적인 절충안으로써 약 0.7%의 오차를 보이는 Test3 격자 크기를 선정하여 수치해석을 수행하였다. Fig.
- 자유수면의 정도 높은 계산을 위해 파가 지나가는 영역에 가장 조밀한 격자를 생성하였으며 앞서 수행한 격자 수렴도 테스트를 바탕으로 한 파장을 기준으로 x방향 약 130개, z방향 12개의 격자를 적용하였다. 또한, 부유체의 운동 응답의 정도를 높이기 위해 부유체 주위에 격자 조밀도를 증가시켰다. 시간영역 계산을 위한 시간간격은 T(주기)/200s이고, 총 사용된 격자의 수는 1,300,000개이다.
- )는 해저면까지의 계류선 길이(Suspend length)에 해당하는 노드(Node)의 x방향과 z방향의 변위, TH와 TV는 수평, 수직 방향의 장력, w는 계류선의 단위 길이 당 무게, K는 계류선의 축방향 강성을 의미한다. 매 시간 간격 마다 계류선의 해저면 접촉지점(Touchdown point)과 Suspend length를 계산하고, 이로부터 현수선 방정식을 이용하여 정적 평형상태의 형상과 페어리드에 작용하는 장력을 계산한다. 여기서, x, z 방향의 정적 평형상태의 형상 및 장력을 나타내는 식 (10), (11)을 만족시키기 위해 비선형 수치해석 방법인 Newton-raphson & Jacobian 방법을 이용하였다.
- 본 연구에서는 오픈소스 CFD 프로그램인 OpenFOAM를 이용하여 부유체와 계류 시스템간의 양방향 연성 모듈을 개발하고 이를 상용프로그램의 수치해석결과와 비교하여 검증하였다. 정적평형 상태에서 계류선의 형상 및 초기장력을 비교했으며, 선수파 조건의 규칙파 조건에서 부유체의 거동 및 페어리드 지점의 장력을 주파수 영역에서의 스펙트럼 및 시간영역에서의 최대값과 평균값을 통해 정량적으로 평가하였다.
- 유한요소법(Finite Element Method, FEM)을 이용하는 상용프로그램인 Orcaflex를 통하여 부유체와 계류 시스템 간의 시간 영역 연성해석을 수행하였다. Orcaflex를 통한 수치해석은 전체 계류선의 세그먼트에 따라 노드점과 스프링의 개수가 정해지고, 이를 통해, 각 노드점에 질량을 부과하는 집중질량법(Lumped mass)을 기반으로 한다.
- 부유체의 운동은 경계요소법(Boundary Element Method, BEM)을 이용하는 ANSYS AQWA를 사용하여 계산하였다. 이상유체를 가정한 선형 포텐셜 이론을 기반으로 단위진폭을 가지는 선형파에 대한 부유체의 운동 응답 및 파랑 하중을 계산 하였다. 주파수 영역에서 부가질량계수(Added mass coefficient), 방사감 쇠계수(Radiation damping coefficient), 파진기력(Wave exciting force) 등을 산출하고 이를 통해 입사파의 주파수에 대한 운동의 크기를 나타내는 진폭 응답 함수(Motion RAO)를 계산하였다.
- 입사파의 각도가 0∘(선수파)인 조건에서 입사파 주기에 따른 부유체 운동 응답 특성 및 계류선의 장력변화에 대한 비교를 수행하였다(Davidson and Ringwood, 2017).
- 6s인 조건에서 OpenFOAM, Orcaflex로부터 계산된 부유체의 운동 응답에 대한 시계열을 나타낸다. 정량적 비교를 위해 두 수치해석 간의 위상(Phase)을 맞추어 비교하였다. Heave와 Pitch의 경우 두 개의 시그널이 정량적으로 잘 일치하는 것을 확인할 수 있다.
- 본 연구에서는 오픈소스 CFD 프로그램인 OpenFOAM를 이용하여 부유체와 계류 시스템간의 양방향 연성 모듈을 개발하고 이를 상용프로그램의 수치해석결과와 비교하여 검증하였다. 정적평형 상태에서 계류선의 형상 및 초기장력을 비교했으며, 선수파 조건의 규칙파 조건에서 부유체의 거동 및 페어리드 지점의 장력을 주파수 영역에서의 스펙트럼 및 시간영역에서의 최대값과 평균값을 통해 정량적으로 평가하였다. 본 연구로부터 다음과 같은 결론을 얻을 수 있다.
- 주파수 영역에서 부가질량계수(Added mass coefficient), 방사감 쇠계수(Radiation damping coefficient), 파진기력(Wave exciting force) 등을 산출하고 이를 통해 입사파의 주파수에 대한 운동의 크기를 나타내는 진폭 응답 함수(Motion RAO)를 계산하였다.
- 입사파는 진폭, 방향, 주기는 Table 1과 같다. 파고는 2cm, 파방향은 선수파 방향이고, 주기는 1.0초부터 2.0초까지 0.2초 간격으로 총 6가지 조건으로 수치해석을 수행하였다.
- 파랑 하중을 적용하기에 앞서 정적 평형상태(Static equilibrium state)에서 계류선의 형상 및 페어리드에 작용하는 장력을 비교하였다. Fig.
- 파랑에 의한 부유체 운동응답의 해석에 앞서 규칙파에 대한 격자 민감도 테스트를 수행하였다. Fig.
- 해석 도메인은 25m(L)⨉ 1m(B) ⨉ 1m(D)로 모델링 하였다. 파랑의 수치 오차를 고려하여 대상 부유체는 조파가 시작되는 경계면으로부터 한 파장 떨어진 곳에 위치하며 유동장의 앞, 뒤, 양옆의 경계면에 파의 감쇠 영역을 설정하여 벽면효과를 최소화하도록 설정하였다. 입사파 조건은 개발된 프로그램의 비교 및 검증을 위해 비선형성을 최소화하고자 하였으며, 이에 따라 입사파의 파경사각을 3∘미만으로 선정하였다.
대상 데이터
- 7은 유동장의 격자분포를 나타낸다. 자유수면의 정도 높은 계산을 위해 파가 지나가는 영역에 가장 조밀한 격자를 생성하였으며 앞서 수행한 격자 수렴도 테스트를 바탕으로 한 파장을 기준으로 x방향 약 130개, z방향 12개의 격자를 적용하였다. 또한, 부유체의 운동 응답의 정도를 높이기 위해 부유체 주위에 격자 조밀도를 증가시켰다.
- 해석 도메인은 25m(L)⨉ 1m(B) ⨉ 1m(D)로 모델링 하였다. 파랑의 수치 오차를 고려하여 대상 부유체는 조파가 시작되는 경계면으로부터 한 파장 떨어진 곳에 위치하며 유동장의 앞, 뒤, 양옆의 경계면에 파의 감쇠 영역을 설정하여 벽면효과를 최소화하도록 설정하였다.
데이터처리
- 여기서 계류선의 관성과 양력 효과는 고려되지 않았다. 개발된 모듈에 대한 비교 및 검증을 위해 부유체 운동 계산을 위한 상용프로그램인 AQWA와 계류선의 장력 계산을 위한 상용프로그램인 Orcaflex를 이용하여 수치해석을 수행하였다. 두 프로그램을 선택한 이유는 AQWA는 포텐셜 기반의 경계요소법을 이용한 파랑 중 부유체의 운동계산 프로그램으로서 유사한 기능을 가진 다른 프로그램에 비해 범용성이 가장 높기 때문이고, Orcaflex는 계류선의 동적 거동을 계산하는 프로그램 중에 비선형 요소들을 추가 할 수 있는 확장성이 가장 높기 때문에 추후 본 연구의 확장에 대응하기 위함이다.
이론/모형
- Newton raphson & Jacobian 방법을 이용하여 매 시간 간격 마다 변화하는 페어리드 좌표로부터 계류선의 장력을 계산하기 위한 수치해석을 수행한다.
- 유한요소법(Finite Element Method, FEM)을 이용하는 상용프로그램인 Orcaflex를 통하여 부유체와 계류 시스템 간의 시간 영역 연성해석을 수행하였다. Orcaflex를 통한 수치해석은 전체 계류선의 세그먼트에 따라 노드점과 스프링의 개수가 정해지고, 이를 통해, 각 노드점에 질량을 부과하는 집중질량법(Lumped mass)을 기반으로 한다. Fig.
- 특히 위치 유지 성능은 구조물의 가동률이나 작업능률에 크게 영향을 미치기 때문에 매우 중요한 요소이다. 본 논문에서는 파랑 중 해양구조물의 위치를 유지하기 위해 계류 시스템을 사용하였으며, 그 중에서 현수선 계류(catenary mooring)를 사용하였다.
- 본 유동장의 유체는 비압축성·점성유체이며 유동장내의 유동에 대한 속도와 압력을 계산하기 위해 연속방정식과 Navier-Stokes 방정식을 사용하였다.
- 부유체의 운동은 경계요소법(Boundary Element Method, BEM)을 이용하는 ANSYS AQWA를 사용하여 계산하였다. 이상유체를 가정한 선형 포텐셜 이론을 기반으로 단위진폭을 가지는 선형파에 대한 부유체의 운동 응답 및 파랑 하중을 계산 하였다.
- 지배방정식 중 시간 항은 1차 정확도의 Euler scheme, 공간 항은 2차 정확도의 Linear upwind scheme을 사용하였다. 속도와 압력의 연성을 위해서 SIMPLE 알고리즘과 PISO 알고리즘을 결합한 PIMPLE 알고리즘을 사용하였다. 수치해석 과정은 Fig.
- 여기서, x, z 방향의 정적 평형상태의 형상 및 장력을 나타내는 식 (10), (11)을 만족시키기 위해 비선형 수치해석 방법인 Newton-raphson & Jacobian 방법을 이용하였다.
- 유체 간의 구분은 격자 내에 차지하는 체적비율인 α를 통해 자유수면을 수치적으로 표현하는 VOF(Volume of fluid) 기법을 사용하였다.
- 유체의 VOF 이송 방정식과 지배방정식은 유한 체적법(Finite Volume Method)에 의하여 이산화 된다. 지배방정식 중 시간 항은 1차 정확도의 Euler scheme, 공간 항은 2차 정확도의 Linear upwind scheme을 사용하였다. 속도와 압력의 연성을 위해서 SIMPLE 알고리즘과 PISO 알고리즘을 결합한 PIMPLE 알고리즘을 사용하였다.
성능/효과
- (1) 개발된 OpenFOAM 코드의 Heave, Pitch 운동은 상용프로그램인 Orcaflex간에 정량적으로 잘 일치하는 것을 확인하였다.
- (2) 개발된 OpenFOAM 코드의 Surge 운동 및 계류선의 장력은 상용프로그램인 Orcaflex와 정량적으로 허용 범위 내의 차이를 보이며 같은 경향을 보임을 확인하였다.
- Fig. 6으로부터 입력한 파 진폭과 계산된 파 진폭을 비교한 결과, 격자 크기가 가장 큰 Test1에서는 약 4%의 오차가 발생하였고, 격자 크기가 가장 작은 Test4에서는 약 0.5%오차가 발생하였다. 계산 정도와 계산시간의 합리적인 절충안으로써 약 0.
- 2rad/s에 위치하는 피크는 계류시스템의 고유진동수와 일치하는 응답을 의미한다. OpenFOAM과 Orcaflex의 결과들 간에 진폭의 차이는 존재하지만 각각의 파주파수와 계류시스템의 고유진동수에 의한 피크위치는 정확히 예상하는 결과와 일치하는 것을 확인 할 수 있다.
- Heave와 Pitch의 경우 두 개의 시그널이 정량적으로 잘 일치하는 것을 확인할 수 있다. Surge 응답에서는 정량적으로 일부 상이한 결과가 나타나지만, 개발된 계류모듈 코드는 Orcaflex의 결과와 정량적으로 유사한 결과를 보여주는 것을 확인하였다. 또한, Surge 응답에서 부유체와 계류시스템 간의 연성 운동으로부터 기인하는 하모닉(Harmonic) 현상이 관측되고 있는데, 이를 해석하고자 각각의 시간이력을 파워스펙트럼으로 Fig.
- 9는 페어리드에서 앵커링 지점까지 총 50개의 세그먼트로 나눈 계류선의 형상을 나타내고 있다. 개발된 OpenFOAM과 상용프로그램인 Orcaflex 간의 전체적인 형상이 잘 일치함을 확인할 수 있다.
- 13, 14는 계류선의 장력에 대해 평균값과 최대값을 나타내었다. 장력의 평균값은 두 방법 간에 큰 차이를 보이지 않았으나 장력의 최대값은 특정 파주기에서 약간의 차이를 보이고 있는데, 이는 개발된 계류시스템용 모듈이 현수선방정식을 사용함으로써 평면내의 운동만 표현이 가능하고, 동역학적인 계류선의 거동을 표현하는데 한계가 있어서 생기는 것으로 판단된다. 이는 추후 모형실험을 통해 추가적인 비교를 수행하거나, 연성모듈의 현수선방정식을 집중질량법으로 개량하는 방법을 통해 검증이 이루어질 예정이다.
- Line1, Line3은 선수, 선미의 대각선 방향에 위치한 계류선을 나타낸다. 현수선 방정식으로부터 계산된 OpenFOAM 계류 시스템의 초기장력은 Orcaflex의 결과와 잘 일치하는 것을 확인할 수 있다.
후속연구
- 개발된 모듈에 대한 비교 및 검증을 위해 부유체 운동 계산을 위한 상용프로그램인 AQWA와 계류선의 장력 계산을 위한 상용프로그램인 Orcaflex를 이용하여 수치해석을 수행하였다. 두 프로그램을 선택한 이유는 AQWA는 포텐셜 기반의 경계요소법을 이용한 파랑 중 부유체의 운동계산 프로그램으로서 유사한 기능을 가진 다른 프로그램에 비해 범용성이 가장 높기 때문이고, Orcaflex는 계류선의 동적 거동을 계산하는 프로그램 중에 비선형 요소들을 추가 할 수 있는 확장성이 가장 높기 때문에 추후 본 연구의 확장에 대응하기 위함이다.
- 장력의 평균값은 두 방법 간에 큰 차이를 보이지 않았으나 장력의 최대값은 특정 파주기에서 약간의 차이를 보이고 있는데, 이는 개발된 계류시스템용 모듈이 현수선방정식을 사용함으로써 평면내의 운동만 표현이 가능하고, 동역학적인 계류선의 거동을 표현하는데 한계가 있어서 생기는 것으로 판단된다. 이는 추후 모형실험을 통해 추가적인 비교를 수행하거나, 연성모듈의 현수선방정식을 집중질량법으로 개량하는 방법을 통해 검증이 이루어질 예정이다.
- (3) 본 논문에서는 점성효과를 고려하는 CFD를 이용하여 부유체와 계류 시스템간의 연성 운동을 구현하였다. 추후 모형 실험 또는 연성모듈의 개선을 통해 추가적인 비교 및 검증이 필요하며, 더 나아가 이를 바탕으로 극한 환경 조건과 계류선의 동적 효과를 고려한 연구가 요구된다.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.