[논문 리뷰] Numerical solutions of some hyperbolic stochastic partial differential equations with mixed derivatives including sine-Gordon equation

Henry C. Tuckwell|arXiv (Cornell University)|2015. 08. 02.

Stochastic processes and financial applications참고 문헌 28인용 수 3

한 줄 요약

이 논문은 공간-시간 가우시안 화이트 노이즈에 의해 구동되는 혼합 도함수를 포함하는 쌍곡형 스토케스틱 편미분방정식(SPDE)을 해결하기 위한 명시적 수치적 방법을 제시한다. 이 방법은 결정론적 및 스토케스틱 영역에서 파동 유사한 구조를 정확하게 포착하며, 특히 작은 노이즈 강도(σ ≤ 0.25)일 경우 초기 조건이 0이어도 노이즈에 의해 유도된 안정적인 수렴성과 응집된 역학적 행동이 나타나는 것으로 나타났다.

ABSTRACT

We consider linear and nonlinear hyperbolic SPDEs with mixed derivatives with additive space-time Gaussian white noise of the form $Y_{xt}=F(Y) + σW_{xt}.$ Such equations, which transform to linear and nonlinear wave equations, including Klein-Gordon, Liouville's and the sine-Gordon equation, are related to what Zimmerman (1972) called a diffusion equation. An explicit numerical scheme is employed in both deterministic and stochastic examples. The scheme is checked for accuracy against known exact analytical solutions for deterministic equations. In the stochastic case with $F=0$, solutions yield sample paths for the Brownian sheet whose statistics match well exact values. Generally the boundary conditions are chosen to be initial values $Y(x,0)$ and boundary values $Y(0,t)$ on the quarter-plane or subsets thereof, which have been shown to lead to existence and uniqueness of solutions. For the linear case solutions are compared at various grid sizes and wave-like solutions were found, with and without noise, for non-zero initial and boundary conditions. Surprisingly, wave-like structures seemed to emerge with zero initial and boundary conditions and purely noise source terms with no signal. Equations considered with nonlinear $F$ included quadratic and cubic together with the sine-Gordon equation. For the latter, wave-like structures were apparent with $σ\le 0.25$ but they tended to be shattered at larger values of $σ$. Previous work on stochastic sine-Gordon equations is briefly reviewed.

연구 동기 및 목표

• 혼합 편도함수와 가감성 공간-시간 화이트 노이즈를 갖는 쌍곡형 SPDE를 해결하기 위한 수치적 방법을 개발하고 검증한다.
• 초기 조건과 경계 조건이 0이지만 노이즈가 존재할 경우 스토케스틱 SPDE에서 파동 유사 행동이 어떻게 나타나는지 조사한다.
• 결정론적 경우의 알려진 해석적 해와의 비교를 통해 수치적 방법의 정확성과 수렴성을 평가한다.
• 다양한 노이즈 강도(σ)에서 비선형 SPDE, 특히 사인-고든 방정식의 거동을 탐구한다.
• 좌표 변환을 통한 혼합 도함수 SPDE 형태의 해와 고전적 파동 방정식 간의 관계를 분석한다.

제안 방법

• 공간-시간 화이트 노이즈 $ W_{xt} $를 갖는 SPDE $ Y_{xt} = F(Y) + \sigma W_{xt} $ 를 해결하기 위해 명시적 유한차분 스킴을 적용한다.
• 선형 및 비선형 결정론적 사례, 즉 파동 방정식과 리우빌 방정식에 대해 정확한 해석적 해와의 비교를 통해 스킴을 검증한다.
• 스토케스틱 케이스에서는 브라운 운동 시트의 표본 경로를 시뮬레이션하며, 통계적 성질을 이론적 기대값과 비교한다.
• 해의 유일성을 보장하기 위해 초기 조건 $ Y(x,0) $ 와 경계 조건 $ Y(0,t) $ 가 일정한 1/4 평면 영역 $[0, T] \times [0, X]$ 에서 수치 실험을 수행한다.
• 수렴성을 평가하기 위해 격자 해상도를 높이는 실험을 수행하며, 다양한 격자 크기(예: 400×400 및 800×800)와 시도 횟수의 해를 비교한다.
• 노이즈 강도 $ \sigma $ 를 체계적으로 변화시키며 제곱, 세제곱, 사인-고든 $ F(Y) $ 를 포함한 비선형 케이스에 스킴을 적용한다.

실험 결과

연구 질문

• RQ1 혼합 도함수와 가감성 공간-시간 화이트 노이즈를 갖는 쌍곡형 SPDE를 명시적 수치 스킴이 정확하게 시뮬레이션할 수 있는가?
• RQ2 초기 조건과 경계 조건이 모두 0이지만 시스템이 순수 노이즈에 의해 구동될 경우, 해에서 파동 유사한 구조가 나타나는가?
• RQ3 노이즈 강도 $ \sigma $ 는 스토케스틱 사인-고든 방정식에서 파동 유사 해의 안정성과 응집성에 어떤 영향을 미치는가?
• RQ4 노이즈 강도 $ \sigma \to 0 $ 일 때 수치적 해가 결정론적 해에 얼마나 수렴하는가? 격자 해상도는 이 과정에 어떤 영향을 미치는가?
• RQ5 혼합 도함수 SPDE $ Y_{xt} = F(Y) + \sigma W_{xt} $ 의 해와 고전적 파동 방정식 $ \phi_{tt} - \phi_{xx} = G(\phi) $ 간의 관계는 어떠한가?

주요 결과

• 명시적 수치 스킴은 결정론적 선형 및 비선형 SPDE에 대해 알려진 해석적 해를 정확히 재현하여 신뢰성을 확인한다.
• 선형 SPDE에서 $ F = 0 $ 인 경우, 시뮬레이션된 표본 경로의 통계적 성질이 이론적 기대값과 일치하여 노이즈 통합의 정당성이 검증되었다.
• 특히 작은 $ \sigma $ 에서는 초기 조건과 경계 조건이 0이어도 노이즈에 의해 파동 유사한 구조가 나타나며, 이는 내재된 노이즈 유도 패턴 형성의 증거이다.
• 사인-고든 방정식의 경우 $ \sigma \leq 0.25 $ 일 때 응집된 파동 유사한 구조가 안정적이었지만, 높은 노이즈 강도에서는 파편화되었다.
• 격자 해상도를 높이는 것(예: 800×800 대비 400×400)은 파동 구조의 정의를 향상시키고 표준편차를 감소시켜 더 미세한 해상도에서의 수렴성을 나타냈다.
• $ \sigma = 0.1 $ 일 때는 평균의 스토케스틱 시뮬레이션과 결정론적 해 사이에 양호한 일치가 관찰되었으며, 특히 800×800 격자와 50회의 시행에서 수치적 안정성이 높은 것으로 나타났다.