The analysis of nonlocal discrete equations driven by fractional powers of the discrete Laplacian on a mesh of size $h>0$ \[ (-Î_h)^su=f, \] for $u,f:\mathbb{Z}_h\to\mathbb{R}$, $0<s<1$, is performed. The pointwise nonlocal formula for $(-Î_h)^su$ and the nonlocal discrete mean value property for discrete $s$-harmonic functions are obtained. We observe that a characterization of $(-Î_h)^s$ as the Dirichlet-to-Neumann operator for a semidiscrete degenerate elliptic local extension problem is valid. Regularity properties and Schauder estimates in discrete Hölder spaces as well as existence and uniqueness of solutions to the nonlocal Dirichlet problem are shown. For the latter, the fractional discrete Sobolev embedding and the fractional discrete Poincaré inequality are proved, which are of independent interest. We introduce the negative power (fundamental solution) \[ u=(-Î_h)^{-s}f, \] which can be seen as the Neumann-to-Dirichlet map for the semidiscrete extension problem. We then prove the discrete Hardy--Littlewood--Sobolev inequality for $(-Î_h)^{-s}$.
As applications, the convergence of our fractional discrete Laplacian to the (continuous) fractional Laplacian as $h\to0$ in Hölder spaces is analyzed. Indeed, uniform estimates for the error of the approximation in terms of $h$ under minimal regularity assumptions are obtained. We finally prove that solutions to the Poisson problem for the fractional Laplacian \[ (-Î)^sU=F, \] in $\mathbb{R}$, can be approximated by solutions to the Dirichlet problem for our fractional discrete Laplacian, with explicit uniform error estimates in terms of~$h$.
논문 ID : 1608.08913제목 : Nonlocal discrete diffusion equations and the fractional discrete Laplacian, regularity and applications저자 : Óscar Ciaurri, Luz Roncal, Pablo Raúl Stinga, José L. Torrea, Juan Luis Varona분류 : math.AP cs.NA math.CA math.FA math.NA math.PR발표 시간 : 2018년 3월 21일 (v2)논문 링크 : https://arxiv.org/abs/1608.08913 본 논문은 이산 라플라시안의 분수 거듭제곱으로 구동되는 비국소 이산 방정식 ( − Δ h ) s u = f (-\Delta_h)^s u = f ( − Δ h ) s u = f 에 대한 심층 분석을 수행하였다. 여기서 u , f : Z h → R u,f: \mathbb{Z}_h \to \mathbb{R} u , f : Z h → R , 0 < s < 1 0 < s < 1 0 < s < 1 , h > 0 h > 0 h > 0 은 격자 크기이다. 본 연구는 ( − Δ h ) s u (-\Delta_h)^s u ( − Δ h ) s u 의 점별 비국소 공식과 이산 s s s -조화 함수의 비국소 이산 평균값 성질을 도출하였다. ( − Δ h ) s (-\Delta_h)^s ( − Δ h ) s 가 반이산 퇴화 타원 국소 확장 문제의 Dirichlet-to-Neumann 연산자로 특성화될 수 있음을 증명하였다. 이산 Hölder 공간에서의 정칙성 성질과 Schauder 추정, 그리고 비국소 Dirichlet 문제 해의 존재성과 유일성을 확립하였다. 음의 거듭제곱(기본해) u = ( − Δ h ) − s f u = (-\Delta_h)^{-s}f u = ( − Δ h ) − s f 를 도입하고 대응하는 이산 Hardy-Littlewood-Sobolev 부등식을 증명하였다. 응용으로서, h → 0 h \to 0 h → 0 일 때 분수 이산 라플라시안이 Hölder 공간에서 연속 분수 라플라시안으로 수렴하는 것을 분석하였으며, 분수 라플라시안 Poisson 문제의 해를 분수 이산 라플라시안 Dirichlet 문제의 해로 근사할 수 있음을 증명하였다.
분수 라플라시안 연산자는 고전 라플라시안의 정수 거듭제곱으로서 위치 이론, 조화 분석, 분수 미적분학, 함수 해석 및 확률론 등 여러 수학 분야에서 오랫동안 존재해왔다. 지난 십 년간 분수 라플라시안 연산자를 포함하는 비국소 확산 방정식은 편미분 방정식 분야에서 가장 활발한 연구 주제 중 하나가 되었다.
이산화 과제 : 분수 라플라시안 연산자 ( − Δ ) s (−\Delta)^s ( − Δ ) s 의 비국소성과 특이성이 수치 근사에 주요 어려움을 야기이론적 공백 : 체계적인 이산 분수 라플라시안 연산자 이론 틀의 부재수렴성 분석 : 이산 모델에서 연속 모델로의 수렴에 대한 엄격한 이론적 보증 필요완전한 이산 분수 라플라시안 연산자 이론 체계 구축 연속 문제의 효과적인 이산 근사 방법 제공 최소 정칙성 가정 하에서 명시적 오차 추정 제시 이산 분수 라플라시안 연산자의 완전한 이론 틀 확립 , 점별 공식, 핵 함수 추정 및 기본 성질 포함이산 Hölder 공간에서의 정칙성 이론 증명 , Schauder 추정 및 최대값 원리 포함분수 이산 Sobolev 임베딩과 Poincaré 부등식 확립 , 비국소 Dirichlet 문제에 대한 이론적 기초 제공이산 Hardy-Littlewood-Sobolev 부등식 제시 , 이산 분수 적분 이론 완성연속 분수 라플라시안으로의 균일 수렴성 증명 , Hölder 공간에서 명시적 오차 추정 제시연속 Poisson 문제 해의 이산 근사 이론 확립 , L ∞ L^\infty L ∞ 노름 하에서의 오차 추정 제공격자 Z h = { h j : j ∈ Z } \mathbb{Z}_h = \{hj : j \in \mathbb{Z}\} Z h = { hj : j ∈ Z } 위의 비국소 이산 확산 방정식 연구:
( − Δ h ) s u = f (-\Delta_h)^s u = f ( − Δ h ) s u = f
여기서 u , f : Z h → R u, f: \mathbb{Z}_h \to \mathbb{R} u , f : Z h → R , 0 < s < 1 0 < s < 1 0 < s < 1 , h > 0 h > 0 h > 0 .
이산 라플라시안의 분수 거듭제곱은 반군 방법으로 정의됨:
( − Δ h ) s u j = 1 Γ ( − s ) ∫ 0 ∞ ( e t Δ h u j − u j ) d t t 1 + s (-\Delta_h)^s u_j = \frac{1}{\Gamma(-s)} \int_0^\infty \left(e^{t\Delta_h}u_j - u_j\right) \frac{dt}{t^{1+s}} ( − Δ h ) s u j = Γ ( − s ) 1 ∫ 0 ∞ ( e t Δ h u j − u j ) t 1 + s d t
여기서 w j ( t ) = e t Δ h u j w_j(t) = e^{t\Delta_h}u_j w j ( t ) = e t Δ h u j 는 반이산 열 방정식의 해:
{ ∂ t w j = Δ h w j , in Z h × ( 0 , ∞ ) w j ( 0 ) = u j , in Z h \begin{cases}
\partial_t w_j = \Delta_h w_j, & \text{in } \mathbb{Z}_h \times (0,\infty) \\
w_j(0) = u_j, & \text{in } \mathbb{Z}_h
\end{cases} { ∂ t w j = Δ h w j , w j ( 0 ) = u j , in Z h × ( 0 , ∞ ) in Z h
정리 1.1 : 0 < s < 1 0 < s < 1 0 < s < 1 과 u ∈ ℓ s u \in \ell_s u ∈ ℓ s 에 대해:
( − Δ h ) s u j = ∑ m ∈ Z , m ≠ j ( u j − u m ) K s h ( j − m ) (-\Delta_h)^s u_j = \sum_{m \in \mathbb{Z}, m \neq j} (u_j - u_m) K_s^h(j-m) ( − Δ h ) s u j = ∑ m ∈ Z , m = j ( u j − u m ) K s h ( j − m )
여기서 이산 핵 함수는:
K s h ( m ) = 4 s Γ ( 1 / 2 + s ) π ∣ Γ ( − s ) ∣ ⋅ Γ ( ∣ m ∣ − s ) h 2 s Γ ( ∣ m ∣ + 1 + s ) K_s^h(m) = \frac{4^s \Gamma(1/2 + s)}{\sqrt{\pi}|\Gamma(-s)|} \cdot \frac{\Gamma(|m| - s)}{h^{2s}\Gamma(|m| + 1 + s)} K s h ( m ) = π ∣Γ ( − s ) ∣ 4 s Γ ( 1/2 + s ) ⋅ h 2 s Γ ( ∣ m ∣ + 1 + s ) Γ ( ∣ m ∣ − s )
상수 0 < c s ≤ C s 0 < c_s \leq C_s 0 < c s ≤ C s 가 존재하여:
c s h 2 s ∣ m ∣ 1 + 2 s ≤ K s h ( m ) ≤ C s h 2 s ∣ m ∣ 1 + 2 s \frac{c_s}{h^{2s}|m|^{1+2s}} \leq K_s^h(m) \leq \frac{C_s}{h^{2s}|m|^{1+2s}} h 2 s ∣ m ∣ 1 + 2 s c s ≤ K s h ( m ) ≤ h 2 s ∣ m ∣ 1 + 2 s C s
수정된 Bessel 함수를 이용한 반이산 열 핵의 정확한 공식:
G ( m , t ) = e − 2 t I m ( 2 t ) G(m,t) = e^{-2t}I_m(2t) G ( m , t ) = e − 2 t I m ( 2 t )
( − Δ h ) s (-\Delta_h)^s ( − Δ h ) s 를 반이산 퇴화 타원 확장 문제의 Dirichlet-to-Neumann 사상으로 특성화하여 국소 및 비국소 연산자 간의 연결 확립.
이산 Hölder 공간 C h k , α C_h^{k,\alpha} C h k , α 의 완전한 이론 확립, 분수 이산 라플라시안 연산자의 정칙성 증명:
정리 1.5 : u ∈ C h k , α u \in C_h^{k,\alpha} u ∈ C h k , α 이고 2 s < α 2s < \alpha 2 s < α 이면, ( − Δ h ) s u ∈ C h k , α − 2 s (-\Delta_h)^s u \in C_h^{k,\alpha-2s} ( − Δ h ) s u ∈ C h k , α − 2 s .
본 논문은 주로 이론 연구로서 다음 방식으로 결과를 검증:
극한 거동 검증 :s → 0 + s \to 0^+ s → 0 + 일 때: ( − Δ h ) s u j → u j (-\Delta_h)^s u_j \to u_j ( − Δ h ) s u j → u j s → 1 − s \to 1^- s → 1 − 일 때: ( − Δ h ) s u j → − Δ h u j (-\Delta_h)^s u_j \to -\Delta_h u_j ( − Δ h ) s u j → − Δ h u j 일관성 검사 : 격자 세분화 시 이산 연산자와 연속 연산자의 일관성 검증핵 함수 성질 : 핵 함수의 대칭성, 양성 및 점근 거동 검증수렴 차수 : Hölder 노름 하에서의 수렴 속도오차 추정 : L ∞ L^\infty L ∞ 노름 하에서의 명시적 오차 한계정칙성 지표 : 이산 Hölder 공간에서의 노름 추정정리 1.7 : U ∈ C 0 , α U \in C^{0,\alpha} U ∈ C 0 , α 이고 2 s < α 2s < \alpha 2 s < α 에 대해:
∥ ( − Δ h ) s ( r h U ) − r h ( ( − Δ ) s U ) ∥ ℓ h ∞ ≤ C [ U ] C 0 , α h α − 2 s \|(-\Delta_h)^s(r_h U) - r_h((-\Delta)^s U)\|_{\ell^\infty_h} \leq C[U]_{C^{0,\alpha}} h^{\alpha-2s} ∥ ( − Δ h ) s ( r h U ) − r h (( − Δ ) s U ) ∥ ℓ h ∞ ≤ C [ U ] C 0 , α h α − 2 s
정리 1.8 : F ∈ C 0 , α F \in C^{0,\alpha} F ∈ C 0 , α 가 컴팩트 지지를 가지고, U U U 가 연속 Poisson 문제 ( − Δ ) s U = F (-\Delta)^s U = F ( − Δ ) s U = F 의 해, u u u 가 대응하는 이산 Dirichlet 문제의 해이면:
∥ u − r h U ∥ ℓ h ∞ ( B R h ) ≤ C ∥ F ∥ C 0 , α R 2 s h α \|u - r_h U\|_{\ell^\infty_h(B_R^h)} \leq C\|F\|_{C^{0,\alpha}} R^{2s} h^\alpha ∥ u − r h U ∥ ℓ h ∞ ( B R h ) ≤ C ∥ F ∥ C 0 , α R 2 s h α
1 < p < q < ∞ 1 < p < q < \infty 1 < p < q < ∞ 이고 1 / q ≤ 1 / p − 2 s 1/q \leq 1/p - 2s 1/ q ≤ 1/ p − 2 s 에 대해:
∥ ( − Δ h ) − s f ∥ ℓ h q ≤ C h 1 / p − 2 s − 1 / q ∥ f ∥ ℓ h p \|(-\Delta_h)^{-s}f\|_{\ell^q_h} \leq \frac{C}{h^{1/p-2s-1/q}} \|f\|_{\ell^p_h} ∥ ( − Δ h ) − s f ∥ ℓ h q ≤ h 1/ p − 2 s − 1/ q C ∥ f ∥ ℓ h p
최적 수렴 차수 : 최소 정칙성 가정 하에서 이론적으로 최적의 수렴 속도 달성균일 추정 : 모든 상수가 격자 크기 h h h 에 무관경계 효과 : 이산 최대값 원리에서 나타나는 R 2 s R^{2s} R 2 s 인수는 비국소 연산자의 본질적 특성 반영연속 분수 라플라시안 연산자 이론 : Caffarelli-Silvestre 확장 문제, Silvestre 정칙성 이론수치 방법 : 유한 차분, 유한 요소, 스펙트럼 방법 등 이산화 기술확률론적 해석 : 점프 과정 및 Lévy 과정과의 연결완전한 이산 분수 라플라시안 연산자 이론 체계 최초 구축 엄격한 오차 분석 및 수렴성 증명 제공 L ∞ L^\infty L ∞ 노름에서 명시적 오차 추정 제시, 이는 해당 분야에서 최초이산 분수 라플라시안 연산자의 완전한 수학 이론 구축 이산 모델에서 연속 모델로의 균일 수렴성 증명 분수 Poisson 문제에 대한 효과적인 수치 근사 방법 제공 차원 제한 : 주요 결과가 1차원에 집중, 고차원 확장에는 새로운 기술 필요매개변수 범위 : 일부 결과가 s < 1 / 2 s < 1/2 s < 1/2 를 요구하여 응용 범위 제한격자 규칙성 : 방법이 균일 격자에 적용 가능, 비균일 격자는 추가 고려 필요고차원 확장 : 다차원 이산 분수 라플라시안 연산자 이론 개발최적 정칙성 : 이산 Dirichlet 문제 해의 최적 Hölder 정칙성 결정적응형 방법 : 오차 추정에 기반한 적응형 격자 방법 개발이론적 완전성 : 기본 정의에서 응용까지 엄격한 논증을 갖춘 체계적이고 완전한 이론 틀 구축기술적 혁신 : 반군 방법과 Bessel 함수 기술의 결합으로 정확한 핵 함수 표현식 제시실용적 가치 : 구현 가능한 수치 방법과 명확한 오차 제어 제공작문 품질 : 논문 구조가 명확하고 증명이 엄밀하며 기술적 세부사항이 적절히 처리됨계산 복잡성 : 핵 함수가 감마 함수와 Bessel 함수를 포함하여 계산 비용이 높음응용 검증 : 이론 예측을 검증하는 구체적인 수치 실험 부재비교 분석 : 다른 이산화 방법과의 비교가 충분하지 않음학술적 기여 : 비국소 연산자의 이산화 이론에 중요한 기초 마련응용 전망 : 분수 미분 방정식의 수치 해법에 이론적 지원 제공재현성 : 이론 결과가 완전하여 후속 연구자의 추진 및 확장 용이분수 확산 방정식의 수치 해법 비국소 모델의 이산화 분석 점프 과정의 수치 시뮬레이션 영상 처리의 비국소 필터 설계 논문은 해당 분야의 중요 문헌을 인용하고 있으며, 다음을 포함:
Caffarelli & Silvestre (2007): 분수 라플라시안 연산자의 확장 문제 Silvestre (2007): 분수 라플라시안 연산자의 정칙성 이론 Savin & Valdinoci (2011): 비국소 변분 모델의 Sobolev 부등식 Stein & Wainger (2000): 조화 분석의 이산 유사 종합 평가 : 이는 이산 분수 라플라시안 연산자의 수학 이론을 체계적으로 구축하고 비국소 미분 방정식의 수치 분석에 중요한 이론적 기초를 제공하는 고품질의 이론 수학 논문이다. 논문의 기술적 깊이와 이론적 완전성이 매우 높은 수준에 도달하였으며, 해당 분야의 발전을 추진하는 데 중요한 의의를 가진다.