[논문 리뷰] Optimized Wakefield Computations Using A Network Model

J.-F. Ostiguy, K.Y. Ng|ArXiv.org|2000. 08. 18.

Electromagnetic Simulation and Numerical Methods인용 수 8

한 줄 요약

이 논문은 선형 충돌기에서 사용되는 감쇠된 딜레이팅 구조(DDS)의 웨이크필드를 계산하기 위한 고도로 최적화된 네트워크 모델을 제시한다. 회로 모델의 계수 행렬의 희소성 특성을 활용하기 위해 마크위츠 재정렬과 희소 선형 해법기를 적용함으로써 저자들은 세 자리 수의 속도 향상을 달성하였으며, 이는 각 구조당 웨이크필드 계산 시간을 수 시간에서 약 14초로 단축시켰다. 이는 처음으로 전체 린액 시뮬레이션을 가능하게 한다.

ABSTRACT

During the course of the last decade, traveling wave accelerating structures for a future Linear Collider have been the object of intense R&D efforts. An important problem is the efficient computation of the long range wakefield with the ability to include small alignment and tuning errors. To that end, SLAC has developed an RF circuit model with a demonstrated ability to reproduce experimentally measured wakefields. The wakefield computation involves the repeated solution of a deterministic system of equations over a range of frequencies. By taking maximum advantage of the sparsity of the equations, we have achieved significant performance improvements. These improvements make it practical to consider simulations involving an entire linac ($\sim 10^3$ structures). One might also contemplate assessing, in real time, the impact of fabrication errors on the wakefield as an integral part of quality control.

연구 동기 및 목표

• 선형 충돌기 가속 구조에서 제작 오차로 인한 비드 범위 품질 열화를 효율적이고 대규모로 시뮬레이션할 수 있도록 하는 것.
• 수천 개의 셀을 포함한 구조를 시뮬레이션할 때 기존 유한 요소법 또는 유한 차분법의 계산 병목 현상을 극복하는 것.
• 고정밀도 실험 측정치와의 일치를 확보하면서도 확장 가능하고 정확한 장거리 웨이크필드 계산 방법을 개발하는 것.
• 가속기 연구 및 개발의 품질 관리에서 제작 오차가 웨이크필드 성능에 미치는 영향을 실시간 평가하는 데 실용적인 방법을 제공하는 것.

제안 방법

• 저자들은 맥스웰 방정식을 저차수 모드 전개를 통해 이산화한 주파수 영역 RF 회로 모델을 사용하여, 시스템을 희소한 복소 대칭 선형 연립방정식으로 표현한다.
• 모드 전압과 전류를 포함하는 블록 구조 행렬 방정식 (1)으로 시스템을 수립하며, TM 및 TE 모드 간의 결합과 전송선을 통한 맨틀 데드닝을 포함한다.
• 임피던스의 허수부로부터 유도된 스펙트럼 함수 S(ω)는 다수의 주파수에서 선형 연립방정식을 반복적으로 해석함으로써 계산된다.
• 희소 선형 시스템을 효율적으로 해석하기 위해 Quasi Minimal Residual (QMR) 반복 해법기를 사용하며, 이전 주파수 단계에서의 좋은 초기 추정치를 활용한다.
• 요소 분해 과정 중에 불필요한 채움을 최소화하기 위해 초기 단계에서 한 번만 마크위츠 재정렬을 적용하여 희소성을 유지하고 계산 비용을 극적으로 감소시킨다.
• 행렬 역행렬을 형성하는 대신 해석 후 루프 전류를 구한 후 직접 행렬 곱셈을 수행함으로써 성능을 더욱 향상시킨다.

실험 결과

연구 질문

• RQ1 네트워크 기반의 회로 모델을 최적화하여 감쇠된 딜레이팅 구조(DDS)의 웨이크필드를 충분한 정확도와 함께 세 자리 수의 속도 향상을 달성할 수 있는가?
• RQ2 마크위츠 제거와 같은 희소성 및 재정렬 전략이 가속기 물리학에서 주파수에 따라 변화하는 큰 크기의 선형 연립방정식을 푸는 데에 계산 비용을 얼마나 줄일 수 있는가?
• RQ3 지금은 실용적인 시간(4시간 이내)으로 계산 가능한 웨이크필드 계산을 통해 전체 린액(1000개 이상의 구조)을 시뮬레이션하는 것이 가능한가?
• RQ4 주파수에 따라 변화하는 행렬의 특성에도 불구하고, 최적화된 코드가 전체 딜레프 모드 주파수 대역에서 수치적 안정성과 정확성을 유지할 수 있는가?

주요 결과

• 최적화된 코드는 550 MHz 펜티엄 III 프로세서에서 단일 RDDS 구조의 웨이크필드를 약 14초 내에 계산한다.
• 이것은 이전에 보고된 계산 시간에 비해 세 자리 수의 성능 향상이며, 전체 린액 시뮬레이션을 4시간 이내로 수행할 수 있게 하였다.
• 스펙트럼 함수와 웨이크필드 출력 결과는 이전에 더 느린 방법을 사용한 SLAC 팀이 생성한 결과와 뛰어난 일치를 보였다.
• 마크위츠 재정렬 전략은 전체 주파수 대역에서 수치적으로 안정적이며, 한 번 계산하여 모든 주파수 단계에 재사용할 수 있다.
• 맨틀 전압 변수 A를 제거하지 않고 유지함으로써 시스템의 희소성을 유지함으로써 불필요한 채움을 방지하고 계산 효율성을 유지한다.
• SPARSE 라이브러리를 통해 구현된 연결 리스트 기반의 희소 행렬 표현 방식은 효율적인 순차적 업데이트와 포인터 기반 접근을 가능하게 하여 고성능 계산을 지원한다.