슬로싱

전산유체역학 유체역학 다상유동 마지막 수정: 2026-07-19 04:12:18

1. 개요[편집]

슬로싱 Sloshing
분야	유체역학 × 자유표면 유동 × 유체-구조 연성
핵심 변수	충전율(filling ratio), 가진 주파수
주력 해석 기법	VOF, SPH, 레벨셋
대표 무대	LNG 운반선 화물창, 로켓 추진제 탱크

컵에 물을 반쯤 담고 걸어봐라. 가득 채웠을 때보다 훨씬 잘 쏟아진다. 그 짜증나는 현상에 산업계는 수십 년째 연구비를 갈아넣고 있다.

슬로싱(sloshing)은 부분적으로 채워진 탱크가 외부에서 흔들릴 때 내부 액체의 자유표면(free surface)이 크게 진동하면서, 벽면에 동적 하중과 국부적 충격압을 가하는 현상이다. 액체가 가득 차 있으면 표면이 없어 슬로싱도 없고, 텅 비어 있어도 없다. 하필 실제 운항 조건에서 가장 흔한 부분 충전 상태에서만 무섭게 살아난다는 점이 이 문제를 성가시게 만든다.¹

슬로싱이 단순한 “출렁임” 이상으로 취급되는 이유는 두 가지다. 첫째, 액체의 자유표면에는 고유진동수가 있고, 선박의 횡동요(roll)나 차량의 조향 주기가 여기에 걸리면 공진이 일어나 진폭이 폭발한다. 둘째, 그렇게 자란 파가 벽이나 천장을 때릴 때 발생하는 충격압(impact pressure)은 지속시간이 밀리초 단위인 대신 크기가 평균 정압의 수십 배까지 튄다. 구조는 평균이 아니라 이 꼬리(tail)에서 깨진다.

2. 고유진동수와 공진[편집]

가장 기본이 되는 결과는 폭 $L$ , 정수(靜水) 수심 $h$ 인 직사각형 탱크의 선형 슬로싱 고유진동수다. 포텐셜 유동 가정 아래 자유표면 경계조건을 풀면 $n$ 차 모드의 각진동수는 다음과 같다.

\omega_n^2 = \frac{n\pi g}{L}\,\tanh\!\left(\frac{n\pi h}{L}\right)

여기서 $g$ 는 중력가속도. 이 식이 말해주는 바는 명확하다. 고유진동수는 탱크의 크기와 채운 깊이만으로 결정되고, 액체의 점성이나 밀도는 (선형 영역에서는) 끼어들지 않는다. 그래서 설계 단계에서 제일 먼저 하는 일이 운항 중 겪을 운동 주기 스펙트럼과 $\omega_1$ 을 겹쳐 그려보는 것이다. 겹치면 그 충전율은 운항 금지 구간이 되거나, 탱크 형상을 바꿔야 한다.

깊은 수심( $h/L \gg 1$ )에서는 $\tanh \to 1$ 이라 $\omega_n^2 \approx n\pi g/L$ 로 수심 무관해지고, 얕은 수심에서는 $\tanh x \approx x$ 근사가 들어가 $\omega_1 \approx (\pi/L)\sqrt{gh}$ 가 된다. 얕은 물 슬로싱이 특히 고약한 건 이때 파가 선형 정상파로 남지 못하고 수력도약(hydraulic jump) 형태의 이동파(traveling bore)로 변해, 벽에 정면으로 달려들기 때문이다. 3차원 탱크에서는 여기에 회전하며 도는 스월(swirl) 모드까지 끼어들어 비선형 난이도가 한 단계 더 올라간다.²

3. 충격압이라는 통계 문제[편집]

슬로싱 해석이 일반적인 전산유체역학 문제와 결정적으로 다른 지점이 여기다. 충격압은 결정론적 값이 아니라 확률변수다. 같은 불규칙 운동 시계열을 조금만 바꿔 다시 돌려도 최대 압력이 두 배씩 왔다 갔다 한다. 이유는 충격 순간의 물리가 지독하게 국소적이기 때문이다.

가스 쿠션(gas cushion): 파가 벽을 때리기 직전 갇힌 기체가 압축되며 완충재 노릇을 한다. 기체가 얼마나 갇히느냐가 압력 피크를 좌우한다.
기포 진동: 갇힌 기포가 압축-팽창을 반복하며 압력 신호에 감쇠 진동 꼬리를 남긴다.
하이드로일래스틱 연성: 벽이 탄성적으로 밀려나면서 압력이 재분배된다. 강체 벽 가정으로 계산한 압력이 실제보다 클 수도 작을 수도 있다. 유체-구조 연성 문제로 넘어가는 지점.

그래서 실무 절차는 “한 케이스 잘 푸는” 것이 아니라 5시간 내외의 불규칙 운동 시뮬레이션(혹은 모형시험)을 수십~수백 회 반복해 압력 피크 표본을 모으고, 극값 통계로 설계 수명 동안의 재현 확률 압력을 외삽하는 쪽으로 굳어졌다. 선급(DNV, ABS, LR 등)의 슬로싱 평가 절차서가 통계 처리 방법을 명시하고 있는 것도 이 때문이다.

4. 축소모형 시험과 상사의 함정[편집]

슬로싱은 아직도 모형시험이 강한 발언권을 갖는 분야다. 문제는 축척을 줄일 때 모든 무차원수를 동시에 맞출 수 없다는 것. 자유표면 유동이니 프루드수를 맞추는 게 국룰인데,

Fr = \frac{U}{\sqrt{gL}}

를 보존하면 속도는 $\sqrt{\lambda}$ 로, 시간은 $\sqrt{\lambda}$ 로 줄어든다. 그런데 여기서 레이놀즈수는 $\lambda^{3/2}$ 배로 어긋나고, 웨버수(표면장력)와 기체-액체 밀도비, 압축성 관련 수는 전부 틀어진다. 무차원수 사이의 이 근본적 충돌 때문에 축소모형 압력을 실선 압력으로 환산하는 스케일링 법칙이 오랫동안 논쟁거리였다. 특히 LNG처럼 액체-기체 밀도비가 중요한 계에서는 상온 물-공기 조합이 실제 조건을 제대로 못 흉내 낸다는 지적이 나와, 무거운 기체나 저온 유체를 쓰는 시험 설비까지 등장했다.³

5. 수치해석 기법[편집]

슬로싱 CFD는 결국 자유표면을 얼마나 잘 붙잡느냐의 싸움이다.

기법	강점	약점
VOF 방법	질량 보존이 본질적으로 정확	계면 번짐, 압축성 기체 처리 부담
레벨셋 방법	매끄러운 계면·곡률	질량 손실
SPH	파 부서짐·비산 자연스러움	압력장 노이즈, 수렴성 확보 어려움
포텐셜/모달 기법	압도적으로 빠름	파 부서짐·충격 표현 불가

상용·오픈소스 진영에서는 유한체적법 기반 VOF가 여전히 주류다. 여기에 요구되는 조건이 만만치 않다. 충격 순간을 잡으려면 시간 간격이 $10^{-5}$ 초 수준까지 내려가야 하고(CFL 조건이 그냥 넘어가 주지 않는다), 압력 피크는 격자 크기에 대해 좀처럼 수렴하지 않기로 악명 높다. “격자를 반으로 줄였더니 압력이 30% 올랐다”가 흔한 일이라, 검증 및 확인 관점에서 슬로싱 충격압은 지금도 어려운 축에 속한다. 기체를 비압축으로 두면 쿠션 효과가 사라져 압력이 과대평가되므로, 최소한 기체상은 압축성으로 두는 것이 요즘 권장 사항이다.

6. 억제와 등가 모델[편집]

슬로싱을 줄이는 고전적 처방은 배플(baffle) 이다. 탱크 내부에 수직·수평 격벽이나 다공판을 넣어 유동을 끊고, 와류 생성으로 에너지를 소산시켜 공진 응답을 낮춘다. 원리는 유동박리에 의한 형상 항력이므로, 배플 근처는 아예 난류 해석 영역이 된다.

반대로 전체 시스템 거동만 필요할 때는 유체를 통째로 풀지 않고 등가 기계 모델로 치환한다. 슬로싱하는 액체를, 탱크에 고정된 강체 질량 + 스프링에 매달린(또는 진자로 매달린) 슬로싱 질량으로 바꾸는 하우스너(Housner) 계열 모델이 대표적이다. 이렇게 하면 저수조 내진 설계나 위성 자세제어 루프에 슬로싱을 몇 개의 상미분방정식으로 끼워 넣을 수 있다. 실제로 우주 발사체와 위성에서 추진제 슬로싱이 자세제어기와 연성되어 불안정을 일으킨 사례가 여러 차례 보고되어, 추진제 탱크의 배플 설계와 슬로싱 모드 식별은 항공우주 쪽의 전통적인 필수 항목이다.⁴

7. 어디서 문제가 되는가[편집]

LNG 운반선 화물창: 슬로싱 연구비의 대부분이 여기서 나온다. 멤브레인형 화물창은 단열 패널이 직접 충격압을 받기 때문에, 부분 충전 운항 허용 여부가 곧 상업적 경쟁력이다.
차량 연료탱크·탱크로리: 급제동·급선회 시 슬로싱 하중이 차량 동역학에 되먹임된다. 탱크로리 전복 사고의 단골 기여 인자.
저수조·수처리조 내진: 지진 시 유체의 대류 성분(convective mode)이 구조에 추가 모멘트를 준다. 모드 해석에 슬로싱 모드를 별도로 포함시킨다.
원자로 및 대형 액체 저장 시설: 자유표면 진동이 내부 구조물과 유체-구조 연성을 일으킨다.

8. 관련 문서[편집]

그래서 실무에서 제일 먼저 던지는 질문이 “충전율 몇 %인가요”다. 보통 위험 구간은 10~~30%(얕은 물, 이동파)와 70~~95%(천장 충격) 양쪽에 몰려 있고, 50% 근처는 의외로 얌전한 경우가 많다. 물론 탱크 형상이 바뀌면 이 감각도 바뀐다. ↩
스월 모드는 정상파가 안정성을 잃고 회전 방향으로 갈라지는 전형적인 분기(bifurcation) 현상이다. 방향이 좌우 어느 쪽으로 도는지는 초기 미세 교란이 정하기 때문에, 같은 조건을 다시 돌리면 반대로 돌기도 한다. 재현성을 요구하는 리뷰어에게 설명하기 가장 곤란한 결과. ↩
스케일링 논쟁이 길어진 근본 원인은 “무엇을 상사시킬 것인가”가 케이스마다 다르기 때문이다. 전체 파 운동은 프루드수가 지배하지만, 국부 충격압은 기체 압축성과 밀도비가 지배한다. 하나의 모형시험으로 둘 다 만족시키는 방법은 아직 없다. ↩
액체 추진제는 발사체 총질량의 압도적 비중을 차지하므로, 그 질량이 탱크 안에서 제멋대로 출렁이면 사실상 거대한 가변 질량-스프링계를 자세제어기에 물려놓은 셈이 된다. 배플은 무게를 늘리는 부품이지만, 이 경우엔 안 달 수가 없다. ↩

슬로싱