[논문 리뷰] On Uncertainty Quantification in Particle Accelerators Modelling

Andreas Adelmann|arXiv (Cornell University)|2015. 09. 27.

Probabilistic and Robust Engineering Design참고 문헌 28인용 수 3

한 줄 요약

이 논문은 입자 가속기 모델링을 위한 다항 혼합 확장(PCE) 기반 비침습적 불확실성 정량화(UQ) 프레임워크를 소개한다. 이는 대체 모델 구축, 소볼' 지수를 통한 국소 감도 분석, 오차 전파를 가능하게 하며, 낮은 수의 고정밀도 시뮬레이션으로도 효율적으로 비드 에미ittance, 에너지 분포, 홀로 형성에 영향을 주는 주요 불확실 매개변수를 특정한다.

ABSTRACT

Using a cyclotron based model problem, we demonstrate for the first time the applicability and usefulness of a uncertainty quantification (UQ) approach in order to construct surrogate models for quantities such as emittance, energy spread but also the halo parameter, and construct a global sensitivity analysis together with error propagation and $L_{2}$ error analysis. The model problem is selected in a way that it represents a template for general high intensity particle accelerator modelling tasks. The presented physics problem has to be seen as hypothetical, with the aim to demonstrate the usefulness and applicability of the presented UQ approach and not solving a particulate problem. The proposed UQ approach is based on sparse polynomial chaos expansions and relies on a small number of high fidelity particle accelerator simulations. Within this UQ framework, the identification of most important uncertainty sources is achieved by performing a global sensitivity analysis via computing the so-called Sobols' indices.

연구 동기 및 목표

• 비침습적 불확실성 정량화(UQ)가 고강도 입자 가속기 모델링에 적용될 수 있음을 시험하기 위해 사이클로트론 기반 모델 문제를 사용한다.
• 다항 혼합 확장(PCE)을 이용해 핵심 비드 품질 지표인 에미ittance, 에너지 분포, 홀로 매개변수에 대한 대체 모델을 구축한다.
• 소볼' 지수를 통한 국소 감도 분석을 수행하여 가속기 시뮬레이션에서 가장 영향력 있는 불확실한 입력 매개변수를 규명한다.
• 최소한의 고정밀도 시뮬레이션을 사용하여 비드 역학 출력의 오차 전파 및 오차 분석을 가능하게 한다.
• 실제 설계 프로젝트인 PSI의 소형 고강도 사이클로트론을 대상으로 프레임워크를 검증한다.

제안 방법

• 기존의 비드 역학 코드(예: OPAL)를 블랙박스로 간주하는 비침습적 PCE 방법을 사용하며, 독립적인 시뮬레이션 실행만 필요로 한다.
• 불확실한 입력 매개변수를 랜덤 변수로 표현하고, 샘플링 기반 방법을 사용하여 관심 있는 출력 양(관심 지표, QoI)을 다항 혼합 기저로 전개한다.
• 각 입력 매개변수의 출력 분산에 미치는 기여도를 정량화하기 위해 첫 번째 순서($S_k$) 및 총 효과($S_k^T$) 소볼' 지수를 활용한다. 이는 국소 감도 분석을 가능하게 한다.
• 고정밀도 시뮬레이션 세트를 사용하여 최소 제곱법 기반의 콜로케이션을 통해 PCE 계수를 계산하며, 시뮬레이션 코드에 대한 침습적 수정을 피한다.
• 균일 분포를 가진 입력에 대해 레전드르 다항식을 사용하고, 입력 확률 분포에 맞추기 위해 위erner-아스키 다항 혼합 확장 프레임워크를 적용한다.
• 주어진 PSI 인젝터 2 유사 사이클로트론 모델의 10턴 중심 영역에 이 방법을 적용하며, 관심 지표로는 비드 품질 지표를 사용한다.

실험 결과

연구 질문

• RQ1 비침습적 다항 혼합 기반 UQ는 고강도 입자 가속기 시뮬레이션에 효과적으로 적용될 수 있는가? 정확한 대체 모델을 구축할 수 있는가?
• RQ2 사이클로트론 설계에서 비드 에미ittance, 에너지 분포, 홀로 형성에 가장 크게 영향을 주는 입력 매개변수는 무엇인가?
• RQ3 소볼' 지수를 통한 국소 감도 분석을 어떻게 활용하여 가속기 모델링에서 가장 중요한 불확실한 매개변수를 식별하고 순위를 매길 수 있는가?
• RQ4 PCE를 통한 대체 모델은 불확실성 전파에서 정확도를 유지하면서 계산 비용을 얼마나 줄일 수 있는가?
• RQ5 이 UQ 프레임워크는 가속기 설계 프로젝트에 통합되어 내구성 및 오차 분석을 지원할 수 있는가?

주요 결과

• 제안된 비침습적 PCE 기반 UQ 프레임워크는 낮은 수의 고정밀도 시뮬레이션으로도 비드 품질 지표(에미ittance, 에너지 분포, 홀로)의 정확한 대체 모델링을 가능하게 한다.
• 소볼' 지수를 통한 국소 감도 분석은 비드 역학에 영향을 주는 가장 영향력 있는 매개변수를 성공적으로 규명하였으며, $S_k^T$ 값은 일부 비드 분포 모멘트가 출력 변동성에 중대한 영향을 준다는 것을 시사한다.
• $S_k^T \text{<<} 1$인 매개변수들은 무시할 수 있으며, 평균값으로 고정함으로써 정확도 손실 없이 모델 단순화가 가능하다.
• 이 방법은 효율적인 오차 전파 및 분산 분해를 보여주어 입력 매개변수의 변동성에 기인한 출력 불확실성의 정량화를 가능하게 하였다.
• 실제로 PSI의 소형 고강도 사이클로트론 설계 프로젝트에 성공적으로 적용되어, 가속기 연구 및 개발에서의 실용적 유용성을 입증하였다.
• 이 방법은 각 시뮬레이션 실행이 상호 독립적이므로, 매우 쉽게 병렬 처리가 가능하며, HPC 자원에서 확장성이 뛰어나다.