[논문 리뷰] Stochastic Solution of Elliptic and Parabolic Boundary Value Problems for the Spectral Fractional Laplacian

Mamikon Gulian, Guofei Pang|arXiv (Cornell University)|2018. 12. 04.

Advanced Mathematical Modeling in Engineering참고 문헌 23인용 수 5

한 줄 요약

이 논문은 비영인 딜리클레 경계 조건을 갖는 스펙트럼 분수 라플라스 연산자 $(-\Delta_{\Omega,g})^{\alpha/2}$를 포함하는 타원형 및 방정형 경계값 문제에 대해 새로운 스토케스틱 해 공식(Feynman-Kac 유형)을 제시한다. 이 방법은 보조 정지 브라운 운동을 사용하며, 몬테카를로 시뮬레이션을 통한 확률적 프레임워크를 제공하여 유한 도메인 내에서 효율적이고 국소적으로 병렬화 가능한 계산이 가능하며, 수렴성과 정(regularity) 분석이 엄밀하게 이루어진다.

ABSTRACT

We prove and implement stochastic solution (or Feynman-Kac) formulas for boundary value problems involving the spectral fractional Laplacian with nonzero Dirichlet boundary condition. The main tools used in the proofs are the abstract Cauchy problem for Feller semigroups together with Balakrishnan's theory of fractional powers. We show the the spectral fractional Laplacian with nonzero Dirichlet boundary conditions is the generator of an appropriate Feller semigroup for subordinate stopped Brownian motion, and obtain a stochastic solution formula for the fractional heat equation in a bounded domain. We then obtain precise regularity and steady-state convergence properties of the parabolic problem using the eigenfunction expansion of the classical solution, which leads to estimates for the survival probability of subordinate stopped Brownian motion. These results allow us to take the long-time limit in the parabolic formula to establish a stochastic solution formula for the Dirichlet boundary value problem. These stochastic solution formulas for the spectral fractional Laplacian with nonzero boundary conditions are novel, and allow for efficient, embarrassingly parallel local solution of the above boundary value problems. We discuss the discretization of these formulas, and verify them in dimensions two and three with benchmark examples. We study the effect of the number of path samples and the time step size for path discretization on the accuracy of the solution.

연구 동기 및 목표

• 비영인 딜리클레 경계 조건을 갖는 스펙트럼 분수 라플라스 연산자에 대한 스토케스틱 해 공식을 개발하고 검증하기.
• 펠러 반군 이론과 보조 이론을 통해 스펙트럼 분수 라플라스 연산자의 확률적 해석을 수립하기.
• 경로 기반 몬테카를로 샘플링을 통한 분수 PDE의 효율적이고 명백하게 병렬화 가능한 수치적 해법을 가능하게 하기.
• 고유함수 전개를 이용해 분수 열 방정식의 정규성과 안정 상태 수렴성을 분석하기.
• 2차원 및 3차원 도메인에서의 수치 벤치마크를 통해 제안된 공식을 검증하기.

제안 방법

• 밸라크리슈난의 분수 거듭제곱 이론을 통해 스펙트럼 분수 라플라스 연산자를 펠러 반군의 생성자로 유도한다.
• 과정 $X^{\Omega,\alpha}_{t} = \sqrt{2}B_{T_{\alpha/2}(t) \wedge \tau_{\Omega}}$를 구성한다. 이는 $\alpha/2$-스테이블 보조자에 의해 시간이 조정된 브라운 운동이며, 도메인 외부로 나가기 전까지 정지된다.
• 비동차 분수 코시 문제에 대한 페인만-카프 공식을 수립한다: $u(t,x) = \mathbb{E}_x\left[ f(X_t) \chi_{\tau_\Omega > T_{\alpha/2}(t)} + g(X_t) \chi_{\tau_\Omega \leq T_{\alpha/2}(t)} + \int_0^t r(t-s,X_s) \chi_{\tau_\Omega > T_{\alpha/2}(s)} ds \right]$.
• 정적 상태 딜리클레 문제의 해를 $u(x) = \mathbb{E}_x\left[ g(X_{T^{-1}_{\alpha/2}(\tau_\Omega)}) \right] + \mathbb{E}_x\left[ \int_0^{T^{-1}_{\alpha/2}(\tau_\Omega)} r(X_s) ds \right]$를 통해 확장한다.
• 고유함수 전개를 이용해 보조 정지 과정의 정규성 및 생존 확률 추정치를 유도한다.
• 경로 수와 시간 간격 크기의 변화에 따른 수렴 분석을 수행하며, 2D 및 3D 도메인에서 경로 샘플링을 통한 수치적 구현과 검증을 수행한다.

실험 결과

연구 질문

• RQ1 비영인 딜리클레 경계 조건을 갖는 스펙트럼 분수 라플라스 연산자는 어떤 스토케스틱 과정을 통해 표현될 수 있는가?
• RQ2 보조 정지 브라운 운동과 분수 편미분방정식의 타원형 및 방정형 해 사이의 관계는 무엇인가?
• RQ3 보조 정지 과정의 생존 확률은 해의 정규성과 수렴성에 어떤 영향을 미치는가?
• RQ4 스토케스틱 해 공식은 보장된 수렴성과 병렬 확장성을 갖는 수치적 구현이 가능한가?
• RQ5 경로 수와 시간 간격 크기의 관점에서 몬테카를로 이산화의 정확도와 수렴 거동은 어떠한가?

주요 결과

• 스펙트럼 분수 라플라스 연산자 $(-\Delta_{\Omega,g})^{\alpha/2}$는 보조 정지 브라운 운동과 관련된 펠러 반군의 생성자로 나타남이 입증되었다.
• 추상 코시 문제와 보조 이론을 이용해 분수 열 방정식의 스토케스틱 해 공식이 도출되고 증명되었다.
• 정적 상태 딜리클레 문제의 해는 파라볼릭 공식에서 $t \to \infty$로 갈 때 얻어지며, 과정의 퇴출 시간을 포함한 경로 기반 기댓값 형태로 표현된다.
• 2D 및 3D에서의 수치 실험 결과 경로 수 증가와 시간 간격 감소에 따라 수렴함을 보여주며, 방법의 정확성을 검증하였다.
• 이 방법은 국소적으로 병렬화 가능하고 격자 없는 분수 PDE의 해법을 가능하게 하며, 고차원에서 우수한 스케일링 성능을 보였다.
• 고유함수 전개를 통해 정규성 및 생존 확률 추정치가 도출되었으며, 이는 해 프레임워크의 잘 정의됨을 확인한다.