Burgersova rovnice

Z testwiki
Skočit na navigaci Skočit na vyhledávání

Burgersova rovnice je jednou ze základních parciálních diferenciálních rovnic mechaniky tekutin. Objevuje se v mnoha partiích aplikované matematiky, jako je například dynamika plynů a modelování dopravního toku. Rovnice je pojmenována po J. M. Burgersovi (1895–1981). Je ekvivalentní Navierově–Stokesově rovnici pro nestlačitelný tok bez tlakového členu.[1]

Pro danou rychlost u and koeficient vazkosti ν je obecný tvar jednorozměrné Burgersovy rovnice (rovněž známé pod pojmem vazká Burgesova rovnice) tvaru:

ut+uux=ν2ux2.

Je-li ν=0, Burgersova rovnice se stává nevazkou Burgersovou rovnicí:

ut+uux=0,

což je jeden z typů rovnic, v jejichž řešení se mohou vyskytnout nespojitosti (rázové vlny). Předešlá rovnice je advekční formou Burgersovy rovnice. Konzervativní forma je tvaru:

ut+12x(u2)=0.

Řešení

Nevazká Burgersova rovnice

Nevazká Burgersova rovnice je parciální diferenciální rovnicí prvního řádu. Její řešení může být zkonstruováno pomocí metody charakteristik. Tato metoda říká, že pokud je X(t) řešením obyčejné diferenciální rovnice

dX(t)dt=u[X(t),t],

pak

U(t):=u[X(t),t]

je konstantní vzhledem k t. Tudíž [X(t),U(t)] je řešením soustavy obyčených diferenciálních rovnic:

dXdt=U,
dUdt=0.

Řešení této soustavy je vyjádřeno pomocí počáteční hodnoty výrazem:

X(t)=X(0)+tU(0),
U(t)=U(0).

Při substituci X(0)=η, kdy platí U(0)=u[X(0),0]=u(η,0), můžeme zapsat soustavu ve tvaru

X(t)=η+tu(η,0)
U(t)=U(0).

Celkově:

u(η,0)=U(0)=U(t)=u[X(t),t]=u[η+tu(η,0),t].

Toto je implicitní vztah určující řešení nevazké Burgesovy rovnice za předpokladu, že se jednotlivé charakteristiky vzájemně neprotínají. Pokud k průniku charakteristik dojde, pak neexistuje klasické řešení rovnice.

Vazká Burgersova rovnice

Vazká Burgersova rovnice může být linearizována Coleovou–Hopfovou transformací [2]

u=2ν1ϕϕx,

z čehož dostáváme rovnici tvaru

x(1ϕϕt)=νx(1ϕ2ϕx2),

která může být přepsána jako

ϕt=ν2ϕx2+f(t)ϕ,

kde f(t) je libovolná funkce. Pokud poslední člen vymizí, obdržíme difuzní rovnici

ϕt=ν2ϕx2.

Můžeme tedy řešit počáteční úlohu:

u(x,t)=2νxln{(4πνt)1/2exp[(xx)24νt12ν0xu(x,0)dx]dx}.

Reference

Šablona:Překlad

Externí odkazy

Šablona:Autoritní data