19.23 수중 유체역학과 부유체 역학

수중 유체역학(underwater hydrodynamics)은 물을 매질로 하는 유동을 다루는 학문 분야이며, 부유체 역학(marine hydrodynamics)은 수면 근처 또는 수면 아래에서 운동하는 물체의 동역학과 주변 유체의 상호작용을 포괄한다. 공기를 매질로 하는 유동과 비교했을 때 물은 밀도가 약 800배 크고, 동점성 계수(kinematic viscosity)가 작으며, 자유 수면(free surface)의 존재로 인해 중력 기반의 파동 현상이 지배적이라는 점에서 해석상의 특징이 뚜렷이 구별된다. 본 절에서는 수중 로봇과 부유체 로봇 플랫폼의 해석에 요구되는 핵심 원리를 체계적으로 정리한다.

1. 해수와 담수의 물리적 특성

해수의 밀도는 표준 조건에서 약 $1025\,\mathrm{kg/m^3}$ 이며, 담수의 밀도는 약 $998\,\mathrm{kg/m^3}$ 이다. 물의 동점성 계수는 $20\,^\circ\mathrm{C}$ 에서 약 $1.0 \times 10^{-6}\,\mathrm{m^2/s}$ 로, 공기의 약 15분의 1 수준이다. 물의 압축성은 매우 작아 일반적인 수중 유동 해석에서는 비압축성 가정이 성립한다. 해수의 상태량은 염분(salinity), 온도, 깊이에 따라 변화하며, 밀도 성층(density stratification)을 형성하여 해양 내부파(internal wave) 현상과 관련된다.

2. 정수압과 깊이에 따른 압력 변화

정지 유체 내부의 압력은 깊이 $h$ 에 대해

$p(h) = p_0 + \rho g h$

의 선형 관계를 따른다. 여기서 $p_0$ 은 자유 수면에서의 기준 압력이며, $g$ 는 중력 가속도이다. 수심이 증가함에 따라 압력이 급격히 상승하므로, 수중 로봇의 내압 설계에서는 목표 작동 수심에 해당하는 정수압을 정확히 산정해야 한다. 예컨대 수심 $100\,\mathrm{m}$ 에서 해수는 약 $1\,\mathrm{MPa}$ 의 정수압을 생성한다.

19.23.3 부력과 자세 안정성

아르키메데스 원리에 따라 물체에 작용하는 부력 $F_b$ 는 배수 체적 $V$ 에 해당하는 유체의 무게와 같다.

$F_b = \rho g V$

수중 로봇의 총 밀도가 유체 밀도와 정확히 일치하는 상태를 중립 부력(neutral buoyancy)이라 하며, 이 조건에서는 외부 힘이 작용하지 않는 한 정지 상태가 유지된다. 부력 중심(center of buoyancy, CB)과 무게 중심(center of gravity, CG)의 상대적 위치는 자세 안정성을 결정한다. CB가 CG보다 위에 있을 때 물체는 복원 모멘트를 통해 수직 자세로 복귀하며, 반대의 경우 불안정해진다. 부유체의 경우 메타센터(metacenter) 개념이 도입되어, 경사 시 메타센터 높이(metacentric height) $GM$ 이 양이어야 복원성이 확보된다.

3. 수중 물체의 항력 특성

수중에서 운동하는 물체에 작용하는 유체 항력은 일반적으로

$F_D = \tfrac{1}{2} \rho V^2 C_D A$

로 표현된다. 저레이놀즈 수 영역에서는 점성 항력이, 고레이놀즈 수 영역에서는 형상 항력이 지배적이다. 수중 로봇은 공기 중 로봇보다 유체 밀도가 훨씬 크기 때문에 동일 속도 조건에서도 항력이 크게 증가하며, 이는 저속·고정밀 운동 제어 전략의 선택을 강제한다. 또한 수중 로봇은 전진, 횡방향, 수직 이동 각각에 대해 서로 다른 항력 계수를 가지며, 비대칭 형상의 경우 6자유도 항력 행렬을 해석 대상에 포함해야 한다.

19.23.5 부가 질량과 유체 관성 효과

수중에서 가속 운동을 수행하는 물체는 주변 유체를 함께 가속시켜야 하며, 이로 인해 추가적인 관성 효과가 발생한다. 이를 부가 질량(added mass) 또는 가상 질량(virtual mass)이라 한다. 6자유도 강체 운동에 대해 부가 질량은 $6 \times 6$ 행렬 $\mathbf{M}_A$ 로 표현되며, 총 유효 관성은 강체 관성 행렬 $\mathbf{M}_{RB}$ 와 부가 질량 행렬의 합으로 주어진다.

$\mathbf{M} = \mathbf{M}_{RB} + \mathbf{M}_A$

예를 들어 반경 $R$ 인 구가 이상 유체 속에서 병진 가속하는 경우 부가 질량은 $\tfrac{2}{3}\pi \rho R^3$ 에 해당하며, 이는 구의 배수 질량의 절반과 같다. 부가 질량 효과는 수중 로봇의 동역학 모델링에서 공기 중 모델에 비해 훨씬 큰 비중을 차지하며, 모델 기반 제어기의 정확성을 좌우한다.

4. 수중 로봇의 6자유도 운동 방정식

해양 로봇 공학에서는 Fossen이 정립한 표준 형식의 수중 로봇 운동 방정식이 널리 사용된다.

$\mathbf{M}\dot{\boldsymbol{\nu}} + \mathbf{C}(\boldsymbol{\nu})\boldsymbol{\nu} + \mathbf{D}(\boldsymbol{\nu})\boldsymbol{\nu} + \mathbf{g}(\boldsymbol{\eta}) = \boldsymbol{\tau}$

여기서 $\boldsymbol{\nu}$ 는 물체 좌표계에서의 선속도와 각속도를 묶은 6차원 속도 벡터, $\boldsymbol{\eta}$ 는 위치와 자세를 기술하는 관성 좌표계 벡터, $\mathbf{C}(\boldsymbol{\nu})$ 는 코리올리 및 구심력 행렬, $\mathbf{D}(\boldsymbol{\nu})$ 는 수력 감쇠 행렬, $\mathbf{g}(\boldsymbol{\eta})$ 는 중력과 부력에 의한 복원 벡터, $\boldsymbol{\tau}$ 는 추력기 및 제어면이 생성하는 일반화 힘이다. 이 표준 식은 원격 조종 무인 잠수정(ROV)과 자율 무인 잠수정(AUV)의 동역학 해석 및 제어기 설계의 기반을 이룬다.

19.23.7 자유 수면과 표면파 현상

부유체가 자유 수면 근처에서 운동하는 경우, 중력 복원력에 의해 수면파(surface wave)가 생성된다. 파동의 분산 관계(dispersion relation)는 심수(deep water) 조건에서

$\omega^2 = g k$

로 주어지며, 천해(shallow water) 조건에서는 $\omega^2 = g k^2 h$ 로 근사된다. 여기서 $\omega$ 는 각주파수, $k$ 는 파수, $h$ 는 수심이다. 부유체에 작용하는 파력(wave force)은 Froude-Krylov 힘, 회절력(diffraction force), 방사력(radiation force)의 합으로 분해되며, 이는 Cummins 방정식이나 주파수 영역 해석을 통해 모델링된다.

5. 부유체의 성능과 무차원 수

부유체의 운동 특성은 Froude 수 $Fr = V / \sqrt{g L}$ , 레이놀즈 수 $Re = V L / \nu$ , 그리고 Strouhal 수 등의 무차원 수로 특징지어진다. 여기서 $L$ 은 대표 길이이다. Froude 수는 중력에 의한 관성력과 자유 수면 파동 형성의 관계를 나타내며, 부유체의 파동 저항(wave-making resistance)이 지배적 항목이 되는 조건을 결정한다. 수중 완전 잠수체의 경우 자유 수면 효과가 배제되므로 레이놀즈 수와 점성 저항이 주요 변수로 작용한다.

6. 캐비테이션

물속에서의 국소 압력이 해당 온도의 포화 증기압(saturation vapor pressure) 이하로 감소하면, 물이 국소적으로 기화하여 기포가 형성되는 캐비테이션(cavitation) 현상이 발생한다. 캐비테이션 수(cavitation number)는 다음과 같이 정의된다.

$\sigma = \dfrac{p_\infty - p_v}{\tfrac{1}{2}\rho V^2}$

여기서 $p_\infty$ 는 자유류 압력, $p_v$ 는 포화 증기압이다. $\sigma$ 가 작을수록 캐비테이션이 발생하기 쉬우며, 캐비테이션이 진행되면 프로펠러 블레이드 표면의 침식, 소음 증가, 추력 저하 등 부정적 결과가 초래된다. 수중 로봇의 추진기 설계에서는 캐비테이션 임계 조건을 회피하는 것이 중요한 설계 지침이 된다.

19.23.10 수중 음향학과 감쇠

물은 공기에 비해 음향 전파 특성이 우수하며, 음속은 약 $1500\,\mathrm{m/s}$ 에 이른다. 음속은 수온, 염분, 압력에 의존하며, UNESCO 표준 식이나 Mackenzie 식과 같은 실험식으로 산정된다. 수중 통신과 소나(sonar)는 전자기파가 물속에서 급격히 감쇠되는 문제를 극복하기 위해 음향파를 주요 전파 수단으로 활용한다. 감쇠 계수는 주파수가 증가할수록 크며, 이는 장거리 수중 통신이 주로 저주파 대역을 사용하는 이유이다.

19.23.11 점성 경계층과 난류

수중 물체 표면의 경계층은 공기 중과 유사하게 층류와 난류로 구분되며, 레이놀즈 수가 증가함에 따라 천이를 겪는다. 물의 점성이 공기보다 크나 밀도 또한 크므로, 동점성 계수가 작은 특성이 지배적으로 작용하여 동일 속도와 길이 조건에서 물속의 레이놀즈 수가 공기 중보다 약 15배 크게 나타난다. 이로 인해 수중 물체는 비교적 낮은 속도에서도 난류 경계층 영역에 진입한다.

19.23.12 로봇 공학에서의 응용

수중 유체역학과 부유체 역학은 해양 로봇 전반의 설계, 해석, 제어에 직접적으로 적용된다. 자율 무인 잠수정의 경우 선체 형상 최적화, 추진기 배치, 항력·부가 질량 모델링, 부력 제어 시스템 설계가 핵심 과제이다. 수상 로봇(USV)과 양서 로봇에서는 자유 수면 효과, 파력 응답, 복원 모멘트가 자세 안정성과 직접적으로 연관된다. 또한 수중 매니퓰레이션 로봇은 말단 부하 운동에 의한 부가 질량 변화와 유체 감쇠를 고려한 제어기를 요구한다. 본 절에서 정립한 이론은 이러한 해양 로봇 응용의 해석적 기반을 제공한다.

출처

Fossen, T. I., Handbook of Marine Craft Hydrodynamics and Motion Control, 2nd ed., Wiley, 2021.
Newman, J. N., Marine Hydrodynamics, 40th anniversary ed., MIT Press, 2018.
Faltinsen, O. M., Sea Loads on Ships and Offshore Structures, Cambridge University Press, 1990.
Brennen, C. E., Cavitation and Bubble Dynamics, Cambridge University Press, 2013.
Triantafyllou, M. S., and Hover, F. S., Maneuvering and Control of Marine Vehicles, MIT, 2003.
Antonelli, G., Underwater Robots, 4th ed., Springer, 2018.

버전

v1.0