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Дата публикации : 21 марта 2018 г. (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 может быть охарактеризован как оператор Дирихле-Неймана полудискретной вырожденной эллиптической задачи расширения. Установлены свойства регулярности в дискретных пространствах Гёльдера и оценки Шаудера, а также существование и единственность решений нелокальной задачи Дирихле. Введены отрицательные степени (фундаментальные решения) u = ( − Δ h ) − s f u = (-\Delta_h)^{-s}f u = ( − Δ h ) − s f и доказаны соответствующие дискретные неравенства Харди-Литтлвуда-Соболева. В качестве приложения проанализирована сходимость дробного дискретного оператора Лапласа к непрерывному дробному оператору Лапласа в пространствах Гёльдера при h → 0 h \to 0 h → 0 , и доказано, что решения задачи Пуассона для дробного оператора Лапласа могут быть аппроксимированы решениями задачи Дирихле для дробного дискретного оператора Лапласа.
Дробный оператор Лапласа как положительная степень классического оператора Лапласа существует уже давно в потенциальной теории, гармоническом анализе, дробном исчислении, функциональном анализе и теории вероятностей. За последнее десятилетие нелокальные уравнения диффузии, включающие дробный оператор Лапласа, стали одной из наиболее актуальных тем исследований в области уравнений в частных производных.
Трудности дискретизации : Нелокальность и сингулярность оператора ( − Δ ) s (−\Delta)^s ( − Δ ) s создают основные сложности при численной аппроксимацииОтсутствие теории : Отсутствует систематическая теоретическая база для дискретного дробного оператора ЛапласаАнализ сходимости : Требуются строгие теоретические гарантии сходимости дискретной модели к непрерывнойПостроить полную систему теории дискретного дробного оператора Лапласа Предоставить эффективные методы дискретной аппроксимации непрерывных задач Дать явные оценки погрешности при минимальных предположениях о регулярности Установлена полная теоретическая база дискретного дробного оператора Лапласа , включая поточечные формулы, оценки ядра и фундаментальные свойстваДоказана теория регулярности в дискретных пространствах Гёльдера , включая оценки Шаудера и принцип максимумаУстановлены дискретные вложения Соболева дробного порядка и неравенства Пуанкаре , обеспечивающие теоретическую основу для нелокальной задачи ДирихлеПолучены дискретные неравенства Харди-Литтлвуда-Соболева , совершенствующие теорию дискретного дробного интегрированияДоказана равномерная сходимость к непрерывному дробному оператору Лапласа с явными оценками погрешности в пространствах ГёльдераУстановлена теория дискретной аппроксимации решений непрерывной задачи Пуассона с оценками погрешности в норме 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 , на Z h × ( 0 , ∞ ) w j ( 0 ) = u j , на Z h \begin{cases}
\partial_t w_j = \Delta_h w_j, & \text{на } \mathbb{Z}_h \times (0,\infty) \\
w_j(0) = u_j, & \text{на } \mathbb{Z}_h
\end{cases} { ∂ t w j = Δ h w j , w j ( 0 ) = u j , на Z h × ( 0 , ∞ ) на 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
Использование модифицированных функций Бесселя для получения точной формулы полудискретного теплового ядра:
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 как оператора Дирихле-Неймана полудискретной вырожденной эллиптической задачи расширения, устанавливающая связь между локальными и нелокальными операторами.
Построена полная теория дискретных пространств Гёльдера 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 Проверка согласованности : Верификация согласованности дискретного оператора с непрерывным при измельчении сеткиСвойства ядра : Проверка симметрии, положительности и асимптотического поведения ядраПорядок сходимости : Скорость сходимости в норме ГёльдераОценки погрешности : Явные границы погрешности в норме L ∞ L^\infty L ∞ Показатели регулярности : Оценки норм в дискретных пространствах ГёльдераТеорема 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 — решение непрерывной задачи Пуассона ( − Δ ) s U = F (-\Delta)^s U = F ( − Δ ) s U = F , u u u — решение соответствующей дискретной задачи Дирихле, тогда:
∥ 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 в дискретном принципе максимума отражает существенные особенности нелокального оператораТеория непрерывного дробного оператора Лапласа : Задача расширения Каффарелли-Сильвестра, теория регулярности СильвестраЧисленные методы : Методы конечных разностей, конечных элементов, спектральные методы и другие техники дискретизацииВероятностная интерпретация : Связь со скачкообразными процессами и процессами ЛевиВпервые построена полная система теории дискретного дробного оператора Лапласа Предоставлены строгие анализ погрешности и доказательства сходимости Получены явные оценки погрешности в норме L ∞ L^\infty L ∞ , что является первым результатом в данной области Построена полная математическая теория дискретного дробного оператора Лапласа Доказана равномерная сходимость дискретной модели к непрерывной Предоставлены эффективные методы численной аппроксимации для задачи Пуассона с дробным оператором Ограничение по размерности : Основные результаты сосредоточены на одномерном случае, обобщение на высшие размерности требует новых техникДиапазон параметров : Некоторые результаты требуют s < 1 / 2 s < 1/2 s < 1/2 , что ограничивает область примененияРегулярность сетки : Методы применимы к равномерным сеткам, неравномерные сетки требуют дополнительного рассмотренияОбобщение на высшие размерности : Развитие теории многомерного дискретного дробного оператора ЛапласаОптимальная регулярность : Определение оптимальной регулярности Гёльдера решений дискретной задачи ДирихлеАдаптивные методы : Разработка адаптивных методов сеток на основе оценок погрешностиПолнота теории : Построена систематическая и полная теоретическая база от фундаментальных определений до приложений с строгими доказательствамиТехнические инновации : Метод полугрупп в сочетании с техникой функций Бесселя дает точные выражения для ядраПрактическая ценность : Предоставлены реализуемые численные методы и явный контроль погрешностиКачество изложения : Статья имеет четкую структуру, строгие доказательства и тщательную обработку технических деталейВычислительная сложность : Ядро включает гамма-функции и функции Бесселя, что требует значительных вычислительных затратОтсутствие численных экспериментов : Не хватает конкретных численных экспериментов для верификации теоретических предсказанийНедостаточное сравнение : Сравнение с другими методами дискретизации проведено недостаточно полноНаучный вклад : Работа закладывает важную основу для теории дискретизации нелокальных операторовПерспективы применения : Предоставляет теоретическую поддержку для численного решения дробных дифференциальных уравненийВоспроизводимость : Полнота теоретических результатов облегчает последующие исследования и расширенияЧисленное решение дробных уравнений диффузии Анализ дискретизации нелокальных моделей Численное моделирование скачкообразных процессов Проектирование нелокальных фильтров в обработке изображений Статья цитирует важные работы в данной области, включая:
Caffarelli & Silvestre (2007): Задача расширения для дробного оператора Лапласа Silvestre (2007): Теория регулярности дробного оператора Лапласа Savin & Valdinoci (2011): Неравенства Соболева для нелокальных вариационных моделей Stein & Wainger (2000): Дискретные аналоги в гармоническом анализе Общая оценка : Это высокачественная теоретическая математическая работа, систематически построившая математическую теорию дискретного дробного оператора Лапласа и предоставившая важную теоретическую основу для численного анализа нелокальных дифференциальных уравнений. Техническая глубина и полнота теории работы достигают очень высокого уровня и имеют важное значение для развития данной области.