[논문 리뷰] A Break of the Complexity of the Numerical Approximation of Nonlinear SPDEs with Multiplicative Noise

Arnulf Jentzen, Michael Roeckner|arXiv (Cornell University)|2010. 01. 15.

Stochastic processes and financial applications참고 문헌 50인용 수 8

한 줄 요약

이 논문은 다중성(trace class) 소음을 가진 비선형 확률 편미분방정식(SPDEs)을 시뮬레이션하기 위한 새로운 무한차원 밀스타인 유형 알고리즘을 제안한다. 이는 이전 방법에 비해 반복적인 확률적 적분의 시뮬레이션을 단순화한다. 이 방법은 계산 복잡도에 있어 뚜렷한 개선을 이룩하였으며, 확률적 열 방정식, 반응-확산 방정식, 버거스 방정식에 대한 수치 실험에서 상당한 효율성 향상을 보였다.

ABSTRACT

A new algorithm for simulating stochastic partial differential equations (SPDEs) of evolutionary type, which is in some sense an infinite dimensional analog of Milstein’s scheme for finite dimensional stochastic ordinary differential equations (SODEs), is introduced and analyzed in this article. The Milstein scheme is known to be impressively efficient for scalar one-dimensional SODEs but only for some special multidimensional SODEs due to difficult simulations of iterated stochastic integrals in the general multidimensional SODE case. It is a key observation of this article that, in contrast to what one may expect, its infinite dimensional counterpart introduced here is very easy to simulate and this, therefore, leads to a break of the complexity (number of computational operations and random variables needed to compute the scheme) in comparison to previously considered algorithms for simulating nonlinear SPDEs with multiplicative trace class noise. The analysis is supported by numerical results for a stochastic heat equation, stochastic reaction diffusion equations and a stochastic Burgers equation showing significant computational savings.

연구 동기 및 목표

• 다중성 흔들림(trace class) 소음을 가진 비선형 SPDEs를 시뮬레이션할 때 발생하는 높은 계산 비용을 해결하기 위해.
• 기존에 SODEs에 대해 잘 알려진 밀스타인 스킴의 효율성을 무한차원 SPDEs의 설정으로 확장하기 위해.
• 일반적인 다차원 SODEs에서 반복적인 확률적 적분을 시뮬레이션하는 데 어려움을 해결하기 위해, 무한차원 사례에서의 구조적 이점들을 활용하기 위해.
• 진화형 SPDEs를 위한 수치적으로 실현 가능하고 계산 효율성이 높은 알고리즘을 개발하기 위해.
• 이론적 복잡도 감소를 수치적 검증을 통해 입증하고, 기존 방법들과 비교하여 계산 복잡도의 상당한 감소를 보여주기 위해.

제안 방법

• 다중성 소음의 구조를 활용하여 SPDEs에 대한 밀스타인 스킴의 무한차원 유사체를 제안한다.
• 스펙트럼 또는 갈레르킨 유사체 공간 근사법을 사용하여 SPDE를 SDE 시스템으로 축소시켜, 밀스타인 유형의 시간 이산화를 적용할 수 있도록 한다.
• 무한차원 설정에서 반복적인 확률적 적분의 시뮬레이션을 단순화하여, 소음의 구조 덕분에 이러한 항목들이 다루기 쉬워지도록 한다.
• 두 번째 차수까지의 확률적 테일러 전개를 적용하여 스킴을 유도함으로써 약한 수준의 수렴 순서를 확보한다.
• 갈레르킨 투영을 사용하여 다양한 SPDEs, 즉 확률적 열 방정식, 반응-확산 방정식, 버거스 방정식에 스킴을 적용한다.
• 수치 실험을 통해 이론적 복잡도 감소를 검증하고 계산 자원 절감 정도를 평가한다.

실험 결과

연구 질문

• RQ1 일반적인 SPDEs에서 반복적인 확률적 적분의 계산 부담을 피하면서도, 밀스타인 스킴의 무한차원 확장판을 구성할 수 있는가?
• RQ2 제안된 알고리즘이 기존의 비선형 SPDEs에 대한 다중성 소음 스킴에 비해 상당한 계산 복잡도 감소를 달성하는가?
• RQ3 스토하스틱 열 방정식과 버거스 방정식과 같은 다양한 비선형 SPDEs에 대해 이 스킴은 수치적으로 안정적이고 효율적인가?
• RQ4 작업 수와 난수 변수의 측면에서, 이 새로운 알고리즘의 복잡도는 표준 오일러 유형 또는 기타 고차수 스킴과 비교해 어떻게 나타나는가?
• RQ5 무한차원에서의 다중성 소음의 구조적 특성을 활용하여 확률적 적분의 시뮬레이션을 단순화할 수 있는가?

주요 결과

• 제안된 무한차원 밀스타인 유형 스킴은 다중성 흔들림(trace class) 소음을 가진 비선형 SPDEs를 위한 기존 알고리즘에 비해 상당한 계산 복잡도 감소를 달성한다.
• 유한차원 SODEs에서는 계산적으로 금기시되는 반복적인 확률적 적분의 시뮬레이션이, 무한차원 설정에서는 소음의 구조 덕분에 실현 가능하고 효율적으로 구현된다.
• 확률적 열 방정식, 확률적 반응-확산 방정식, 확률적 버거스 방정식에 대한 수치 실험에서 상당한 계산 자원 절감이 이루어졌다.
• 알고리즘은 요구되는 난수 변수 수와 계산 작업 수를 크게 줄이면서도 높은 정확도를 유지한다.
• 특히 고정밀도가 요구되는 문제에서는 표준 오일러 유형 스킴보다 효율성이 뛰어나다.
• 핵심 통찰은 무한차원 설정이 유한차원 사례에서는 존재하지 않는 단순화를 가능하게 하며, 이는 고차수 스킴의 실용적 구현을 가능하게 한다.