[논문 리뷰] Numerical approximation of the stochastic heat equation with a distributional reaction term

Ludovic Goudenège, El Mehdi Haress외 1명|arXiv (Cornell University)|2024. 05. 13.

Advanced Thermodynamics and Statistical Mechanics인용 수 1

한 줄 요약

이 논문은 일차원 공간에서 분포적 반응 항을 가진 확률적 열 방정식에 적용된 테임드 오일러 유한차분 스킴에 대한 수렴 속도를 확립한다. 스토하스틱 실링 기법과 비소르 정규성의 드리프트를 활용하여, 드리프트의 정규성에 따라 의존하는 L^m 수렴 속도를 증명한다. 드리프트가 유계 가측 함수일 경우 공간에서 거의 최적의 (1/2 − ε) 및 시간에서 (1/4 − ε) 수렴 속도를 달성하며, 드리프트가 유한 측도일 경우조차도 양의 수렴 속도를 확보한다.

ABSTRACT

We study the numerical approximation of the stochastic heat equation with a distributional reaction term. Under a condition on the Besov regularity of the reaction term, it was proven recently that a strong solution exists and is unique in the pathwise sense, in a class of H\"older continuous processes. For a suitable choice of sequence $(b^k)_{k\in \mathbb{N}}$ approximating $b$, we prove that the error between the solution $u$ of the SPDE with reaction term $b$ and its tamed Euler finite-difference scheme with mollified drift $b^k$, converges to $0$ in $L^m(\Omega)$ with a rate that depends on the Besov regularity of $b$. In particular, one can consider two interesting cases: first, even when $b$ is only a (finite) measure, a rate of convergence is obtained. On the other hand, when $b$ is a bounded measurable function, the (almost) optimal rate of convergence $(\frac{1}{2}-\varepsilon)$-in space and $(\frac{1}{4}-\varepsilon)$-in time is achieved. Stochastic sewing techniques are used in the proofs, in particular to deduce new regularising properties of the discrete Ornstein-Uhlenbeck process.

연구 동기 및 목표

• 스토하스틱 열 방정식의 수치 수렴 결과를 유계 가측 드리프트를 초월하여 비소르 공간에 속하는 분포적 드리프트로 확장하는 것.
• 경로적으로 병합된 드리프트를 가진 테임드 오일러 유한차분 스킴의 강수렴 속도를 확립하는 것.
• 반응 항 b의 비소르 정규성에 기반하여 L^m(Ω)에서의 수렴 속도를 정량화하는 것.
• b가 유한 측도일 경우에도 수렴이 가능함을 보여주는 것, 특히 함수가 아니더라도.

제안 방법

• 공간에 중심 차분, 시간에 명시적 스텝을 사용한 테임드 오일러 유한차분 스킴을 구성하며, 이를 위해 모리피케이션 드리프트 수열 (bk)k∈N을 사용한다.
• 스토하스틱 실링 기법을 통해 분석을 수행하며, 특히 이산 오르누이즈-울렌벡 과정과 그 정규화 성질을 다룰 수 있도록 적응시켰다.
• 진짜 해 u와 수치 해 un,k 사이의 오차는 모멘트 추정과 드리프트 차이 ∥b − bk∥B^{γ−1}_p의 비소르 노름 제어를 통해 유계화된다.
• 약한 정규성 가정 하에 수치 오차의 성장을 제어하기 위해 로그 인자와 함께 비슷한 그로엔발의 부등식이 사용된다.
• 특수한 지수를 가진 수정된 스토하스틱 실링 보조정리의 응용에 기반하여, 특이 드리프트가 존재하는 상황에서도 오차 제어가 가능하다.
• 모멘트 유계, 이산 세미군트 추정, 해의 경로적 정규성의 조합을 통해 수렴을 확립한다.

실험 결과

연구 질문

• RQ1 반응 항 b가 비소르 공간 B^γ_{p,∞}에 속하는 분포일 경우, 스토하스틱 열 방정식의 수치 해에 대해 달성 가능한 수렴 속도는 무엇인가?
• RQ2 b가 유한 측도일 뿐만 아니라 함수가 아닐 경우, 테임드 오일러 유한차분 스킴이 L^m(Ω)에서 수렴할 수 있는가?
• RQ3 b가 유계 가측 함수일 경우 시간과 공간에서의 최적 수렴 속도는 무엇이며, 기존에 알려진 정확한 수렴 속도인 시간 1/4 및 공간 1/2에 얼마나 가까이 갈 수 있는가?
• RQ4 스토하스틱 실링 기법과 이산 오르누이즈-울렌벡 과정의 정규화 성질은 특이 드리프트 존재 하에서 오차 추정에 어떻게 기여하는가?
• RQ5 로그 항을 포함한 비슷한 그로엔발 보조정리를 사용하여, 저정규성 드리프트를 가진 수치 스킴에서 오차 전파를 제어할 수 있는가?

주요 결과

• 수치 스킴은 비소르 정규성의 드리프트 b에 따라 의존하는 L^m(Ω) 수렴 속도를 보이며, 특히 ∥b − bk∥B^{γ−1}_p에 따라 결정된다.
• b가 유한 측도일 경우(즉, γ − 1/p ≤ −1)에도 불구하고 양의 수렴 속도를 확보한다.
• 유계 가측 드리프트의 경우(γ = 0, p = ∞) 스킴은 공간에서 (1/2 − ε) 및 시간에서 (1/4 − ε)의 거의 최적 수렴 속도를 달성한다.
• 약한 정규성 하에서도 이산 해 오차를 제어할 수 있도록 특수 지수를 가진 새로운 스토하스틱 실링 보조정리의 응용에 기반한다.
• 이산 오르누이즈-울렌벡 과정의 새로운 정규화 성질이 도출되었으며, 이는 수치 스킴의 오차 유계화에 핵심적인 역할을 한다.
• 로그 항을 포함한 비슷한 그로엔발 보조정리의 사용은 드리프트의 비립시츠 성질을 다루고 스킴 내 오차 성장을 제어할 수 있게 한다.