Burgersova rovnice

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 {\displaystyle u} and koeficient vazkosti ν {\displaystyle \nu } je obecný tvar jednorozměrné Burgersovy rovnice (rovněž známé pod pojmem vazká Burgesova rovnice) tvaru:

u t + u u x = ν 2 u x 2 {\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}} .

Je-li ν = 0 {\displaystyle \nu =0} , Burgersova rovnice se stává nevazkou Burgersovou rovnicí:

u t + u u x = 0 , {\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=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:

u t + 1 2 x ( u 2 ) = 0. {\displaystyle {\frac {\partial u}{\partial t}}+{\frac {1}{2}}{\frac {\partial }{\partial x}}{\big (}u^{2}{\big )}=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 ) {\displaystyle X(t)} řešením obyčejné diferenciální rovnice

d X ( t ) d t = u [ X ( t ) , t ] , {\displaystyle {\frac {\mathrm {d} X(t)}{\mathrm {d} t}}=u[X(t),t],}

pak

U ( t ) := u [ X ( t ) , t ] {\displaystyle U(t):=u[X(t),t]}

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

d X d t = U , {\displaystyle {\frac {\mathrm {d} X}{\mathrm {d} t}}=U,}
d U d t = 0. {\displaystyle {\frac {\mathrm {d} U}{\mathrm {d} t}}=0.}

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

X ( t ) = X ( 0 ) + t U ( 0 ) , {\displaystyle \displaystyle X(t)=X(0)+tU(0),}
U ( t ) = U ( 0 ) . {\displaystyle \displaystyle U(t)=U(0).}

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

X ( t ) = η + t u ( η , 0 ) {\displaystyle \displaystyle X(t)=\eta +tu(\eta ,0)}
U ( t ) = U ( 0 ) . {\displaystyle \displaystyle U(t)=U(0).}

Celkově:

u ( η , 0 ) = U ( 0 ) = U ( t ) = u [ X ( t ) , t ] = u [ η + t u ( η , 0 ) , t ] . {\displaystyle \displaystyle u(\eta ,0)=U(0)=U(t)=u[X(t),t]=u[\eta +tu(\eta ,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 , {\displaystyle u=-2\nu {\frac {1}{\phi }}{\frac {\partial \phi }{\partial x}},}

z čehož dostáváme rovnici tvaru

x ( 1 ϕ ϕ t ) = ν x ( 1 ϕ 2 ϕ x 2 ) , {\displaystyle {\frac {\partial }{\partial x}}{\Bigl (}{\frac {1}{\phi }}{\frac {\partial \phi }{\partial t}}{\Bigr )}=\nu {\frac {\partial }{\partial x}}{\Bigl (}{\frac {1}{\phi }}{\frac {\partial ^{2}\phi }{\partial x^{2}}}{\Bigr )},}

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

ϕ t = ν 2 ϕ x 2 + f ( t ) ϕ , {\displaystyle {\frac {\partial \phi }{\partial t}}=\nu {\frac {\partial ^{2}\phi }{\partial x^{2}}}+f(t)\phi ,}

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

ϕ t = ν 2 ϕ x 2 . {\displaystyle {\frac {\partial \phi }{\partial t}}=\nu {\frac {\partial ^{2}\phi }{\partial x^{2}}}.}

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

u ( x , t ) = 2 ν x ln { ( 4 π ν t ) 1 / 2 exp [ ( x x ) 2 4 ν t 1 2 ν 0 x u ( x , 0 ) d x ] d x } . {\displaystyle u(x,t)=-2\nu {\frac {\partial }{\partial x}}\ln {\Bigl \{}(4\pi \nu t)^{-1/2}\int _{-\infty }^{\infty }\exp {\Bigl [}-{\frac {(x-x')^{2}}{4\nu t}}-{\frac {1}{2\nu }}\int _{0}^{x'}u(x'',0)\mathrm {d} x''{\Bigr ]}\mathrm {d} x'{\Bigr \}}.}

Reference

V tomto článku byl použit překlad textu z článku Burgers' equation na anglické Wikipedii.

  1. Burgers Equation. [s.l.]: [s.n.] Dostupné online. Je zde použita šablona {{Citation}} označená jako k „pouze dočasnému použití“.
  2. HOPF, Eberhard. The partial differential equation ut + uux = μxx. Communications on Pure and Applied Mathematics. September 1950, s. 201–230. DOI 10.1002/cpa.3160030302. Je zde použita šablona {{Cite journal}} označená jako k „pouze dočasnému použití“.

Externí odkazy

  • (anglicky) Burgers' Equation at EqWorld: The World of Mathematical Equations.
  • (anglicky) Burgers' Equation Archivováno 26. 6. 2010 na Wayback Machine. at NEQwiki, the nonlinear equations encyclopedia.
  • (anglicky) Burgers shock-waves and sound in a 2D microfluidic droplets ensemble Phys. Rev. Lett. 103, 114502 (2009).