[논문 리뷰] A spectral, quasi-cylindrical and dispersion-free Particle-In-Cell algorithm

Remi Lehé, Manuel Kirchen외 3명|2015. 07. 16.

Spectroscopy and Laser Applications인용 수 3

한 줄 요약

이 논문은 원통대칭 시스템에서 전자기파 전파에 대한 분산 없음을 보장하는 새로운 스펙트럼적, 준원통좌표계 기반의 입자-장(PI) 알고리즘을 제안한다. 원통좌표계에서의 라디얼 및 축방향에 대해 스펙트럼 방법을 활용함으로써, 특히 제로 오더 수치 체렌코프 효과를 포함한 비물리적 수치 분산을 제거하면서도 레이저-플라즈마 및 비드-플라즈마 시뮬레이션의 계산 효율성을 유지한다.

ABSTRACT

We propose a spectral Particle-In-Cell (PIC) algorithm that is based on the combination of a Hankel transform and a Fourier transform. For physical problems that have close-to-cylindrical symmetry, this algorithm can be much faster than full 3D PIC algorithms. In addition, unlike standard finite-difference PIC codes, the proposed algorithm is free of numerical dispersion. This algorithm is benchmarked in several situations that are of interest for laser-plasma interactions. These benchmarks show that it avoids a number of numerical artifacts, that would otherwise affect the physics in a standard PIC algorithm - including the zero-order numerical Cherenkov effect.

연구 동기 및 목표

• 근원통대칭 시스템에 대해 전체 3D 카르테시안 PIC 시뮬레이션의 높은 계산 비용을 해결하기 위해.
• 표준 유한차분 PIC 코드에서 주요 오차 원인이 되는 진공 내 비물리적 수치 분산을 제거하기 위해, 특히 상대론적 입자 및 레이저 펄스에 대해.
• 준원통좌표계 알고리즘의 효율성과 스펙트럼(PSATD 유사) 방법의 분산 없음 정확성을 통합한 프레임워크를 구축하기 위해.
• 수치적 오염(예: 수치 체렌코프 효과) 없이 레이저 웨이브필드 가속 및 기타 상대론적 플라즈마 현상의 고정밀 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 라디얼 성분에 대해 헤이켈 변환, 축방향 성분에 대해 푸리에 변환을 사용하는 스펙트럼적 PIC 알고리즘을 수립한다.
• 미분 연산자의 유한차분 근사 없이 스펙트럼 영역에서 맥스웰 방정식을 해석하는 의사-스펙트럼적 접근을 적용한다.
• 실공간 격점과 라디얼 스펙트럼 모드 간의 변환을 위해 이산 헤이켈 변환(DHT) 및 역이산 헤이켈 변환(IDHT)을 사용한다.
• 특히 비영인 원위상 모드에 대해 정확한 재구성과 직교성을 보장하기 위해 변환 행렬에 제약 조건을 부과한다.
• DHT 및 IDHT 행렬을 사전 계산하고 역행렬화하여 정확성과 안정성을 유지하는 행렬 기반 형식을 채택한다.
• 특히 상대론적 속도에서도 비물리적 위상 속도 오차 없이 필드 방정식을 해석하는 방식으로 적용한다.

실험 결과

연구 질문

• RQ1 진공 내에서 수치 분산이 없는 동시에 준원통좌표계 기반의 스펙트럼적 PIC 알고리즘을 구축할 수 있는가?
• RQ2 헤이켈 및 푸리에 변환과 같은 스펙트럼 방법을 어떻게 PIC 입자 이동기와 융합하여 원통대칭 시스템에서 정확성을 유지할 수 있는가?
• RQ3 제안된 알고리즘이 표준 유한차분 PIC 코드에 비해 제로 오더 수치 체렌코프 효과를 어느 정도 억제하는가?
• RQ4 비원위상 대칭 효과를 고려한 전체 3D 정확성을 유지하면서도 2D 시뮬레이션 수준의 계산 효율성을 달성할 수 있는가?
• RQ5 비영인 원위상 모드에 대해 이산 헤이켈 변환을 구현할 때 발생하는 수치적 과제는 무엇이며, 이를 어떻게 완화할 수 있는가?

주요 결과

• 해석적 해와의 벤치마킹 결과, 진공 내에서 비물리적 수치 분산이 완전히 제거됨을 확인하였다.
• 표준 PIC 시뮬레이션에서 잘 알려진 문제인 제로 오더 수치 체렌코프 효과를 성공적으로 억제하였다.
• 레이저 웨이브필드 가속 시나리오의 경우, 표준 유한차분 PIC 코드가 수치 오차로 인해 잘못된 탈진 길이를 예측하는 데 반해, 본 알고리즘은 정확한 탈진 길이를 예측하였다.
• 적절한 행렬 제약 조건을 적용한 이산 헤이켈 변환을 통해 라디얼 모드 전반에 걸쳐 정확하고 안정적인 필드 재구성 가능함을 입증하였다.
• 근원통대칭 시스템에 대해 전체 3D 카르테시안 PIC 시뮬레이션 대비 상당한 속도 향상을 달성하였다.
• empirical 제약 조건을 통해 $n=0$ 및 $n\neq0$ 원위상 모드 모두에서 수치적 안정성과 정확성이 보장되며, 알고리즘이 강건함을 입증하였다.