Multigrid Preconditioning for FEEC using Mass-Lumping and Transforming Smoothers

기술 요약: 질량 럼핑(Mass-Lumping) 및 변환 스무더(Transforming Smoothers)를 이용한 FEEC용 멀티그리드 프리컨디셔닝

문제 정의
본 논문은 드 람 복합체(de Rham complex) 상에 정의된 편미분 방정식(PDE)의 유한 요소 외사적 계산(FEEC) 이산화로부터 발생하는 선형 시스템을 해결하는 과제를 다룹다. 구체적으로, 호지-디락 연산자(Hodge-Dirac operator), 호지-라플라시안(Hodge-Laplacian), 그리고 마그네토스태틱 시스템(코울롱 게이지에서의 맥스웰 방정식)에 대한 혼합 및 사들 포인트(saddle-point) 정식화(formulations)를 대상으로 한다. 이러한 이산화는 일반적으로 부정 definite(indefinite) 성질을 가지며 사들 포인트 구조를 갖는 선형 시스템을 생성한다. 가우스-사이델(Gauss-Seidel)이나 자코비(Jacobi)와 같은 표준 반복 스무더는 이러한 부정 definite 시스템에 효과적이지 않은데, 이는 해당 연산자들이 릴랙세이션 스킴에 필요한 양의 확정성(positive definiteness)을 결여하고 있기 때문이다. 기존 연구들이 슈어 보완(Schur complement)에 기반한 멀티그리드 접근 방식을 제안해 왔으나, 본 논문은 슈어 보완 역행렬 계산의 복잡성을 피하는 대안적인 전략을 모색한다.

방법론
제안된 프레임워크는 일관된(consistent) FEEC 시스템을 위한 효과적인 프리컨디셔너를 구축하기 위해 **질량 럼핑(mass-lumping)**과 **변환 스무더(transforming smoothers)**를 결합한다.

1. 질량 럼핑 및 스펙트럼 등가성(Spectral Equivalence):
   저자들은 일관된 질량 행렬이 대각 행렬(또는 명시적으로 역행렬 계산이 가능한 행렬)로 대체된 질량 럼핑 버전의 FEEC 연산자를 도입한다. 이를 통해 밀집된(dense) 일관된 질량 행렬을 사용하면 비용이 매우 많이 드는 리츠 동형 사상(Riesz isomorphism)의 명시적 역행렬 계산을 가능하게 한다.

• 일관된 내적과 질량 럼핑된 내적 사이의 완만한 $h$-균등 노름-등가성(norm-equivalence) 가정 하에, 저자들은 질량 럼핑된 노름에 대한 이산 포앵카레 부등식(discrete Poincaré inequality)을 증명한다.
• 이 안정성 결과는 일관된 FEEC 연산자와 질량 럼핑된 연산자 사이의 스펙트럼 등가성에 대한 증명으로 이어진다. 이 등가성은 질량 럼핑된 시스템을 원래 시스템의 프리컨디셔닝을 위한 대리물(proxy)로 사용하는 것을 정당화한다.

2. 변환 스무더:
   질량 럼핑된 연산자의 부정 definite 성질을 처리하기 위해, 방법론은 변환 스무더를 채택한다. 핵심 아이디어는 연산자를 대각 블록이 양의 확정적인 형태의 블록 형태로 대수적으로 변환하여, 가우스-사이델과 같은 표준 스무더를 사용할 수 있게 하는 것이다.

• 호지-디락, 호지-라플라시안, 그리고 마그네토스태틱 문제를 위해 특정 좌측 및 우측 변환 행렬($T_L, T_R$)이 구성된다.
• 이러한 변환은 원래의 부정 definite 시스템을 하삼각 블록 시스템(또는 양의 확정 대각 블록을 가진 시스템)으로 매핑한다.
• 가우스-사이델 릴랙세이션이 이 변환된 시스템에 적용되며, 그 결과는 다시 원래의 공간으로 매핑된다.

3. 멀티그리드 프리컨디셔닝:
   멀티그리드 알고리즘은 질량 럼핑된, 변환된 시스템에 특화되어 설계되었다. 이 멀티그리드 사이클($M_{MG}$)은 질량 럼핑된 연산자의 근사 역행렬 역할을 한다.

• 원래의 FEEC 시스템 $G_h u = g$를 풀기 위해, 프리컨디셔너 $B$는 $B \approx J_h^{-1} M_{MG} \tilde{H}_h (J_h^*)^{-1}$로 구성된다 (여기서 $J_h$는 표준 공간과 질량 럼핑된 공간 사이의 동형 사상이며, 실제로는 항등 행렬인 경우가 많고, $\tilde{H}_h$는 질량 럼핑된 리츠 맵이다).
• 결과적으로 생성된 프리컨디셔닝 시스템이 여전히 부정 definite이므로, 저자들은 PCG나 MINRES와 같은 표준 방법 대신 GMRES(Generalized Minimal Residual method)를 외부 솔버로 활용한다.

주요 기여

• 알고리즘 설계: 본 논문은 FEEC 문제를 위해 질량 럼핑과 변환 스무더를 통합하는 통일된 멀티그리드 프리컨디셔닝 프레임워크를 제안한다.
• 이론적 안정성: 질량 럼핑된 시스템의 안정성을 확립하고, 완만한 가정 하에 원래의 FEEC 연산자와의 스펙트럼 등가성을 증명함으로써 프리컨디셔닝 전략의 이론적 토대를 제공한다.
• 부정 definite 처리: 변환 스무더가 어떻게 부정 definite 블록 연산자를 슈어 보완 근사 없이도 고전적인 릴랙세이션 기법에 적합한 형태로 변환할 수 있는지를 보여준다.
• 복잡한 위상에 대한 강건성: 본 방법은 베티 수(Betti numbers)에 의해 시스템의 랭크 결핍(rank deficiency)이 결정되는 비자명한 위상을 가진 도메인에서도 효과적으로 작동함을 보여준다.

수치 결과
재진입 코너(reentrant corners) 및 구멍(holes)이 있는 다양한 메쉬(비자명한 위상 포함)를 포함하여 2D 및 3D에서 광범 위한 수치 실험이 수행되었다.

• 수렴성: 변환 스무더와 결합된 멀티그리드 사이클(V-cycle 및 W-cycle)은 견고한 수렴성을 보였다. 오차 전파 연산자의 스펙트럴 반경(spectral radii)은 메쉬 크기 $h$와 무관하게 1에서 떨어진 상태로 유지되었다.
• 질량 럼핑 전략: 세 가지 질량 럼핑 전략(scaled identity, row sums, barycentric duals)이 테스트되었다. 결과는 이러한 서로 다른 전략들이 수렴 거동에 있어 뚜렷한 질적 차이를 유발하지 않음을 나타냈다.
• 솔버 성능: GMRES의 프리컨디셔너로 사용되었을 때, 이 방법은 메쉬 세분화 수준과 무관한 반복 횟수를 산출했다. W-cycle이 V-cycle보다 일반적으로 더 견고한 수렴성을 제공했다.
• 조화 형식(Harmonic Forms): 실험을 통해 본 방법이 조화 형식과 관련된 랭크 결핍을 올바르게 처리함을 확인하였다. 즉, 수렴하지 않는 고유 모드의 수가 조화 공간의 차원과 일치하였다.

의의 및 주장
저자들은 본 연구를 공식적인 수렴 이론 증명이라기보다는 알고리즘적 기여로 규정한다. 저자들은 스펙트럼 등가성과 안정성을 증명했음에도 불구하고, GMRES 반복 횟수에 대한 균등한 경계(uniform bounds)를 설정하는 엄격한 수렴 이론은 향후 연구 과제로 남겨두었다고 명시적으로 밝히고 있다.

이 접근 방식의 의의는 명시적이고 저렴한 연산(질량 럼핑)과 대수적 변환을 사용하여 부정 definite FEEC 시스템을 위한 실용적이고 견고한 프리컨디셔너를 제공할 수 있다는 점에 있다. 수치적 증거는 이 접근 방식이 전자기학 및 관련 분야의 혼합 유한 요소 문제를 해결하는 데 유효하고 유망한 방향임을 시사하며, 향후 스톡스 흐름(Stokes flow)이나 계수 점프(coefficient jumps)가 있는 시스템으로의 이론적 연구 및 확장을 필요로 한다.