그레이드 메시와 SOE 기반 템퍼드 시간 분수 대류‑확산 방정식 2차 정확도 해법

초록

본 논문은 템퍼드 시간‑분수 대류‑확산 방정식에 대해, 합지수(SOE) 근사와 시간‑그레이드 메시를 결합한 2차 정확도 시간스텝 스킴을 제안한다. 반시간 레벨에서의 차분과 2차 중심 차분을 이용해 공간을 이산화하고, 저장·연산 복잡도를 크게 낮추면서 안정성과 수렴성을 엄밀히 증명한다.

상세 분석

본 연구는 템퍼드 시간‑분수 미분 연산자의 핵심인 약한 특이성을 효율적으로 처리하기 위해 두 가지 핵심 기술을 도입한다. 첫째, 커널 (t‑s)^{‑α}e^{‑λ(t‑s)}를 합지수(SOE) 형태 ∑{ℓ=1}^{N{exp}} ω_ℓ e^{‑s_ℓ(t‑s)} 로 근사함으로써 과거 시간 단계에 대한 누적 연산을 재귀식으로 전개한다. 이때 N_{exp}는 지정 오차 ε에 따라 적당히 작은 값으로 선택될 수 있어, 저장 요구량이 O(M N) → O(M N_{exp}) 로, 연산 복잡도는 O(M N^2) → O(M N N_{exp}) 로 감소한다. 둘째, 초기 시점의 강한 특이성을 완화하기 위해 시간 격자를 t_n = T (n/N)^r (r≥1) 로 그레이드한다. r값을 적절히 조정하면 t≈0 근처에서 작은 시간 간격을 확보하면서 전체 오류를 (N^{‑r(δ+1‑α)}) 수준으로 제어한다. 또한, 기존 L1 스킴이 1차 정확도에 머무는 반면, 저자들은 반시간 레벨 \bar t_n = (t_n + t_{n+1})/2 에서 차분을 수행하고, 히스토리 부분과 로컬 부분을 각각 SOE와 선형 보간(L1)으로 처리함으로써 전역 2차 정확도를 달성한다. 오류 분석에서는 먼저 L1 근사의 truncation error를 상세히 추정하고, 그 후 SOE 근사에 의한 추가 오차를 ε와 결합해 최종적으로 |F C R_u(\bar t_n)| ≤ C·(N^{‑r(δ+1‑α)} + ε) (n=0) 및 C·((n+1)^{‑min{2‑α, r(1+δ)}} + ε) (n≥1) 를 얻는다. 안정성 증명은 에너지 추정법을 활용해 비대칭 행렬의 양의 정부호성을 확보하고, 격자 비율 조건 r≥1 하에서 무조건적인 안정성을 보인다. 이러한 이론적 결과는 수치 실험에서 확인되며, 특히 기존 2차 L2‑1σ 스킴 대비 메모리 사용량이 10배 이상 절감되는 효과가 입증된다. 전반적으로 SOE와 그레이드 메시의 결합은 템퍼드 분수 미분 방정식의 고정밀·고효율 수치 해법을 제공한다는 점에서 의미가 크다.