Equazione di Burgers

Da testwiki.
Versione del 13 mar 2024 alle 12:13 di imported>Lulo
(diff) ← Versione meno recente | Versione attuale (diff) | Versione più recente → (diff)
Vai alla navigazione Vai alla ricerca

In matematica, l'equazione di Burgers, il cui nome si deve a Johannes Martinus Burgers, è un'equazione differenziale alle derivate parziali fondamentale per la meccanica dei fluidi, e utile anche in numerose aree della matematica applicata, quali la modellazione della gasdinamica e del flusso del traffico.

Per una data funzione y(x1,x2) di due variabili, la forma generale dell'equazione di Burgers è:

yx2+yyx1a2yx12=0

Quando a=0, l'equazione diventa inviscida:

yx2+yyx1=0

che è un prototipo per equazioni per le quali la soluzione può sviluppare discontinuità a funzione gradino (onde d'urto). La precedente equazione è la "forma avvettiva" dell'equazione di Burgers, mentre la "forma conservativa" è:

yx2+12x1(y2)=0

Soluzione

L'equazione di Burgers inviscida è un'equazione differenziale alle derivate parziali del primo ordine. La sua soluzione può essere costruita con il metodo delle caratteristiche. Seguendo questo metodo, se X(t) è una soluzione dell'equazione differenziale ordinaria:

dx1dx2=f[x1,x2]

la funzione f[x1,x2] è costante come funzione di x2. Allora [x1(x2),y(x2)] è una soluzione del sistema di equazioni ordinarie:

dx1dx2=y
dydx2=0

La soluzione di questo sistema, in termini dei valori iniziali, è:

x1(x2)=x1(0)+x2y(0)
y(t)=y(0)

Sostituendo x1(0)=a, e y(0)=y[x1(0),0]=y(a,0), il sistema diventa:

x1(x2)=a+x2y(a,0)
y(x2)=y(0)

In conclusione:

y(a,0)=y(0)=y(t)=y[x1(x2),x2]=y[a+x2y(a,0),x2]

Questa è una relazione implicita che determina la soluzione dell'equazione di Burgers inviscida, solo se le caratteristiche non si intersecano. Se le caratteristiche si intersecano, non esiste una soluzione classica all'equazione.

L'equazione di Burgers nel caso viscoso può essere linearizzata con la sostituzione di Cole-Hopf:

y=2ν1ϕϕx1

che la trasforma nell'equazione del calore:

ϕx2=ν2ϕx12

Questo permette di risolverla come un problema ai valori iniziali:

y(x1,x2)=2νx1ln{(4πνx2)1/2exp[(x1x1)24νx212ν0x1y(x1,0)dx1]dx1}

Bibliografia

Voci correlate

Collegamenti esterni

Template:Portale