Efficient Krylov solvers for inverse source problem in 2D space-time fractional diffusion equation

기술 요약: 2차원 시공간 분수 확산 방정식의 역원천 문제를 위한 효율적인 크릴로프 솔버(Krylov Solvers)

문제 정의
본 연구는 가변 계수를 갖는 2차원 시공간 분수 확산 방정식(FDE)과 관련된 역원천 문제(inverse source problem)를 다룬다. 지배 방정식은 시간 측면에서의 템퍼드 분수 미분($T D^{\xi,\rho}_t$)과 두 개의 공간 차원에서의 카푸토(Caputo) 분수 미분($CD^{\eta_1}_x, CD^{\eta_2}_y$)을 포함한다. 목표는 초기 조건과 노이즈가 포함된 최종 시간 데이터를 이용하여 미지 원천 항 $f(x, y)$를 재구성하는 것이다.

이 문제는 하다마르(Hadamard) 의미에서 본질적으로 부적정(ill-posed)하다. 즉, 입력 데이터의 미세한 섭동(노이즈 수준 $\epsilon$)이 해의 상당한 오차를 유발할 수 있다. 이를 완화하기 위해 저자들은 준경계값 방법(quasi-boundary value method)을 채택하였는데, 이는 비국소적 정규화 기법으로 부적정한 문제를 매개변수 의존적인 적정(well-posed) 정규화 문제로 변환한다.

방법론
제안된 접근 방식은 수치 이산화, 일반화된 국소 토플리츠(Generalized Locally Toeplitz, GLT) 이론을 통한 스펙트럼 분석, 그리고 효율적인 반복 솔버를 위한 프리컨디셔너(preconditioner) 설계를 결려한다.

1. 이산화:

   • 정규화된 문제는 유한 차분 스킴을 사용하여 이산화된다. 템퍼드 분수 시간 미분은 L1 공식을 사용하여 근사화되며, 공간 카푸토 미분 또한 L1 근사를 통해 이산화된다.
   • 이 과정은 $\tilde{N} = (m+1)\nu(n)$인 대규모 "올-앳-원스(all-at-once)" 선형 시스템 $A_{\tilde{N}} v_\lambda = z_\lambda$를 생성한다 (여기서 $\tilde{N}$은 시간 단계 $m$과 공간 격자점 $\nu(n) = n_1 n_2$를 결합한 총 자유도이다).
   • 결과적인 계수 행렬 $A_{\tilde{N}}$은 복잡한 다층 블록 구조를 나타낸다. 가변 계수 $a(x, y)$로 인해 이 행렬은 엄격한 토플리츠 구조는 아니지만, 대각 샘플링 행렬에 의해 섭동된 다층 블록 토플리츠 유사 구조를 가진다.

2. 스펙트럼 분석 (GLT 이론):

   • 저자들은 GLT 이론을 사용하여 행렬 시퀀스의 점근적 스펙트럼 분석을 수행한다. 이 프레임워크를 통해 이산화된 미분 연산자로부터 발생하는 대규모 행렬 시퀀스의 특이값 및 고윳값 분포를 규명할 수 있다.
   • 분석을 통해 계수 행렬 시퀀스의 **GLT 심볼(symbol)**을 식별한다. 심볼은 시간 단계 블록(템퍼드 미분과 관련된)과 공간 연산자 블록(가변 계수 및 공간 분수 미분을 포함하는)의 심볼의 합으로 도출된다.
   • 분석은 행렬 시퀀스가 이 심볼에 따라 분포됨을 확인하며, 이는 프리컨디셔닝의 이론적 토대를 제공한다.

3. 프리컨디셔닝 전략:

   • GLT 분석을 바탕으로, 저자들은 블록 삼각 프리컨디셔너($P_{\tilde{N}}$)를 제안한다. 이 프리컨디셔너는 블록 하삼각 구조를 유지하면서 대각 블록을 단순화하여 원래의 행렬 $A_{\tilde{N}}$을 근사한다.
   • 이론적 결과(정리 5.2)는 프리컨디셔닝된 행렬 시퀀스 $\{P_{\tilde{N}}^{-1} A_{\tilde{N}}\}$가 GLT 의미에서 항등 함수(심볼 = 1)로 분포됨을 입증한다. 결과적으로, 프리컨디셔닝된 시스템의 고윳값들은 1 주변에 **약하게 클러스터링(weakly clustered)**되는데, 이는 GMRES와 같은 크릴로프 부공간 방법의 수렴에 바람직한 특성이다.
   • 프리컨디셔너의 대각 블록을 포함하는 선형 시스템의 해는 표준적인 기술(예: 순환 근사 또는 가변 계수의 적분 평균 사용)을 통해 처리된다.

주요 기여

• 정식화: 본 논문은 가변 계수를 갖는 2D 시공간 FDE에 대한 역원천 문제를 정식화하고, 노이즈에 대한 안정성을 보장하기 위해 준경계값 방법을 적용한다.
• 구조적 분석: GLT 이론을 사용하여 결과적인 "올-앳-원스" 이산화 행렬에 대한 엄밀한 스펙트럼 분석을 수행하고, 다층 블록 구조에 대한 심볼을 명시적으로 도출한다.
• 프리컨디셔너 설계: 이산화의 다층 구조를 보존하는 새로운 블록 삼로 프리컨디셔너를 구축한다. 저자들은 이 프리컨디셔너가 1 주변의 고윳값 클러스터링을 유도하여 이론적으로 효율적인 수렴을 보장함을 증명한다.
• 계산 효율성: 제안된 방법은 행렬-벡터 곱 및 프리컨디셔너 적용의 계산 비용을 $O(m^2 \nu(n))$ (토플리츠 해법을 위한 로그 인자 포함)으로 줄이는 행렬 구조를 활용한다.

수치 결과
수치 실험은 영역 $\Omega = (0, 1)^2$에서 다양한 분수 차수($\xi, \eta_1, \eta_2$)와 격자 크기에 대해 수행되었다.

• 수렴성: 프리컨디셔닝된 GMRES 방법은 프리컨디셔너가 없는 버전을 크게 능가했다. 예를 들어, 격자 크기가 $n_1=n_2=m=25$일 때, 특정 분수 매개변수에 대해 반복 횟수가 프리컨디셔너 미적용 시 186회에서 적용 시 39회로 감소하여 약 80%의 감소를 보였다.
• 강건성: 프리컨디셔닝된 시스템의 반복 횟수는 시간 이산화($m$)가 증가함에 따라 거의 일정하게 유지되었으며, 공간 이산화($n_1, n_2$)가 증가함에 따라 완만하게만 증가했다. 반면, 프리컨디셔너가 없는 시스템은 격자 세분화에 따라 반복 횟수가 급격히 증가했다.
• 비용: 프리컨디셔너 적용에 따른 오버헤드가 발생하지만, 수렴을 위해 필요한 GMRES 반복 횟수의 극적인 감소 덕분에 전체 CPU 시간은 실질적으로 단축되었다.

의의 및 주장
본 논문은 제안된 프레임워크가 분수 확산 모델의 역원천 문제를 해결하기 위한 효율적이고 확장 가능한 접근 방식을 제공한다고 주장한다. 정규화, 구조 보존 이산화, GLT 기반 스펙트럼 분석을 통합함으로써, 이 방법은 GMRES 솔버를 성공적으로 가속화한다.

저자들은 프리컨디셔너가 분수 문제와 관련된 격자 세분화 시 발생하는 스펙트럼 특성의 악화를 효과적으로 완화한다고 강조한다. 그들은 클러스터링이 (다층 문제에서 일반적이듯) "약한" 수준이지만, 결과적인 수렴 동작이 다양한 분수 차수와 이산화 수준에 대해 강건하다고 언급한다. 이 연구는 기존의 1차원 연구를 더 복잡한 2차원 사례로 확장하여, 고급 행렬 구조의 활용이 비국소적 분수 연산자와 부적정 역문제에 의한 계산적 난관을 극복할 수 있음을 보여준다.

저자들은 향-후 연구로서 본 방법론을 더 복잡한 기하학적 구조, 고차 유한 요소, 등각 이소제메트릭 분석(isogeometric analysis) 및 고차원 공간으로 확장할 수 있음을 제시하며, 제시된 이론적 프레임워크가 이러한 시나리오에도 적응 가능하다는 점을 명시하였다.