[논문 리뷰] Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise

Rikard Anton, David Cohen외 2명|arXiv (Cornell University)|2015. 02. 28.

Stochastic processes and financial applications참고 문헌 17인용 수 10

한 줄 요약

이 논문은 다중성 노이즈를 가진 반선형 스토크라스틱 웨이브 방정식을 위한 완전 이산 수치적 방법을 제안한다. 이는 선형 유한요소 공간 이산화와 시간 적분을 위한 스토크라스틱 삼각함수 방법(STM)을 조합한 것이다. 주요 기여는 공간 메쉬 크기에 영향을 받지 않는 오차 한계를 가진 평균 제곱수렴을 달성하는 것으로, 기존의 에일러-마리야무 방법과 같은 표준적 방법보다 장기적인 에너지 행동에서 뛰어나다. 또한 예상 에너지의 선형 드리프트(선형 드리프트)를 가지는 트레이스 공식을 유지한다.

ABSTRACT

A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size restriction as in the often used Störmer-Verlet-leap-frog scheme. Furthermore, it satisfies an almost trace formula (i.e., a linear drift of the expected value of the energy of the problem). Numerical experiments are presented and confirm the theoretical results.

연구 동기 및 목표

• 다중성 노이즈에 의해 영향을 받는 반선형 스토크라스틱 웨이브 방정식을 위한 완전 이산 수치적 방법을 개발한다.
• 공간 이산화에 영향을 받지 않는 평균 제곱수렴을 달성하며, 오차 한계를 확보한다.
• 예상 에너지의 선형 드리프트를 가지는 트레이스 공식과 같은 기하학적 성질을 수치적 해에서 유지한다.
• 기존의 스토름어-버렛-리프로그와 같은 전통적 방법에서 발생하는 단계 크기 제약을 극복한다.

제안 방법

• 유한요소 방법을 사용하여 경계가 다각형인 유계 볼록 영역에서 표준 선형 유한요소 방법을 적용해 공간 이산화를 수행한다.
• 시간 적분은 스토크라스틱 웨이브 방정식에 특화된 명시적 시간 적분기인 스토크라스틱 삼각함수 방법(STM)을 사용한다.
• 공간 반이산화 이후 얻어진 반이산계를 대상으로 STM을 적용하여 안정성과 정확도를 확보한다.
• 특히 가중치가 부여된 소볼레프 노름에서 힐베르트 공간 내의 에너지 추정과 보간 기법을 사용하여 방법을 분석한다.
• 정확한 해, 유한요소 반이산해, 완전 이산해에 대해 트레이스 공식을 유도하였으며, 예상 에너지의 선형 드리프트를 보여준다.
• 수치 실험을 통해 이론적 수렴률과 다양한 시간 단계 및 메쉬 크기에서의 트레이스 공식 유지 여부를 검증한다.

실험 결과

연구 질문

• RQ1 다중성 노이즈를 가진 반선형 스토크라스틱 웨이브 방정식을 위한 완전 이산 수치적 방법을 구성할 수 있는가? 이 방법은 공간 메쉬 크기에 영향을 받지 않는 평균 제곱수렴을 달성할 수 있는가?
• RQ2 제안된 스토크라스틱 삼각함수 방법은 정확한 해의 트레이스 공식(예상 에너지의 선형 드리프트)을 유지하는가?
• RQ3 장기적인 에너지 진동 시뮬레이션에서 STM은 에일러-마리야무나 크랭크-니콜슨-마리야무와 같은 전통적 방법보다 성능이 뛰어나게 되는가?
• RQ4 완전 이산 수치적 방법의 평균 제곱수렴 속도는 무엇이며, 이는 공간 및 시간 이산화 매개변수에 어떻게 의존하는가?
• RQ5 큰 시간 단계를 사용할 경우에도 제안된 방법은 에너지 드리프트 행동과 같은 기하학적 성질을 유지할 수 있는가?

주요 결과

• 제안된 완전 이산 수치적 방법은 공간 메쉬 크기에 영향을 받지 않는 평균 제곱수렴을 달성하며, 스토름어-버렛-리프로그 방법에서 발생하는 단계 크기 제약을 피한다.
• 스토크라스틱 삼각함수 방법은 트레이스 공식을 유지하여 수치적 해의 예상 에너지가 정확한 해와 일치하는 선형 드리프트를 보인다.
• 수치 실험 결과 이론적 수렴률이 확인되었으며, 관측된 기울기들이 시간과 공간에서 예측된 수렴 차수와 일치한다.
• STM은 장기적인 시뮬레이션에서 에일러-마리야무 유형의 방법보다 뛰어나며, 예상 에너지의 선형 드리프트를 정확히 포착하는 반면, 에일러 방법은 잘못된 에너지 증가를 보인다.
• 큰 시간 단계(예: k=0.1)와 긴 시간 간격(예: [0,250])에서도 STM은 예상 에너지 변화에 대해 뛰어난 행동을 유지하여 강건성을 확인한다.
• 수치 실험에서 몬테카를로 오차는 극히 미미하다(예: <0.01), 이는 관측된 수렴성과 트레이스 공식 행동의 신뢰성을 검증한다.