Space-time fractional equations
European Congress of Mathematics, Sevilla
Aim
\[ \text{``}{\color{blue}\partial_t^\alpha}\text{''} u + \text{``}{\color{red} (-\Delta)^s} \text{''} u = f \]
Fractional in time: “\({\color{blue} \partial_t}\)”
Fractional in time
Let \({\color{blue} I}: C([0,T]) \to C([0,T])\) the Riemann integral \[ {\color{blue} I} f(t) = \int_0^t f(s) ds. \] Then \[\begin{align*} {\color{blue} I}^2 f (t) &= \int_0^t I_1f (s) ds = \int_0^t \int_0^s f(\sigma) d\sigma ds = \int_0^t f(\sigma) \int_s^t ds d \sigma \\ &= \int_0^t f(\sigma) (t-s) d \sigma, \\ {\color{blue} I}^n f (t) &= \int_0^t f(\sigma) K_n(t-s) d \sigma, \qquad K_n(t) = \frac{t^{n-1}}{(n-1)!}. \end{align*}\]
We define the Riemann-Liouville integral for \(\alpha > 0\) \[ \iRL f (t) = \int_0^t f(\sigma) K_\alpha(t-s) d \sigma, \qquad {\color{blue} K_\alpha}(t) = \frac{t^{\alpha-1}}{\Gamma(\alpha)}. \]
Key properties
\[ \iRL f (t) = \int_0^t f(\sigma) K_\alpha(t-s) d \sigma, \qquad K_\alpha(t) = \frac{t^{\alpha-1}}{\Gamma(\alpha)}. \]
Semigroup: \(\iRL {\color{blue} I_\beta} = {\color{blue} I_{\alpha + \beta}}\)
Monomial, for \(\beta > -1\) \[ \iRL t^\beta =\frac{\Gamma(\beta+1)}{\Gamma(\beta+1+\alpha)} t^{\beta+\alpha}. \]
Laplace transform symbol \[ {\mathfrak L}[\iRL f](s) = s^{-\alpha} {\mathfrak L}[f](s). \]
Fractional derivatives
A fractional derivative should have symbol \(s^\alpha\). Thus we get two choices:
Riemann-Liouville derivative \[ \dRL u = \partial_t {\color{blue} I_{1-\alpha}} u \]
Caputo derivative \[ \dC u = {\color{blue} I_{1-\alpha}} \partial_t u \]
Using the properties of \(I_\alpha\) and \(\partial_t\) it follows directly that \[\begin{align*} \dC t^\beta &= \dRL t^\beta = C(\beta, \alpha) t^{\beta - \alpha} \qquad \text{if } \beta > 0. \\ \dC 1 &= 0 \\ \dRL 1 &= \frac{t^{\alpha - 1}}{\Gamma(1-\alpha)} \ne 0 \\ \dRL t^{\alpha - 1} &= 0. \end{align*}\]
Integral formulation
\[\begin{align*} \dC u = f &\iff u = u(0) + \iRL f \\ \dRL v = f &\iff v = C_\alpha {\color{blue} I_{1-\alpha}} v(0) t^{\alpha-1} + \iRL f \end{align*}\]
Who cares? Stochastic volatility model for asset pricing.
Let \(S_t\) be the discounted spot price of an asset. We model \[ dS_t = S_t \sqrt{V_t} dW_t, % \qquad \text{i.e. } S_t = S_0 + \int_0^t S_s \sqrt{V_s} dW_s, \] where \(W_t\) be a Brownian motion in risk-neutral probability measure.
Consider \(B_t\) a second Brownian motion with \(\langle dW_t, dB_t \rangle = \rho\).
- Black-Scholes (1973): \(V_t = \sigma^2\) constant.
Heston (1993): \[ V_t = V_0 + \int_0^t \kappa( \theta - V_s ) ds + \int_0^t \sigma \sqrt{V_t} d B_s, \]
Rough Heston (2014) model: \[ V_t = V_0 + \int_0^t \kappa( \theta - V_s ) {\color{blue} K_\alpha (t-s)} ds + \int_0^t \sigma \sqrt{V_t}{\color{blue} K_\alpha (t-s)} d B_s. \] Several problem related to RH involve fractional ODEs.
Does \(\alpha\) matter?
Arbitrage-free price of European call option: \(e^{-rT} \mathbb E[(S_T - K)_+]\)
Fractional in space: “\({\color{red} (-\Delta)^s}\)”
Fractional Laplacian
\[ {\color{red} (-\Delta)^s} u (x) = {C(d,s)} \,\, \underbrace{\lim_{\varepsilon \to 0} \int_{\mathbb R^d \setminus B_\varepsilon (x) }}_{\text{ P.V. }\int _{\mathbb R^d}} \frac{u(x) - u(y)}{|x-y|^{d+2s}} d y. \]
Equivalent definitions (up-to ten [Kwasnicki, 2017])
The unique operator such that
\[ \mathfrak F[ {\color{red} (-\Delta)^s} u ] = |\xi|^{2s} \mathfrak F[u] \]
i.e. the spectral fractional power of \(-\Delta\).
The infinitesimal generator of a Lévy process \(X_h\), characterised by “long jumps”
\[ {\color{red} (-\Delta)^s} u (x) = \lim_{h \to 0} \frac{ \mathbb E[ f(x) - f(x - X_h) ] }{ h } \]
Fractional Elliptic/Heat/Porous Medium in \(\mathbb R^d\)
Analysis of the operator: [Caffarelli & Silvestre, 2007] extension
Obstacle: [Silvestre, 2007], [Barrios, Figalli & Ros-Oton, 2018]
Heat: We can compute estimates for the heat kernel \(K(t,\cdot) = \mathfrak F^{-1} [e^{-t|\xi|^{2s}}]\).
Porous Medium: [De Pablo, Quirós, Rodríguez & Vázquez, 2011], [Vázquez & Volzone, 2014]
Connection to geometry: [Chang & González, 2011]
Semilinear elliptic: [Nguyen & Véron, 2018],
[Chan, González, Huang, Mainini & Volzone, 2020]
There is a different Porous-Medium-type fractional equation \[ \partial_t u = \mathrm{div} (u^{m-1} \nabla (-\Delta)^{-s} u) \]
We will not discuss it:
\(m = 2\) [Caffarelli & Vazquez, 2011], \(m \ne 2\) [Stan, Del Teso & Vázquez, 2016]
Fractional LaplacianS
in bounded domains
Restricted fractional Laplacian (RFL). Singular integral operator: \[ {\color{red} (-\Delta)^s_{\mathrm {RFL}}} u (x)= C(d,s) \, \mathrm{P.V.}\int_{ \mathbb R^d } \frac{u(x) - u(y)}{|x-y|^{d+2s}} \; dy. \] If we work only in \(\Omega\), we must prescribe \(u\) in \(\Omega^c = \mathbb R^d \setminus \Omega\).
Censored fractional Laplacian (CFL). For \(s > \frac 1 2\) \[\begin{equation} {\color{red} (-\Delta)^s_{\mathrm{CFL}}} u (x) = C(d,s) \, \mathrm{P.V.} \int_{ \Omega } \frac{u(x) - u(y)}{|x-y|^{n+2s}} \; dy, \end{equation}\]
Spectral fractional Laplacian (SFL). Operational power. \(-\Delta \varphi_m = \lambda_m \varphi_m \textrm{ in } \Omega\), \(\varphi_m = 0 \textrm{ on } \partial \Omega.\) one defines \[ u(x) = \sum_{m=1}^{+\infty} {u_m} \varphi_m(x) \qquad \longmapsto \qquad {\color{red} (-\Delta)_{\mathrm{SFL}}^s} u (x) = \sum_{m=1}^{+\infty} \lambda_m^s u_m \varphi_m (x). \] The “boundary condition” is \(u = 0\) on \(\partial \Omega\).
Green kernels
The solution of the Laplace equation \[ \begin{dcases} {\color{red} \mathcal L} u = f & \Omega \\ u = 0 & \partial \Omega \text{ or } \mathbb R^d \setminus \overline \Omega \\ \end{dcases} \]
can be written via the Green function \[ u(x) = {\color{red} \mathrm G}[f] (x) = \int_\Omega {\color{red} \mathbb G}(x,y) f(y) dy. \]
The probabilistic approach provides in each of our cases a similar shape \[ {\color{red} \mathbb G}(x,y) \asymp \frac{1}{|x-y|^{d-2s}} \left(1 \wedge \frac{\delta(x)}{|x-y|} \right)^\gamma\left(1 \wedge \frac{\delta(y)}{|x-y|} \right)^\gamma \]
Since we deal only with self-adjoint \({\color{red} \mathcal L}\), then \({\color{red} \mathbb G}(x,y) = {\color{red} \mathbb G}(y,x)\).
The probabilistic approach provides in each of our cases a similar shape \[ {\color{red} \mathbb G}(x,y) \asymp \frac{1}{|x-y|^{d-2s}} \left(1 \wedge \frac{\delta(x)}{|x-y|} \right)^\gamma\left(1 \wedge \frac{\delta(y)}{|x-y|} \right)^\gamma \]
And \(\gamma \in (0,1]\) depends on the setting
RFL: \(\gamma = s\).
SFL: \(\gamma = 1\)
CFL: \(\gamma = 2s-1\), only \(s \in (\frac 1 2, 1)\).
Parabolic problems with homogeneous Dirichlet data
Linear fractional ODEs
Let \(E_{\alpha, \beta}\) be the Mittag-Leffler function
We define
\(E_\alpha(z) = E_{\alpha,\alpha}(z)\) and
\(P_{\alpha}(t;\lambda) = t^{\alpha-1} E_{\alpha,\alpha}(-\lambda t^\alpha)\).
Using Laplace transform \[\begin{align*} &\dC u = -\lambda u + f \\ &\qquad \implies u (t) = E_\alpha (-\lambda t^\alpha) u(0) + \int_0^t P_{\alpha}(t-\tau;\lambda) f(\tau) d \tau.\\ &\dRL v = -\lambda v + g \\ &\qquad \implies v(t) = P_{\alpha}(t;\lambda) I_{1-\alpha}v (0) + \int_0^t P_{\alpha} (t-\tau;\lambda) g(\tau) d \tau \end{align*}\]
Duhamel’s formulas for space-time problems
(homogeneous Dirichlet conditions)
Using the spectral decomposition, we write \[\begin{align*} &\dC u = -{\color{red} \mathcal L} u + f \\ &\qquad \implies u (t) = E_\alpha (-t^\alpha {\color{red} \mathcal L}) u(0) + \int_0^t P_{\alpha}(t-\tau;{\color{red} \mathcal L}) f(\tau) d \tau.\\ &\dRL v = -{\color{red} \mathcal L} v + g \\ &\qquad \implies v(t) = P_{\alpha}(t;{\color{red} \mathcal L}) I_{1-\alpha}v (0) + \int_0^t P_{\alpha} (t-\tau;{\color{red} \mathcal L}) g(\tau) d \tau \end{align*}\]
Here we are using the notation for \(L\varphi_i = \lambda_i \varphi_i\) and \(F:\mathbb R_+\to \mathbb R\) \[ F(L) \sum_i a_i \varphi_i = \sum_i a_i F(\lambda_i) \varphi_i. \]
Singular boundary data for
Laplace equation
Estimates for interior data for \((-\Delta)^s u = f\)
Given that \({\color{red} \mathbb G}(x,y) \asymp \frac{1}{|x-y|^{d-2s}} \left(1 \wedge \frac{\delta(x)}{|x-y|} \right)^\gamma\left(1 \wedge \frac{\delta(y)}{|x-y|} \right)^\gamma\)
For \(K \Subset \Omega\) \[ \int_K |{\color{red} \mathrm G}(f)| \le C_K \int_\Omega |f| \delta^\gamma. \]
So \({\color{red} \mathrm G}: L^1 (\Omega, \delta^\gamma) \to L^1_{loc} (\Omega)\)
If \(f \ge 0\) then, using the kernel estimates \[ \tag{Hopf} {\color{red} \mathrm G}[f] (x) \ge c \delta^\gamma(x) \int_\Omega f(y) \delta(y)^\gamma dy. \]
If \(0\le f \notin L^1 (\Omega, \delta^\gamma)\), then \({\color{red} \mathrm G}[f] = +\infty\).
Large solutions with boundary blow-up for \({\color{red} \mathcal L}\)
[Abatangelo, GC & Vázquez, 2022] denoting \(\delta(x) = \mathrm{dist}(x, \partial \Omega)\)
We can compute \({\color{red} \mathrm G}[\delta^\beta] \asymp \delta^\alpha\) with
\[ {\color{red} \mathrm G}[\delta^\beta] \asymp \begin{dcases} \delta^\gamma & \beta +2s > \gamma \\ \delta^{\gamma} \log|\delta| & \beta +2s = \gamma \\ \delta^{\beta + 2s} & \beta +2s < \gamma \\ &\text{ and } \beta > -1-\gamma \\ \infty & \beta \le -1-\gamma \end{dcases} \]
\(\infty \not \equiv {\color{red} \mathrm G}[\delta^\beta] \notin L^\infty \iff -1-\gamma < \beta < -2s\).
For \(-\Delta\) we have \(\gamma = 1 = s\). Either \({\color{red} \mathrm G}[\delta^\beta] \in L^\infty\) or \({\color{red} \mathrm G}[\delta^\beta] \equiv + \infty\)
Singular solutions of the elliptic “Dirichlet” problem
The harmonic singular solution
Take \(s \in (0,1)\)
\[ u(x) = \begin{cases} (1-|x|^2)^{s-1} & \text{if } |x|<1 \\ 0 & \text{if } |x| \ge 1 \\ \end{cases} \]
satisfies
\[ (-\Delta)_{\mathrm{RFL}}^s u(x) = 0 \qquad \text{if } |x| < 1. \]
No right-hand side? There is an extra “boundary” condition
Martin problem
Because of the scaling we consider \[ f_m (x) = \frac{|\partial \Omega|}{ |A_m| } \frac{\chi_{A_m} (x)}{\delta (x)^\gamma} \qquad \text{ where } A_m = \left\{ x: \frac 1 m \delta(x) < \frac 2 m \right\} \]
[Abatangelo, GC & Vázquez, 2022] We show that \(\mathrm{supp}(f_m) \to \partial \Omega\) and \({\color{red} \mathrm G}[f_m] \to u^\star \ne 0\).
In fact, for any \(h \in L^1 (\partial \Omega)\) the problem
\[ \begin{dcases} {\color{red} \mathcal L} u = 0 & \Omega \\ u = 0 & \mathbb R^d \setminus \overline \Omega \\ \lim_{x \to z} \frac{u(x)}{u^\star (x)} = h(z) & \text{for all } z \in \partial \Omega \end{dcases} \]
admits a solution
(constructed in similar fashion).
We show that \(u^\star \asymp \delta^{2s-\gamma -1 }.\)
RFL: \(s-1\)
SFL: \(2(s-1)\)
CFL: \(0\), i.e. bounded.
\(u_j \to u^\star\)
Singular boundary data
in other problems
Schrödinger-type problems
[Chan, GC & Vázquez, 2021] : we study
\[ \begin{cases} {\color{red} \mathcal L} u = \lambda u + f & \text{in }\Omega, \\ \frac{u}{u^\star} = h & \text{on } \partial \Omega. \end{cases} \]
We use the eigenfunction expansion to understand the blow-up as
\[\lambda \to \mathrm{spectrum} ({\color{red} \mathcal L}).\]
Parabolic problem
\[ \begin{dcases} \frac{\partial u}{\partial t} + {\color{red} \mathcal L} u = f & \text{in } (0,\infty) \times \Omega, \\ \frac{u}{u^\star} = h &\text{on } (0,\infty) \times \partial \Omega \\ u = u_0 & \text{at } t = 0. \end{dcases} \]
We construct the heat kernel, and localise \(f\) to the boundary.
Fractional in space and time
[Chan, Gómez-Castro & Vázquez, 2024]
\[ \begin{dcases} \partial^\alpha_t u + {\color{red} \mathcal L} u = f & \text{in } (0,\infty) \times \Omega, \\ \frac{u}{u^\star} = h &\text{on } (0,\infty) \times \partial \Omega \\ + \text{with suitable initial condition} \end{dcases} \]
where \(\partial^\alpha_t\) is either in the Caputo or Riemann-Liouville setting.
Questions?
Now,
coffee,
Bibliography
Appendix
Laplace equation
Laplace equation in \(\mathbb R^d\)
\[ \begin{cases} {\color{red} (-\Delta)^s} u (x) = f(x) & \textrm{for all } x \in \mathbb R^d \\ u(x) \to 0 & \textrm{as } |x| \to \infty \end{cases} \]
Transforming the equation \(|\xi|^{2s} \widehat u(\xi) = \widehat f(\xi)\).
Formally, we can recover \(u = \mathfrak F^{-1}[|\xi|^{-2s}] * f\).
For \(d > 2s\), \(\mathfrak F^{-1}[|\xi|^{-2s}]\) is the Riesz kernel so we get
\[ u(x) = C(d,-2s) \int_{\mathbb R^d} \frac{f(y)}{|x-y|^{d-2s}} dy. \]
For any \(d, 2s\) we can also look for this solution by the minisation of the energy functional
\[ J(u) = \frac {C(d,s)} 2 \int_{\mathbb R^d} \int_{\mathbb R^d } \frac{|u(x) - u(y)|^2}{|x-y|^{d+2s}} dx \, dy - \int_\Omega f u \]
in homogeneous fractional Sobolev spaces [Brasco, Gómez-Castro & Vázquez, 2021].
Laplace equation in a bounded domain
\[ \begin{dcases} {\color{red} \mathcal L} u = f & \Omega \\ u = 0 & \partial \Omega \text{ or } \mathbb R^d \setminus \overline \Omega \\ \end{dcases} \]
We observe
\({\color{red} \mathcal L} = (-\Delta)^s_{\mathrm{RFL}}, (-\Delta)^s_{\mathrm{CFL}}\) are sub-differentials of energies \[ J(u) = \int_A \int_A \frac{|u(x) - u(y)|^2 }{|x-y|^{d+2s}} dx \, dy. \]
\({\color{red} \mathcal L} = (-\Delta)^s_{\mathrm{SFL}}\) is just a power.
The inverse is naturally \((-\Delta_\Omega)^{-s}\), and in works between “powers” of \(H_0^1(\Omega)\).
Self-adjoint compact operators. Furthermore \(\lambda_1 > 0\).
Nice theory of energy solutions.
Higher regularity: RFL [Ros-Oton & Serra, 2014], CFL [Chen, 2018], [Fall & Ros-Oton, 2021]
Numerics
Finite Elements
For the RFL and CFL integrate by parts: [Acosta & Borthagaray, 2017], [Borthagaray, Nochetto & Salgado, 2019]
For the SFL use semigroup formula: [Cusimano, del Teso & Gerardo-Giorda, 2020]
Finite Difference schemes \[ (-\Delta)^s_h u(x) = \sum_{i \in \mathbb Z^d} \omega_{i} (u(x) - u(x+ih)) \]
[Jakobsen, Karlsen & La Chioma, 2008], [Biswas, Jakobsen & Karlsen, 2010],
\(d = 1\): [Huang & Oberman, 2014].
See also: [Del Teso, Endal & Jakobsen, 2018]
Discontinuous Galerkin: [Cifani, Jakobsen & Karlsen, 2011]
Singular boundary condition
Satisfying the boundary equation.
Continuous data
When \(h \in \mathcal C(\partial \Omega)\) it is a direct computation that
\(\displaystyle \frac{\mathcal M[h](x) }{u^\star (x)} - h (\zeta_0)\) \(\displaystyle = \frac{\int_{\partial \Omega} \mathbb M(x,\zeta) h(\zeta) d \zeta }{ \int_{\partial \Omega} \mathbb M (x, \zeta) d \zeta } - \frac{\int_{\partial \Omega} \mathbb M(x,\zeta) d \zeta }{ \int_{\partial \Omega} \mathbb M (x, \zeta) d \zeta } h(\zeta_0)\)
Thus
\[ \left| \frac{\mathcal M[h] (x) }{u^\star (x)} - h (\zeta_0) \right| \le \int_{\partial \Omega} \frac{ \mathbb M(x,\zeta) }{\int_{\partial \Omega} \mathbb M (x, z) d z} | h(\zeta) - h(\zeta_0)| d \zeta \]
For any \(\varepsilon > 0\), due the kernel estimates
\[ \limsup_{x \to \zeta_0} \left| \frac{\mathcal M[h] (x) }{u^\star (x)} - h (\zeta_0) \right| \le \int_{|\zeta - \zeta_0| \le \varepsilon} \frac{ \mathbb M(x,\zeta) }{\int_{\partial \Omega} \mathbb M (x, z) d z} | h(\zeta) - h(\zeta_0)| d \zeta \]
As \(\varepsilon \to 0\) the RHS vanishes due to the continuity of \(h\).
Satisfying the boundary equation.
\(L^1\) data
A more involved argument works for \(h \in L^1 (\partial \Omega)\), and gives integral convergence.
For \(\psi \in \mathcal C (\overline \Omega)\) we have
\[ \frac{1}{\eta} \int_{\delta < \eta} \frac{\mathrm M[h]}{u^\star} \phi \longrightarrow \int_{\partial \Omega} h \phi, \qquad \text{ as } \eta \searrow 0. \]
Very weak formulation of singular boundary value problem
Let \(u_m = {\color{red} \mathrm G}[f_m]\). Then, for \(\psi \in L^\infty_c (\Omega)\)
\[ \int_\Omega u_m \psi = \frac{|\partial \Omega|}{ |A_m| } \int_{A_m} \widetilde h \frac{{\color{red} \mathrm G}[\psi]}{\delta^\gamma}. \]
Notice that as \(x \to \zeta \in \partial \Omega\)
\[ \lim_{x \to \zeta} \frac{{\color{red} \mathrm G}[\psi] (x)}{\delta(x)^\gamma} = \lim_{x \to \zeta} \int_\Omega \frac{{\color{red} \mathbb G}(y,x)}{\delta(x)^\gamma} \psi(y) dy \]
by the previous hypothesis
\[ D_\gamma {\color{red} \mathrm G}[\psi] (\zeta) = \int_\Omega \mathbb M (y, \zeta) \psi(y) dy \]
Given that
\[ \int_\Omega u_m \psi = \frac{|\partial \Omega|}{ |A_m| } \int_{A_m} \widetilde h \frac{{\color{red} \mathrm G}[\psi]}{\delta^\gamma}. \]
by compactness we know \(u_m \to u\) in \(L^1 (\mathrm{supp} \, \psi)\).
Large solutions
Large solutions
By large solution we mean solutions that satisfy \[ {\color{red} \mathcal L} u = F(x,u) \text{ in } \Omega \]
By large solution we mean \(\| u \|_{L^\infty} = \infty\).
We are mostly here in boundary blow-up: \[ u(x) \to \infty \text{ as } \mathrm{dist}(x, \partial \Omega) \to 0. \]
From now on
\[ \delta(x) = \mathrm{dist}(x, \partial \Omega) . \]
For the usual Laplacian
For the problems \[ -\Delta u + F(u) = 0 \]
Keller-Osserman condition:
Canonical example: \(-\Delta u + u^p = 0\) with certain \(p\).
This kind of singular behaviour exists also in the fractional case: [Chen & Véron, 2014], [Abatangelo, 2015], [Nguyen & Véron, 2018], [Chen & Véron, 2022]
Laplace equation
Weak dual formulation
Multiplying by \(\varphi\) and integrating \[ \int_\Omega ( {\color{red} \mathcal L} u ) \varphi = \int_\Omega f \varphi \]
If \(\varphi\) is in the suitable class of homogeneous boundary conditions \[ \int_\Omega u ( {\color{red} \mathcal L} \varphi )= \int_\Omega f \varphi \]
Choosing the set of test functions \(\varphi\) depends on the problem.
In particular, if we take \(\varphi = \mathrm{G} (\psi)\) we a weak dual solution if \(u \in L^1_{loc} (\Omega)\) and
\[\begin{equation} \tag{WDF} \int_\Omega u \psi = \int_\Omega f \mathrm{G} (\psi) \qquad \forall \psi \in L^\infty_c(\Omega) \end{equation}\]
See [Bonforte, Sire & Vázquez, 2015], [Bonforte, Figalli & Vázquez, 2018]
From the interior to the boundary.
Usual Laplacian
Laplace equation
\[ \begin{dcases} -\Delta u = f & \Omega \\ u = 0 & \partial \Omega \\ \end{dcases} \]
Poisson
\[ \begin{dcases} -\Delta v = 0 & \Omega \\ v = h & \partial \Omega \\ \end{dcases} \]
For \(\varphi \in W_0^{1,\infty} (\Omega)\) we have that
\[ \int_\Omega u (-\Delta\varphi) = \int_\Omega f \varphi \]
\[ \int_\Omega v (-\Delta\varphi) = - \int_{\partial \Omega} h \frac{\partial \varphi}{\partial n} \]
We aim to prove \(u_m = {\color{red} \mathrm G}[f_m] \to v\). Can we make \[ \int_\Omega f_m \varphi \to - \int_{\partial \Omega} h \frac{\partial \varphi}{\partial n} ? \]
Tubular neighbourhood of \(\partial \Omega\)
The map \[ \begin{aligned} (-R,R) \times \partial \Omega &\to U_R \\ (r,z) &\mapsto z - r \vec n(z) \end{aligned} \] is smooth and invertible for \(R\) small enough. This defines a tubular neighbourhood of \(\partial \Omega\).
A function in \(h \in L^1(\partial \Omega)\) can be extend to \(L^1(U_R)\) by \[\widetilde h(x) = h(z).\]
For integration, there is a Jacobian such that \(J(z,0) = 1\) and \[ \int_{a < \delta < b} f(x) = \int_{a}^b \int_{\partial \Omega} f(z - r n(z)) \, J(z,r) d z \, d r \]
From the kernel side
\[ \begin{aligned} {\color{red} \mathrm G}[f_m] (x) &= \frac{|\partial \Omega|}{|A_m|}\int_{A_m} \frac{{\color{red} \mathbb G}(x,y)}{\delta(y)} \widetilde h(y) dy \\ & \longrightarrow -\int_{\partial \Omega} \frac{\partial {\color{red} \mathbb G}}{\partial n_y} (x,\zeta) h(\zeta) d\zeta \\ &= v(x) \end{aligned} \]
This is the Poisson kernel.
Localising to the boundary. Usual Laplacian
Can we make \(\int_\Omega f_m \varphi \to -\int_{\partial \Omega} h \frac{\partial \varphi}{\partial n} ?\)
Extend \(h\) towards the interior by the tubular neighbourhood mapping and \[ f_m (x) = \widetilde h (x) \frac{|\partial \Omega|}{ |A_m| } \frac{\chi_{A_m} (x)}{\delta (x)} \qquad \text{ where } A_m = \left\{ x: \delta(x) < \frac 1 m \right\} \]
We get
\(\displaystyle\int_\Omega f_m \varphi = \frac{|\partial \Omega|}{|\{ \delta < \frac 1 m \}|} \int_{ \{ \delta < \frac 1 m \} } \widetilde h \frac{\varphi}{\delta}\)
\(\displaystyle \phantom{\int_\Omega f_m \varphi} = \frac 1 m \int_0^{\frac 1 m} \int_{\partial \Omega} h(z) \frac { \varphi( z - r n(z) ) }{ r } J(z,r) \, dz \, dr\)
\(\displaystyle \phantom{\int_\Omega f_m \varphi} \longrightarrow -\int_{\partial \Omega} h \frac{\partial \varphi}{\partial n} \qquad \text{ as } m \to \infty.\)
It works!
Localisation to the boundary
The kernel point of view
As we did before \[ \begin{align} {\color{red} \mathrm G}[f_m] &= \frac{|\partial \Omega|}{|A_j|} \int_{A_j} \frac{{\color{red} \mathbb G}(x,y)}{\delta(y)^\gamma} \widetilde h (y) d y \\ &{{} \to \int_{\partial \Omega} \mathbb M(x,\zeta) h(\zeta) d\zeta = \mathrm M[h]} \end{align} \]
where we have to assume that the following limit happens suitably (it does in the examples)
\[\begin{equation} \tag{H} \mathbb M(x,\zeta) = \lim_{y \to \zeta} \frac{{\color{red} \mathbb G}(x,y)}{\delta(y)^\gamma} = D_\gamma {\color{red} \mathbb G}(x,\zeta). \end{equation}\]
using the Green kernel estimates, we deduce that
\[ \mathbb M(x,\zeta) \asymp \frac{\delta(x)^\gamma}{|x-\zeta|^{d+\gamma - (2s-\gamma)}}. \]
We get define que representative \(u^\star (x) = \mathrm M[1].\)
[Abatangelo, GC & Vázquez, 2022]
Let \(h \in L^1 (\partial \Omega)\), then \(G(f_m) \to u\) in \(L^1_{loc} (\Omega)\).
It satisfies \[\begin{equation*} \tag{WDF$_s$} \label{eq:WDFs} \int_\Omega u \psi = \int_{\partial \Omega} h D_\gamma {{\color{red} \mathrm G}[\psi]}, \qquad \forall \psi \in L^\infty_c(\Omega). \end{equation*}\] Let \(u \in L^1_{loc} (\Omega)\) be such that \(\eqref{eq:WDFs}\). Then, \(u = \mathrm M [h]\).
If \(h \in C(\partial \Omega)\) then \(\frac{u(x)}{u^*(x)} \to h(z)\) as \(x \in \Omega \to z \in \partial \Omega\).
Martin problem. Literature
[Abatangelo, 2015] for the RFL
[Abatangelo & Dupaigne, 2017] for the SFL
[Chen, 2018] for the CFL
[Abatangelo, GC & Vázquez, 2022] unified theory for the RFL, CFL, and SFL.