Umění najít řešení diferenciální rovnice je sympatické a naučíme se v úvodním příkladě. Není to však nic proti umění sestavit model (naučili jsme se již ve druhém týdnu a připomeneme si v následujícím příkladě s tloušťkou ledu), umění posoudit jednoznačnost řešení (většina modelů se řeší numericky a musíme být přesvědčeni o smysluplnosti takové činnosti) a stabilitu řešení (řešení, která nejsou stabilní, jsou sice v souladu s přírodními zákony, ale pravděpodobnost jejich spontánního výskytu je nulová). Jednoznačnost a zjednodušenou verzi stability řešení (stabilita konstantních řešení) jsme viděli na přednášce a připomeneme v dalších příkladech.
Řešení
Takzvaný Stefanův zákon (J. Stefan, "Uber die Theorie der Eisbildung, insbesondere "uber die Eisbildung im Polarmeere, 1891) vyjadřuje že tloušťka ledu na hladině moře roste ve stabilních podmínkách rychlostí nepřímo úměrnou této tloušťce. Zapište tento fakt pomocí vhodného matematického modelu a najděte řešení vzniklé diferenciální rovnice.
Řešení
\[ \begin{aligned} \frac{\mathrm dh}{\mathrm dt}&=\frac kh\\ h\,\mathrm dh&=k\, \mathrm dt\\ \int h\,\mathrm dh&=\int k\, \mathrm dt\\ \frac {h^2}{2}&=kt+C\\ \end{aligned} \]
Z fyziky je známo, že rychlost s jakou vytéká tekutina otvorem u dna nádoby je úměrná odmocnině výšky hladiny (protože se mění potenciální energie úměrná výšce na kinetickou energii úměrnou druhé mocnině rychlosti). Proto je i rychlost s jakou se zmenšuje objem vody v nádrži úměrná odmocnině výšky hladiny.
Ukažte, že matematickým popisem procesu je diferenciální rovnice. Napište rovnici pro výšku hladiny vody v nádrži jako funkci času. Uvažujte tři případy: nádrž cylindrického tvaru (válec postavený na podstavu), nádrž ve tvaru kvádru a nádrž ve tvaru kužele otočeného vrcholem dolů (trychtýř).
V tomto příkladě vystupuje derivace jak rychlost, ale po přepisu zadání do modelu máme v rovnici dvě různé veličiny, které se mění: objem vody a výšku hladiny. Musíme ještě najít a použít vztah mezi rychlostmi změn těchto veličin. Fyzikální zákon je formulován pro derivaci objemu a nás zajímá derivace výšky.
Řešení
Buď \(\displaystyle V\) objem vody a \(\displaystyle h\) výška hladiny od dna. Podle zadání ve všech případech platí \[\frac {\mathrm dV}{\mathrm dt}=-k_1\sqrt h\] a musíme derivaci \(\displaystyle \frac {\mathrm dV}{\mathrm dt}\) vyjádřit pomocí \(\displaystyle \frac {\mathrm dh}{\mathrm dt}\).
Pro cylindr, kvádr nebo jakoukoliv nádrž se svislými stěnami je objem úměrný výšce hladiny, \(\displaystyle V=k_2 h\), a proto \(\displaystyle \frac {\mathrm dV}{\mathrm dt}=k_2\frac {\mathrm dh}{\mathrm dt}\). Odsud \[k_2\frac {\mathrm dh}{\mathrm dt}=\frac {\mathrm dV}{\mathrm dt}=-k_1\sqrt h,\] tj. \[\frac {\mathrm dh}{\mathrm dt}=-\frac{k_1}{k_2}\sqrt h\] a pro \(\displaystyle k=\frac{k_1}{k_2}\) má model tvar \[\frac {\mathrm dh}{\mathrm dt}=-k\sqrt h.\]
Pro kužel platí díky konstantnímu úhlu u vrcholu vztah \(\displaystyle V=k_3h^3\) (díky podobnosti je objem přímo úměrný třetí mocnině libovolného délkového parametru) a proto \(\displaystyle \frac {\mathrm dV}{\mathrm dt}=k_3 \times 3h^2 \frac {\mathrm dh}{\mathrm dt}\). Odsud \[3k_3 h^2 \frac {\mathrm dh}{\mathrm dt}=\frac {\mathrm dV}{\mathrm dt}=-k_1\sqrt h,\] tj. \[\frac {\mathrm dh}{\mathrm dt}=-\frac{k_1}{3k_3}h^{-3/2}\] a po přeznačení konstanty má model pro kuželovou nádrž tvar \[\frac {\mathrm dh}{\mathrm dt}=-kh^{-3/2}.\]
Dříve jsme odvodili rovnici \[\frac{\mathrm dh}{\mathrm dt}=-k\sqrt h\] popisující úbytek hladiny vody v nádrži tvaru kvádru, ze které vypouštíme vodu.
Řešení
Ad 1: Nabídneme dvě varianty, pro argumentaci je možno použít kteroukoliv z nich.
Pokud je tedy v nádrži nějaká voda, je jednoznačně dáno, jak bude vytékat a je možné vypočítat, jaká bude v libovolném okamžiku hladina.
Ad 2: Pro \(\displaystyle h=\frac 14 k^2 t^2\) a \(\displaystyle t<0\) dostáváme \[\begin{aligned} \frac{\mathrm dh}{\mathrm dt}&=\frac 14 k^2 \cdot 2t = \frac 12 k^2 t\\ -k\sqrt h&=-k\sqrt{\frac 14 k^2 t^2} = - k \frac 12 |k| \cdot |t| = - k \frac 12 k (-t) = \frac 12 k^2 t \end{aligned} \] a obě strany rovnice jsou stejné. Pro \(\displaystyle h=0\) je dosazení triviální.
Je-li \(\displaystyle h(t_0)=0\), může to být proto, že voda v čase \(\displaystyle t_0\) právě vytekla, nebo proto, že vytekla před hodinou nebo proto, že v nádrži nikdy voda nebyla. Proto je nejednoznačnost přirozená. Například \(\displaystyle h(t)=0\) je řešení odpovídající tomu, že voda v nádrži nikdy nebyla. Funkce \(\displaystyle h(t)=\frac 14 k^2t^2\) pro \(\displaystyle t<0\) odpovídá tomu, že pro \(\displaystyle t<0\) v nádrži voda byla a vytekla v čase \(\displaystyle t=0\).
Hromada sypkého materiálu má tvar kužele. Úhel u vrcholu je konstantní, daný mechanickými vlastnostmi materiálu a je nezávislý na objemu. Předpokládejme, že personál stavebnin přisypává na hromadu materiál konstantní rychlostí (v jednotkách objemu za jednotku času). Tato hromada je však v poměrně otevřené krajině a vítr rozfoukává materiál po okolí. Je rozumné předpokládat, že rozfoukávání (opět v jednotkách objemu za jednotku času) se děje rychlostí úměrnou povrchu návětrné strany pláště. Vyjádřete proces kvantitativně pomocí derivací. Napište rovnici pro derivaci objemu hromady podle času.
Toto je podobný model jako model vypouštění nádrže, ale kratší. Opět máme po přepisu zadání do matematického modelu dvě veličiny měnící se s časem v jedné rovnici. Derivace objemu, která nás zajímá, již v rovnici přítomna naštěstí je. Stačí vyjádřit obsah pomocí objemu, nejlépe pomocí rozměrové analýzy.
Řešení
Rychlost s jakou se mění objem je \(\displaystyle \frac{\mathrm dV}{\mathrm dt}\), rychlost přisypávání označme \(\displaystyle R\), povrch návětrné strany \(\displaystyle S\). Podle zadání platí \[ \frac{\mathrm dV}{\mathrm dt} = R - k_0S.\] Protože kužel má stále stejný tvar, objem jednoznačně determinuje rozměry, povrch kužele, nebo i povrch poloviny pláště, tj. povrch návětrné strany. Z rozměrové analýzy na základě Buckinghamova Pi-teorému z přednášky je zřejmé, že musí platit úměrnost mezi takovými mocninami těchto veličin, pro které jednotky “pasují”, Existuje tedy konstanta taková, že \[S=k_1V^{\frac 23}.\] Spojením těchto dvou vztahů dostáváme \[ \frac{\mathrm dV}{\mathrm dt} = R - k V^{\frac 23},\] kde \(\displaystyle r\) a \(\displaystyle k=k_0k_1\) jsou konstanty.
Hromada sypkého materiálu má tvar kužele. Úhel u vrcholu je konstantní, daný mechanickými vlastnostmi materiálu a je nezávislý na objemu. V předchozím příkladě jsme sestavili diferenciální rovnici popisující růst hromady ve tvaru \[\frac{\mathrm dV}{\mathrm dt}=R-kV^{\frac 23},\] kde \(\displaystyle R\) je rychlost přisypávání a \(\displaystyle k\) konstanta.
Řešení
Označme \(\displaystyle f(V)=R-kV^{\frac 23}\). Konstantní řešení je řešením rovnice \(\displaystyle f(V)=0\), tj. \[R-kV^{\frac 23}=0.\] Odsud \[V_0=\left(\frac{R}{k}\right)^{3/2}.\] Protože \(\displaystyle f\) klesá v bodě \(\displaystyle V_0\), je toto řešení stabilní.
Protože \(\displaystyle f(0)>0\), malá hromada vždy roste a proto nemůže skončit celá rozfoukaná. Pro malý objem je přisypávání intenzivnější než rozfoukávání.
Protože \(\displaystyle f\) je pro velké \(\displaystyle V\) záporná, pro velkou hromadu objem ubývá (více se rozfouká než přisype) a hromadu není možné navršit libovolně velkou.