비선형 LTE 선 복사전달 문제 해결을 위한 사전조건 비공액 그라디언트 방법

초록

본 연구는 비국부열평형(non‑LTE) 조건에서 1차원 평면 평행 구성을 갖는 두 레벨 원자 모델의 선 복사전달 방정식을 풀기 위해, 기존의 ALI‑Jacobi, Gauss‑Seidel 및 SOR 방식보다 빠른 수렴을 보이는 사전조건 비공액 그라디언트(BiCG‑P) 알고리즘을 제안하고, 그 성능을 수치 실험을 통해 검증한다.

상세 분석

논문은 먼저 전통적인 ALI(Accelerated Λ‑Iteration) 방법이 Λ 연산자의 대각 성분을 근사함으로써 Jacobi 형태의 반복을 수행한다는 점을 상기한다. 이후 Trujillo Bueno와 Fabiani Bendicho가 제안한 Gauss‑Seidel(GS) 및 Successive Over‑Relaxation(SOR) 방식이 비대칭 행렬에도 적용 가능하도록 확장되었으나, 여전히 수렴 속도와 메모리 요구량에서 한계가 있다. 저자들은 이러한 한계를 극복하기 위해, 방사전달 방정식을 선형 시스템 A x = b 형태로 재구성하고, A = I − (1−ε)Λ 로 정의한다. 여기서 ε는 충돌 파괴 확률이며, Λ는 전통적인 Λ 연산자이다.

비대칭이면서 양의 정부호가 아닌 A에 대해, 일반적인 CG(Conjugate Gradient) 방법은 적용이 불가능하므로, 비공액 그라디언트(BiCG) 방법을 선택한다. 그러나 BiCG는 행렬 전치 Aᵀ를 필요로 하며, ε가 매우 작아 알베도가 1에 가까워질 경우 조건수가 급격히 악화된다. 이를 완화하기 위해 저자들은 M을 A의 대각 성분으로 설정한 사전조건(preconditioner) M⁻¹을 도입한다. 사전조건을 적용하면 실제 연산은 M⁻¹A x = M⁻¹b 형태가 되며, 대각 행렬 M의 역은 매우 저렴하게 계산된다.

알고리즘은 초기 추정 x₀(보통 평면 평행 슬랩의 초기 소스 함수)에서 시작해, 잔차 r₀ = b − A x₀와 그 전치 버전 ˜r₀을 정의한다. 이후 반복마다 M z = r, Mᵀ ˜z = ˜r 를 풀어 z, ˜z를 얻고, 스칼라 ρ, β, α 등을 계산해 탐색 방향 p, ˜p와 새로운 근사 x, 잔차 r, ˜r을 업데이트한다. 수렴 판단은 ‖r‖₂/‖b‖₂ < 10⁻² 로 설정한다.

수치 실험에서는 광학 깊이 τ_max = 2 × 10⁸, ε = 10⁻⁴ ~ 10⁻¹² 범위, 1 decade당 5 ~ 20 격점(공간 그리드) 등을 사용했다. 결과는 다음과 같다. 첫째, BiCG‑P는 GS와 SOR에 비해 반복 횟수가 2~3배 정도 감소했으며, 전체 실행 시간은 ALI과 비슷하거나 약간 빠른 수준을 유지했다. 둘째, 격자 해상도를 높이면 true error(정확도) 플래토가 감소하고, 수렴 속도는 격자 수에 거의 선형적으로 비례한다는 점을 확인했다. 셋째, 초기 소스 함수에 고주파 오류를 인위적으로 삽입했을 때, BiCG‑P는 첫 번째 반복에서 이미 오류를 크게 감소시키고, 이후 반복에서 완전히 스무딩되는 특성을 보였다. 이는 다중 격자(Multigrid) 방법에서 평활 연산자로 활용될 가능성을 시사한다.

마지막으로, BiCG‑P는 Aᵀ 연산이 필요하다는 점이 구현 복잡성을 증가시킨다. 저자들은 현재 이 추가 연산 비용과 수렴 이득 사이의 트레이드‑오프를 정량화하는 작업을 진행 중이며, 다중 레벨 및 다차원(2‑D, 3‑D) 문제에 대한 확장 가능성도 검토하고 있다.