[논문 리뷰] A Dispersion Minimized Mimetic Method for Cold Plasma

Vrushali A. Bokil, Vitaliy Gyrya외 1명|arXiv (Cornell University)|2016. 04. 05.

Electromagnetic Simulation and Numerical Methods참고 문헌 1인용 수 3

한 줄 요약

이 논문은 냉성 플라즈마에서 맥스웰 방정식을 위한 산란 최소화 미메틱 유한차분(MFD) 방법을 제안한다. 이 방법은 지수 시간 차분(ETD)과 새로운 m-적응 최적화 기법을 사용한다. 공간에서 일반화된 질량 러핑 MFD 스킴에 m-적응을 적용하고 시간 영역에서 ETD를 사용함으로써, 4차 수준의 수치적 산란 오차를 달성한다. 이는 기존의 표준 2차 스킴 대비 오차를 두 개의 주자 정도 감소시킨다. 이는 수치적 및 이론적 방법으로 검증되었다.

ABSTRACT

In this paper we consider the lowest edge-based mimetic finite difference (MFD) discretization in space for Maxwell's equations in cold plasma on rectangular meshes. The method uses a generalized form of mass lumping that, on one hand, eliminates a need for linear solves at every iteration while, on the other hand, retains a set of free parameters of the MFD discretization. We perform an optimization procedure, called m-adaptation, that identified a set of free parameters that lead to the smallest numerical dispersion. The choice of the time stepping proved to be critical for successful optimization. Using exponential time differencing we were able to reduce the numerical dispersion error from second to fourth order of accuracy in mesh size. It was not possible to achieve this order of magnitude reduction in numerical dispersion error using the standard leapfrog time stepping. Numerical simulations independently verify our theoretical findings.

연구 동기 및 목표

• 냉성 플라즈마 매질 내 전자기파 시뮬레이션에서 수치적 산란 오차를 줄이기 위해, 장시간 및 전기적으로 큰 영역에서의 시뮬레이션에 있어 이 오차가 매우 중요하기 때문이다.
• 이전에 진공에서 성공적으로 적용된 m-적응 기법을, 추가적인 역학이 존재하는 선형 분산 매질인 냉성 플라즈마로 확장하기 위해 필요하다. 이는 문제의 비선형성으로 인해 비틀림이 발생하기 때문이다.
• 자유 매개변수를 최적화하여 수치적 산란을 최소화하면서도 명시적 시간 스텝을 유지하는, 미메틱 유한차분 스킴의 가족 내 최적의 자유 매개변수 집합을 규명하기 위해 필요하다.
• 지수 시간 차분(ETD)이 표준 레프로그 시간 스텝과 달리 고차수 산란 감소를 가능하게 하므로, 이 방법이 고차수 산란 감소를 가능하게 함을 보여주기 위해 필요하다.
• 기존의 4차 수준 방법보다 더 작은 스텐실 크기를 갖는 방법을 개발하여 계산 효율성을 향상시키기 위해 필요하다.

제안 방법

• 직각 격자에서 냉성 플라즈마 내 맥스웰 방정식의 공간 이산화를 위해 최저차수의 모서리 기반 미메틱 유한차분(MFD) 방법을 사용한다.
• 선형 해법이 매 시간 단계에서 필요 없도록 일반화된 질량 러핑을 적용함으로써, 최적화를 위한 자유 매개변수를 유지한다.
• m-적응 절차는 푸리에 도메인에서의 산란 관계 분석을 사용하여 이 자유 매개변수를 최적화하여 수치적 산란을 최소화한다.
• 시간 적분을 위해 지수 시간 차분(ETD)을 사용하여, 강한 선형 항을 정확하게 통합하고 고차수 정밀도를 달성한다.
• 결과적으로, 지수 시간 미메틱 유한차분(ETMFD) 스킴은 표준 ETMFD 스킴의 기초가 되는 2차 오차와는 달리 4차 수준의 수치적 산란 오차를 달성한다.
• 이 방법은 E 및 J 필드의 H에 대한 명시적 스타깅을 유지하여 반의미적 시간 평균을 피한다.

실험 결과

연구 질문

• RQ1 진공에서 성공적으로 적용된 m-적응 기법이, 분극 전류를 위한 보조 미분 방정식을 포함하는 냉성 플라즈마 모델로 성공적으로 확장될 수 있는가?
• RQ2 왜 레프로그 시간 스텝은 냉성 플라즈마의 MFD 가족에서 2차 수준을 초월한 산란 감소를 지원하지 못하는가? 반면 ETD는 이를 성공적으로 수행하는가?
• RQ3 일반화된 질량 러핑 MFD 스킴 내에서 수치적 산란을 최소화하는 데 최적의 자유 매개변수 조합은 무엇인가?
• RQ4 ETMFD 방법은 4차 수준의 수치적 산란 오차를 달성하는가? 그리고 이는 시뮬레이션을 통해 독립적으로 확인될 수 있는가?
• RQ5 최적화된 ETMFD 방법의 스텐실 크기는 기존 문헌의 다른 4차 수준 방법들과 비교해 어떻게 되는가?

주요 결과

• ETD를 사용할 경우 m-적응 절차는 메쉬 크기의 제곱에 따라 수치적 산란 오차를 2차에서 4차로 감소시킨다.
• 이론적 분석과 수치 시뮬레이션을 통해 ETMFD 방법이 수치적 산란 오차에서 4차 수렴을 보임을 확인하였다.
• ETMFD 스킴 내 최적의 자유 매개변수와 국소 질량 행렬은 진공의 경우와 동일하다. 이는 다양한 모델 간의 구조적 일관성을 시사한다.
• ETMFD 방법은 컴팩트한 스텐실 크기를 유지하여, 고차수 기반 이산화에 의존하는 다른 4차 수준 방법들보다 계산적 이점을 제공한다.
• 표준 레프로그 시간 스텝은 냉성 플라즈마 모델에서 m-적응을 통한 고차수 산란 감소를 지원하지 못함을 보여주며, 시간 적분 선택의 중요성을 강조한다.
• 특수 해의 수치 시뮬레이션 결과, 전기장과 전류 밀도의 L² 오차에서 4차 수렴이 관찰되어 이론적 결과를 검증하였다.