유동박리

1. 개요[편집]

유동박리(flow separation)는 물체 표면을 따라 흐르던 경계층이 표면에서 떨어져 나가면서 벽 근처의 유동이 역류하기 시작하는 현상이다. 박리가 시작되는 지점을 박리점(separation point)이라 부르며, 그 뒤로는 소용돌이가 가득한 무질서한 후류(wake)가 형성된다.

골프공을 봤다면 이미 유동박리를 봤다. 표면의 딤플(dimple)은 바로 이 박리를 지연시켜 공을 더 멀리 날리기 위한 장치다.¹ 비행기가 갑자기 양력을 잃고 떨어지는 실속(stall), 자동차 뒤꽁무니에 생기는 저압 후류로 인한 항력, 굴뚝 뒤에 규칙적으로 떨어지는 카르만 와열 — 이 모든 게 유동박리가 주범이다. 공력해석에서 박리를 잡느냐 못 잡느냐가 곧 해석의 성패를 가른다.

2. 왜 떨어지는가: 역압력구배[편집]

유동박리의 근본 원인은 역압력구배(adverse pressure gradient), 즉 유동 방향으로 갈수록 압력이 높아지는 상황이다.

유선을 따라가면 유체는 압력이 낮아지는 쪽(순압력구배)으로는 신나게 가속하지만, 압력이 높아지는 쪽으로는 언덕을 오르듯 감속해야 한다. 문제는 경계층 안쪽, 벽에 가까운 유체다. 이 녀석들은 이미 벽면 마찰로 운동량을 거의 다 잃어 기진맥진한 상태라, 압력 언덕을 끝까지 오를 힘이 없다. 결국 어느 지점에서 벽 근처 유체가 멈췄다가 거꾸로 흐르기 시작하고, 바로 그 순간 경계층이 벽에서 떨어져 나간다.

수학적으로 박리점은 벽면에서의 속도 기울기, 즉 벽면 전단응력이 0이 되는 지점으로 정의된다.

$$\left. \frac{\partial u}{\partial y} \right|_{y=0} = 0$$

이 값이 음수로 바뀌면(즉 벽 근처 속도가 역방향이 되면) 완전히 박리된 상태다. 순압력구배 구간에서는 박리가 절대 일어나지 않는다는 점이 중요하다 — 물체 형상 설계가 곧 압력구배 설계인 이유다.

3. 층류 박리 vs 난류 박리[편집]

여기서 유체역학의 유명한 역설이 하나 등장한다. 난류 경계층이 층류 경계층보다 박리에 더 강하다.

• 층류 경계층은 얌전하지만 벽 근처 운동량이 빈약해서, 역압력구배를 만나면 비교적 일찍 박리된다.
• 난류 경계층은 무질서한 혼합 덕분에 바깥의 빠른 유체가 벽 근처로 계속 운동량을 실어 나른다. 그래서 압력 언덕을 더 오래 버티고, 박리점이 뒤로 밀린다.

골프공 딤플의 정체가 바로 이것. 딤플이 경계층을 일부러 난류로 만들어 박리를 지연시키고, 후류를 좁혀 압력 항력을 줄인다. 매끈한 공보다 딤플 공이 훨씬 멀리 나가는 이유다. 레이놀즈수가 특정 임계값을 넘으면 항력계수가 갑자기 뚝 떨어지는 “항력 위기(drag crisis)” 현상도 층류→난류 전이로 박리점이 이동하며 생긴다.²

4. 후류와 항력[편집]

박리가 일어나면 물체 뒤편에 넓은 저압 영역, 즉 후류가 형성된다. 이 후류가 공학적으로 두 가지 골칫거리를 만든다.

첫째는 **형상 항력(form drag)**이다. 물체 앞면은 정체점에서 압력이 높은데 뒷면은 후류라 압력이 낮으니, 앞뒤 압력 차이가 물체를 뒤로 끌어당긴다. 유선형(streamlined) 형상이 각진 형상보다 저항이 훨씬 작은 것은 박리를 최대한 늦춰 후류를 좁게 만들기 때문이다. 물방울 모양이 자연이 찾아낸 최적해인 셈.

둘째는 **비정상성(unsteadiness)**이다. 원기둥처럼 뭉툭한 물체 뒤에서는 후류의 소용돌이가 위아래로 번갈아 떨어져 나가며 카르만 와열을 만든다. 이 주기적 와도 방출이 물체에 진동 하중을 걸어, 다리·굴뚝·송전선을 흔들고 심하면 파괴한다.³ 유체-구조 연성 해석이 필요해지는 대표적 상황.

5. 수치해석에서의 박리[편집]

박리는 CFD가 가장 애를 먹는 현상 중 하나다. 이유는 명확하다. 박리점의 위치와 후류의 크기가 난류 모델링에 극도로 민감하기 때문이다.

• RANS 모델은 박리를 곧잘 놓친다. 특히 $k$-$\varepsilon$ 계열은 역압력구배 하에서 박리를 너무 늦게 예측하는 경향이 악명 높다. 그래서 박리가 중요한 문제에는 벽 근처를 잘 다루는 $k$-$\omega$ SST 모델이 준(準)표준으로 쓰인다.⁴
• 벽면 격자 해상도가 결정적이다. 박리는 벽면 전단응력이 0이 되는 지점을 잡는 문제라, 경계층을 충분히 얇은 격자로 분해해야 한다($y^+ \approx 1$). 벽함수로 대충 넘기면 박리점이 엉뚱한 데 찍힌다.
• 비정상 해석이 필요할 때가 많다. 박리 후류는 본질적으로 시간에 따라 요동치므로, 정상상태(steady) 해석으로는 수렴 자체가 안 되고 잔차가 끝없이 진동한다. 이럴 땐 LES나 DES 같은 스케일 분해 접근으로 넘어가야 물리를 제대로 잡는다.

요컨대 “박리 잘 맞추면 CFD 좀 하는 사람”이라는 말이 괜히 있는 게 아니다. 격자, 난류 모델, 시간 이산화가 삼위일체로 맞아떨어져야 비로소 실험과 비슷한 박리점이 나온다.

6. 관련 문서[편집]

• 경계층
• 와도 · 카르만 와열
• 공력해석
• 나비에-스토크스 방정식
• 난류 모델링 · 벽함수
• 레이놀즈수
• 유체-구조 연성
• 전산유체역학 · 포텐셜 유동

7. Footnotes[편집]

1. 그래서 포텐셜 유동만 알던 시절의 이론은 골프공이 왜 딤플로 더 멀리 나가는지 설명하지 못했다. 박리는 점성과 경계층이 빚어내는 순수 실전형 현상이라, 비점성 이론의 사각지대에 정확히 놓여 있다. ↩

2. 매끈한 원기둥의 항력계수가 레이놀즈수 약 $3 \times 10^5$ 부근에서 1.2에서 0.3 언저리로 급락한다. 처음 데이터를 본 사람은 측정 실수인 줄 안다. 하지만 물리다. ↩

3. 1940년 타코마 내로스 다리 붕괴가 종종 카르만 와열의 사례로 인용되지만, 실제로는 비틀림 플러터가 주원인이었다는 게 정설이다. 그래도 유동박리가 구조물을 흔든다는 교훈 자체는 유효하다. ↩

4. 멘터(Menter)의 SST 모델은 벽 근처에서는 $k$-$\omega$, 자유 유동에서는 $k$-$\varepsilon$을 쓰도록 혼합한 하이브리드다. “둘 다 애매하면 섞어 쓴다”는 공학의 위대한 전통을 계승한 모델. ↩