복잡한 형상 플라즈마 시뮬레이션을 위한 임의 곡선좌표 PIC 기법

초록

본 논문은 2차원 전기정적 플라즈마를 대상으로, 물리적 비균일 격자를 논리적 단위 정사각형 격자로 매핑하여 PIC 연산을 수행하는 ACC‑PIC(Arbitrary Curvilinear Coordinate PIC) 방법을 제안한다. 논리 격자에서 입자 이동을 비분리 형태의 방정식으로 기술하고, 이를 보존적이며 심볼릭한 특성을 유지하도록 수정된 반암시적 Leapfrog(Modified Leapfrog) 알고리즘으로 통합한다. 또한, 논리 격자에 맞춘 곡선좌표 형태의 포아송 방정식을 도입해 전기장을 계산하고, 전하 밀도와 전기장을 모두 논리 격자 상에서 가중·보간함으로써 물리 격자의 복잡성을 제거한다.

상세 분석

ACC‑PIC 방법은 전통적인 PIC이 물리적 격자에서 입자와 필드를 직접 연결해야 하는 복잡성을 해소한다는 점에서 혁신적이다. 물리 도메인의 복잡한 경계와 비균일 격자는 논리 도메인인 단위 정사각형으로 매핑되며, 이때 좌표 변환 Jacobian과 메트릭 텐서를 이용해 물리량을 변환한다. 핵심은 논리 격자 상에서 입자 운동 방정식이 일반적인 직교 좌표계에서와 달리 비분리 형태가 된다는 점이다. 이를 해결하기 위해 저자들은 반암시적 Modified Leapfrog(M‑L) 스킴을 확장했으며, 이 스킴은 시간 중심화된 전기장을 사용해 입자 위치와 속도를 동시에 업데이트한다. 중요한 것은 M‑L이 symplectic 구조를 보존한다는 점이다. 즉, 장기 시뮬레이션에서 에너지와 위상공간 부피가 인위적으로 증감하지 않아 물리적 정확도가 유지된다.

전기장 계산 부분에서는 곡선좌표 형태의 포아송 방정식을 도출한다. 물리적 전하 밀도 ρ(x,y)를 논리 좌표 (ξ,η) 로 변환할 때 Jacobian J와 메트릭 g^{ij}를 포함한 변환식이 필요하며, 이는 논리 격자 상에서 균일한 라플라시안 연산자를 적용할 수 있게 만든다. 전하 가중은 전통적인 Cloud‑in‑Cell(CIC) 방식과 동일하게 논리 격자에 수행되며, 물리 격자에서 발생할 수 있는 비균일 셀 크기에 따른 가중 오류를 회피한다. 전기장은 논리 격자에서 푸리에 변환 혹은 다중 그리드 해법 등 효율적인 선형 솔버로 풀고, 다시 메트릭을 이용해 물리 입자 위치에 보간된다.

이러한 구조는 복잡한 경계(예: 곡선형 전극, 비정형 진공 챔버)와 다중 스케일 격자(예: 고해상도 영역과 저해상도 영역을 동시에 포함)에서 PIC을 적용할 때 큰 장점을 제공한다. 또한, 기존 물리 격자 기반 PIC에서 발생하는 “가중·보간 오류 누적” 문제를 논리 격자에서 일관되게 처리함으로써 수치 확산과 전하 보존성을 크게 개선한다. 다만, 좌표 변환 과정에서 Jacobian 계산과 메트릭 텐서 저장이 추가 메모리를 요구하고, 비분리 방정식의 해석적 특성이 복잡해져서 구현 난이도가 상승한다는 단점도 존재한다.