Nonlinear parabolic problem with time fractional derivative

Titel: Niet-lineair parabolisch probleem met tijdsfractionele afgeleide

Auteurs: Nikolai Kutev en Tsviatko Rangelov
Instituut: Instituut voor Wiskunde en Informatica, Bulgaarse Academie van Wetenschappen

---

1. Probleemstelling

Het artikel onderzoekt een tijdsfractioneel parabolisch probleem dat de $p$-Laplace-operator combineert met een dubbel singuliere Hardy-type potentiaal. Het centrale probleem wordt beschreven door de volgende partiële differentiaalvergelijking (PDV) in het domein $Q = (0, T) \times \Omega$:

$$\begin{cases} D_t^\alpha u - \text{div}(|\nabla u|^{p-2}\nabla u) - \mu W(x)|u|^{p-2}u = 0, & \text{in } Q, \\ u(t, x) = 0, & \text{voor } (t, x) \in [0, T] \times \partial\Omega, \\ u(0, x) = u_0(x), & \text{voor } x \in \Omega. \end{cases}$$

Belangrijke parameters en definities:

• $D_t^\alpha$: De Liouville-Caputo fractionele afgeleide van orde $\alpha \in (0, 1)$ ten opzichte van de tijd.
• $\Delta_p u = \text{div}(|\nabla u|^{p-2}\nabla u)$: De $p$-Laplace-operator, met $p > 2$ en $n > p$ (waarbij $n$ de ruimtedimensie is).
• $\mu$: Een positieve parameter die de sterkte van de potentiaal bepaalt.
• $W(x)$: Een potentiaal die singulier is op de rand $\partial\Omega$ en in de oorsprong. Deze wordt gegeven door:
  $$W(x) = |x|^{-p} \left[ 1 - \left( \frac{|x|}{\phi(x)} \right)^{\frac{n-p}{p}} \right]^{-p}$$
  waarbij $\phi(x)$ een positieve homogene functie van orde 0 is die het domein $\Omega$ definieert (een ster-vormig domein).
• Optimale Hardy-constante: $\Lambda(n, p) = \left(\frac{n-p}{p}\right)^p$.

Het doel is om het gedrag van zwakke oplossingen te analyseren, specifiek gericht op het bestaan van globale oplossingen versus het optreden van "finite-time blow-up" (exploderen in eindige tijd), afhankelijk van de verhouding tussen $\mu$ en de Hardy-constante $\Lambda(n, p)$.

---

2. Methodologie

De auteurs hanteren een gestructureerde aanpak die bestaat uit de volgende stappen:

1. Truncatie van het probleem:
   Omdat de potentiaal $W(x)$ singulier is, wordt eerst een afgesneden (truncated) probleem geïntroduceerd met $W_N(x) = \min\{N, W(x)\}$. Dit zorgt voor een goed gedefinieerd probleem waarvoor klassieke technieken kunnen worden toegepast.

2. Vergelijkingsprincipes (Comparison Principles):

   • Voor fractionele ODE's: Er wordt eerst een vergelijkingsprincipe bewezen voor gewone fractionele differentiaalvergelijkingen (FODE) van de vorm $D_t^\alpha w = f(t, w)$. Dit is cruciaal om de monotonie van oplossingen te garanderen.
   • Voor fractionele parabolische PDV's: Met behulp van een "fundamentele identiteit" (gebaseerd op werk van Zacher) wordt een vergelijkingsprincipe voor de tijdsfractionele $p$-Laplace-vergelijking bewezen. Dit stelt de auteurs in staat om sub- en superoplossingen te vergelijken.

3. A-priori schattingen:
   Voor het afgesneden probleem worden energiestrategieën en Hardy-ongelijkheden gebruikt om uniforme schattingen af te leiden voor de oplossingen $u_N$. Deze schattingen zijn essentieel om de convergentie van de rij van oplossingen $u_N$ naar een limietoplossing $u$ te garanderen.

4. Analyse van Blown-up:
   Voor het geval van explosie wordt een subsolutie geconstrueerd die afhangt van de eerste eigenfunctie van een gerelateerd eigenwaardeprobleem en een oplossing van een gewone fractionele differentiaalvergelijking die in eindige tijd explodeert.

---

3. Belangrijkste Resultaten

Het artikel presenteert twee hoofdstellingen die het gedrag van de oplossingen volledig karakteriseren op basis van de parameter $\mu$:

Stelling 1.1 (Globale Existentie):
Als $\mu < \Lambda(n, p)$ (d.w.z. de potentiaal is zwakker dan de kritieke Hardy-grens) en de beginwaarde $u_0 \in W_0^{1,p}(\Omega)$, dan bestaat er een globale zwakke oplossing.

• De oplossing blijft begrensd in de normen $L^2(\Omega)$ en $W_0^{1,p}(\Omega)$ voor alle tijdstippen.
• Er treedt geen "blow-up" op; de oplossing bestaat voor alle $t > 0$.

Stelling 1.2 (Finite-Time Blow-up):
Als $\mu > \Lambda(n, p)$ (d.w.z. de potentiaal is sterker dan de kritieke Hardy-grens) en de beginwaarde $u_0$ is positief en begrensd ($u_0 \in L^\infty \cap W_0^{1,p}$), dan bestaat er geen globale oplossing.

• De oplossing explodeert in een eindige tijd $T_m < \infty$.
• Dit wordt bewezen door een subsolutie te construeren die in eindige tijd divergeert en het vergelijkingsprincipe toe te passen.

Kritieke Drempel:
De overgang tussen globale existentie en explosie wordt bepaald door de optimale Hardy-constante $\Lambda(n, p)$. Dit bevestigt dat de Hardy-ongelijkheid de fundamentele limiet vormt voor de stabiliteit van dit type niet-lineaire systemen, zelfs in de context van tijdsfractionele dynamica.

---

4. Bijdragen en Significantie

• Uitbreiding van bestaande theorie: Het artikel breidt bekende resultaten voor klassieke parabolische problemen (zonder fractionele afgeleiden) uit naar het domein van tijdsfractionele evolutionaire PDV's met de $p$-Laplace-operator.
• Dubbele singulariteit: Het behandelt een complex scenario met singulariteiten zowel in de oorsprong als op de rand van het domein, wat een veralgemening is van eerdere werken die zich vaak beperkten tot één type singulariteit.
• Methodologische bijdrage: De combinatie van het vergelijkingsprincipe voor fractionele ODE's met de "fundamentele identiteit" voor $p$-Laplace-systemen biedt een robuust raamwerk voor het analyseren van niet-lineaire fractionele systemen met singuliere potentiaaltermen.
• Toepassingsgebied: De resultaten zijn relevant voor fysica en engineering waar memory-effecten (via fractionele afgeleiden) en niet-lineaire diffusie (via $p$-Laplacian) samenkomen, zoals in modellen voor anomale diffusie in poreuze media, visco-elasticiteit en kwantummechanica met singuliere potentiaaltermen.

Samenvattend biedt dit artikel een rigoureuze wiskundige analyse die de stabiliteitsgrenzen van een belangrijke klasse van niet-lineaire fractionele PDV's vastlegt, waarbij de Hardy-ongelijkheid fungeert als de doorslaggevende factor voor het bestaan van oplossingen.