코로나 방전 시뮬레이션을 위한 효율적인 반라그랑지안 알고리즘

초록

본 논문은 전하 운반을 지배하는 대류‑확산 연속 방정식을 풀기 위해 위치‑상태 분리 방법(POSS)을 제안한다. Eulerian 스킴으로 확산·반응·가속을, Lagrangian 스킴으로 선형 대류를 처리해 플럭스 보정 없이 높은 해상도와 낮은 연산 비용을 달성한다. 특히 전기장이 천천히 변하는 코로나 방전 상황에서 CFL 제한보다 큰 시간 단계가 가능해 장거리·장시간 시뮬레이션에 유리하다.

상세 분석

POSS는 연속 방정식의 대류 항을 ‘위치’와 ‘상태’ 두 부분으로 분리하는 혁신적인 접근법이다. 위치 파트는 입자들의 이동 경로를 Lagrangian 방식으로 추적해 선형 대류를 정확히 해결하고, 상태 파트는 전하 농도·전위·반응 등을 Eulerian 격자상에서 전통적인 유한차분·유한체적 스킴으로 계산한다. 이때 플럭스 보정 절차가 필요 없으며, 대류‑확산 결합 문제에서 흔히 발생하는 수치 확산과 비물리적 진동을 최소화한다. 시간 복잡도는 대류 단계에서 입자 이동을 단순히 좌표 변환으로 처리하므로 O(N)이며, N은 자유도(격자점) 수와 동일하다. 또한 전기장이 크게 변하지 않는 경우(코로나 방전의 전형적 상황)에는 대류 속도가 느리게 변하므로 CFL 조건에 얽매이지 않고 시간 간격을 크게 늘릴 수 있다. 이는 기존 고정밀 Eulerian 스킴이 요구하는 작은 Δt와 비교해 시뮬레이션 전체 시간을 수십 배 단축한다. 실험에서는 1‑D 평면, 2‑D 직교, 원통형·구형 좌표계 등 다양한 기하학적 설정에서 POSS가 정확도와 효율성 면에서 기존 고해상도 스킴을 능가함을 보였다. 특히 전하 구름이 급격히 이동하거나 경계층이 얇은 경우에도 안정적인 결과를 제공한다. 이러한 특성은 장거리 전극 간격(수십 cm 이상)과 수백 µs 이상의 장시간 전압 펄스를 다루는 실제 코로나 방전 모델링에 직접 적용 가능하게 만든다.