예인수조

유체역학 전산유체역학 시뮬레이션 마지막 수정: 2026-08-16 04:51:07

1. 개요[편집]

예인수조 Towing Tank
최초 시설	윌리엄 프루드, 토키 (1872)
규모	길이 수십~수백 m, 예인 전차 + 조파기 + 소파빔
상사 전략	프루드 수를 맞추고 마찰은 계산으로 보정
마찰선	ITTC-1957: C_F = 0.075/(log₁₀Re − 2)²
형상계수	(1+k), 프로하스카 저속 외삽
확장법	ITTC-1978 성능예측법

배 한 척 짓기 전에 미니어처를 물에 띄워 끌어 보는 곳. 150년째 같은 짓을 하는데 아직 아무도 대체하지 못했다.

예인수조(towing tank)는 가늘고 긴 수조 위를 달리는 예인 전차에 축소 모형을 매달아 일정 속도로 끌면서 저항·자세·유동을 계측하는 시설이다. 폭이 좁고 길이가 압도적으로 긴 특이한 형상은 “직선 등속 구간을 최대한 길게 확보한다”는 단 하나의 목적에서 나왔다.

수조가 오늘날까지 살아남은 이유는 계측이 정확해서만이 아니다. 모형에서 잰 값을 실선으로 옮기는 절차가 국제 표준으로 합의되어 있고, 그 절차가 규제와 계약의 근거로 쓰이기 때문이다. 수조는 물리 실험실이면서 동시에 인증 기관에 가깝다.

2. 물리적 구성[편집]

수조 본체. 길이가 성능을 결정한다. 가속·감속 구간을 빼고 정상상태 계측 구간을 확보하려면 수십 m로는 부족해서, 주요 시설은 100~300 m대이고 큰 곳은 500 m를 넘는다. 폭과 수심은 뒤에 나올 봉쇄효과를 억제할 만큼 넉넉해야 한다.
예인 전차. 수조 양옆 레일 위를 달리는 구조물로, 계측 장비와 사람이 함께 탄다. 요구 사양의 핵심은 최고 속도가 아니라 속도 안정도다. 저항 계측의 반복 정밀도가 전차 속도 변동에 직결되므로 0.1% 수준의 정속 제어가 상식이다.
조파기와 소파빔. 한쪽 끝의 플랩식·피스톤식 조파기로 규칙파·불규칙파를 만들고, 반대편의 소파 사면이나 다공성 흡수 구조가 반사파를 잡는다. 다방향 파를 만들려면 조파판을 폭 방향으로 잘게 나눈 분절형이 필요하다. 반사가 남으면 다음 런까지 물이 안 가라앉아 대기 시간이 실험 일정을 지배한다.
모형. 축척비 $\lambda = L_s/L_m$ 는 보통 15~40이고 모형 길이는 몇 m급이다. 예전에는 파라핀 왁스를 깎았고 지금은 목재·강화플라스틱·고밀도 폼을 NC 가공한다. 표면은 실선의 조도를 흉내 내지 않고 매끈하게 만드는 것이 원칙이다 — 조도 효과는 나중에 식으로 더한다.

3. 상사의 딜레마 — 이 시설이 존재하는 이유[편집]

모형과 실선의 유동이 완전히 닮으려면 관련된 무차원수가 모두 같아야 한다. 자유표면 유동에서 지배적인 것은 둘이다.

\mathrm{Fr} = \frac{U}{\sqrt{gL}}, \qquad \mathrm{Re} = \frac{UL}{\nu}

프루드 수를 맞추려면 $U_m = U_s/\sqrt{\lambda}$ 로 느리게 끌어야 한다. 그러면 레이놀즈수는 $\lambda^{3/2}$ 배 작아진다. $\lambda=25$ 면 무려 125배다. 둘을 동시에 맞추려면 $\nu_m = \nu_s/\lambda^{3/2}$ 인 유체, 즉 물보다 백 배 넘게 묽은 액체가 필요한데 그런 것은 없다.¹

윌리엄 프루드가 1870년대에 내놓은 타협이 지금까지 쓰인다. 전저항을 두 조각으로 쪼개되, 하나는 $\mathrm{Re}$ 만의 함수, 다른 하나는 $\mathrm{Fr}$ 만의 함수라고 가정한다.

C_T(\mathrm{Re},\mathrm{Fr}) = C_F(\mathrm{Re}) + C_R(\mathrm{Fr})

즉 맞출 수 있는 것(프루드)은 실험으로 맞추고, 못 맞추는 것(레이놀즈)은 계산으로 갈아 끼운다. 이 분해는 엄밀하지 않다 — 점성과 조파는 실제로 상호작용한다 — 지만, 이 가정 하나가 조선공학을 150년간 굴러가게 했다.

4. 마찰선과 형상계수[편집]

$C_F$ 를 계산하려면 평판 마찰 공식이 필요하다. 국제적으로 합의된 것이 ITTC-1957 모형-실선 상관선이다.

C_F = \frac{0.075}{(\log_{10}\mathrm{Re}-2)^{2}}

이름이 “마찰선”이 아니라 상관선(correlation line)인 것이 중요하다. 이 곡선은 순수한 평판 마찰의 최선 근사가 아니라, 당시 관행과 실선 실적을 맞추기 위해 3차원 효과를 조금 얹어 놓은 공학적 규약이다. 순수 평판에 더 가까운 것은 휴즈 선( $C_{F0}=0.066/(\log_{10}\mathrm{Re}-2.03)^2$ )이나 쇤헤르 선이고, 어느 선을 쓰느냐에 따라 뒤에 나오는 $(1+k)$ 값이 통째로 달라진다. 마찰선과 형상계수는 짝으로만 의미가 있다.

프루드의 원래 분해에서 $C_R$ 은 조파 성분과 점성 압력(형상)저항이 뒤섞인 얼룩이었다. 이를 정리하려고 휴즈가 도입한 것이 형상계수 $(1+k)$ 다.

C_T = (1+k)\,C_F + C_W

선체가 평판이 아니라 3차원 물체이기 때문에 생기는 점성저항 증가분을 배율로 흡수하고, 남은 $C_W$ 만 프루드 수의 함수로 취급하는 것이다. $(1+k)$ 는 실무에서 대략 1.1~1.4 범위이며 비대한 선형일수록 크다.

값은 프로하스카 방법으로 뽑는다. 조파저항이 거의 없는 저속 구간( $\mathrm{Fr}\lesssim 0.2$ )에서 조파 성분이 $\mathrm{Fr}^4$ 에 비례한다고 보고

\frac{C_T}{C_F} = (1+k) + c\,\frac{\mathrm{Fr}^{4}}{C_F}

를 직선으로 회귀해 $y$ 절편에서 $(1+k)$ 를 읽는다. 이 플롯은 조선 수조의 상징 같은 그림인데, 실제로 그려 보면 저속에서 저항이 작아 계측 잡음이 크고 층류 잔존 가능성까지 겹쳐 점들이 예쁘게 안 놓인다. $(1+k)$ 가 몇 % 흔들리면 실선 마력 예측이 그만큼 흔들리므로, 여기가 수조 결과 불확실성의 주요 발원지다.

5. 실선 확장 — ITTC-1978[편집]

모형에서 잰 $C_{TM}$ 을 실선 $C_{TS}$ 로 옮기는 표준 절차가 ITTC-1978 성능예측법이다. 구조는 단순하다.

모형에서 $C_{WM} = C_{TM} - (1+k)C_{FM}$ 로 조파 성분을 분리한다.
프루드 상사에 따라 $C_{WS} = C_{WM}$ 로 그대로 옮긴다.
실선 레이놀즈수에서 $C_{FS}$ 를 다시 계산하고, 여기에 보정항들을 더한다.

C_{TS} = (1+k)\,C_{FS} + \Delta C_F + C_A + C_{AAS} + C_{WS}

$\Delta C_F$ — 조도 허용치. 실선 표면은 도장·용접비드·부착생물로 거칠다. 표준 조도 $k_s = 150\ \mu\mathrm{m}$ 를 가정한 경험식으로 더한다. 취항 후 오손이 진행되면 이 항이 커지고, 그것이 “선저 청소로 연비가 개선된다”는 이야기의 정량적 근거다.
$C_A$ — 상관 허용치. 남은 모든 차이를 흡수하는 항으로, 사실상 각 수조가 자기 시설의 과거 실선 실적으로 보정해 온 경험 상수다. 같은 모형을 다른 수조에서 시험하면 결과가 조금씩 다른 이유의 상당 부분이 여기 있다.
$C_{AAS}$ — 수면 위 선체의 공기저항.

자항 시험(모형에 프로펠러를 달고 미는 시험)까지 하면 추력감소계수, 반류계수, 상대회전효율을 얻어 실선 소요 마력과 회전수까지 예측한다. 프로펠러 단독 시험은 별도로 하고, 이 조합이 속도-마력 곡선이라는 최종 산출물을 만든다.

6. 수조가 스스로 만드는 오차[편집]

수조 실험의 어려움은 대부분 “실물이 아니라서” 생긴다.

층류 잔존과 난류 촉진. 모형 레이놀즈수가 $10^6{\sim}10^7$ 대라 선수 부근이 층류로 남을 수 있다. 그러면 마찰이 실제보다 작게 나오고 $C_R$ 이 오염된다. 그래서 선수에서 선체 길이의 5% 정도 되는 위치에 스터드(작은 원기둥 열)·트립와이어·모래띠를 붙여 강제로 난류로 만든다. 문제는 이 촉진장치 자체가 저항을 만든다는 것이라, 촉진장치의 항력을 별도로 추정해 빼 주는 절차가 표준에 포함되어 있다. 오차를 없애려고 넣은 장치가 또 다른 보정 대상이 되는 전형적인 실험 구조.
봉쇄효과. 수조는 무한 수역이 아니다. 모형이 단면적을 차지하면 주변 유속이 빨라지고(고체 봉쇄), 벽과 바닥이 파를 반사하며, 수심이 얕으면 파의 분산 관계 자체가 바뀐다. 대응은 두 갈래다 — 모형 중앙 단면적을 수조 단면적의 1% 이하로 억제하는 설계 지침, 그리고 슈스터·다무라 계열의 봉쇄 보정식으로 유효 속도를 수정하는 사후 보정. 수심 프루드 수 $\mathrm{Fr}_h = U/\sqrt{gh}$ 가 0.7을 넘어가면 천수 효과로 저항 곡선이 통째로 왜곡되므로 시험 속도 상한이 여기서 정해진다.
런 간 대기. 앞 런이 만든 파와 순환류가 가라앉기 전에 다음 런을 시작하면 데이터가 오염된다. 큰 수조에서 하루에 소화할 수 있는 런 수가 생각보다 적은 이유이고, 수조 사용료가 비싼 이유이기도 하다.
계측계. 저항 동력계(로드셀)로 예인력을, LVDT나 광학 추적으로 트림·침하를, 파고계로 파형을 잰다. 자유항주 모형과 PMM(평면운동장치) 시험은 조종성 미계수를 얻기 위한 별도 계열이며, 이쪽은 예인수조가 아니라 넓은 조종·내항 수조에서 한다.

이 모든 항목에 대해 국제예인수조회의(ITTC)가 절차서와 불확실성 해석 지침을 발행한다. 1933년 첫 회의 이래 3년 주기로 모여 마찰선·확장법·시험 절차를 개정해 온 조직이고, 수조 데이터가 시설을 넘어 비교 가능한 것은 이 표준화 덕분이다.

7. CFD와의 관계 — 대체재가 아니라 심급[편집]

자유표면 RANS가 성숙하면서 “수조는 곧 없어진다”는 말이 30년쯤 돌았는데, 실제로 벌어진 일은 역할 분담이었다.

수조는 기준 데이터 공급자다. KCS 컨테이너선, KVLCC2 유조선, JBC 벌크선, DTMB 5415 함정 같은 공개 선형의 저항·자세·국부 유동장·파형 절단 데이터가 국제 CFD 워크숍의 표준 문제로 쓰인다. 검증 및 확인에서 코드 검증(verification)은 계산 안에서 닫히지만, 확인(validation)은 실험값 없이는 정의 자체가 성립하지 않는다. 그 실험값의 출처가 수조다.
CFD가 수조를 보완한다. 축척 효과 — 특히 $(1+k)$ 가 정말 레이놀즈수에 무관한가 — 를 따지려면 모형 스케일과 실선 스케일을 같은 방법으로 계산해 비교해야 하는데, 이건 수조가 원리적으로 못 하는 일이다. 반대로 CFD의 난류 모델링과 격자 의존성은 수조 데이터 없이 못 잡는다. 서로의 사각지대를 메우는 관계.
인증의 최종심급은 여전히 수조다. 선박 에너지효율 규제에서 요구하는 속도-마력 곡선 검증은 표준 절차에 따른 모형시험을 근거로 삼고, 조선소-선주 계약의 속력 보증도 수조 예측과 실선 시운전 결과로 정산된다. 계산으로 대체하려면 수치해석의 불확실성이 수조의 불확실성보다 작다는 것을 먼저 증명해야 하는데, 자유표면 문제에서 아직 그 단계는 아니다.²

대표 시설로는 미 해군의 데이비드 테일러 모형시험수조(카더록), 네덜란드 MARIN, 독일 HSVA, 스웨덴 SSPA, 일본 NMRI, 한국의 선박해양플랜트연구소(KRISO)와 여러 대학 수조가 있다. 신설이 뜸한 대신 기존 수조는 계측계·조파기 현대화와 자율운항선 시험 같은 새 수요로 살아남고 있다.³

8. 관련 문서[편집]

그렇다고 시도가 없었던 건 아니다. 압력을 걸어 캐비테이션 상사까지 맞추는 감압 수조나, 회전 수조·순환 수조처럼 다른 축을 포기하고 특정 상사를 얻는 시설들이 존재한다. 다만 “물보다 훨씬 묽은 액체를 몇백 톤 채운다”는 방향은 여전히 아무도 못 갔다. ↩
실무자들 사이에는 “CFD로 경향을 보고, 수조로 숫자를 받고, 시운전에서 진실을 안다”는 3단 격언이 있다. 셋의 값이 다 다를 때 누구를 믿느냐는 질문에는 대체로 계약서가 답한다. ↩
예인수조 견학에서 가장 인기 있는 순간은 언제나 전차가 출발하는 30초다. 그 뒤로 이어지는 “물결이 가라앉을 때까지 20분 대기”는 아무도 사진을 찍지 않는다. 실험 과학의 시간 배분이란 대개 이런 모양이다. ↩