[논문 리뷰] Adaptive Multi-Dimensional Particle In Cell

Giovanni Lapenta|ArXiv.org|2008. 06. 04.

Electrohydrodynamics and Fluid Dynamics인용 수 3

한 줄 요약

이 논문은 다차원 입자-격자(Particle-In-Cell, PIC) 방법을 적응형으로 개선하여 시간 및 길이 스케일이 광범위한 다스케일 플라즈마 현상을 효율적으로 시뮬레이션할 수 있도록 한다. 격자 적응과 입자 재조정을 통합함으로써 충돌 없는 충격파 및 먼지 전하화와 같은 관심 영역에서 정확도를 유지하면서 계산 비용을 줄였다. 비적응형 시뮬레이션에 비해 수치적 안정성과 통계 정확도 향상이 입증되었다.

ABSTRACT

Kinetic Particle In Cell (PIC) methods can extend greatly their range of applicability if implicit time differencing and spatial adaption are used to address the wide range of time and length scales typical of plasmas. For implicit differencing, we refer the reader to our recent summary of the implicit moment PIC method implemented in our CELESTE3D code [G. Lapenta, Phys. Plasmas, 13, 055904 (2006)]. Instead, the present document deals with the issue of PIC spatial adaptation. Adapting a kinetic PIC code requires two tasks: adapting the grid description of the fields and moments and adapting the particle description of the distribution function. Below we address both issues. First, we describe how grid adaptation can be guided by appropriate measures of the local accuracy of the solution. Based on such information, grid adaptation can be obtained by moving grid points from regions of lesser interest to regions of higher interest or by adding and removing points. We discuss both strategies. Second, we describe how to adapt the local number of particles to reach the required statistical variance in the description of the particle population. Finally two typical applications of adaptive PIC are shown: collisionless shocks and charging of small bodies immersed in a plasma.

연구 동기 및 목표

• 플라즈마에서 넓은 시간 및 길이 스케일을 다루는 데 있어 표준 명시적 PIC 방법의 계산 비효율성을 해결하기 위해.
• 국소적 기울기와 동적 입자 분포를 포착하는 데 한계가 있는 고정 격자, 고정 입자 PIC 시뮬레이션의 한계를 극복하기 위해.
• 암시적 PIC 방법에 격자 적응과 입자 재조정을 통합하여 정확도를 향상시키고 노이즈를 줄이는 프레임워크를 개발하기 위해.
• 충돌 없는 충격파 및 우주 환경에서의 먼지 전하화와 같은 복잡한 플라즈마 시스템을 정확하게 시뮬레이션할 수 있도록 하기 위해.
• 적응형 해상도와 입자 제어를 통해 병렬 구현에서 에너지 보존과 로드 밸런싱을 확보하기 위해.

제안 방법

• 지역 해의 정확도를 기반으로 격자 적응을 이끄는 사후 오차 지표를 사용하여 고해상도가 필요한 영역을 식별한다.
• 격자 적응 전략을 두 가지 적용: 격자 이동(고해상도 영역으로 점을 이동) 및 격자 정밀화(필요한 곳에 점을 추가).
• 통계 정확도를 유지하고 노이즈를 줄이기 위해 세포당 입자 수를 동적으로 조정함으로써 입자 재조정을 적용한다.
• 먼지 입자 경계 조건을 플라즈마 시뮬레이션에서 모델링하기 위해 유동 경계 방법에 감수성 및 감쇠 항목을 통합한다.
• 정확도와 계산 비용의 균형을 이루기 위해 이동 격자 방법(Moving Mesh Approach, MMA)과 적응형 격자 정밀화(Adaptive Mesh Refinement, AMR)를 구현한다.
• 안정성 제약을 완화하고 적절한 경우 더 큰 시간 간격과 더 흐린 격자를 허용하기 위해 암시적 시간 차분을 통합한다.

실험 결과

연구 질문

• RQ1 다차원 PIC 시뮬레이션에서 지역 정확도 측정치에 기반해 격자 적응을 어떻게 이끌 수 있는가?
• RQ2 적응형 PIC 방법에서 격자 이동과 격자 정밀화 중 어떤 전략이 더 효과적인가?
• RQ3 변화하는 입자 밀도 영역에서 입자 재조정은 어떻게 통계 정확도를 유지할 수 있는가?
• RQ4 적응형 격자 및 입자 제어는 병렬 PIC 코드에서 에너지 보존과 로드 밸런싱에 어떤 영향을 미치는가?
• RQ5 적응형 PIC는 충돌 없는 충격파 및 먼지 전하화와 같은 복잡한 플라즈마 현상의 시뮬레이션을 어떻게 향상시키는가?

주요 결과

• 사후 오차 지표에 기반한 격자 적응은 충격파 전면 및 셰일드와 같은 고기울기 영역에서 해상도를 크게 향상시킨다.
• 입자 재조정은 세포당 일관된 수의 입자를 유지함으로써 저밀도 영역에서 통계 노이즈가 증가하는 것을 방지한다.
• 적응형 격자와 입자 제어의 조합은 먼지 전하화 시뮬레이션을 안정적이고 정확하게 수행할 수 있게 하여 평형 상태에서 전하 축적을 달성한다.
• 입자 재조정 없이 시뮬레이션할 경우 먼지 주변의 세포당 입자 수가 시간이 지남에 따라 감소하여 정확도가 떨어지고 평형 상태에 수렴하지 못한다.
• 적응형 PIC 방법은 특히 강한 공간적 및 시간적 스케일 분리가 있는 시스템에서 고정 격자, 고정 입자 시뮬레이션에 비해 더 높은 계산 효율성과 정확도를 달성한다.