시각화로 이해하는 Reynolds 수와 유동 천이

Reynolds 수란 무엇인가#

유체 유동을 이해할 때 가장 중요한 무차원수 중 하나가 **Reynolds 수(Re)**다. 이 수는 관성력과 점성력의 비를 나타낸다.

$Re = \frac{\rho U L}{\mu} = \frac{U L}{\nu}$

여기서:

• $\rho$: 유체 밀도 $[\text{kg/m}^3]$
• $U$: 특성 속도 $[\text{m/s}]$
• $L$: 특성 길이 $[\text{m}]$
• $\mu$: 동점성 계수 $[\text{Pa·s}]$
• $\nu = \mu/\rho$: 동점도 $[\text{m}^2/\text{s}]$

유동 체계의 구분#

유동 천이의 물리적 메커니즘#

층류에서 난류로의 천이는 켈빈-헬름홀츠(Kelvin-Helmholtz) 불안정성에서 시작된다. 속도 구배가 존재하는 경계에서 작은 교란이 증폭되면서 와류(vortex)가 형성되고, 이것이 연쇄적으로 에너지 캐스케이드를 일으킨다.

속도 프로파일의 변화를 수식으로 표현하면:

$\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla)\mathbf{u} = -\frac{1}{\rho}\nabla p + \nu \nabla^2 \mathbf{u}$

Navier-Stokes 방정식에서 좌변의 비선형 항 $(\mathbf{u} \cdot \nabla)\mathbf{u}$가 관성을 담당하고, 우변의 $\nu \nabla^2 \mathbf{u}$가 점성 감쇠를 담당한다. Reynolds 수는 이 두 항의 상대적 크기를 나타낸다.

속도장 시각화: Reynolds 수에 따른 유동 변화#

아래 시뮬레이션을 통해 Reynolds 수 변화에 따른 속도장의 변화를 직접 확인해보세요:

관찰 포인트:

• 저 Re (Re ≈ 1~50): 벡터가 부드럽고 규칙적으로 배열됨 — 점성이 교란을 즉시 감쇠
• 중간 Re (Re ≈ 100~500): 장애물 후류에서 비대칭성 출현 — 관성력이 점성력에 필적
• 고 Re (Re > 1000): 벡터 방향이 불규칙해지고 와류 구조가 발달

원기둥 주위 유선: 카르만 소용돌이길#

원기둥(cylinder)을 지나는 유동에서 Reynolds 수가 증가하면 주기적인 와류 방출 현상인 **카르만 와열(Kármán vortex street)**이 나타난다.

방출 주파수는 Strouhal 수 $St$로 정규화된다:

$St = \frac{f D}{U_\infty}$

여기서 $f$는 와류 방출 주파수, $D$는 원기둥 직경, $U_\infty$는 자유류 속도다. 원기둥의 경우 $St \approx 0.2$ ($100 < Re < 10^5$)로 거의 일정하게 유지된다.

아래 시뮬레이션에서 원기둥 후류의 유선 패턴을 확인하세요:

관찰 포인트:

• 원기둥 전면의 **정체점(stagnation point)**에서 유선이 갈라지는 것을 확인하세요
• 후류 영역에서 대칭이 깨지며 주기적 와류 방출이 발생하는 과정을 관찰하세요
• speed 값을 높이면 (빠른 유속) 후류 불안정성이 더 강해짐을 볼 수 있습니다

수치해석에서의 Reynolds 수: 격자 요구사항#

난류 시뮬레이션(DNS, Direct Numerical Simulation)에서 가장 작은 길이 스케일인 **콜모고로프 미세 스케일(Kolmogorov microscale)**까지 해상도가 필요하다:

$\eta = \left(\frac{\nu^3}{\varepsilon}\right)^{1/4}$

전체 격자 수는 Reynolds 수에 대해 다음과 같이 스케일된다:

$N \sim Re^{9/4}$

즉 $Re$가 10배 증가하면 격자 수는 약 178배 증가해야 한다. 이것이 실용적인 난류 시뮬레이션에 RANS, LES 같은 난류 모델이 필요한 이유다.

정리#

• Reynolds 수는 관성력/점성력의 비로, 유동 특성을 결정하는 핵심 무차원수
• 임계 $Re$를 넘으면 층류 → 천이 → 난류의 경로를 따름
• 원기둥 후류에서는 카르만 와열이 발생하며 Strouhal 수로 특성화
• DNS 난류 시뮬레이션의 격자 비용은 $Re^{9/4}$로 스케일되므로, 고 Re 유동에는 난류 모델이 필수

다음 주에는 **유한체적법(FVM)**을 이용해 이 Navier-Stokes 방정식을 어떻게 이산화하는지, 그리고 업윈드 vs 중앙차분 스킴이 정확도에 미치는 영향을 다룬다.