[논문 리뷰] A Hybrid Nodal-Staggered Pseudo-Spectral Electromagnetic Particle-In-Cell Method with Finite-Order Centering

Edoardo Zoni, Remi Lehé외 6명|arXiv (Cornell University)|2021. 06. 23.

Laser-Plasma Interactions and Diagnostics참고 문헌 61인용 수 4

한 줄 요약

이 논문은 유한차수 중심화를 사용하여 노드형 및 스태그어드 그리드를 조합하는 새로운 하이브리드 가짜 스펙트럼 전자기 PIC 방법을 제안한다. 이 방법은 수치적 안정성과 정확도를 향상시키기 위해 맥스웰 방정식을 스태그어드 Yee 그리드에서 풀고, 고차수 보간을 통해 전류를 노드 그리드에 기판하고 힘을 수집함으로써 수치적 체렌코프 불안정성(NCI)을 감소시킨다. 이는 레이저-플라즈마 가속 및 쌍 생성과 같은 다양한 플라즈마 물리 시나리오에서 순수한 노드형 또는 스태그어드 형식보다 뛰어난 성능을 보인다.

ABSTRACT

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

연구 동기 및 목표

• 상대론적 및 저밀도 플라즈마에서 전자기 PIC 시뮬레이션의 수치적 체렌코프 불안정성(NCI)을 해결한다.
• 그리드 스태그어딩으로 인한 고분산 또는 불안정성으로 인해 악화되는 순수한 노드형 또는 스태그어드 PIC 방법의 한계를 극복한다.
• 장치 해석에 안정성이 뛰어난 스태그어드 그리드를 사용해 장력을 풀고, 입자 푸시어의 힘 수집 정확도를 높이기 위해 노드 그리드의 장점을 활용하는 하이브리드 접근법을 개발한다.
• 노드 그리드와 스태그어드 그리드 간의 필드 및 전류 데이터를 중심화하기 위해 Fornberg 계수를 기반으로 한 유한차수 보간을 도입한다.
• 레이저 구동형 및 빔 구동형 플라즈마 웨이브필드 가속 및 전자-양전자 쌍 생성과 같은 다양한 물리 영역에서 향상된 성능을 입증한다.

제안 방법

• 전하 및 전류 기판이 이루어지는 노드 그리드와 맥스웰 방정식이 해결되는 스태그어드 Yee 그리드를 결합한 하이브리드 PIC 사이클을 사용한다.
• 노드 그리드와 스태그어드 그리드 간의 필드 및 전류 데이터를 고정밀도로 이동하기 위해 Fornberg 계수에 기반한 유한차수 보간을 적용한다.
• 시간 영역 정확도를 향상시키기 위해 세 가지 유형의 PSATD 방법—표준형, 갈릴레이형, 평균 갈릴레이형—을 스태그어드 그리드에 적응시킨다.
• 맥스웰 방정식의 푸리에 공간 해석적 적분을 사용하여 수치적 분산을 최소화하고 안정성을 유지한다.
• Fornberg 보간 계수를 사용해 그리드 양상의 푸리에 변환을 조절하는 전달 연산자를 정의함으로써, 보간 계수 ζF가 1에서 벗어나는 정도를 효과적으로 제어한다.
• 더 높은 차수의 보간을 통해 ζF가 1에서 벗어나는 영역(0 < ζF < 1)을 최소화함으로써 동역량 보존을 확보하고 NCI를 감소시킨다.

실험 결과

연구 질문

• RQ1 노드형 및 스태그어드 그리드를 조합한 하이브리드 PSATD-PIC 방법이 순수한 노드형 또는 스태그어드 형식보다 수치적 체렌코프 불안정성(NCI)을 더 효과적으로 줄일 수 있는가?
• RQ2 노드형과 스태그어드 그리드 간의 유한차수 중심화가 상대론적 플라즈마 시뮬레이션에서 전자기장 해석의 안정성과 정확도에 어떤 영향을 미치는가?
• RQ3 표준형, 갈릴레이형, 평균 갈릴레이형 등의 PSATD 변종 선택이 다양한 물리 영역에서 하이브리드 방법의 성능에 어떤 영향을 미치는가?
• RQ4 보간 계수 ζF는 NCI 억제에 어떤 역할을 하는가? 그리고 이 값이 1에서 벗어나는 정도는 어떤 유한차수 중심화의 차수에 따라 어떻게 달라지는가?
• RQ5 전자-양전자 쌍 생성 및 고강도 레이저-플라즈마 상호작용과 같은 극한 조건에서도 하이브리드 방법은 정확도와 안정성을 유지하는가?

주요 결과

• 하이브리드 방법은 순수한 노드형 또는 스태그어드 PIC 방법보다 수치적 체렌코프 불안정성(NCI)을 크게 감소시키며, 특히 저밀도 및 상대론적 영역에서 뚜렷한 효과를 보인다.
• 높은 차수의 유한차수 중심화(최대 2m = 8)는 보간 계수 ζF가 1에서 벗어나는 정도를 줄여, 불안정성이 발생하기 쉬운 영역(0 < ζF < 1)을 최소화한다.
• 레이저 구동형 및 입자 빔 구동형 플라즈마 웨이브필드 가속과 같은 다양한 물리적 시나리오에서 높은 정확도와 안정성을 유지한다.
• 수치적 시험 결과, 하이브리드 방법은 표준 노드형 및 스태그어드 PSATD 방법보다 필드 정확도와 장기 시뮬레이션 안정성 측면에서 뛰어난 성능을 보였다.
• Fornberg 기반 보간을 사용함으로써 노드 그리드와 스태그어드 그리드 간의 고차수 데이터 전달이 가능해져 정확도와 계산 효율성을 효과적으로 균형 잡았다.
• 분산 분석 결과, ζF가 1에 가까울수록 NCI가 억제됨을 확인하였으며, 이는 특히 하이브리드 구성에서 더 높은 차수의 중심화를 통해 더욱 효과적으로 달성된다.