대한조선학회논문집 Journal of the Society of Naval Architects of Korea pissn:5-3, Vol. 5, No., pp. 338-38, August 5 eissn:87-7355, http://dx.doi.org/.37/snak.5.5..338 단파장 영역에서의 부가저항 해석 양경규 서민국 김용환 서울대학교 조선해양공학과 Analysis of Added Resistance in Short Waves Kyung-Kyu Yang Min-Guk Seo Yonghwan Kim Department of Naval Architecture and Ocean Engineering, Seoul National University This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. In this study, the added resistance of ships in short waves is systematically studied by using two different numerical methods - Rankine panel method and Cartesian grid method and existing asymptotic and empirical formulae. Analysis of added resistance in short waves has been preconceived as a shortcoming of numerical computation. This study aims to observe such preconception by comparing the computational results, particularly based on two representative three-dimensional methods, and with the existing formulae and experimental data. In the, a near-field method based on direct pressure integration is adopted. In the Cartesian grid method, the wave-body interaction problem is considered as a multiphase problem, and volume fraction functions are defined in order to identify each phase in a Cartesian grid. The computational results of added resistance in short waves using the two methods are systematically compared with experimental data for several ship models, including S75 containership, KVLCC and Series hulls (C B =.7,.8). The present study includes the comparison with the established asymptotic and empirical formulae in short waves. Keywords : Added resistance(부가저항), (랜킨패널법), Cartesian grid method(직교격자법), Short-wave approximation(단파장 근사법), EEDI(에너지효율설계지수). 서 론 전통적으로 선박의 저항 문제는 입사파가 없는 정수상태에서 전진하는 경우에 대해 해석되었다. 하지만, 파랑에 의한 부가저 항 문제는 최근 대두되고 있는 선박의 EEDI 도입과 파랑 중 저항 감소를 통한 고효율 선형 설계의 필요성과 함께 중요한 문제로 부각되고 있다 (Kim & Kim, ). 파랑에 의한 부가저항 문제 는 운동과 저항을 동시에 해석해야 하는 문제이기 때문에 해석에 많은 어려움이 있다. 실험의 경우, 계측장비나 실험방법, 모형선 의 크기 등에 따라 계측결과에 차이가 있고, 수치해석의 경우, 격자 및 수치기법 등에 따라 파랑 중 부가저항 값이 민감하게 변 화한다. 파랑 중 부가저항 문제는 실험에 기반을 두어 많은 연구가 수 행되었다. Fujii and Takahashi (975)와 Nakamura and Naito (977)는 S75 컨테이너선에 대해 실험을 진행하였고, 최근에는 Sadat-Hosseini, et al. (3)가 KVLCC 선형에 대해서 파랑 중 부가저항 실험을 수행한 바 있다. 수치적 방법을 이용한 파랑 중 부가저항 계산의 경우 이상유 동(ideal flow)에 근거한 포텐셜 기반 수치 기법들을 활용한 연구 가 활발하게 진행되어왔다. 이상유동에 근거한 파랑 중 부가저항 계산 기법은 크게 두 가지 방법으로 구분되며, 이는 모멘텀 보존 방법(far-field method)과 압력직접 적분방법(near-field method) 등이다. 모멘텀 보존방법은 Maruo (9)에 의해서 처음으로 제 안되었으며, 물체 표면에서의 압력을 직접 계산할 필요가 없기 때문에 상대적으로 간단하게 계산할 수 있다. 따라서 현재까지도 세장체 이론(slender-body theory)에 근거한 운동해석에 추가되 어 실제 공학적 문제에 많이 사용된다. 예를 들어, Kashiwagi, et al. (9)은 개선된 통합이론(enhanced unified theory)을 적용 한 모멘텀 보존방법을 사용하여 파랑 중 부가저항을 계산한 바 있다. 압력직접적분법의 경우에는 물체 표면의 차 동유체 압력 을 직접 적분하는 방법으로, 모멘텀 보존방법에 비해서 계산 과 정이 좀 더 복잡하지만 파랑 중 부가저항 성분에 대한 물리적 해 석이 용이하며, 다중 물체 문제나 비선형 문제로의 확장이 수월 하다는 장점이 있다. 예를 들어, Faltinsen, et al. (98)은 스트 립 방법을 적용한 운동해석과 이 방법을 결합하여 파랑 중 부가 저항을 계산한 바 있다. 이외에도 Choi, et al. () 이 주파수 영역에서의 그린함수법(green function method)에 기반을 둔 압 Received : 8 November 3 Revised : August Accepted : 3 August 5 Corresponding author : Yonghwan Kim, yhwankim@snu.ac.kr
양경규 서민국 김용환 력직접 적분방법을 사용하여 파랑 중 부가저항 해석을 수행하기 도 하였다. 최근에 Joncquez (9) 및 Kim and Kim (; )은 시 간 영역에서 고차의 랜킨패널법을 이용하여 규칙 파에서 부가저 항을 계산한 바 있다. 두 연구에서는 운동량 보존방법과 압력직 접 적분방법을 모두 사용하였으며, 경계 조건 문제를 선형화하기 위한 두 가지 대표적인 방법인 Neumann-Kelvin 및 double-body 선형화 방법을 적용하였다. 특히, Kim and Kim (; )의 연구에서는 규칙 파에서뿐만 아니라 시간영역에 서 불규칙 파로 인한 부가저항을 해석하여 불규칙 파에서의 부가 저항 해석을 위한 기준을 제시한 바 있다. 이에 반해 CFD를 이용한 연구는 아직 초기 단계이며, 주로 선 박의 운동에 초점을 맞추어 진행되어 왔다. Orihara and Miyata (3)는 선수규칙파 중 SR-8 컨테이너선에 대한 운동응답과 부가저항을 WISDAM-X라는 자체 개발 프로그램을 이용하여 해 석한 바 있다. Visonneau, et al. (8)은 ISIS-CFD라는 프로그 램을 이용하여 선수규칙파 중 선박의 상하운동과 부가저항을 계 산하였는데, 비정렬격자(unstructured mesh)를 사용하고 Analytical Weighting Mesh Deformation 기법을 적용하여 선박의 움직임을 구현하였다. Hu and Kashiwagi (7)는 Constrained Interpolation Profile(CIP) 기법을 적용하여 직교격자에서 선박의 운동을 계산하는 프로그램을 개발하였고, 물체의 부피 비를 이용 하여 직교격자 내부에서 물체를 표현하였다. 이러한 여러 기법들로부터 예측되는 부가저항 해석에서 선박 의 선형설계에서 중요하게 다루어져야 하는 것은 단파장에서의 부가저항 예측이다. 이는 실제 상선들이 운항하는 일반적인 해양 환경이 대부분 대파고를 가지는 파랑상태가 아닌, 선박의 길이에 비해 파장이 짧은 영역에 속하기 때문이다. 즉, 단파장에서의 부 가저항 성능이 선박의 연료효율과 관련하여 중요성이 높은 설계 인자이다. 특히 최근 선박설계에서 반드시 고려가 되어야 하는 EEDI(에너지효율설계지수)에서 부가저항은 매우 중요한 역할을 하는 f w 의 결정적 인자이기 때문에, 최근 들어 단파장에서의 부 가저항 문제는 점차 중요성을 더해가고 있다. 본 연구에서는 최근 중요성이 높아지는 단파장에서의 부가저 항을 예측하는 기법에 대하여 체계적으로 살펴보고자 한다. 특 히, 대표적인 3차원 수치기법인 랜킨패널법과 CFD 기반의 두 수 치기법을 적용하여 단파장에서의 부가저항을 계산하고, 이를 기 존의 연구에서 제시된 점근적(asymptotic) 방법과 경험적 식들, 그리고 실험값들과 비교함으로써, 현재 사용되고 있는 수치기법 들의 정확성과 활용성을 관찰하고자 한다. 본 연구에서 적용된 수치기법으로서는, 본 연구에 앞서 개발 된 랜킨패널법 (Seo, et al., 3) 과 직교격자에 기반을 둔 CFD 코드 (Yang, et al., )를 사용하였다. 랜킨패널법에서 는 비스플라인(B-spline) 함수를 도입한 고차 패널법을 적용하였 으며, 경계면을 이산화 하여 3차원 랜킨소스를 선박의 표면과 자 유표면에 분포하였다. 파랑 중 부가저항을 계산하기 위해서 물체 표면에서의 차 동유체 압력을 구하여 적분하는 압력직접 적분 방법을 적용하였다. 직교격자법에서는 유한체적법에 근간하여 이산화를 하였고, 자유표면과 물체를 구분하기 위하여 주어진 셀에서 특정 물질이 차지하는 부피 비를 나타내는 밀도 함수를 도입하였다. 정수 중 과 파랑 중 전진하는 선박에 대한 계산을 통해 선체에 작용하는 힘을 구하고, 각 경우의 전후방향 힘에 대한 시간 평균값 차이를 통해 파랑 중 부가저항을 계산하였다. 이 때, 상하동요(heave)와 종동요(pitch) 방향 운동만 고려되었으며, 선체에 작용하는 힘은 압력을 선체표면에서 적분하여 계산하였다. 기존 연구가 (Seo, et al., 3; Yang, et al., 3) 각 수치 기법에 대한 개발과 이를 이용한 넓은 영역의 파장에 대한 부가 저항 해석에 초점을 맞추었다면, 본 연구에서는 단파장 영역에서 부가저항 해석에 초점을 두었다. 특히, 여러 선형에 따른 개발된 수치해석기법과 기존의 이론식 및 경험식이 갖는 정확도와 한계 등을 고찰하고자 한다.. 해석기법. 랜킨패널법 입사파가 존재하는 무한수심의 자유표면상을 일정한 전진속 도 U로 진행하고 있는 선박에 대해서 Fig. 과 같이 정의되는 물체 고정 좌표계를 생각할 수 있다. 선박으로 입사하는 파는 규 칙파이며, A, ω는 입사파의 진폭과 주파수를 나타내며, β는 입 사각을 의미한다. 선박의 운동은 평균위치에서 정의되며, 병진 운동 ξ T = (ξ, ξ, ξ 3 )과 회전운동 ξ R = (ξ, ξ 5, ξ )의 자유도 운동을 한다. Fig. Coordinate system for 선박운동을 해석하기 위해서 일반적으로 많이 사용되는 속도 포텐셜을 도입하였으며, 이상유체의 비회전 유동을 가정하면 라 플라스 방정식과 다음의 경계조건을 만족하는 속도포텐셜, ϕ를 도입할 수 있다. in fluid domain () JSNAK, Vol. 5, No., August 5 339
단파장 영역에서의 부가저항 해석 () (3) () 는 임의의 지점에서의 선박의 변위를 나타내며, ζ는 파고(wave elevation)를 의미한다. 은 선박의 표 면에서의 법선벡터를 의미하며, g는 중력가속도를 나타낸다. 식 (3)은 운동학적 자유표면 경계조건(kinematic free-surface boundary condition) 이며, 식 ()는 동역학적 자유표면 경계조 건 (dynamic free-surface boundary condition)이다. 위의 경계 조건 문제를 선형화하기 위해서 속도포텐셜과 파고 를 다음과 같이 정의하였다. (5) () 여기서 Φ는 기본 포텐셜을 나타내며, 아래첨자 I와 d는 각각 입 사파 및 교란파와 관련된 성분을 의미한다. Neumann-Kelvin 선형 화 방법을 적용하는 경우 기본 포텐셜 Φ는 균일류 포텐셜 (-Ux)이 되며, double body 선형화 방법을 적용하는 경우에는 z=을 대칭 으로 이미지 물체를 두어 무한 유체장에서 을 만 족하는 경계조건 문제를 풀어 Φ을 계산한다. 위의 문제에 대한 선형화된 경계치 문제는 다음과 같이 나타 낼 수 있다. (8) (7) 여기서 는 평균위치에서의 물체 표면을 의미하며, m는 j m- 항이라고 불리며, 속도에 의해서 발생하는 정상해와 선박의 운동 에 의한 비정상해의 연성효과를 포함하고 있다. 선박이 강체인 경우에 운동방정식은 아래와 같이 표현될 수 있다. () 여기서 M jk 는 질량을 나타내고, F F.K., F Res. 는 Froude-Krylov 힘 및 복원력을 의미한다. F H.D. 는 Froude-Krylov 힘과 복원력을 제외한 선박의 운동과 파의 산란으로 발생하는 유체동역학적 힘 이다. 본 연구에서 적용된 랜킨패널법에서는 파랑 중 부가저항을 해석 하기 위해서 다음과 같은 수식을 도입하여 압력을 직접 적분한다. () () (9) (3) 은 선박의 법선 벡터이며,, 는 각각 차, 차 법선 벡 터를 의미한다. WL은 선박의 수선면을 나타낸다. 위 식은 다른 비선형 하중은 고려하지 않고 선형해로부터 구해지는 차 유체 3 대한조선학회논문집 제 5 권 제 호 5년 8월
양경규 서민국 김용환 력 만을 정리한 것으로서, 부가저항은 이 식으로 구해지는 유체 력의 평균값으로 정의된다. 운동 및 부가저항의 수치해석에 대한 자세한 내용은 Seo, et al. (3)에서 확인할 수 있다.. 직교격자법 비정상, 비점성, 비압축성 유체에 대한 연속방정식과 Euler 방 정식은 다음과 같다. () (5) 여기서 Ω는 전체 영역을 의미하고 Γ는 영역의 경계를 나타낸 다. 은 경계면에서 영역 바깥쪽을 향하는 법선벡터이고, 는 외력항을 의미한다. 속도, 와 압력, p의 연성은 비정상 유동계 산에 효율적인 다단계 방법(fractional step method)을 사용하였 다. 본 수치기법에서는 대류 항에 대한 시간 차분에 대해서는 차 의 명시적 오일러(explicit Euler) 기법을, 압력-포아송 방정식에 대 한 수치계산에서는 Bi-conjugate gradient stabilized(bi-cgstab) 기법을 이용하였다. 또한, 대류 항의 공간에 대한 이산화는 각 셀 경계에서의 플럭스를 통해 이루어지며, monotonized central 제 한자(limiter)를 이용하여 각 셀 경계에서의 속도를 계산하였다. 그리고 3차원 효과를 고려하기 위하여 directional splitting 방법 을 사용하였고, 압력 항에 대해서는 통상적으로 사용되는 중심차 분법(central difference scheme)을 사용하여 이산화 하였다. 직교 격자계(cartesian grid)를 사용하여 선박의 운동과 파랑 중 부가저항을 계산하기 위하여 Fig. 와 같이 각 격자에서 서로 다른 물질이 차지하는 부피 비를 나타내는 밀도 함수(ϕ m, m=: 물, m=: 공기, m=3: 물체)를 정의하였다. ϕ 대한 이송 방정식 은 VOF 계열 기법 중 하나인 tangent of hyperbola for interface capturing(thinc, Xiao, et al., 5) 기법과 다차원 효과를 고 려하는 weighed line interface calculation(wlic, Yokoi, 7) 기법을 사용하여 계산하였다. 각 셀에서 물체에 대한 밀도 함수 값을 알고 있는 경우 다음 식을 이용하여 물체경계 조건을 만족 시키게 된다 (Hu & Kashiwagi, 7). () 여기서 는 속도가 정의된 검사체적 중심에서 물체의 속 도를 나타내고, 는 각 속도 검사체적에 대한 물체의 부피 비 를 의미한다. 위 식을 사용하는 경우 물체경계 조건을 정확한 물 체표면에서 만족하지 않지만, 관성력이 지배적인 문제에는 정확 한 결과를 제공한다. 선체에 작용하는 힘과 모멘트는 다음 식과 같이 면적분을 이용 하여 계산되며, 본 연구에서는 압력에 의한 성분만 고려하였다. (7) (8) 여기서 'nface'는 물체 표면을 표현하는 삼각형 요소 전체 개 수를 의미하고, p, ΔS, n은 각각 직교 격자로부터 구해진 선체표 면에서의 압력과 선체표면 요소의 면적 및 법선벡터를 의미한다. 아래첨자 l은 각 선체표면 요소의 색인(index)을 의미하고 x c 와 x cg 는 각각 물체를 표현한 삼각형 요소의 중심 및 물체의 무게중 심 좌표를 뜻한다. 수치기법에 대한 보다 자세한 내용은 Yang, et al. ()에서 찾을 수 있다. Fig. Coordinate system for Cartesian grid method.3 이론식 및 경험식 단파장 영역에서 부가저항을 계산할 때, 선수부에서 산란되는 파의 비선형성을 고려해야 하기 때문에 수치적으로 정확한 계산 에 어려움이 있다. 이를 보완하기 위하여 단순화된 가정에 근거 한 이론식과 실험에 근거한 경험식이 제안된 바 있다. Faltinsen, et al. (98)은 파랑 중 부가저항에 대한 다음의 이론적 점근식(asymptotic equation)을 유도하였다. 파장이 짧은 파가 입사하는 경우, 파에 노출된 영역에서는 완전반사가 발생한 다고 가정하였다. 또한, 국부 정상 유동(local steady flow)에 대 한 속도를 단순히 물체 표면에 대한 선박 전진속도의 접선 성분 으로 가정하였고, 물체를 정수면에서 단면을 갖는 실린더 형상으 로 생각하였다. JSNAK, Vol. 5, No., August 5 3
단파장 영역에서의 부가저항 해석 coscos (9) sin sin cos cos sin 식 (9)에서 정의되는 적분영역은 Fig. 3에서 파에 직접적으 로 노출된 영역, A-F-B 구간이다. 이 경험식의 보완점이라고 할 수 있는 부분은 반사계수(α d )에 대한 식에서 전진속도를 고려한 파수(wave number)를 사용한 점, 전진속도의 영향을 고려한 보정계수(+α U )에서 실험값을 반 영한식을 사용한 점이며, 이를 통해 날씬한 형상을 가진 선박에 서도 좋은 결과를 기대할 수 있다. 3. 수치해석 3. 계산 모델 및 격자 본 연구에서 사용한 선박 모델은 Series (C B =.7,.8), S75 컨테이너선 및 KVLCC 선형이다. 랜킨패널법의 경우는 네 가지 경우에 대해 모두 계산을 수행하였고, 직교격자법의 경 우 계산 시간의 제약으로 S75컨테이너선과 KVLCC선형에 대 해서만 계산이 수행되었다. 입사파는 모든 계산에 대해 선수파(β =8 ) 경우만 고려하였다. 각 선형에 대한 제원은 Table 에 정 리되어 있다. Fig. 3 Coordinate system for the added resistance calculation in short waves (Faltinsen, et al., 98) 단파장에서의 부가저항을 계산하는 다른 방법으로, 고정된 수 직의 원기둥에 가해지는 표류력에 대한 식을 선박에 적용할 수 있도록 몇 가지 가정을 통하여 보정계수(α d, +α U )를 적용한 경 험식을 Fujji and Takahashi (975)가 제시한 바 있다. 최근에 일본의 National Maritime Research Institute(NMRI) 에서는 Fujji and Takahashi (975)의 식을 기본으로 하되, 여러 실험값을 사 용하여 보정계수의 보완한 식을 다음과 같이 제시하였다 (Tsujimoto, et al., 8; Kuroda, et al., 8). sin sin sin sin Table Coefficients of empirical formula Fujii and Takahashi (975) NMRI Fn max[., -3B (β)+8] f ω e d/g () Table Main dimension of test models Model L(m) B(m) D(m) C B..9 5.7.7 Series. 5.39.5.8 S75 75. 5. 9.5.5 KVLCC 3. 58..8.898 Fig. 는 S75 컨테이너선과 KVLCC 선형에 대해 직교격자 법에서 사용한 삼각형 표면 격자와 이를 통해 계산된 물체 표면 (ϕ 3 =.5)을 나타낸다. 직교격자 내부에서 선박의 형태가 비교적 정확하게 표현되는 것을 확인할 수 있으며, 정량적인 비교는 Yang, et al. ()에서 찾을 수 있다. Fig. Triangle surface meshes and iso-surface (ϕ 3 =.5) of calculated body volume fraction for S75 containership (upper) and KVLCC (lower) Fig. 5는 선박의 운동 계산에 사용되는 패널과 직교격자를 나 타낸다. 랜킨패널법의 경우 원형 격자를 사용하였으며, 선박으로 부터 다섯 파장의 계산 영역을 가진다. 이 때, 방사조건을 만족 하기 위하여 끝 쪽 두 파장에 대하여 감쇠영역을 사용하였다. 직 교격자법의 경우 역시 끝 쪽 두 파장에 감쇠영역을 사용하였고, 3 대한조선학회논문집 제 5 권 제 호 5년 8월
양경규 서민국 김용환 입사파 영역은 파장 길이의 두 배로 설정하였다. 이 때, 입사파 의 파장 (λ/l=.3~.)에 따라 계산 영역이 달라지며, 각 파장에 따라 서로 다른 패널과 직교격자가 사용되었다. 이에 대한 자세한 내용은 Yang, et al. (3)의 기존 연구에 잘 설명이 되어 있다..5.5 Exp. (H/=/, Fonseca, ) Exp. (H/=/, Fonseca, ) Cartesian grid method (H/=/) 3 /A.5 (a) Panel model for.5.5.5.5 /L (a) Heave Exp. (H/=/, Fonseca, ) Exp. (H/=/, Fonseca, ) Cartesian grid method (H/=/).5 5 /ka.5 (b) Grid distribution for Cartesian grid method Fig. 5 Examples of solution grid for two different numerical methods 3. 운동해석 Fig. 와 Fig. 7은 S75 컨테이너선과 KVLCC 선형에 대하 여 선수파 중 상하동요와 종동요의 파장에 따른 운동응답을 계산 한 결과이다. 각각의 운동응답 결과는 실험 결과와 (Fonseca & Soares, ; Lee, et al., 3) 비교하였고, 상하동요 및 종 동요를 제외한 나머지 운동은 구속된 상태로 계산하였다. 직교격 자법의 경우 실험 조건 중 파고가 높은 H/λ=/을 적용하여 계 산을 수행하였고, 랜킨패널법의 경우 선형화 가정이 만족하는 파 고 (H/λ</)를 적용하였다. 운동응답의 경우 특히 공진 주파수 부근에서 파고에 의한 영 향이 두드러지게 나타나는 것을 확인 할 수 있다. 상하동요 및 종동요의 경우 큰 파고에서는 운동응답이 작은 파고의 경우에 비 하여 조금 작게 나타나는 경향을 보이고 있는데, 이는 비선형 운 동에서 흔히 관찰되는 현상이다. 쉽게 이해할 수 있는 바와 같이 랜킨패널법의 경우 H/λ=/의 실험 결과와 유사한 경향을 보 이며, 직교격자법의 경우는 H/λ=/에 대한 실험결과와 유사한 것을 확인 할 수 있다..5.5.5 /L (b) Pitch Fig. Vertical motion responses of S75 containership (Fn =.5) 직교격자법의 경우 파장과 파고에 따라 각각 x방향과 z방향 격자가 결정된다. 때문에 파고를 일정하게 하면 장파장 영역에서 격자의 종횡비가 지나치게 증가하여 계산이 불안정하거나, 입사 파가 부정확하게 생성되는 문제가 발생한다. 이러한 이유로 KVLCC의 경우, 실험보다 큰 파고에 대해 파 기울기를 일정하 게 하여 계산이 수행되었다. Fig. 7은 KVLCC 선형에 대한 운동응답 계산결과를 실험과 비 교한 것이다. 파장이 비교적 긴 영역에서 종동요 운동 응답의 계 산 결과는 서로 유사하지만 계산결과들이 실험에서 관찰된 값들 보다 조금 크게 계산되었다. 실험의 경우 전후동요를 고정하지 않 은 반면, 계산에서는 전후동요를 고정하였다. 이러한 차이가 장파 에서 종동요 운동 응답에 영향을 미치는 것으로 일부 계산과 실험 (Sadat-Hosseini, et al., 3) 에서 관찰되었으나, 추후 보다 자 세한 비교연구가 필요한 부분이다. 이를 제외한 나머지 영역에서 는 실험과 두 계산 결과가 잘 일치하는 것으로 판단된다. JSNAK, Vol. 5, No., August 5 33
단파장 영역에서의 부가저항 해석 3 /A 5 /ka.5.5.5 Exp. (H/L=/, Lee et al., 3) Cartesiangridmethod(H/=/).5.5.5.5.5.5 /L (a) Heave Exp. (H/L=/, Lee et al., 3) Cartesiangridmethod(H/=/).5.5.5 /L (b) Pitch Fig. 7 Vertical motion responses of KVLCC (Fn =.) 3.3 파랑 중 부가저항해석 단파장 영역에서 부가저항을 계산하기 위해서는 선수부 주변 의 국부적인 유동을 정확하게 해석해야 한다. 이를 위해 물체 주 변의 패널 혹은 격자수가 입사파와 산란되는 파의 특성을 충분히 표현할 만큼 많아야 하며, 이에 대한 수렴도 테스트를 반드시 수 행해야 한다. Fig. 8은 KVLCC 선형에 대하여 패널법과 직교격자법의 격자 수렴도 테스트를 수행한 결과이다. 랜킨패널법의 경우 원형 격자 를 사용하였기 때문에 선체의 길이방향에 대한 패널수는 패널의 반경방향에 대한 패널수로 나타낼 수 있다. 서로 다른 파장에 대 하여 반경방향 패널수를 증가시킴에 따라 부가저항 값이 수렴하 는 것을 확인 할 수 있다. 즉, 단파장 영역에서 선체에 의해 산란 되는 파 성분을 랜킨패널법을 사용하여 정확하게 해석하기에는 어려움이 있지만, 많은 수의 패널을 사용하여 수렴된 부가저항 값을 얻을 수 있음을 확인하였다. 본 연구에서는 수렴성 테스트 결과를 바탕으로 부가저항 값의 차이가 약 % 내외가 되는 패널 을 기준으로 삼았으며, 입사파 파장에 따라 반경방향으로 패널을 8~개의 분포시켰을 때 계산된 결과를 택하였다. 직교격자법의 경우 기존 연구 (Yang, et al., 3) 결과를 바 탕으로 격자를 구성하였다. 가장 작은 격자 크기 기준으로, B/ y min =3, A/ z min = 이며 길이 방향에 대해서는 L/ x bow =8,, 8의 격자에 대해 결과를 비교하여 Fig. 8에 나타내었다. KVLCC 선형의 경우 파랑 중 부가저항이 선체의 길이방향에 대 한 의존도가 큰 것으로 관찰되었으며, 선미 부근의 격자수 보다 는 선수 부근의 격자수에 민감하게 변하는 것을 확인하였다. R/gA B /L R/gA B /L 5 3 8 5 3 L/x stern = L/x (a) L/x stern =8 L/x bow /L =.3 /L =.5 /L=.5 8 3 (b) Cartesian grid method Fig. 8 Grid convergence test of added resistance for KVLCC (Fn =.) 3 대한조선학회논문집 제 5 권 제 호 5년 8월
양경규 서민국 김용환 Fig. 9는 직교격자법에 대해 선수 부근 격자크기에 따른 선수 주위의 파형을 보여준다. 격자의 조밀도에 따라 선수 부근에서 약간의 파고차이가 있음을 알 수 있고, 이는 결국 선수면에 가해 지는 압력의 차이, 그리고 궁극적으로 파랑 중 부가저항의 차이 를 가져온다. 다만, 직교격자법의 경우 길이 방향 격자를 조밀하 게 하면, Courant-Friedrichs-Lewy(CFL) 조건에 따른 시간 간격 (time step)에 대한 제약이 강해져 계산 시간이 증가하게 된다. 때문에 랜킨패널법에 비해 수렴된 결과를 얻는 데 어려움이 있으 며, 정량적인 값을 실험과 비교하는 데 주의가 요구된다. 대해서도 유사하게 나타난다. 반면 직교격자법의 경우, 선형에 관계없이 랜킨패널법 보다 실험에 근접한 결과를 보여주었다. 다 만, 직교격자법의 경우 격자에 대한 민감도가 크기 때문에 정량 적인 값을 비교하는 데 주의가 요구된다. R/gA B /L 8 Exp. (Fn=.7, Strom-Tejsen et al.,973) Exp. (Fn=., Strom-Tejsen et al.,973) Short wave (Fujii,975) Short wave (Faltinsen et al.,98) Short wave (NMRI) -....8 /L (a) C B =.7, Fn =. Fig. 9 Instantaneous wave contours computed with three different mesh resolutions for KVLCC (Fn =., Cartesian grid method) R/gA B /L 8 Exp. (Fn=.7, Strom-Tejsen et al.,973) Exp. (Fn=.5, Strom-Tejsen et al.,973) Short wave (Fujii,975) Short wave (Faltinsen et al.,98) Short wave (NMRI) 앞서 소개 되었던 네 가지 선형에 대하여 단파장 영역에서 부 가저항을 계산하여 실험 (Storm-Tejsen, et al., 973; Fujii & Takahashi, 975; Nakamura & Naito, 977; Lee, et al., 3), 이론식 및 경험식 (Faltinsen, et al., 98; Fujii & Takahashi, 975; Kuroda, et al., 8)과 비교하여 Fig. 과 에 나타내었다. Series 선형의 경우 두 가지 서로 다른 CB 값을 갖는 선형 을 대상으로 파랑 중 부가저항을 계산하였는데, C B =.7인 경우 NMRI에서 제안된 식이 실험과 가장 유사한 값을 주었고 다른 경 험식과 이론식, 랜킨패널법에 의한 계산 결과들은 실험결과보다 작은 값들이 예측되었다. 반면 C B =.8인 경우에는 랜킨패널법의 결과가 이론식 및 경험식과 유사하게 나타났다. 이는 뭉툭한 선 형인 경우 랜킨패널법도 NMRI의 방법과 비슷한 정도의 단파장 영역에서 부가저항 값을 계산할 수 있다는 것을 의미한다. 이러 한 경향은 Fig. 에 나타낸 S75컨테이너선과 KVLCC 선형에 -....8 /L (b) C B =.8, Fn=.5 Fig. Added resistance in short waves for Series 이러한 차이를 해석하기 위하여 이론식에 적용된 가정을 살펴 볼 필요가 있다. 이론식을 유도하는 데 있어 두 가지 중요한 가 정 은, 국부적인 정상유동 속도에 대한 가정과 물체형상을 정수 면 까지만 고려하는 것이다. 두 번째 가정은 랜킨패널법에서도 적용되는데, 선형(linear) 가정 하에 문제를 해석하기 때문에 정 수면 위쪽의 형상이 고려되지 않고 정수면 까지 존재하는 물체만 고려된다. 반면, 직교격자법이나 실험의 경우 정수면 위쪽의 형 상 변화에 대한 영향도 고려할 수 있다. 즉, 정수면 아래 위의 선 형변화에 따른 비선형 특성이 파랑 중 부가저항에 영향을 미친 것으로 판단된다. JSNAK, Vol. 5, No., August 5 35
단파장 영역에서의 부가저항 해석 R/gA B /L 8 Experiment (Fujii,975) Experiment (Nakamura,977) Cartesiangridmethod(H/=/) Short wave (Fujii,975) Short wave (Faltinsen et al.,98) Short wave (NMRI) -....8 /L (a) S75 containership, Fn=. (a) S75 containership R/gA B /L 8 Exp. (H/L=., Lee et al.,3) Exp. (H/L=.5, Lee et al.,3) Cartesiangridmethod(H/=/) Short wave (Fujii,975) Short wave (Faltinsen et al.,98) Short wave (NMRI) (b) KVLCC Fig. Comparison of wave contour around ship bow region: (upper) and Cartesian grid method (lower) -....8 /L (b) KVLCC, Fn=. Fig. Added resistance in short waves for S75 and KVLCC 랜킨패널법과 직교격자법을 이용하여 얻은 S75 컨테이너선 과 KVLCC에 대한 선수 근처 파형을 비교하여 Fig. 에 나타 내었다. 랜킨패널법의 경우 정수면 기준 선박의 단면 형상과 파 형을 나타내지만, 직교격자법의 경우 파형을 구하기 위해 3차원 밀도 함수값(ϕ, ϕ 3 )을 이용하여 차원 결과로 변환하였다. 때문 에, 구상선수가 존재하는 선수부 쪽 형상 표현에 차이가 발생한 것을 알 수 있다. KVLCC의 경우, 산란되는 파의 형태나 등고선 의 밀집도 등이 두 수치방법에서 유사하게 나타났지만, S75 선 형의 경우 선수 근처 파형에서 차이가 나타나는 것을 확인 할 수 있다. 특히, KVLCC에 비해 S75 선형의 경우 선수부터 중앙평 행부가 나타나기 이전까지 선체 표면에서 파고의 높이가 차이 나 는 것을 알 수 있고, 이러한 차이가 결과적으로 파랑 중 부가저 항의 차이로 나타났다고 판단된다. 본 계산 모델 중 KVLCC 선형의 경우에는 정수면 근처 위아 래의 단면의 차이가 크지 않지만, S75 선형의 경우 위쪽으로 갈수록 단면의 폭이 넓어지는 선형이다. 이러한 물체 및 자유표 면 상에서의 비선형성 효과가 파랑 중 부가저항에 미치는 영향에 대해서는 추후 더욱 체계적인 연구가 필요하다.. 결 론 본 연구에서는 랜킨패널법과 직교격자 기반의 수치기법 및 경 험식과 이론식을 이용하여 전진하는 선박에 대한 단파장에서의 부가저항을 해석하였다. 각 기법을 통해 얻어진 결과를 모형시험 결과와 비교하였고, 이상의 연구를 통하여 다음과 같은 결론을 얻을 수 있다. ㆍ랜킨패널법이나 직교격자법과 같은 수치기법을 이용하여 단파 장에서의 부가저항을 해석함에 있어 격자 수렴도 테스트는 필 수적으로 수행되어야 하며, 랜킨패널법의 경우 반경 방향에 패 널이 8개 이상 사용되는 경우 파랑 중 부가저항 값은 실험과 약 5%내외의 차이를 보였다. 직교격자법은 선수 근처의 길이 방향 격자에 따라 부가저항 값이 크게 차이 났으며, 가장 조밀 3 대한조선학회논문집 제 5 권 제 호 5년 8월
양경규 서민국 김용환 한 격자에서 실험과 약 8% 정도 차이가 나타났다. 다만, 직교 격자법의 경우 격자에 대한 민감도가 크고 계산 시간의 제약이 있기 때문에 정량적인 값을 비교하는 데 주의가 요구된다. ㆍ단파장 영역에서 부가저항은 선수부의 형상, 특히 선수부의 뭉 툭함(bluntness)과 정수면 위쪽의 변화에 민감한 결과를 보였 다. 이론식과 랜킨패널법의 경우 상대적으로 얇은(slender) Series(C B =.7)와 S75 선형에 대해서는 실험보다 약 ~% 작은 단파장에서의 부가저항을 예측하였고, 뭉툭한 (blunt) 선형인 Series(C B =.8)과 KVLCC에 대해서는 실험 과 유사한 결과를 얻을 수 있었다. 반면, 경험식의 경우 실험에 의해 전진속도에 대한 보정계수를 수정한 NMRI의 식이 모든 선 형에 대해 상대적으로 실험 결과에 근사한 예측을 하였다. ㆍ선수 부근에서 정수면 위쪽의 형상 변화가 심하여 물체의 비선 형성이 고려되어야 하는 경우, 직교격자법이나 실험을 통한 경 험식이 단파장 영역에서 부가저항을 해석하는 데 사용되어야 한다. 반면, 선수부가 상대적으로 비대하거나 정수면 위쪽의 형상 변화가 심하지 않아 선형(linear) 가정에 부합하는 선박에 대해서는 랜킨패널법이나 이론식만으로도 단파장 영역에서 부 가저항을 비교적 정확하게 예측 할 수 있다. 후 기 본 연구는 한국해양과학기술원 주요사업인 해양구조물의 이송 설치를 위한 작업성능 평가기술 개발 및 지식경제부의 에너지 절감을 위한 선박 저항감소 및 추진성능 향상 핵심기술 개발 (과 제번호: 3), 그리고 로이드기금 연구센터 (LRF*-Funded Research Center)에서 지원하는 연구과제의 일환으로 진행되었습 니다. 모든 지원기관들에게 감사를 드립니다. (*LRF: Lloyd s Register Foundation) 또한 서울대학교 공학연구소와 해양시스템 공학연구소(RIMSE)의 행정적 지원에도 감사를 표합니다. References Choi, Y.R. Hong, S.Y. & Choi, H.S.,. An Analysis of Second-Order Wave Forces on Floating Bodies by using a Higher-Order Boundary Element Method. Ocean Engineering, 8, pp.7-38. Faltinsen, O.M. Minsaas, K.J. Liapis, N. & Skjørdal, S.O., 98. Prediction of Resistance and propulsion of a ship in a seaway. Proceeding of 3th Symposium on Naval Hydrodynamics, Tokyo, Japan, - October 98. Fonseca, N. & Soares, C.G.,. Experimental Investigation of the Nonlinear Effects on the Vertical Motions and Loads of a Containership in Regular Waves. Journal of Ship Research, 8(), pp.8-7. Fujii, H. & Takahashi, T., 975. Experimental study on the resistance increase of a ship in regular oblique waves. Proceeding of the th International Towing Tank Conference, Ottawa, September 975. Hu, C. & Kashiwagi, M., 7. Numerical and experimental studies on three-dimensional water on deck with a modified wigley model. 9th International Conference on Numerical Ship Hydrodynamics, Ann Arbor, Michigan, 5-8 August 7. Joncquez, S.A.G., 9. Second-Order Forces and Moments Acting on Ships in Waves. PhD Thesis, Technical University of Denmark. Kashiwagi, M. Takehiro, I. & Takuma, S., 9. Effect of forward speed of a ship on added resistance in waves. Proceeding of 9th International Offshore and Polar Engineering Conference, Osaka, Japan, - June 9. Kim, K.H. & Kim, Y.,. Numerical Analysis of Added Resistance on Ships by a Time-Domain Rankine Panel Method. Journal of the Society of Naval Architects of Korea, 7(3), pp.398-9. Kim, K.H. & Kim, Y.,. Numerical Study on Added Resistance of Ships by using a Time-Domain Rankine Panel Method. Ocean Engineering, 38, pp.357-37. Kuroda, M. Tsujimoto, M. Fujiwara, T. Ohmatsu, S. & Takagi, K., 8. Investigation on Components of Added Resistance in Short Waves. Journal of the Japan Society of Naval Architects and Ocean Engineers, 8, pp.7-7. Lee, J.H. Seo, M.G. Park, D.M. Yang, K.K. Kim, K.H. & Kim, Y., 3. Study on the effects of hull form on added resistance. Proceeding of the th international symposium on practical design of ships and other floating structures, Changwon, Korea, -5 October 3. Maruo, H., 9. The Drift of a Body Floating on Waves. Journal of Ship Research, (3), pp.-. Nakamura, S. & Naito, S., 977. Propulsive Performance of a Containership in Waves. Journal of the Society of Naval Architects of Japan, 5, pp.-8. Orihara, H. & Miyata, H., 3. Evaluation of Added Resistance in Regular Incident Waves by Computational Fluid Dynamics Motion Simulation using an Overlapping Grid System. Journal of Marine Science and Technology, 8(), pp.7-. Sadat-Hosseini, H. Wu, P.C. Carrica, P.M. Kim, H. Toda, Y. & Stern, F., 3. CFD Verification and JSNAK, Vol. 5, No., August 5 37
단파장 영역에서의 부가저항 해석 Validation of Added Resistance and Motions of KVLCC with Fixed and Free Surge in Short and Long Head Waves. Ocean Engineering, 59, pp.-73. Seo, M.G. Kim, K.H. Park, D.M. & Kim, Y., 3. Comparative Study on Added Resistance for Different Hull Forms by using Weakly-Nonlinear Seakeeping Formulations. Journal of the Society of Naval Architects of Korea, 5(), pp.9-58. Storm-Tejsen, J. Yeh, H.Y.H. & Moran, D.D., 973. Added Resistance in Waves. Transactions The Society of Naval Architects and Marine Engineers, 8, pp.5-79. Tsujimoto, M. Shibata, K. Kuroda, M. & Takagi, K., 8. A Practical Correction Method for Added Resistance in Waves. Journal of the Japan Society of Naval Architects and Ocean Engineers, 8, pp.-. Visonneau, M. Queutey, P. Leroyer, A. Deng, G.B. & Guilmineau, E., 8. Ship motions in moderate and steep waves with an interface capturing method. 8th International Conference on Hydrodynamics, Nantes, France, 3 September-3 October 8. Xiao, F. Honma, Y. & Kono, T., 5. A Simple Algebraic Interface Capturing Scheme Using Hyperbolic Tangent Function. International Journal for Numerical Methods in Fluids, 8(9), pp.3-. Yang, K.K. Nam, B.W. Lee, J.H. & Kim, Y.,. Analysis of Large-Amplitude Ship Motions using a Cartesian-Grid-based Computational Method. Journal of the Society of Naval Architects of Korea, 9(), pp.-8. Yang, K.K. Lee, J.H. Nam, B.W. & Kim, Y., 3. Analysis of Added Resistance using a Cartesian-Grid-based Computational Method. Journal of the Society of Naval Architects of Korea, 5(), pp.79-87. Yokoi, K., 7. Efficient Implementation of THINC Scheme: A Simple and Practical Smoothed VOF Algorithm. Journal of Computational Physics, (), pp.985-. 양 경 규 서 민 국 김 용 환 38 대한조선학회논문집 제 5 권 제 호 5년 8월