플라즈마 시뮬레이션

계산물리 물리 마지막 수정: 2026-07-21 04:31:58

1. 개요[편집]

플라즈마 시뮬레이션 Plasma Simulation
대표 기법	PIC (Particle-In-Cell)
기술 스케일	유체(MHD) ↔ 운동론(kinetic)
응용	핵융합, 반도체 식각, 우주추진, 레이저-플라즈마

기체를 너무 뜨겁게 데우면 전자가 도망가서 이온만 남는다. 이 무법지대를 시뮬레이션하는 게 플라즈마 시뮬레이션이다.

플라즈마 시뮬레이션은 원자가 전자를 잃거나 얻어 이온화된 기체, 즉 플라즈마의 거동을 컴퓨터로 수치 모사하는 기법을 통칭한다. 플라즈마는 흔히 “물질의 제4상태”로 불리는데, 전기적으로 중성인 유체와 달리 하전 입자들이 전자기력으로 서로 얽혀 있어서 일반 유체역학의 직관이 곧잘 배신당한다. 거시적 규모에서는 자기유체역학(MHD)의 유체 근사로 뭉뚱그려 풀 수 있지만, 입자 하나하나의 속도 분포가 중요해지는 순간 그 근사가 무너지고 훨씬 미시적인 운동론적(kinetic) 기술이 필요해진다.

핵융합로 내부의 초고온 플라즈마 가둠, 반도체 웨이퍼를 식각하는 저온 플라즈마, 인공위성의 홀 추진기, 레이저로 물질을 순간 이온화하는 관성 봉입 핵융합까지 — 스케일도 온도도 천차만별인 이 현상들을 하나의 은탄환 기법으로 다 풀 수는 없다. 그래서 플라즈마 시뮬레이션은 유체·운동론·혼성(hybrid) 기법을 상황에 맞춰 골라 쓰는, 방법론 자체가 하나의 스펙트럼인 분야다.

2. 플라즈마의 특징 스케일[편집]

플라즈마를 “그냥 뜨거운 기체”가 아니라 별도의 물리로 취급하게 만드는 두 가지 핵심 척도가 있다. 하나는 데바이 길이 $\lambda_D$ 로, 하전 입자가 만든 전기장이 다른 입자들에 의해 가려지는(screening) 거리다. 플라즈마 내부에서는 $\lambda_D$ 보다 먼 거리에서 전기적으로 거의 중성처럼 보이는데, 이를 준중성(quasi-neutrality)이라 한다. 다른 하나는 플라즈마 진동수 $\omega_{pe} = \sqrt{n_e e^2/(\epsilon_0 m_e)}$ 로, 전자가 평형에서 살짝 밀려났을 때 되돌아오며 진동하는 고유 주파수다. 이 두 척도가 시뮬레이션에 필요한 격자 크기와 시간 간격의 하한선을 사실상 결정해버린다.¹

3. 유체 기술과 운동론적 기술 사이[편집]

플라즈마를 기술하는 방법은 크게 두 갈래로 나뉜다.

유체(fluid) 기술: 입자 분포를 밀도·평균속도·온도 같은 거시 모멘트로 뭉개서 나비에-스토크스 방정식 꼴의 보존법칙으로 다룬다. 자기유체역학이 대표적이며 계산 비용이 싸지만, 입자 속도 분포가 맥스웰 분포에서 크게 벗어나는 상황(무충돌 충격파, 빔-플라즈마 불안정성 등)에서는 정확도가 무너진다.
운동론적(kinetic) 기술: 위상공간에서 입자 분포함수 $f(\mathbf{x}, \mathbf{v}, t)$ 가 시간에 따라 어떻게 진화하는지를 직접 추적하는 블라소프 방정식을 푼다. 정확하지만 위상공간 차원이 6차원(3차원 위치 + 3차원 속도)이라 격자로 직접 풀기엔 계산량이 감당이 안 된다.

이 딜레마를 절충한 것이 다음 절의 PIC 기법이다.

4. PIC (Particle-In-Cell) 기법[편집]

PIC는 분포함수를 격자로 직접 이산화하는 대신, 수많은 “슈퍼입자(super-particle)” — 실제 입자 수백만 개를 대표하는 가상 입자 — 를 공간에 뿌려놓고 이들의 궤적을 추적해 분포함수를 통계적으로 근사한다. 한 스텝은 대략 이런 순환으로 이뤄진다.

가중(weighting): 입자들의 위치·전하를 격자점에 뿌려 전하밀도 $\rho$ , 전류밀도 $\mathbf{J}$ 를 계산.
장 풀이(field solve): $\rho$ , $\mathbf{J}$ 를 소스로 맥스웰 방정식(또는 정전근사의 포아송 방정식)을 격자 위에서 풀어 $\mathbf{E}$ , $\mathbf{B}$ 를 구한다. 이 단계에 흔히 FDTD 기법이 동원된다.
밀기(push): 격자의 장을 각 입자 위치로 보간해 로렌츠 힘을 구하고, 뉴턴 운동방정식으로 입자를 한 스텝 전진시킨다. 여기서 표준으로 쓰이는 것이 시간 대칭적이고 에너지 보존이 좋은 **보리스 알고리즘(Boris algorithm)**이다.²
가중 되돌리기: 새 입자 위치로 다시 1번으로 돌아간다.

이 격자-입자-격자를 오가는 구조 때문에 “Particle-In-Cell”이라는 이름이 붙었다. 입자를 개별적으로 추적하니 비맥스웰 분포나 극단적 비평형 현상도 자연스럽게 잡아내는 것이 최대 장점. 정전기적 상호작용만 고려하는 **정전 PIC(electrostatic PIC)**는 포아송 방정식만 풀면 되어 가볍고, 자기장까지 스스로 발생시키는 강한 전류가 있는 문제는 맥스웰 방정식을 통째로 푸는 **전자기 PIC(electromagnetic PIC)**로 넘어가야 한다.

플라즈마와 별개로 중성 기체 입자의 희박 기체 유동을 다루는 직접모사 몬테카를로법(DSMC)도 구조상 PIC와 사촌뻘이다. 격자-입자 순환 대신 입자 간 충돌을 확률적으로 처리한다는 점이 다르지만, “입자를 통계적으로 추적해 연속체 근사가 무너지는 영역을 메운다”는 철학은 똑같다.

5. 수치적 난제[편집]

PIC의 정확도는 슈퍼입자 수, 격자 해상도, 시간 간격 세 가지에 인질로 잡혀 있다. 격자 간격이 데바이 길이보다 크면 전자가 자기 자신의 그림자를 못 알아보고 인위적으로 가열되는 격자 가열(grid heating) 현상이 생긴다. 시간 간격도 플라즈마 진동수의 역수보다 훨씬 작아야 한다는 CFL 조건 유사 제약이 붙는다 ( $\omega_{pe}\Delta t \ll 1$ ). 여기에 슈퍼입자 수가 유한한 탓에 생기는 통계적 노이즈까지 겹쳐, 실무자는 “격자를 촘촘히 할수록 정확해지지만 그럴수록 실행 시간이 우주적으로 늘어난다”는 트레이드오프와 매일 씨름한다.³ 이 비용 문제 때문에 이온은 유체로, 전자만 입자로 다루는 하이브리드 기법이나, 몇 개 스텝을 묶어 넘기는 암시적 PIC 변형도 활발히 연구된다.

6. 응용 분야[편집]

핵융합: 토카막 가장자리의 난류 수송, 자기 재결합 등 유체 근사가 못 잡는 미시적 불안정성 해석.
반도체 공정: 웨이퍼 식각·증착에 쓰이는 저온 플라즈마 방전 시뮬레이션. 이온 에너지 분포가 식각 프로파일을 좌우해 반도체 소자 시뮬레이션 공정 최적화와 맞닿아 있다.⁴
우주 추진: 홀 추진기 같은 전기추진 장치 내부 이온 가속 과정.
관성 봉입 핵융합: 레이저가 표적을 순간 이온화하며 만드는 초고밀도 플라즈마의 몬테카를로-PIC 혼성 해석.

7. 관련 문서[편집]

지구 전리층 플라즈마는 데바이 길이가 센티미터 단위지만, 핵융합로 코어 플라즈마는 마이크로미터 단위다. 같은 “플라즈마”라는 이름 아래 격자 설계 난이도는 하늘과 땅 차이라는 뜻. ↩
보리스(1970)가 제안한 이 알고리즘은 로렌츠 힘의 전기장·자기장 항을 분리해 자기장 회전을 기하학적으로 정확히 처리하는 트릭을 쓴다. 겉보기엔 소박한 벡터 연산 몇 줄인데, 반세기가 지나도 이보다 나은 표준 대체재가 잘 안 나온다 — 플라즈마 시뮬레이션계의 스테디셀러. ↩
실무자들 사이에서는 “PIC 시뮬레이션은 컴퓨터 클러스터가 다 한다”는 말이 있다. 입자 수를 늘리는 순간 시간은 기하급수로, 전기요금은 산술급수로(사실 이것도 기하급수로) 늘어나기 때문. ↩
웨이퍼 식각 공정 엔지니어에게 “PIC 시뮬레이션 결과가 이렇습니다”라고 보고하면 열에 아홉은 “그래서 수율이 오르나요?”라고 되묻는다. 물리는 우아해도 공장은 숫자로만 말한다. ↩