Regolarità Hölder di equazioni paraboliche con condizioni al contorno di Dirichlet e applicazione a sistemi di reazione-diffusione e reazione-cross-diffusione

Informazioni Fondamentali

• ID Articolo: 2511.18872
• Titolo: Hölder regularity of parabolic equations with Dirichlet boundary conditions and application to reaction-diffusion and reaction-cross-diffusion systems
• Autori: Hector Bouton, Laurent Desvillettes, Helge Dietert (Université Paris Cité and Sorbonne Université, CNRS, IMJ-PRG)
• Classificazione: math.AP (Analisi di Equazioni Differenziali Parziali)
• Data di Sottomissione: 24 novembre 2025
• Link Articolo: https://arxiv.org/abs/2511.18872v1

Riassunto

Il presente articolo generalizza i recenti risultati degli autori riguardanti condizioni al contorno di Neumann al caso delle condizioni di Dirichlet. L'oggetto principale di studio è l'equazione parabolica con coefficienti ruvidi $a$ della forma $a\partial_t w - \Delta w = f$, con condizioni al contorno di Dirichlet omogenee e l'ipotesi speciale $\partial_t w \geq 0$. Basandosi su questo risultato, si dimostra l'esistenza di soluzioni forti globali per il sistema di cross-diffusione triangolare di Shigesada-Kawasaki-Teramoto (SKT) tridimensionale con termini di reazione di Lotka-Volterra, nonché stime per le soluzioni di sistemi di reazione-diffusione chimici reversibili, tutti sotto condizioni al contorno di Dirichlet.

Contesto e Motivazione della Ricerca

Sfondo del Problema

1. Problema Centrale: Studio della regolarità Hölder di equazioni paraboliche con coefficienti ruvidi non costanti sotto condizioni al contorno di Dirichlet. In particolare, si considera l'equazione: $a(t,x)\partial_t w(t,x) - \Delta w(t,x) = f(t,x), \quad \partial_t w \geq 0$ dove $a$ soddisfa la condizione di ellitticità $0 < a_0 \leq a \leq c_0 a_0 < \infty$.
2. Importanza:
   • La regolarità Hölder è un problema fondamentale nella teoria delle equazioni differenziali parziali, cruciale per comprendere le proprietà delle soluzioni
   • Questo risultato trova applicazioni in matematica biologica, cinetica delle reazioni chimiche e altri campi delle scienze applicate
   • Fornisce strumenti tecnici chiave per l'esistenza globale di sistemi di cross-diffusione e sistemi di reazione-diffusione
3. Limitazioni dei Metodi Esistenti:
   • Il risultato classico di Krylov-Safonov (1981), sebbene copra situazioni più generali, presenta costanti non esplicite e richiede ipotesi spaziali più forti
   • I lavori precedenti degli autori boutondesvillettes2025 trattano solo condizioni al contorno di Neumann
   • Sotto condizioni di Dirichlet, il comportamento della soluzione vicino al contorno è più complesso, richiedendo nuove tecniche
4. Motivazione della Ricerca:
   • Sfruttare l'ipotesi di monotonia $\partial_t w \geq 0$ per semplificare la dimostrazione e ottenere costanti esplicite
   • Richiedere solo il controllo di $f$ in spazi di Lebesgue critici
   • Trattare completamente le condizioni al contorno, fornendo supporto teorico esplicito per le applicazioni

Contributi Principali

1. Risultato Teorico Principale (Teorema 1): Si dimostra che sotto la condizione $\gamma := 2 - \frac{2}{p} - \frac{d}{q} > 0$, la soluzione non negativa $w$ dell'equazione (1) appartiene a $C^{0,\alpha}([0,T]\times\Omega)$ con la stima: $\|w\|_{C^{0,\alpha}([0,T]\times\Omega)} \leq C_*\left(\|f\|_{L^p((0,T];L^q(\Omega))} + \|w^{in}\|_{C^1(\Omega)}\right)$
2. Applicazione a Sistemi Chimici Reversibili (Proposizione 1): Si dimostra l'esistenza globale di soluzioni forti per il sistema di reazione chimica reversibile $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,\ldots,4$ in spazi di dimensione al più quattro sotto condizioni al contorno di Dirichlet
3. Applicazione al Sistema SKT di Cross-Diffusione (Proposizione 2): Si dimostra l'esistenza globale di soluzioni forti per il sistema SKT triangolare\partial_t u - \Delta[(d_1+\sigma v)u] = f_u(u,v)\\ \partial_t v - d_2\Delta v = f_v(u,v) \end{cases}$$ in spazi di dimensione al più quattro sotto condizioni al contorno di Dirichlet
4. Innovazione Tecnica: Sviluppo di tecniche specializzate di decadimento dell'oscillazione adatte alle condizioni al contorno di Dirichlet, che differiscono essenzialmente dal caso di Neumann

Spiegazione Dettagliata dei Metodi

Definizione del Compito

Input:

• Dominio limitato $C^2$ $\Omega \subset \mathbb{R}^d$
• Coefficiente $a(t,x)$ soddisfacente la condizione di ellitticità
• Termine forzante $f \in L^p((0,T];L^q(\Omega))$, dove $\gamma = 2-\frac{2}{p}-\frac{d}{q} > 0$
• Dato iniziale $w^{in} \in C^1(\Omega)$ soddisfacente $w^{in}|_{\partial\Omega} = 0$

Output:

• Continuità Hölder della soluzione: $w \in C^{0,\alpha}([0,T]\times\Omega)$
• Stima esplicita della norma Hölder

Vincoli:

• Condizione al contorno di Dirichlet omogenea: $w = 0$ su $(0,T]\times\partial\Omega$
• Monotonia: $\partial_t w \geq 0$
• Non negatività: $w \geq 0$

Quadro Tecnico Principale

1. Stime del Nucleo di Calore di Dirichlet (Lemmi 1-2)

Lemma 2 (Stima di Limite Inferiore Cruciale): Per il nucleo di calore di Dirichlet $\Gamma_{B(0,R),D}$ nella palla $B(0,R)$, prendendo $T_R := a_0\frac{9R^2}{32d}$, si ha $\inf_{x,y\in B(0,R/4)} \inf_{t\in[\frac{T_R}{2c_0a_0}, \frac{T_R}{a_0}]} \Gamma_{B(0,R),D}(t,x,y) \gtrsim_{c_0,d} R^{-d}$

Idea della Dimostrazione:

• Costruzione di una funzione di confronto $\Psi(t,x,y) := p(t,|x-y|) - \sup_{0\leq s\leq t}p(s,3R/4)$
• Applicazione del principio del massimo per ottenere $\Gamma_{B(0,R),D}(t,x,y) > \Psi(t,x,y)$
• Per $x,y\in B(0,R/4)$ e intervalli di tempo appropriati, calcolo esatto che fornisce un limite inferiore positivo

2. Decadimento dell'Oscillazione (Proposizione 3)

Idea Centrale: Dimostrazione del decadimento geometrico dell'oscillazione in regioni lontane dal contorno.

Proposizione 3: Per $R>0$ soddisfacente $\beta R^2 \leq T$ (dove $\beta = \frac{9a_0}{32d}$), la soluzione $w:[−\beta R^2, 0]\times B(0,R)\to[0,1]$ soddisfa $\text{osc}_{(-\beta R^2/16,0]\times B(0,R/4)} w \leq 1-\delta + C_f R^\gamma \|f\|_{L^p((−\beta R^2,0);L^q(B(0,R)))}$

Strategia della Dimostrazione:

• Divisione in due casi: $|\{w(-T_R,\cdot)\geq\frac{1}{2}\}\cap B(0,R/4)| \geq \frac{1}{2}|B(0,R/4)|$ oppure il suo complementare
• Costruzione di una funzione di confronto $v$ soddisfacente $(a_0c_0\partial_t - \Delta)v = f$
• Applicazione del principio di confronto e del Lemma 2 per il limite inferiore, ottenendo $w(t,x) > \delta - C_2R^\gamma\|f\|$
• Iterazione per ottenere il tasso di decadimento geometrico

3. Stime Vicino al Contorno (Lemma 3)

Lemma Chiave 3: Per $z\in\Omega$, sia $d_z := d(z,\partial\Omega)$ la distanza dal contorno, allora $w(t,x) \leq C(\|w^{in}\|_{\text{Lip}(\Omega)} + \|f\|_{L^p([0,T],L^q(\Omega))}) d_x^{\tilde{\gamma}}$ dove $\tilde{\gamma} := \min(\gamma, 1-2\varepsilon)$

Punti Chiave della Dimostrazione:

• Caso 1 (tempo piccolo $t^{1/2}\leq d_x$): Utilizzo della rappresentazione integrale del nucleo di calore e del Lemma 1
• Caso 2 (tempo grande $t > d_x^2$): Utilizzo della stima Gaussiana $\Gamma_{\Omega,D}(t,x,y) \lesssim \frac{d_x d_y}{t^{\frac{d+2}{2}}}e^{-c\frac{|x-y|^2}{t}}$ Decomposizione dell'integrale in tre parti con stima dettagliata

4. Stima per Tempi Piccoli (Lemma 4)

Lemma 4: $|w(t,x) - w(0,x)| \leq C(\|w^{in}\|_{C^1(\Omega)} + \|f\|_{L^p([0,T],L^q(\Omega))}) t^{\frac{\min(1,\gamma)}{2}}$

Metodo della Dimostrazione:

• Utilizzo del principio di confronto: $v_{a_0c_0} \leq w \leq v_{a_0}$
• Decomposizione di $v_{a_0}$ in $v_{a_0,f} + \tilde{v}_{a_0}$
• Applicazione delle stime di Schauder a $\tilde{v}_{a_0}$ per ottenere la continuità Hölder parabolica
• Applicazione del Lemma 1 a $v_{a_0,f}$

Architettura Completa della Dimostrazione del Teorema 1

Trattamento Unificato di Tre Casi:

1. Vicino al Contorno ($\min(d_x,d_y) \leq 2|x-y|$):
   • Applicazione diretta del Lemma 3
   • Ottenimento di $|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
2. Tempi Piccoli ($t'^{1/2}, t^{1/2} \leq |x-y|$):
   • Applicazione del Lemma 4 e della proprietà di Lipschitz del dato iniziale
   • Ottenimento di $|w(t,x)-w(t',y)| \lesssim |x-y|^{\tilde{\gamma}}$
3. Lontano dal Contorno e Tempi Grandi ($|x-y| \leq \min(t'^{1/2}, t^{1/2}, \frac{\min(d_x,d_y)}{2})$):
   • Scelta della scala $R := \min(t^{1/2}, \min(d_x,d_y))$
   • Se $\min(d_x,d_y) \leq t^{1/2}$: utilizzo del Lemma 3
   • Se $\min(d_x,d_y) > t^{1/2}$: utilizzo del Lemma 4
   • Applicazione della Proposizione 3 sul decadimento dell'oscillazione, mediante iterazione si ottiene: $o_k \lesssim R^{\tilde{\gamma}}(\|f\| + \|w^{in}\|)\Lambda^k$ dove $\Lambda > \max(\frac{1}{4^\gamma}, 1-\delta)$
   • Scelta di $k_0$ tale che $R4^{-k_0} \leq \max(|x-y|, (\frac{|t-t'|}{\beta})^{1/2})$
   • Ottenimento di $|w(t,x)-w(t',y)| \lesssim \max(|x-y|, |t-t'|^{1/2})^{\alpha_\Lambda}$

Spiegazione Dettagliata dei Sistemi di Applicazione

Applicazione 1: Sistema Chimico Reversibile (Proposizione 1)

Modello del Sistema: $\partial_t u_i - d_i\Delta u_i = (-1)^i(u_1u_3 - u_2u_4), \quad i=1,2,3,4$ che rappresenta la reazione reversibile $A_1 + A_3 \rightleftharpoons A_2 + A_4$

Strategia della Dimostrazione:

1. Stime a Priori Standard (Equazione 12): $\sum_{i=1}^4 \int_\Omega \frac{u_i^{p+1}}{p+1}(T) + \frac{4p}{(p+1)^2}\sum_{i=1}^4 d_i\int_0^T\int_\Omega |\nabla(u_i^{\frac{p+1}{2}})|^2 \leq \text{dato iniziale} + C_p\sum_{i=1}^4\int_0^T\int_\Omega u_i^{p+2}$
2. Definizione di Funzione Ausiliaria: $w := \int_0^t(\sum_{i=1}^4 d_i u_i)$
   • Verifica: $w\geq 0$, $\partial_t w \geq 0$, $w|_{\partial\Omega}=0$
   • Calcolo: $\Delta w = \sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}$
   • Definizione di $a := \frac{\sum_{i=1}^4 u_i}{\sum_{i=1}^4 d_i u_i}$ soddisfacente $\frac{1}{\max d_i} \leq a \leq \frac{1}{\min d_i}$
   • Ottenimento di: $a\partial_t w - \Delta w = \sum_{i=1}^4 u_i^{in}$
3. Applicazione del Teorema 1: Ottenimento di $\sum_{i=1}^4 u_i \in C^{0,\alpha}([0,T]\times\Omega)$
4. Disuguaglianza di Interpolazione (Equazione 13): $\left\|\sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}\right\|_{L^{\frac{2(3-\alpha)}{2-\alpha}}(\Omega)}^3 \leq C\|w\|_{C^{0,\alpha}(\Omega)}^{\frac{3}{3-\alpha}} \left\|\nabla\left[\sum_{i=1}^4 u_i - \sum_{i=1}^4 u_i^{in}\right]\right\|_{L^2(\Omega)}^{\frac{3(2-\alpha)}{3-\alpha}}$
5. Argomento di Bootstrap: Combinazione delle stime (12) e (13), ottenimento di $u_i \in L^{3+\delta}([0,T]\times\Omega)$, da cui seguono tutte le stime $L^p$

Applicazione 2: Sistema SKT di Cross-Diffusione (Proposizione 2)

Modello del Sistema: