유체역학

편집 역사 토론
유체역학 물리 마지막 수정: 2026-07-10 05:12:44

유체역학
Fluid Mechanics
분야고전역학 × 연속체역학
지배 방정식나비에-스토크스 방정식
핵심 무차원수레이놀즈수, 마하수, 프란틀수
하위 분야유체정역학, 유체동역학
전산 도구전산유체역학

1. 개요[편집]

유체역학(Fluid Mechanics)은 정지하거나 운동하는 유체(액체·기체)의 거동과, 그 유체가 경계나 물체에 작용하는 힘을 다루는 연속체역학의 한 분야이다. 여기서 “유체”란 아무리 작은 전단응력(shear stress)이라도 받으면 끊임없이 변형하는(계속 흐르는) 물질을 말한다. 고체는 전단응력을 받으면 어느 정도 변형한 뒤 버티며 멈추지만, 유체는 응력이 사라질 때까지 멈추지 않고 흐른다. 이 한 줄이 유체와 고체를 가르는 정의다.1

유체역학은 항공기 날개의 양력부터 수도관의 압력강하, 혈관 속 혈류, 은하 규모의 성간 가스까지 스케일을 가리지 않고 등장한다. 그리고 이 방정식들이 대부분 손으로 안 풀리기 때문에 컴퓨터를 갈아 넣는 전산유체역학이 태어났다. 즉 유체역학은 CFD의 상위 개념이자 부모 학문이다.

2. 연속체 가정[편집]

유체는 사실 어마어마하게 많은 분자로 이루어져 있지만, 유체역학은 분자 하나하나를 추적하지 않는다. 대신 유체를 연속적으로 분포한 물질로 취급하는 연속체 가정(continuum hypothesis)을 세운다. 밀도·속도·압력 같은 물리량을 공간의 각 점에서 정의된 매끄러운 장(field)으로 다루는 것이다.

이 가정이 성립하려면 우리가 관심 갖는 길이 스케일이 분자 평균자유행로(mean free path)보다 훨씬 커야 한다. 이 비율을 크누센수(Knudsen number) Kn=λ/LKn = \lambda/L로 정의하는데, Kn1Kn \ll 1이면 연속체 가정이 안전하다. 반대로 초고고도 희박기체나 미소유동(microflow)에서는 KnKn이 커져서 연속체 가정이 깨지고, 분자동역학이나 볼츠만 방정식 기반 접근이 필요해진다.2

3. 기술 방식 — 오일러 vs 라그랑주[편집]

유동을 수학적으로 기술하는 관점은 크게 둘로 나뉜다. 이 둘의 차이를 이해하는 게 유체역학 입문의 첫 관문이다.

  • 라그랑주 기술(Lagrangian): 개별 유체 입자를 따라다니며 그 입자의 위치와 속도를 시간의 함수로 추적한다. 강물 위에 나뭇잎을 띄우고 따라가는 관점. 입자 기반 기법인 SPHFLIP-PIC 유체가 이 관점을 쓴다.
  • 오일러 기술(Eulerian): 공간에 고정된 관측점을 두고, 그 지점을 지나가는 유동을 관찰한다. 다리 위에서 강물이 흘러가는 걸 내려다보는 관점. 대부분의 전산유체역학 코드가 격자에 고정된 오일러 관점을 채택한다.

두 관점을 잇는 다리가 물질미분(material derivative)이다.

DϕDt=ϕt+(u)ϕ\frac{D\phi}{Dt} = \frac{\partial \phi}{\partial t} + (\mathbf{u} \cdot \nabla)\phi

좌변은 유체 입자를 따라가며 본 물리량 ϕ\phi의 변화율(라그랑주), 우변 첫 항은 고정점에서의 시간변화, 둘째 항은 이류(advection)에 의한 변화(오일러)다. 나비에-스토크스 방정식의 그 악명 높은 비선형성이 바로 이 (u)u(\mathbf{u}\cdot\nabla)\mathbf{u} 항에서 나온다.

4. 지배 방정식[편집]

유체역학이 푸는 방정식은 결국 세 가지 보존 법칙으로 요약된다. 질량 보존(연속 방정식):

ρt+(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0

운동량 보존, 즉 뉴턴 제2법칙을 유체에 적용한 것이 그 유명한 나비에-스토크스 방정식이다. 뉴턴 유체(전단응력이 변형률에 비례하는 유체)에 대해서는 다음과 같이 쓴다.

ρDuDt=p+μ2u+ρg\rho\frac{D\mathbf{u}}{Dt} = -\nabla p + \mu \nabla^2 \mathbf{u} + \rho \mathbf{g}

여기에 열이 오가면 에너지 보존 방정식이 추가된다. 점성을 무시하는 이상유체 극한에서는 이 식이 오일러 방정식으로 단순해지고, 정상·비회전·비점성 유동에서는 그 유명한 베르누이 방정식 — “빠른 곳은 압력이 낮다” — 이 유도된다.

5. 층류와 난류[편집]

같은 유체라도 흐르는 양상은 천차만별이다. 조용히 층층이 미끄러지는 층류(laminar flow)와, 소용돌이가 뒤엉켜 카오스적으로 요동치는 난류(turbulent flow)를 가르는 기준이 레이놀즈수다.

Re=ρULμ=관성력점성력Re = \frac{\rho U L}{\mu} = \frac{\text{관성력}}{\text{점성력}}

레이놀즈수가 작으면 점성이 교란을 잡아먹어 층류가 유지되고, 크면 관성이 이겨 난류로 천이한다. 원관 내부 유동에서는 대략 Re2300Re \approx 2300이 임계값으로 통한다. 문제는 공학적으로 의미 있는 유동이 거의 전부 난류라는 것. 그래서 난류 모델링이 CFD의 절반쯤을 차지하는 거대한 하위 분야가 되었다.3

6. 무차원수 — 유체역학의 언어[편집]

유체역학은 무차원수(dimensionless number)로 말한다. 방정식을 무차원화하면 문제를 지배하는 소수의 무차원 그룹만 남고, 이 값들이 같으면 크기가 달라도 물리적으로 닮은 유동이 된다(상사 법칙). 대표 선수들:

무차원수정의 비율지배하는 물리
레이놀즈수관성/점성층류-난류 천이
마하수유속/음속압축성 유동
프란틀수운동량 확산/열 확산열전달 해석
프루드수관성/중력자유표면·조파 저항

이 덕분에 실물 항공기 대신 작은 축소 모형을 풍동에서 시험할 수 있다. 자세한 목록과 유도는 무차원수 문서 참고.

7. 왜 컴퓨터가 필요한가[편집]

나비에-스토크스 방정식은 비선형 연립 편미분방정식이라, 극히 단순한 형상을 빼면 해석해가 존재하지 않는다.4 그래서 공간을 격자로 쪼개고 유한체적법 등으로 이산화해 컴퓨터로 근사해를 구한다. 이 접근 전체가 바로 전산유체역학이며, 유체역학이라는 부모 학문의 방정식을 실용적으로 부려먹는 도구다. 유체역학이 “왜”와 “무엇을”을 정의하면, CFD가 “어떻게 푸느냐”를 맡는 셈이다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 물론 현실은 늘 애매하다. 케첩이나 치약처럼 응력이 어느 문턱을 넘어야 흐르기 시작하는 빙햄 유체(Bingham fluid), 젓는 속도에 따라 점도가 변하는 비뉴턴 유체들이 이 깔끔한 정의의 경계를 흐린다. 실리 퍼티(silly putty)는 천천히 누르면 흐르고 세게 던지면 튄다. 유체냐 고체냐, 답은 “시간 스케일 나름”.

  2. 크누센수가 커지는 영역에서는 나비에-스토크스가 아니라 격자 볼츠만 방법이나 DSMC(직접 시뮬레이션 몬테카를로) 같은 기법이 등판한다. 인공위성 표면의 희박기체 항력 계산이 대표적인 무대다.

  3. 난류는 리처드 파인만이 “고전물리학에 남은 가장 중요한 미해결 문제”라고 부른 걕물이다. 방정식은 1845년부터 있었는데 아직도 인류는 난류를 “이해”하지 못했다. 그런데도 매일 여객기를 잘만 띄우고 있으니, 공학이란 참 뻔뻔한 학문이다.

  4. 3차원 나비에-스토크스 방정식 매끄러운 해의 존재성·유일성은 클레이 밀레니엄 문제 중 하나로 상금이 100만 달러다. 유체역학은 인류가 매일 쓰면서도 그 수학적 바닥을 모르는 몇 안 되는 학문이다.