극대형 열전도 문제를 위한 새로운 안정적 명시적 해법

초록

본 논문은 균일한 열방정식의 대규모 이산화된 시스템에 대해, 기존 명시적·암시적 방법보다 높은 안정성과 효율성을 보이는 새로운 명시적 수치 알고리즘을 제시한다. 2차원 불균일 매개변수를 가진 두 사례에 적용해 시간‑스텝 제한을 크게 완화하고, 계산 비용을 크게 절감함을 실험적으로 입증한다.

상세 분석

열전도 방정식의 공간적 이산화는 일반적으로 수십만~수백만 차원의 강성 ODE 시스템을 만든다. 전통적인 명시적 방법은 CFL 조건에 의해 허용되는 시간 스텝 Δt가 격자 크기 h²에 비례해 매우 작아져, 전체 시뮬레이션 시간이 비현실적으로 길어진다. 반면 암시적 방법(예: Crank‑Nicolson, 후방 Euler)은 무조건적인 안정성을 제공하지만, 매 시간 스텝마다 대규모 선형 시스템을 풀어야 하므로 메모리와 연산량이 급증한다. 특히 매개변수가 공간적으로 무작위로 변동하는 경우, 계수 행렬이 비대칭·비정형이 되어 전형적인 전처리(preconditioner) 설계가 어려워진다.

저자들은 이러한 문제점을 해결하기 위해 “안정화된 명시적 스키마(stabilized explicit scheme)”를 도입한다. 핵심 아이디어는 다중 단계·다중 스텝 Runge‑Kutta‑Chebyshev(RKC) 혹은 ROCK2와 유사한 고차 안정 영역을 갖는 다항식 근사법을 사용해, 허용 가능한 Δt를 h²가 아닌 O(h) 수준까지 확대하는 것이다. 구체적으로, 시간 전진 연산을 다음과 같이 구성한다.

1. 초기값에서 시작해 Chebyshev 다항식의 근을 이용해 여러 중간 스텝을 계산한다.
2. 각 중간 스텝에서는 계수 행렬 A의 행벡터에 대한 행렬‑벡터 곱만 수행하므로, 메모리 요구량은 O(N)이며, 병렬화가 용이하다.
3. 최종 스텝에서는 가중 평균을 통해 2차 정확도를 유지한다.

안정성 분석에서는 다항식 근의 실수 구간이 λ·Δt∈