자기유체역학

편집 역사 토론
전산유체역학 계산물리 물리 마지막 수정: 2026-07-21 04:27:19

1. 개요[편집]

자기유체역학
Magnetohydrodynamics
약칭MHD
분야유체역학 × 전자기학 × 계산물리
지배 방정식유체 방정식 + 맥스웰 방정식
응용핵융합, 태양물리, 액체금속 냉각재

유체도 어렵고 전자기학도 어려운데, 이 둘을 한 방정식계에 우겨넣은 게 자기유체역학이다.

자기유체역학(magnetohydrodynamics, MHD)은 전기를 통하는 유체 — 플라즈마, 액체금속, 전해질 등 — 의 운동과 그 유체가 만들거나 겪는 전자기장의 상호작용을 하나의 연속체 이론으로 기술하는 학문이다. 유체가 자기장 속에서 움직이면 유도 전류가 흐르고, 그 전류가 다시 로렌츠 힘으로 유체를 밀어내는 되먹임(feedback) 구조를 갖는다. 즉 나비에-스토크스 방정식과 맥스웰 방정식이 서로의 소스항으로 끼어드는, 두 분야의 불행한 결혼이라 할 수 있다.

핵융합로 속 플라즈마를 자기장으로 가두는 문제, 태양 흑점과 코로나 질량방출, 지구 핵의 대류가 지구 자기장을 만드는 다이나모 작용, 반도체 공정의 액체금속 교반까지 — 전도성 유체가 있는 곳엔 어김없이 MHD가 등장한다. 문제는 이 방정식계가 순수 유체역학보다 훨씬 뻣뻣하고(stiff) 다루기 까다롭다는 것.1

2. 지배 방정식[편집]

이상 MHD(ideal MHD)의 뼈대는 질량·운동량·에너지 보존에 로렌츠 힘 항과 유도 방정식(induction equation)을 더한 것이다. 운동량 방정식은 압력 경도항 외에 자기력 J×B\mathbf{J} \times \mathbf{B}가 추가된다.

ρ(∂u∂t+u⋅∇u)=−∇p+J×B+μ∇2u\rho \left( \frac{\partial \mathbf{u}}{\partial t} + \mathbf{u}\cdot\nabla\mathbf{u} \right) = -\nabla p + \mathbf{J}\times\mathbf{B} + \mu\nabla^2\mathbf{u}

자기장이 유체와 함께 어떻게 진화하는지는 유도 방정식이 결정한다.

∂B∂t=∇×(u×B)+η∇2B\frac{\partial \mathbf{B}}{\partial t} = \nabla \times (\mathbf{u} \times \mathbf{B}) + \eta \nabla^2 \mathbf{B}

우변 첫 항은 자기장이 유체 운동에 실려 이류(advect)되는 항이고, 둘째 항은 저항(resistivity) η\eta에 의한 확산항이다. 여기에 자기장은 발산이 없어야 한다는 ∇⋅B=0\nabla \cdot \mathbf{B} = 0 제약이 항상 따라붙는다 — 자석을 아무리 쪼개도 N극과 S극만 나오는 단극자 없음의 수학적 표현이다.

3. 무차원수: 하르트만 수와 자기 레이놀즈수[편집]

MHD 유동의 성격은 몇 가지 무차원수로 요약된다. 하르트만 수(Hartmann number) Ha=BLσ/μHa = BL\sqrt{\sigma/\mu}는 자기력과 점성력의 비를 나타내며, HaHa가 클수록 자기장이 유동 단면을 납작하게 눌러 하르트만 층(Hartmann layer)이라는 특유의 경계층을 만든다. 자기 레이놀즈수 Rm=μ0σuLRm = \mu_0 \sigma u L는 자기장의 이류와 확산 중 어느 쪽이 우세한지를 가른다. Rm≫1Rm \gg 1이면 자기력선이 유체에 “얼어붙어(frozen-in)” 함께 움직이는 이상 MHD 영역이고, Rm≪1Rm \ll 1이면 자기장이 그냥 확산해 흩어져버린다. 실험실 액체금속에서는 RmRm이 대개 1보다 훨씬 작지만, 천체 규모 플라즈마에서는 거리 LL이 압도적으로 커서 RmRm이 천문학적으로 커진다 — 그래서 우주의 자기장은 잘 안 사라지고, 실험실 다이나모는 재현하기가 지독히 어렵다.2

4. 이상 MHD와 저항 MHD[편집]

저항 η=0\eta = 0으로 놓은 이상 MHD에서는 앞서 말한 동결(frozen-in flux) 정리가 정확히 성립해, 자기력선이 유체 입자에 꼬리표처럼 붙어 다닌다. 이 근사는 태양풍이나 핵융합 플라즈마의 거시적 거동을 잘 설명하지만, 자기력선이 끊어졌다 다시 이어지는 자기 재결합 같은 국소적 사건은 설명하지 못한다. 자기 재결합은 태양 플레어의 폭발적 에너지 방출 메커니즘으로 지목되는데3, 이걸 제대로 다루려면 저항이나 운동론적(kinetic) 효과를 되살린 저항 MHD 혹은 확장 MHD가 필요하다. 유체 근사 자체가 무너지는 더 미시적인 스케일은 플라즈마 시뮬레이션 쪽의 영역이다.

5. 수치기법: div B = 0을 지키는 싸움[편집]

MHD 코드의 가장 큰 골칫거리는 수치 오차로 ∇⋅B\nabla \cdot \mathbf{B}가 조금씩 쌓여 물리적으로 존재해선 안 될 “자기 단극자”가 격자 위에 생겨나는 것이다. 이걸 억누르는 접근법은 크게 세 갈래다.

  • 제약 수송법(Constrained Transport): 자기장을 셀 면(face)에 정의하고 패러데이 법칙을 이산적으로 정확히 만족시켜 ∇⋅B=0\nabla \cdot \mathbf{B} = 0을 기계적으로 보장한다.
  • 발산 세정(Divergence Cleaning): 여분의 스칼라장을 도입해 쌓인 발산 오차를 파동처럼 도메인 밖으로 흘려보낸다.
  • 투영법(Projection method): 유한체적법 기반 코드에서 매 시간 스텝 후 포아송 방정식을 풀어 발산 성분을 빼낸다.

MHD의 지배방정식 자체는 압축성 압축성 유동과 마찬가지로 쌍곡형이라 리만 솔버 기반 Godunov류 기법이 표준이며, 여기에 위 세 기법 중 하나를 얹어 ∇⋅B\nabla \cdot \mathbf{B} 오차를 관리하는 것이 현대 MHD 코드의 공통 설계다.

6. 응용 분야[편집]

  • 핵융합: 토카막이나 스텔러레이터에서 고온 플라즈마를 자기장으로 벽에 안 닿게 가두는 자기 가둠(magnetic confinement)의 안정성 해석에 MHD가 핵심이다. 플라즈마가 자기장을 뚫고 도망가려는 불안정성(kink, ballooning mode)을 예측하는 것이 로 설계의 사활이 걸린 문제.
  • 천체물리: 태양 코로나, 항성풍, 강착원반, 은하 자기장 — 우주의 거의 모든 뜨거운 물질은 전도성이라 MHD 없이는 설명이 안 된다. 지구·태양 자기장 자체의 기원을 설명하는 다이나모 이론도 결국 MHD 방정식을 도는 대류 시뮬레이션의 산물이다.
  • 산업: 액체금속 냉각 원자로, 전자기 교반(EM stirring)을 이용한 제철 공정, 전자기 펌프 등.4

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 유체 방정식과 맥스웰 방정식의 시간 스케일이 크게 다른 경우가 많아, 명시적 시간 전진 시 CFL 조건이 순수 유체 문제보다 훨씬 빡빡해진다. MHD 솔버 개발자들이 암시적 기법에 목매는 이유. ↩

  2. 지구 자기장을 만드는 핵의 다이나모 작용을 실험실 규모로 재현하려는 시도(리가·카를스루에 다이나모 실험 등)는 있지만, 자기 레이놀즈수를 자연 규모로 끌어올리기가 워낙 어려워 성공 사례가 손에 꼽힌다. ↩

  3. 재결합 자체는 초당 몇 미터 수준으로 느릴 것 같은데 실제 플레어는 수 분 만에 태양 표면을 뒤흔든다. 이 속도 미스매치를 메우는 “빠른 재결합” 메커니즘은 아직도 완전히 정리되지 않은 현역 연구 주제다. ↩

  4. 제철소의 전자기 교반 설비 앞에서 “이게 자기유체역학입니다”라고 말하면 대부분의 현장 기술자는 그냥 모터 취급을 한다. 방정식이 아무리 우아해도 현장에선 그저 “휘젓는 기계”일 뿐. ↩