Towards joint optimization of stellarator coils and support structures

기술 요약: 스텔라레이터 코일 및 지지 구조의 공동 최적화를 향하여

문제 정의
스텔라레이터 핵융합로 설계는 플라즈마 가둠을 위한 정확한 자기장이 필요한 요구사항과, 거대한 로런츠 힘 및 열 응력으로부터 취성이 있는 고온 초전도(HTS) 도체를 보호해야 하는 기계적 요구사항 사이의 균형을 맞추어야 하는 중대한 과제에 직면해 있다. 현재의 설계 프로세스는 순차적이며 노동 집약적이다. 코일 기하학적 구조는 자기 성능 측면에서 단독으로 최적화되며, 그 후 공학적 응력 한계를 충족하기 위해 유한 요소 분석(FEA)을 통한 수동 반복 과정을 거쳐 지지 구조가 추가된다. 이러한 "후처리" 방식은 지지 구조의 위치가 코일 응력에 미치는 영향과, 하중 하에서의 코일 변형이 자기장 오차를 유발하는 방식을 고려하지 못한다. 기존의 최적화 도구들은 지지 구조를 무시하는 단순화된 해석적 대리 모델(proxy)에 의존하거나, 기울기 기반(gradient-based) 코일 최적화 루프와 인터페이스할 수 없는 폐쇄형 상용 FEA 도구를 사용한다.

방법론
저자들은 스텔라레이터 코일 최적화 루프에 통합된 오픈 소스 미분 가능 구조 역학 도구인 coil-fem을 제시한다. JAX-FEM 라이 라이브러리를 기반으로 구축된 coil-fem은 기울기 기반 방법을 사용하여 코일 기하학적 구조와 지지 구조 위치를 동시에 최적화할 수 있게 한다.

주요 방법론적 구성 요소는 다음과 같다:

• 미분 가능한 FEA: 이 도구는 Robin(스프링 기초) 경계 조건을 사용하여 선형 탄성 방정식 $\nabla \cdot \sigma(u) + F_{body} = 0$을 해결한다. 이 근사는 고정된 디리클레(Dirichlet) 경계 조건을 대체하며, 이를 통해 지지 클램프 위치를 FEA 시스템의 희소 행렬 구조를 변경하지 않고도 최적화 가능한 변수로 취급할 수 있게 한다.
• 물리 모델: 이 모델은 중력, 상호 자기력(1D 필라멘트로 모델링됨), 자가 힘(직사각형 단면에 대한 Landreman-Hurwitz 공식을 통해 계산됨)을 포함한 체적력을 고려한다. 또한 응력이 없는 상태에서 운전 상태까지의 약 200K 온도 강하로 인한 열 변형을 포함한다.
• 메싱 및 기하학: 코일은 회전 최소 프레임(rotation-minimizing frame)을 사용하여 중심선을 따라 직사각형 단면을 스윕(sweep)하여 생성된 유한 부피체로 모델링된다. 최적화 과정 중 미분 가능성과 계산 효율성을 보장하기 위해 메쉬 토폴로지는 고정된다.
• 최적화 루프: 저자들은 자기장 오차, 기계적 하중(폰 미세스 응력의 제곱의 체적 적분), 그리고 표준 기하학적 제약 조건(코일 길이, 간격, 곡률)을 포함하는 다목적 페널티 함수를 구성한다. 코일 푸리에 계수와 지지 클레의 파라메트릭 위치 모두에 대해 이 목적 함수를 최소화하기 위해 L-BFGS 알고리즘이 사용된다.

주요 기여

1. 최초의 결합 최적화: 본 논문은 스텔라레이터 코일 기하학적 구조와 지지 구조 위치의 기울기 기반 공동 최적화를 다룬 문헌 중 최초의 사례를 제시한다.
2. coil-fem 도구: 지지 미분 가능한 FEA를 코일 최적화 루프에 직접 통합하여, 지지 클램프 위치를 최적화 가능한 변수로 취급하는 오픈 소스 도구를 공개한다.
3. 지지의 중요성 입증: 본 연구는 지지 구조가 고정되어 있거나 수동으로 배치될 때 흔히 간과되는 요소인, 기계적 응력을 줄이는 데 있어 지지 구조의 배치가 결정적인 요인이라는 수치적 증거를 제공한다.

결과
저자들은 coil-fem을 DOLFINx FEA 제품군 및 상용 스타일의 대리 모델 접근법과 벤치마킹하여, 응력 및 변위 지표에서 1~5% 이내의 일치도를 입증하였다. 이후 이 도구를 단순화된 W7-X 코일 세트에 네 가지 시나리오에 대해 적용하였다:

• Baseline (기준점): 고정된 상/하단 클램프를 가진 원래의 W7-X 코일.
• Case A (최적화된 지지 구조): 고정된 코일 기하학, 최적화된 클램프 위치.
• Case B (최적화된 코일): 최적화된 코일 기하학, 고정된 클램프 위치.
• Case C (공동 최적화): 코일 기하학 및 클램프 위치 모두 최적화.
• Case D (힘 최소화): 고정된 클램프를 가진 상태에서 로런츠 힘(표준 대리 모델)을 최소화하도록 최적화된 코일 기하학.

주요 결과:

• 응력 감소: Case C는 기준점 대비 RMS 폰 미세스 응력을 2.4배 감소시켰다.
• 지지 구조의 역할: 고정된 지지 구조와 함께 코일 기하학만을 최적화했을 경우(Case B) 응력 개선은 미미했다. 반대로, 지지 위치만을 최적화했을 경우(Case A) 응력이 크게 감소하였으며, 이는 지지 구조 배치가 구조 역학에서 지배적인 요인임을 확인시켜 준다.
• 대리 모델의 한계: Case D는 로런츠 힘을 최소화하는 것(흔히 쓰이는 대리 모델)이 반드시 폰 미세스 응력의 최소화와 상관관계가 있지는 않음을 보여주었으며, 이는 통합된 구조 FEA의 필요성을 정당화한다.
• 열 효과: 로런츠 힘이 지배적인 응력원이지만, 열 수축이 최대 응력에 상당한 영향(피크 위치를 약 17% 이동시킴)을 미친다는 점은 설계 루프에 열 하중을 포함하는 것이 필수적임을 강조한다.

의의 및 주장
본 논문은 coil-fem이 지지 구조를 배치하는 수동적인 시행착오 과정을 제거함으로써 스텔라레이터 설계의 반복 속도를 크게 가속화한다고 주장한다. 자기장 오차와 기계적 응력을 동시에 줄임으로써, 이 도구는 더 경제적이고 견고한 핵융합로 설계를 위한 경로를 제공한다. 저자들은 이 작업이 고차원적 열-탄성 모델 및 복잡한 지지 구조(예: 빔 네트워크 또는 쉘)로의 확장이 가능한 "개념 증명(proof of concept)"임을 명시하며, 미분 가능한 구조 역학이 핵융합 코일 설계에서 실현 가능하다는 것을 입증하였다. 저자들은 현재 모델이 단순화된 가정(예: 등방성 재료, 스프링으로 모델링된 케이지형 지지 구조)을 사용하고 있지만, 이 프레임워크가 향러 고충실도 모델로 확장될 수 있음을 밝히며, 본 연구가 모든 스텔라레이터 기계적 문제를 해결했다고 주장하는 것이 아니라 코일과 그 지지 구조를 공동으로 최적화할 수 있는 최초의 계산 프레임워크를 제공하는 데 목적이 있음을 분명히 하였다.