국제해사기구인 IMO(International Maritime Organization)는 지구온난화와 같은 환경문제를 개선하고자 해양환경에 영향을 미치는 선박 운항에 대한 규제를 강화해왔다. IMO의 주요 협약 중 하나인 MARPOL(International Convention for the Prevention of Pollution from Ships)은 해양환경의 오염방지를 위한 선박 설계‧건조 단계의 규제를 다루고 있다. 특히 MARPOL Annex VI는 선박에서 배출되는 SOx, NOx 등의 GHG(Green House Gas) 배출량을 감소시킴으로써 대기오염을 방지하도록 규제하고 있다.

또한 MEPC(Marine Environment Protection Committee)에서는 MARPOL 협약을 준수하여 선박의 GHG 배출량 감축을 목표로 하고 있으며, IMO GHG Study는 2008년부터 2018년까지의 선종 별 GHG 배출량 통계 및 규제 적용에 따른 미래의 GHG 배출량 시나리오를 예측하고 있다 (IMO, 2020a). 운항하는 선박에 대한 대표적인 적용 규제로는 EEOI(Energy Efficiency Operational Indicator), SEEMP(Ship Energy Efficiency Management Plan)를 들 수 있으며, 최근에는 EEXI(Energy Efficiency Existing Ship Index)와 CII(Carbon Intensity Indicator)에 대한 적용을 MEPC 75차에서 승인하였고, MEPC 76차에서 채택을 검토 중에 있다 (IMO, 2020b).

위의 규제 및 조약 모두 선박의 GHG 배출량 감소와 에너지효율 증가를 의무적으로 요구하고 있으며, 이를 충족시키기 위한 연구들이 활발하게 수행되고 있다. 대표적인 방안으로는 선박에 직접적으로 적용하는 기술적 방안(technical measures), 대체연료 및 에너지원(alternative fuels or energy sources), 유지관리전략(maintenance strategies), 운항적 방안(operational measure), 그리고 선박의 GHG 배출량을 경제적으로 고려하여 적용하는 시장기반방안(MBM, Market Based Measures)이 있으며, Table 1에 각 방안 별 세부적인 사항들에 대해 정리하였다 (GloMEEP, 2018; Joung et al., 2020).

선박에 직접적으로 적용하는 대부분의 방안들은 비용적, 시간적인 경제성 검토 후 나온 제안들이며, 그 중 트림 최적화(Trim optimization)는 발라스트수의 조절 또는 화물의 배치를 통해 현실적으로 선박에 비교적 쉽게 적용할 수 있는 방법이다. 또한 모든 선박에 적용이 가능하며, 선박의 효율을 높여 연료소비량을 줄임으로써 비용적인 부분을 감소시킬 수 있다 (Kishev et al., 2014).

선박의 저항 및 추진성능 분석에 수치해석의 활용이 점차 증가하고 있다. 더불어, 지속적인 컴퓨팅 파워와 병렬처리 기술의 발달로 수치해석을 이용한 선박의 성능평가에 대한 신뢰도는 모형시험 수준에 이르고 있다. 수치해석은 모형시험을 통해 얻는데 한계가 있는 국부 영역에 대한 추가 정보 획득과 함께 모사대상 영역 전역에 걸친 유동장의 가시화가 가능한 장점이 있다. 그러나 추진성능 평가에 있어 수치해석은 프로펠러를 직접 모사(direct numerical simulation)하거나 동체력(body force)를 활용하여 모사하는 등의 기법이 존재하나, 직접 모사하는 방법은 현업에서 활용하기에 막대한 계산 시간(computing time)이 필요하고, 동체력을 활용하는 방법은 트림 변화와 같이 미소한 계산 조건 변화에 따른 모사 결과의 변별력을 기대하기 어렵다. 이에 따라 모형시험과 수치해석의 장점을 활용하여 속도마력 성능을 추정하는 연구를 수행하였다.

선박의 트림변화는 선체의 저항에 영향을 미치며, 선체의 저항은 선박의 소요동력과 높은 상관관계가 있다 (Park et al., 2014). 따라서 저항성능 분석에 대하여 모형시험 (Iakovatos et al., 2013; Liarokapis et al., 2017)과 수치해석 (Seo et al., 2015; Islam & Soares, 2019; Gao et al., 2019)을 이용한 많은 연구가 수행되었다. 이 연구들은 저항성능 분석을 통해 트림변화에 따른 선체의 저항, 자세, 파형의 변화를 언급하였으며, 더 나아가 수치해석을 이용하여 저항을 성분 별로 분류하여 검토한 바 있다 (Park et al., 2013; Han et al., 2015). 선박의 소요마력을 분석하기 위한 연구도 많이 수행되었는데, 이는 모형시험을 활용하거나 (Sogihara et al., 2018; Alexandrina, 2015), 수치해석의 저항결과와 모형시험의 자항결과를 함께 활용하여 분석을 수행하였다 (Kim et al., 2015).

본 연구에 사용된 대상선은 11,000TEU 컨테이너선으로 실 운항에서의 연비효율 등을 고려하여 설계 속도에 비하여 저속으로 운항 중이다. 대상선을 보유 및 관리하는 선사의 제공 자료를 바탕으로 실제 운항 조건을 반영한 흘수 및 속도를 연구 수행 조건으로 선정하였다.

각 흘수에 대해 초기 트림상태의 배수량을 유지한 상태에서 선수 및 선미트림을 고려하였으며, 실제 운항 조건인 저속 범위에서 연구를 수행하였다. 대상선의 소요마력에 미치는 영향 분석에 있어 저항성능 분석은 수치해석을, 추진성능 분석은 모형시험과 수치해석을 이용하였다.

대상선은 2018년 건조하여 현재까지 운항을 하고 있는 선박으로 주요제원은 Table 2와 같다.

대상선의 Stowage plan, Trim and Stability Booklet 등의 선사에서 제공한 자료를 바탕으로 연구 수행 조건을 구체화하였다. T1과 T2의 두 개 흘수를 선정하였으며, Fig. 1에 각 흘수에서의 Stowage plan을 나타내었다. 대상선의 선수 및 선미에서의 각 흘수의 높이를 Fig. 2에 나타내었으며, T2 흘수의 경우 Even 상태에서 대상선의 선수벌브 상단이 수면위로 돌출됨을 확인할 수 있다.

Fig. 1 
Stowage plan of target vessel

Fig. 2 
Schematic configuration of T1 and T2 at bow and stern

실제운항조건을 반영한 T1과 T2의 두 개의 흘수에 대하여 트림변화에 따른 저항추진성능 변화를 분석하였다. 트림은 FP 흘수에서 AP 흘수를 뺀 값(T_FP-T_AP)으로 정의하였으며, Fig. 3과 같이 선수트림은 양의 값(+), 선미트림은 음의 값(-)으로 정의하였다. 각 흘수에서의 배수량은 트림에 관계없이 Even 상태의 배수량과 동일하게 적용하였다.

Fig. 3 
Definition of trim

또한, 대상선은 실제 설계 속도보다 낮은 속도로 운항하기 때문에, 저속운항 조건을 포함한 범위에서 연구를 수행하였다. Table 3에 본 연구의 모형시험 및 수치해석을 위한 연구 수행 조건을 요약하여 나타내었다. 모형시험 및 수치해석 결과를 활용하여 트림변화에 따른 저항추진 성능을 Even 상태에서의 성능과 비교하여 분석하였다.

본 연구에서는 2차원 ITTC 1978 해석법을 사용하여 모형선 결과를 실선 결과로 확장하였으며, Table 4에 저항 및 추진성능 분석에 사용되는 대표적인 식과 각 식에 사용되는 물리량에 대하여 파란색(모형시험)과 빨간색(수치해석)으로 구분하여 정리하였다. 한정된 자원 등의 현실적인 상황 및 해당 방법의 신뢰도를 고려하여 저항성능 분석에는 수치해석 결과를 이용하였으며, 추진성능 분석에는 수치해석 및 모형시험 결과를 이용하였다.

11,000TEU 컨테이너선을 대상으로 트림변화에 따른 소요마력의 변화를 분석하였다. 선사의 제공 자료를 바탕으로 상대적으로 빈도가 높은 실제운항조건을 반영하였으며, 선수트림을 양의 값으로하여 선수 및 선미트림을 고려하였다. 위의 조건들을 적용하여 대상선의 트림변화에 따른 저항성능 및 추진성능을 분석하였으며, 2차원 ITTC 1978 해석법을 사용하여 실선의 속도-마력성능을 분석하였다. 이 때, 저항성능 분석은 상용 수치해석 코드인 STAR-CCM+ v13.06을 이용한 수치해석 결과를 바탕으로 수행하였다. 그리고 추진성능 분석에는 수치해석 결과와 KRISO 선형시험수조에서 수행한 자항시험 결과를 바탕으로 수행하였다. 다음은 본 연구의 전반적인 내용을 요약하여 나타내었다.

1) 선사의 제공 자료를 바탕으로 실제 운항 조건을 반영한 흘수(T1, T2) 및 속도(Fn=0.11, 0.14, 0.16)를 선정하였다. 트림은 선수트림을 양의 값으로하여 선수 및 선미트림을 고려하였으며, 배수량은 Even 상태와 동일하게 적용하였다.

2) 한정된 자원 등의 현실적인 상황 및 해당 방법의 신뢰도를 고려하여 저항성능 분석에는 수치해석 결과를 이용하였으며, 추진성능 분석에는 수치해석 및 모형시험 결과를 이용하였다. 결과적으로 수치해석을 통해 유효동력을, 모형 자항시험을 활용하여 추진효율을 얻어 소요마력을 계산하였다.

3) 저항성능의 경우, 전반적으로 선수트림이 증가할수록 유효동력이 작아지는 경향을 확인하였다. 단, 설계흘수인 T1 흘수에서는 even 상태가 작은 경향을 보이는 속도가 Fn=0.14, 0.16에서 나타났다. Off-design 조건이었던 T2 흘수에서는 트림변화에 따른 선수벌브 상단의 돌출정도가 선수 파형 분포의 변화를 가져오며, 이는 결과적으로 유효동력에 영향을 미쳤음을 확인하였다.

4) 준추진효율의 경우, T1 흘수에서는 Fn=0.16 & Trim +2m 조건의 준추진효율이 가장 크게 나타났으나, Even 상태와의 변별력이 없었다. 이 외의 속도조건에서는 Even 상태의 준추진효율이 가장 높음을 확인하였다. T2 흘수에서는 Fn=0.16 & 선수트림일 때 준추진효율의 미소한 증가를 보이지만, 이 외의 속도 조건에서는 Even 상태의 준추진효율이 가장 높음을 확인하였다.

5) 전달동력의 경우, T1 흘수에서는 Even 상태에서 전달동력이 가장 작음을 확인하였고, T2 흘수에서는 동일속도에서 선수트림이 증가할수록 전달동력이 작아지는 경향을 확인하였다.

6) T2 흘수는 선수벌브 상단이 돌출될 정도로 낮은 흘수이므로 수면과 프로펠러 사이의 여유가 충분하지 않다. 분석결과에 따르면 선수트림이 증가할수록 전달동력이 작아지는 경향을 보였으나, 과도한 선수트림에서는 프로펠러가 수면으로 돌출되어 프로펠러 성능이 급격히 저하될 우려가 있다.

7) 전반적으로 T1과 T2 두 흘수 모두 유효동력과 전달동력이 유사한 경향을 보였고, 두 동력 간에는 높은 상관관계를 가지고 있음을 확인하였다.