구면 위 반선형 확률 파동 방정식의 전산적 전시간·공간 전산화

초록

본 논문은 구면 S² 위에 정의된 반선형 확률 파동 방정식을, 구면 조화함수 기반의 스펙트럴 갈루아 방법과 명시적 확률 삼각 적분기법을 결합한 전시간·공간 전산화 스킴으로 근사한다. 전역적인 강수렴률과 거의 확실한 수렴률을 이론적으로 증명하고, 수치 실험을 통해 이론적 속도를 확인한다.

상세 분석

논문은 먼저 구면 S² 위의 등거리 측정과 라플라스‑벨트라미 연산자 Δ_{S²} 의 고유값·고유함수인 구면 조화함수 Y_{ℓ,m} 를 도입한다. 이를 통해 Q‑Wiener 과정의 카루노‑로베 확장식을 구면 조화함수 기반으로 표현하고, 잡음의 공분산 연산자 Q 의 스펙트럼 A_ℓ 을 명시한다. 반선형 확률 파동 방정식
∂{tt}u = Δ{S²}u + f(u) + g(u)·Ẇ(t)
을 추상적인 2차 연산자 A 와 비선형 연산자 F, G 로 재작성하고, 힐베르트 공간 H^s(S²) 위에서 약해(또는 약) 해의 존재·유일성을 보인다. 주요 가정은 f, g가 전역 리프시츠이며, g·Q^{1/2}가 H^{δ-1} 에 대한 Hilbert‑Schmidt 연산자를 형성한다(δ≥1이면 trace‑class 잡음).

시간 전산화는 “stochastic trigonometric integrator”(STI)라 불리는 명시적 스킴을 사용한다. 이 스킴은 정확히 cos와 sin 연산자를 이용해 선형 부분을 정확히 통합하고, 비선형·잡음 항을 Euler‑Maruyama 형태로 처리한다. 구면 조화함수 기반의 스펙트럴 Galerkin 근사는 차원 N 까지의 고유모드만 보존함으로써 공간 오차를 제어한다.

핵심 정리는 두 가지이다.

1. Theorem 9: 스펙트럴 근사에 대한 강수렴률 ‑                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                       

Theorem 9 (연속‑공간 수렴) : N →∞ 에 대해 ‖X(t)−X_N(t)‖{L^2(Ω;H)} ≤ C N^{‑α}, α 는 초기 데이터와 δ (노이즈 정규성) 에 의존.
Theorem 11 (시간 수렴) : Δt →0 에 대해 ‖X_N(t_n)−X{N,Δt}^n‖_{L^2(Ω;H)} ≤ C (Δt)^{β}, β ∈(0,½] (노이즈와 비선형성에 따라).

마지막으로 섹션 5에서는 2차원 구면 위에 정의된 테스트 문제(예: f(u)=−u³, g(u)=u)와 Q‑Wiener 잡음(ℓ‑스펙트럼 A_ℓ∼ℓ^{‑2})를 사용해 수치 실험을 수행한다. 실험 결과는 이론적 수렴률(공간 α≈1, 시간 β≈½)을 정확히 재현하며, 경로별(Almost‑sure) 수렴도 관측된다.