[논문 리뷰] Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation

William McLean, Kassem Mustapha|arXiv (Cornell University)|2020. 12. 26.

Fractional Differential Equations Solutions참고 문헌 29인용 수 6

한 줄 요약

이 논문은 시간 분수 Fokker-Planck 방정식에 대한 공간 이산화 갈레르킨 근사에 대해 α ∈ (0,1]인 분수 확산 지수를 갖는 경우에 대해 균일 안정성을 확립한다. 저자들은 α에 대해 일관되게 유계인 상수를 갖는 안정성 추정을 유도하고, 비동차 소스 항을 고려하며, 기울기 안정성도 제공한다. 이러한 결과들은 이전 연구를 향상시켜 안정성 상수의 α에 의존하는 폭주 문제를 제거하고, 유한요소 방법의 오차 분석을 단순화한다.

ABSTRACT

We prove stability estimates for the spatially discrete, Galerkin solution of a fractional Fokker-Planck equation, improving on previous results in several respects. Our main goal is to establish that the stability constants are bounded uniformly in the fractional diffusion exponent $\alpha\in(0,1]$. In addition, we account for the presence of an inhomogeneous term and show a stability estimate for the gradient of the Galerkin solution. As a by-product, the proofs of error bounds for a standard finite element approximation are simplified.

연구 동기 및 목표

• 모든 α ∈ (0,1]에 대해 시간 분수 Fokker-Planck 방정식의 반연속 갈레르킨 해에 대해 균일 안정성 경계를 확립하는 것, 특히 고전적 극한 α = 1을 포함한다.
• 이전 연구가 g ≡ 0에 국한되었던 바에 비해, 안정성 분석에서 비동차 소스 항 g를 고려하는 것.
• 갈레르킨 해의 기울기에 대한 안정성 추정을 도출하여 공간 정(regularity) 분석을 향상시키는 것.
• 새로운 안정성 프레임워크를 활용하여 표준 유한요소 근사에 대한 오차 경계를 단순화하는 것.
• 불연속 갈레르킨 시간 이산화에 대한 분석을 확장하여, 적절한 조건 하에 α에 대해 균일한 안정성을 입증하는 것.

제안 방법

• 저차원 정규성 자료를 다루기 위해 시간 적분된 변분 형태 (4)를 사용하여 약한 해를 정의함으로써 잘 정의됨을 보장한다.
• 모든 시험 함수 χ ∈ X에 대해 변분 방정식 (6)을 만족하는 반연속 갈레르킨 해 uX ∈ X ⊂ H¹₀(Ω)를 정의한다.
• 분수 적분 추정과 에너지 기법을 활용하여 α에 의존하지 않는 유계 경계를 도출함으로써, α → 1일 때 폭주를 피한다.
• Poincaré 부등식과 Ψ(α)를 포함하는 분수 적분 부등식 (38)을 사용하여 Riemann–Liouville 분수 도함수 항을 제어한다.
• 조각별 다항식 근사를 사용한 시간에 대한 불연속 갈레르킨 방법을 적용하여, α에 의존하지만 유계인 상수를 갖는 완전 이산 안정성을 도출한다.
• 안정성 프레임워크를 활용하여, α와 정규성 매개수 r에 대한 명시적 의존성을 갖는 조각별 선형 유한요소 해에 대해 최적의 오차 경계를 유도한다.

실험 결과

연구 질문

• RQ1 시간 분수 Fokker-Planck 방정식의 반연속 갈레르킨 해에 대한 안정성 상수가 분수 지수 α ∈ (0,1]에 대해 독립적으로 유지될 수 있는가?
• RQ2 비영인 소스 항 g가 갈레르킨 해의 안정성에 어떤 영향을 미치며, 여전히 균일한 경계를 확보할 수 있는가?
• RQ3 갈레르킨 해의 기울기 행동은 α에 대해 균일한 안정성 측면에서 어떻게 되는가?
• RQ4 새로운 안정성 프레임워크는 기존 표준 유한요소 근사에 대한 오차 추정을 단순화하거나 향상시킬 수 있는가?
• RQ5 불연속 갈레르킨 시간 이산화가 분수 확산 방정식에 대해 α에 대해 균일한 안정성을 유지하는가?

주요 결과

• 논문은 반연속 갈레르킨 해에 대해 균일 안정성을 확립한다: 모든 t ∈ (0,T] 및 α ∈ (0,1]에 대해 ‖uX(t)‖ + t^{α/2}‖∇uX(t)‖ ≤ C(‖u₀X‖ + ∫₀ᵗ ‖g(s)‖ds + (1/t ∫₀ᵗ ‖sg(g)‖² ds)^{1/2})이며, C는 α에 독립적이다.
• α → 1일 때 안정성 상수가 유계로 유지되며, 이는 이전 연구에서 Cα ∼ (1−α)^{-1}로 α→1일 때 폭주하던 문제를 해결한다.
• 갈레르킨 해의 기울기에 대해 유사한 α-균일 추정이 성립하여, 개선된 정규성 분석이 가능해진다.
• 조각별 선형 유한요소 해 uh에 대해 오차 경계는 ∥uh − u∥ = O(t^{-α(2−r)/2} h²) 및 ∥∇uh − ∇u∥ = O(t^{-α(2−r)/2} h)이며, r ∈ [0,2]는 해의 정규성 정도를 정량화한다.
• r = 2 이고 u₀h = Rₕu₀일 경우 오차 경계는 시간에 대해 균일해진다: ∥uh(t) − u(t)∥ + h∥∇uh(t) − ∇u(t)∥ ≤ Cα h² (∥u₀∥₂ + M).
• 불연속 갈레르킨 시간 이산화의 경우, 해는 ∥Un⁻∥² + ∑ⱼ∥JUKⱼ∥² + ∫₀ᵗ ⟨∂¹⁻ᵃᵗ κ∇U, ∇U⟩ dt ≤ ∥u₀X∥² + CΩ Ψ(α)/κmin tₙ^{1−α} ∫₀ᵗ ‖g(s)‖² ds를 만족하며, Ψ(α)는 α→1일 때 유계이지만 α→0일 때 폭주한다.