블라소프 방정식

1. 개요[편집]

블라소프 방정식(Vlasov equation)은 입자 간 충돌을 무시할 수 있는 무충돌 플라즈마(collisionless plasma)에서, 위상공간(phase space) 상의 입자 분포함수 $f(\mathbf{x}, \mathbf{v}, t)$의 시간 진화를 지배하는 운동론적(kinetic) 방정식이다. 6차원 위상공간(위치 3 + 속도 3)에서 정의된 분포함수가 자기 자신이 만드는 장(field) 속을 흐른다는, 자기일관(self-consistent) 구조가 핵심이다. 러시아 물리학자 아나톨리 블라소프(Anatoly Vlasov)가 1938년에 제안했다.

플라즈마 시뮬레이션의 이론적 토대이며, 유체 근사(자기유체역학)가 놓치는 속도공간 구조 — 대표적으로 란다우 감쇠 같은 무충돌 감쇠 — 를 담아낼 수 있다는 점에서 플라즈마 물리의 정중앙에 있는 방정식이다.

2. 방정식의 형태[편집]

블라소프 방정식은 사실상 위상공간에서의 연속 방정식이다. 분포함수 $f$가 위상공간 흐름을 따라 보존된다($Df/Dt = 0$)는 진술이다.

$$\frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{x}} f + \frac{\mathbf{F}}{m} \cdot \nabla_{\mathbf{v}} f = 0$$

전자기 플라즈마에서 힘은 로런츠 힘 $\mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B})$이다. 결정적인 부분은 이 $\mathbf{E}$, $\mathbf{B}$가 외부에서 주어지는 게 아니라 분포함수 자신이 만든다는 점. 전하밀도와 전류밀도를 분포함수의 속도 모멘트로 얻고

$$\rho = q\int f\, d^3v, \qquad \mathbf{J} = q\int \mathbf{v} f\, d^3v$$

이걸 맥스웰 방정식(정전 근사에서는 푸아송 방정식)에 넣어 장을 구한 뒤, 그 장을 다시 블라소프 방정식에 돌려주는 폐루프를 이룬다. 이 비선형 자기일관 결합이 방정식을 어렵고도 풍부하게 만든다. 우변이 0인 이유는 충돌항을 버렸기 때문이며, 충돌을 넣으면 우변에 충돌 적분이 붙어 볼츠만 방정식이 된다.

3. 란다우 감쇠[편집]

블라소프 방정식이 유명해진 결정적 사건은 1946년 레프 란다우의 란다우 감쇠(Landau damping) 예측이다. 충돌이 전혀 없는데도 플라즈마 파동이 감쇠한다는, 당시로선 충격적인 결과였다. 물리적으로는 파동의 위상속도와 비슷한 속도로 달리는 입자들이 파동과 공명(resonance)하며 에너지를 주고받는데, 파동보다 조금 느린 입자가 조금 빠른 입자보다 많으면(맥스웰 분포의 꼬리에서 기울기가 음수) 알짜로 파동이 입자에게 에너지를 빼앗겨 감쇠한다. 엔트로피 증가 없이(가역적으로) 일어나는 감쇠라는 점이 물리적으로 미묘하며, 이 발견 이후 위상공간에서만 보이는 미세 구조(filamentation)의 중요성이 부각되었다.¹

4. 수치 해법: PIC와 오일러 방법[편집]

6차원 위상공간을 격자로 이산화하면 저장량이 폭발한다. 각 차원을 100점만 잡아도 $100^6 = 10^{12}$점이라 메모리가 순식간에 증발한다.² 그래서 두 갈래 접근이 발전했다.

• PIC(Particle-In-Cell): 분포함수를 다수의 대표 입자(macro-particle)로 표본화하고, 입자를 장 속에서 움직이며 장은 격자에서 푼다. 라그랑주적 입자와 오일러적 격자를 오가는 하이브리드다. 저장량이 위상공간 차원에 지수적으로 늘지 않아 사실상 표준이 되었다. 자세한 내용은 플라즈마 시뮬레이션 참고.
• 오일러/세미-라그랑주 블라소프 솔버: 위상공간에 직접 격자를 깔고 분포함수를 이류시킨다. 통계 잡음이 없어 저밀도 꼬리나 미세 구조를 깨끗하게 보지만, 고차원의 저주로 계산량이 막대하다.

PIC의 최대 약점은 유한 입자 수에서 오는 통계 잡음이며, 이는 입자 수의 제곱근에 반비례해 줄어든다. 그래서 “입자를 더 갈아 넣으면 조용해진다”는 것이 플라즈마 시뮬레이션의 국룰.

위 시뮬레이션은 다수 입자가 각자 만드는 장 속에서 스스로 궤도를 그리는 자기일관 다체계를 보여준다. 블라소프-PIC의 발상 — “분포함수를 입자로 표본화하고, 그 입자들이 만든 장이 다시 입자를 움직인다” — 과 구조가 같다. 다만 이 데모는 중력 커널을 쓰는 근사 시각화이지 전자기 플라즈마 솔버가 아니다.³

5. 성질과 보존량[편집]

블라소프 방정식은 위상공간에서 비압축 흐름(리우빌 정리)을 이루므로 위상공간 부피와 분포함수 값이 흐름을 따라 보존된다. 그 결과 $\int G(f)\, d^3x\, d^3v$ 꼴의 모든 카시미르 불변량(총 입자 수, 운동에너지, 엔트로피, $L^2$ 노름 등)이 보존된다. 수치 솔버가 이 보존량들을 얼마나 잘 지키는지가 품질의 척도이며, 특히 에너지와 위상공간 부피를 지키는 심플렉틱(symplectic) 시간 적분이 선호된다. 해밀토니안 역학 구조를 그대로 물려받기 때문이다.

한 가지 반직관적인 결과는 분포함수가 시간이 지날수록 속도공간에서 점점 가늘게 접히는 필라멘테이션(filamentation)이다. 위상공간 부피가 보존되므로 구조가 뭉개질 수 없고 대신 무한히 가늘어진다. 이 미세 구조는 어떤 유한한 격자로도 결국 해상되지 않아, 오일러 블라소프 솔버에서 수치적 필터링이나 재격자화가 불가피해진다. 란다우 감쇠 역시 이 필라멘테이션의 거시적 그림자로 이해할 수 있다. 실용적 규모의 문제에서는 완전한 6차원 대신 자이로 평균으로 차원을 줄인 자이로운동론(gyrokinetics)이 토카막 난류 해석의 표준 도구로 쓰인다.

6. 관련 문서[편집]

• 플라즈마 시뮬레이션
• 자기유체역학
• 데바이 길이
• 해밀토니안 역학
• 전자기학
• 계산물리
• 나비에-스토크스 방정식

7. Footnotes[편집]

1. 란다우 감쇠는 오랫동안 “정말 물리적으로 일어나는가, 수학적 인공물인가” 논쟁의 대상이었다. 2011년 무쇼-비야니(Mouhot & Villani)가 비선형 란다우 감쇠를 엄밀히 증명해 빌라니는 필즈상을 받았다. 감쇠 하나 증명하는 데 60여 년이 걸린 셈. ↩

2. 이걸 “차원의 저주(curse of dimensionality)“라 부른다. 위상공간이 6차원이라는 사실 하나가 무충돌 플라즈마 시뮬레이션 전체를 어렵게 만든다. 그래서 다들 PIC로 도망친다. ↩

3. 중력은 항상 인력이고 전자기력은 인력·척력이 공존하므로 물리가 근본적으로 다르다(중력에는 데바이 차폐가 없다). 이 데모는 어디까지나 “자기일관 다체 궤도”라는 발상을 눈에 보이게 하려는 근사일 뿐이다. ↩