[논문 리뷰] Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration

Olga V. Shapoval, Remi Lehé외 4명|arXiv (Cornell University)|2021. 04. 28.

Laser-Plasma Interactions and Diagnostics참고 문헌 29인용 수 7

한 줄 요약

이 논문은 플라즈마 기반 가속기의 부스팅 프레임 입자-격자(PI-Cell) 시뮬레이션에서 시간단계 제한을 극복하기 위해 평균화된 갈릴레이 PSATD 알고리즘을 제안한다. 시간 주기 동안 전자기장의 분석적 평균을 푸리에 공간에서 수행함으로써, 수치적 불안정성을 유발하지 않으면서도 큰 시간단계(c∆t ≫ ∆x, ∆y)를 허용하게 되어, 고해상도 횡방향 시뮬레이션에서 부스팅 프레임 기법의 전체 계산 속도 향상 잠재력을 해방시킨다.

ABSTRACT

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

연구 동기 및 목표

• 횡방향 격자 크기 제약에도 불구하고 종방향 해상도가 완화된 상황에서 부스팅 프레임 PIC 시뮬레이션의 시간단계 제한 문제를 해결하기 위해.
• PSATD 기반 부스팅 프레임 시뮬레이션에서 큰 시간단계에서 재발하는 수치적 체렌코프 불안정성(NCI)을 해결하기 위해.
• PSATD 시뮬레이션에서 c∆t ≫ ∆x, ∆y 수준의 큰 시간단계를 사용하면서도 시뮬레이션 정확도와 안정성을 유지하기 위해.
• 고해상도 횡방향 해상도가 필요한 다양한 플라즈마 가속기 구성에 대해 부스팅 프레임 기법의 계산 속도 향상을 확장하기 위해.

제안 방법

• 시간 간격 [(n−1/2)∆t, (n+1/2)∆t] 동안 전자기장을 푸리에 공간에서 분석적으로 평균화하는 평균화된 갈릴레이 PSATD 알고리즘을 도입한다.
• 시간 평균 전기장 ⟨ˆE⟩n 과 자기장 ⟨ˆB⟩n 을 계산하기 위해 PSATD 프레임워크를 수정하여, 시간단계 동안 전류가 일정한 상황에서 맥스웰 방정식의 해석적 해를 활용한다.
• 갈릴레이 이동과 스펙트럼 정확도를 고려한 복소계수(Ψ1, Ψ2, A1, A2, Cρ)를 포함하는 새로운 필드 갱신 식을 유도한다.
• 입자 추진을 위해 시간 평균 필드를 사용하도록 표준 PIC 루프를 수정하여, 필드 수집 및 운동량 갱신의 일관성을 확보한다.
• 다양한 격자 해상도에서 높은 정확도와 안정성을 유지하기 위해 푸리에 공간에서 p차수 이산화를 적용한다.
• CPU 및 GPU 플랫폼을 지원하는 오픈소스 WarpX 코드에 알고리즘을 구현한다.

실험 결과

연구 질문

• RQ1 큰 시간단계(c∆t ≫ ∆x, ∆y)를 부스팅 프레임 PSATD 시뮬레이션에서 안정적으로 사용할 수 있는가, 수치적 불안정성이 발생하지 않는가?
• RQ2 CFL 자유인 표준 갈릴레이 PSATD 알고리즘이 c∆t > ∆x, ∆y 인 경우에도 불안정해지는 이유는 무엇인가?
• RQ3 푸리에 공간에서 시간 평균을 통해 필드 수집 과정을 수정함으로써 이 불안정성을 완화시킬 수 있는가?
• RQ4 새로운 알고리즘이 실제 플라즈마 가속기 시나리오에서 물리적 정확성과 계산 속도 향상을 어느 정도 유지하는가?
• RQ5 이러한 방법은 고해상도 횡방향 해상도가 필요하지만 고정밀 시간 해상도가 필요하지 않은 다른 PIC 시뮬레이션에 일반화 가능한가?

주요 결과

• 평균화된 갈릴레이 PSATD 알고리즘은 큰 시간단계(c∆t ≫ ∆x, ∆y)에서 재발하는 수치적 체렌코프 불안정성(NCI)을 성공적으로 억제하여 안정적인 시뮬레이션을 가능하게 한다.
• c∆t = ∆z 수준(횡방향 격자 크기보다 훨씬 큰)의 시뮬레이션에서 NCI 모드의 성장이 관찰되지 않으며, 이는 표준 갈릴레이 PSATD 알고리즘과 대조된다.
• 2D 균일 플라즈마 시뮬레이션에서 신규 알고리즘은 c∆t가 ∆z 수준까지 안정성을 유지하지만, 표준 방법은 c∆t > ∆x일 때 불안정해진다.
• 2D 레이저 웨이브필드 가속(LWFA) 시뮬레이션에서 평균화된 PSATD 방법은 큰 시간단계에서도 정확한 웨이브필드 구조와 전자 빔 역학을 재현하며, 기준 시뮬레이션과 일치한다.
• 3D 플라즈마 웨이브필드 가속(PWFA) 시뮬레이션에서 이 방법은 c∆t ≫ ∆x, ∆y 수준에서 안정적이고 고속의 시뮬레이션을 가능하게 하여 뚜렷한 계산 속도 향상을 달성한다.
• 알고리즘은 WarpX 코드에 구현되고 검증되었으며, 적응형 메시 리피팅 및 다중 플랫폼(CPU/GPU) HPC 환경과의 호환성을 입증했다.