Formulazione debole
# Formulazione debole

Nell'ambito delle equazione differenziali, in particolare delle equazioni alle derivate parziali, è di grande importanza lo studio della formulazione debole dei problemi differenziali classici, che per dualità vengono anche chiamati problemi in forma forte o classica. Risolvere un problema in forma debole significa trovare una soluzione, detta soluzione debole, le cui derivate possono non esistere, ma che è comunque soluzione dell'equazione in qualche modo ben preciso. Molto spesso si tratta delle uniche soluzioni che è possibile trovare.

Il concetto di soluzione debole è legato a quello di derivata debole: si tratta di definire la nozione di derivata anche per funzioni integrabili ma non necessariamente differenziabili.

## Introduzione

Un problema relativo ad un'equazione differenziale si dice ben posto se possiede una soluzione, se tale soluzione è unica e se dipende in modo continuo dai dati forniti dal problema. Un problema ben posto contiene tutte le caratteristiche ideali al fine di studiarne la risolubilità. La soluzione di un'equazione alle derivate parziali di ordine $k$ si definisce informalmente soluzione classica o soluzione forte se è una funzione differenziabile fino all'ordine $k$-esimo e tutte le derivate esistono e sono continue: per risolvere una PDE in senso classico bisogna dunque cercare una funzione liscia o almeno di classe $C^{k}$. La maggior parte delle equazioni differenziali alle derivate parziali non ammette soluzioni classiche, come ad esempio le equazioni di continuità. Se si ammette una funzione non differenziabile come soluzione di un problema ben posto tale soluzione è una soluzione debole, anche detta "soluzione generalizzata" o "soluzione integrale". La formulazione debole di un problema deriva da quella forte, e una soluzione del problema forte è anche soluzione del problema debole.

## Descrizione generale

L'idea di fondo delle formulazioni deboli è quella che portò anche all'introduzione in matematica delle distribuzioni, o "funzioni generalizzate": si tratta di funzionali lineari definiti sullo spazio di funzioni costituito dalle funzioni dette funzioni di test. Lo spazio delle distribuzioni è lo spazio duale di quello delle funzioni di test. Si tratta di funzioni in un senso più generale: alcune distribuzioni, se viste come funzioni, possono anche non avere alcun corrispettivo nell'analisi tradizionale (si veda ad esempio la delta di Dirac). Intuitivamente, se questo spazio "di test" è abbastanza grande e se ha certe proprietà, è ragionevole pensare di ricostruire la funzione (generalizzata) sapendo come essa agisce su ogni funzione test dello spazio.

Presa un'equazione, per trovare una soluzione debole si procede generalmente col moltiplicare ambo i termini per una funzione test $\varphi$, e di integrare poi entrambi i membri su tutto il dominio di interesse. Dopodiché si "scaricano" le derivate (integrando per parti) dalla funzione $u$ sulla funzione test $\varphi$ quanto basta per poter richiedere la minor regolarità possibile sia a $u$ che a $\varphi$. Per poter effettuare le integrazioni è necessario che sia $u$ che $\varphi$ stiano almeno in $L^{2}$ (altrimenti l'integrale non ha senso); inoltre per poter integrare anche i prodotti tra le derivate occorre che stiano anche nello spazio di Sobolev $H^{k}$, dove $k$ indica il massimo ordine di derivazione che compare dopo aver scaricato le derivate di $u$ su $\varphi$. Si consideri dunque un operatore differenziale lineare in un insieme aperto $W$ in $\mathbb {R} ^{n}$:

    $P(x,\partial )u(x)=\sum a_{\alpha _{1},\alpha _{2},\dots ,\alpha _{n}}(x)\partial ^{\alpha _{1}}\partial ^{\alpha _{2}}\cdots \partial ^{\alpha _{n}}u(x)$

in cui il multi-indice $(\alpha _{1},\alpha _{2},\dots ,\alpha _{n})$ spazia in un sottoinsieme finito di $\mathbb {N} ^{n}$ ed i coefficienti $a_{\alpha _{1},\alpha _{2},\dots ,\alpha _{n}}$ sono funzioni sufficientemente lisce di $x$. L'equazione $P(x,\partial )u(x)=0$, dopo essere stata moltiplicata per una funzione di test $\varphi$ liscia e avente un supporto compatto in $W$, può essere integrata per parti $\alpha _{1}+\alpha _{2}+\dots +\alpha _{n}$ volte in modo che viene ad essere scritta come:

    $\int _{W}u(x)Q(x,\partial )\varphi (x)\,\mathrm {d} x=0,$

dove l'operatore differenziale $Q(x,\partial )$ è dato da:

    $Q(x,\partial )\varphi (x)=\sum (-1)^{|\alpha |}\partial ^{\alpha _{1}}\partial ^{\alpha _{2}}\cdots \partial ^{\alpha _{n}}\left[a_{\alpha _{1},\alpha _{2},\dots ,\alpha _{n}}(x)\varphi (x)\right].$

Il numero

    $(-1)^{|\alpha |}=(-1)^{\alpha _{1}+\alpha _{2}+\cdots +\alpha _{n}}$

appare poiché ogni integrazione per parti produce una moltiplicazione per -1. L'operatore $Q(x,\partial )$ è l'operatore aggiunto di $P(x,\partial )$.

Si vede quindi che se l'originale formulazione (formulazione forte) richiede di trovare una funzione $u$ (soluzione forte) definita su $W$, differenziabile $|\alpha |$-volte e tale che:

    $P(x,\partial )u(x)=0,\qquad \forall x\in W,$

allora una funzione integrabile $u$ è una soluzione debole se:

    $\int _{W}u(x)Q(x,\partial )\varphi (x)\,\mathrm {d} x=0,$

per ogni funzione liscia a supporto compatto $\varphi$.

Su domini limitati una soluzione forte è anche soluzione debole, in quanto le procedure di integrazione per parti sono lecite. Se ci si pone il problema inverso, cioè se una soluzione $u$ del problema debole soddisfa anche il problema forte, si vede che $u$ non può essere soluzione forte se si interpretano le derivate in senso classico per due motivi:

* La funzione $u$ appartiene a $H^{1}$ e dunque non può avere in generale una derivata seconda continua (altrimenti sarebbe anche in $H^{2}$), come invece richiesto dalla soluzione forte.
* Nella formulazione debole non è nemmeno chiesto che $u$ sia definita ovunque. Affinché ogni integrale di Lebesgue abbia senso, $u$ può assumere valori arbitrari anche in un'infinità numerabile di punti del dominio (più precisamente in un insieme con misura di Lebesgue nulla, o quasi ovunque).

Si spiega quindi il motivo di considerare $u$ non più come una funzione, ma come una distribuzione. Assumendo ciò e interpretando le derivate nel senso delle distribuzioni, si può dire che $u$ soddisfa il problema forte (nel senso delle distribuzioni). Anche l'assunzione dei dati al bordo è problematica: per quanto detto sopra, considerando che il bordo del dominio ha sempre misura nulla, parlare del valore di $u$ sul bordo non ha senso classicamente. La soluzione a questo problema si ha considerando il dato al bordo come limite (nel senso di $L^{2}$) di funzioni di classe $C^{\infty }$ a supporto compatto che approssimano $u$ nel senso di $H^{1}$.

### Esempio

Per illustrare il concetto, si consideri l'equazione del trasporto:

    ${\frac {\partial u}{\partial t}}+{\frac {\partial u}{\partial x}}=0$

in cui $u(t,x)$ è differenziabile con continuità su $R^{2}$. Moltiplicando l'equazione per una funzione liscia e a supporto compatto $\varphi$, e integrando si ottiene:

    $\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }{\frac {\partial u(t,x)}{\partial t}}\varphi (t,x)\,\mathrm {d} t\mathrm {d} x+\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }{\frac {\partial u(t,x)}{\partial x}}\varphi (t,x)\,\mathrm {d} t\mathrm {d} x=0.$

Grazie al teorema di Fubini è possibile scambiare l'ordine di integrazione, in modo che integrando per parti in $t$ il primo termine e in $x$ il secondo:

    $-\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }u(t,x){\frac {\partial \varphi (t,x)}{\partial t}}\,\mathrm {d} t\mathrm {d} x-\int _{-\infty }^{\infty }\int _{-\infty }^{\infty }u(t,x){\frac {\partial \varphi (t,x)}{\partial x}}\,\mathrm {d} t\mathrm {d} x=0.$

Si nota che gli integrali vanno da −∞ a ∞, ma sono sostanzialmente valutati su un dominio chiuso in quanto $\varphi$ ha supporto compatto. Esiste quindi una funzione $u$, che può non essere differenziabile, che soddisfa quest'ultima equazione per ogni $\varphi$ ma che non è una soluzione dell'equazione del trasporto: si tratta di una soluzione debole.

Ad esempio:

    $u(t,x)=|t-x|,\qquad \forall t,x,$

è una soluzione debole, come si mostra integrando per parti ai lati della retta $x=t$.

## Lemma di Lax-Milgram

Sia $V$ uno spazio di Banach. Si vuole trovare una soluzione $u\in V$ dell'equazione:

    $Au=f,$

dove $A\colon V\to V'$ e $f\in V'$, con $V'$ lo spazio duale di $V$.

Il calcolo delle variazioni mostra come questo sia equivalente a trovare $u\in V$ tale che per tutti i $v\in V$ vale:

    $[Au](v)=f(v).$

Si può considerare $v$ un vettore o funzione di test.

La formulazione debole del problema significa trovare $u\in V$ tale che:

    $a(u,v)=f(v),\quad \forall v\in V,$

definendo la forma bilineare:

    $a(u,v):=[Au](v).$

### Enunciato

Il lemma di Lax-Milgram può essere applicato alle forme bilineari, anche se non ne è la versione più generale. Sia $V$ uno spazio di Hilbert e $a(\cdot ,\cdot )$ una forma bilineare su $V$ che è limitata

    $|a(u,v)|\leq C\|u\|\|v\|$

e coercitiva

    $|a(u,u)|\geq c\|u\|^{2}.$

Allora, per ogni $f\in V'$ esiste una soluzione unica $u\in V$ per l'equazione:

    $a(u,v)=f(v)$

e si ha:

    $\|u\|\leq {\frac {1}{c}}\|f\|_{V'}.$

### Sistema di equazioni lineari

Ad esempio, nel caso di un sistema di equazioni lineari si ha $V=\mathbb {R} ^{n}$, e $A:V\to V$ è una trasformazione lineare. La formulazione debole dell'equazione:

    $Au=f$

consiste nel trovare $u\in V$ tale che per ogni $v\in V$ vale l'equazione:

    $\langle Au,v\rangle =\langle f,v\rangle ,$

dove $\langle \cdot ,\cdot \rangle$ denota il prodotto interno.

Dato che $A$ è una mappa lineare è sufficiente provare i vettori di base $e_{i}$:

    $\langle Au,e_{i}\rangle =\langle f,e_{i}\rangle \quad i=1,\ldots ,n.$

Utilizzando l'espansione come combinazione lineare dei vettori di base

    $u=\sum _{j=1}^{n}u_{j}e_{j},$

si ottiene la forma matriciale dell'equazione:

    $\mathbf {A} \mathbf {u} =\mathbf {f} ,$

dove $a_{ij}=\langle Ae_{j},e_{i}\rangle$ e $f_{i}=\langle f,e_{i}\rangle$.

La forma bilineare associata a tale formulazione debole è:

    $a(u,v)=\mathbf {v} ^{T}\mathbf {A} \mathbf {u} .$

Si nota che tutte le forme bilineari su $\mathbb {R} ^{n}$ sono limitate e in particolare:

    $|a(u,v)|\leq \|A\|\,\|u\|\,\|v\|.$

Per quanto riguarda la coercitività, significa che la parte reale degli autovalori di $A$ non deve essere più piccola di $c$. Questo implica che nessun autovalore può essere nullo, e il sistema è quindi risolvibile. Inoltre, si può stimare:

    $\|u\|\leq {\frac {1}{c}}\|f\|,$

dove $c$ è la più piccola parte reale assunta dagli autovalori di $A$.

## Esempio monodimensionale

Si consideri il seguente problema di Poisson con condizioni al bordo miste omogenee:

    $-u''(x)=f,\qquad x\in (-1,1),$
    $u(-1)=0,\qquad u'(1)=0.$

Moltiplicando a destra e a sinistra per una funzione test $v$, per il momento senza specificare a quale spazio appartiene, e integrando per parti tra $-1$ e $1$ si ha:

    $\int \limits _{-1}^{1}u'v'-[u'v]_{-1}^{1}=\int \limits _{-1}^{1}fv.$

Sfruttando quindi le condizioni al bordo per $u$ si può scrivere:

    $\int \limits _{-1}^{1}u'v'+u'(-1)v(-1)=\int \limits _{-1}^{1}fv,$

dove sia $u$ che $v$ devono stare in $H^{1}(-1,1)$ affinché gli integrali abbiano senso. Spesso, soprattutto in analisi numerica, si preferisce effettuare il cambio di incognita ponendo:

    $u={\tilde {u}}+Rg,$

dove $Rg$ è detta "rilevamento" di $u$ sul bordo. La funzione $Rg$, infatti, assume al bordo gli stessi valori di $u$, in modo che ${\tilde {u}}$ sia nulla sul bordo. Inoltre $Rg$ deve appartenere anch'essa a $H^{1}$, di modo che sostituendo ${\tilde {u}}$ nell'equazione si ottenga:

    $\int \limits _{-1}^{1}{\tilde {u}}'v'=\int \limits _{-1}^{1}fv-\int \limits _{-1}^{1}Rg'v'-(u'-Rg')(-1)v(-1).$

Se ora si sceglie come spazio delle funzioni test lo spazio:

    $V=\left\{v{\text{ di }}H^{1}(-1,1):v(-1)=0\right\},$

allora ${\tilde {u}}$ e $v$ stanno nello stesso spazio. Questo è molto utile poiché risulta possibile applicare il lemma di Lax-Milgram per verificare se il problema è ben posto, cioè se ammette un'unica soluzione e se questa dipende con continuità dai dati.

## Formulazione per equazioni ellittiche del secondo ordine

Un'equazione differenziale lineare alle derivate parziali ellittica del secondo ordine in $n$ variabili indipendenti $\mathbf {z} =\left(z_{1},z_{2},\dots ,z_{n}\right)^{T}$ definita su insieme aperto $\Omega \subset \mathbb {R} ^{n}$ può essere scritta in modo generale come:

    $-\sum _{i,j=1}^{n}{\frac {\partial }{\partial x_{i}}}\cdot \left(a_{ij}{\frac {\partial u}{\partial x_{j}}}\right)+\sum _{i=1}^{n}\left({\frac {\partial }{\partial x_{i}}}\left(b_{i}u\right)+c_{i}{\frac {\partial u}{\partial x_{i}}}\right)+a_{0}u=f,$

dove le variabili sono tutte funzioni di $\mathbf {z}$.

È possibile scrivere tale equazione anche nella forma:

    $-\nabla \cdot \left(a_{1}\nabla u\right)+a_{0}u=f$

assumendo $a_{ij}\left(\mathbf {x} \right)=a_{1}\left(\mathbf {x} \right)\delta _{ij}$ e $\mathbf {b} \left(\mathbf {x} \right)=\mathbf {c} \left(\mathbf {x} \right)=0$ in $\Omega$.

La soluzione classica di tale problema consiste nella determinazione di una funzione $u\in C^{2}\left(\Omega \right)$ che soddisfi l'equazione nella sua forma generale per tutti i vettori $\mathbf {x} \in \Omega$ e che soddisfi inoltre le condizioni al bordo per tutti i vettori $\mathbf {x} \in \partial \Omega$. Tale problema non risulta risolvibile in generale, e per questo motivo si introduce la formulazione debole del problema.

La sua derivazione consiste in quattro passi:

* Moltiplicazione ad entrambi i membri per una funzione di test $v\in C^{+\infty }\left(\Omega \right)$:

       $-\nabla \cdot \left(a_{1}\nabla u\right)\cdot v+a_{0}u\cdot v=f\cdot v.$
* Integrazione su $\Omega$:

       $-\int _{\Omega }\nabla \cdot \left(a_{1}\nabla u\right)\cdot vd\mathbf {x} +\int _{\Omega }a_{0}u\cdot vd\mathbf {x} =\int _{\Omega }f\cdot vd\mathbf {x} .$
* Utilizzo del lemma di Green per la riduzione del grado massimo delle derivate:

       $\int _{\Omega }\left(a_{1}\nabla u\right)\cdot \nabla vd\mathbf {x} +\int _{\Omega }a_{0}u\cdot vd\mathbf {x} =\int _{\Omega }f\cdot vd\mathbf {x} +\oint _{\partial \Omega }v\cdot a_{1}\cdot {\frac {\partial u}{\partial n}}d\mathbf {x} ,$
    con $n$ normale alla frontiera di $\Omega$. È possibile altresì dividere il bordo a seconda delle condizioni che vengono fornite per esso. Assumendo $\partial \Omega =\Gamma _{D}\cup \Gamma _{N}$, dove $\Gamma _{D}$ indica i punti del bordo dove vengono fornite condizioni di Dirichlet e $\Gamma _{N}$ i punti del bordo dove vongono fornite condizioni di Neumann. L'equazione precedente si può quindi sviluppare come:
       ${\begin{aligned}\int _{\Omega }\left(a_{1}\nabla u\right)\cdot \nabla vd\mathbf {x} +\int _{\Omega }a_{0}u\cdot vd\mathbf {x} &=\int _{\Omega }f\cdot vd\mathbf {x} +\oint _{\partial \Omega }v\cdot a_{1}\cdot {\frac {\partial u}{\partial n}}d\mathbf {x} \\&=\int _{\Omega }f\cdot vd\mathbf {x} +\oint _{\Gamma _{N}}v\cdot a_{1}\cdot {\frac {\partial u}{\partial n}}d\mathbf {x} +\oint _{\Gamma _{D}}v\cdot a_{1}\cdot {\frac {\partial u}{\partial n}}d\mathbf {x} \\&=\int _{\Omega }f\cdot vd\mathbf {x} +\oint _{\Gamma _{N}}v\cdot a_{1}\cdot {\frac {\partial u}{\partial n}}d\mathbf {x} \end{aligned}}$
* Determinazione dei più ampi spazi funzionali tali per cui $u$ e $v$ siano funzioni con integrale finito:

       $\left\{{\begin{array}{l}u,v\in H_{0}^{1}\left(\Omega \right)\\f\in L^{2}\left(\Omega \right)\end{array}}\right.$
    con $H$ indicante lo spazio di Sobolev.

La formulazione debole richiede quindi a questo punto la determinazione della funzione $u\in H_{0}^{1}\left(\Omega \right)$ che verifica l'equazione all'ultimo punto. Chiaramente la formulazione classica determina una funzione che soddisfa anche la formulazione debole.

