Diferenciální rovnice

Obecné řešení

Najděme obecné  řešení diferenciální rovnice y= MA-DiferencRov_1.gifx

Definujme funkci dvou proměnných f(x,y) = MA-DiferencRov_2.gifx

MA-DiferencRov_3.gif

MA-DiferencRov_4.gif

Nyní  najdeme obecné řešení diferenciální rovnice.

MA-DiferencRov_5.gif

MA-DiferencRov_6.gif

Obecné řešení  MA-DiferencRov_7.gif, x R, MA-DiferencRov_8.gifR

Partikulární řešení

Najděme partikukární  řešení diferenciální rovnice y= MA-DiferencRov_9.gifx
splňující počáteční podmínku y(1) = -MA-DiferencRov_10.gif

MA-DiferencRov_11.gif

MA-DiferencRov_12.gif

Partikulární řešení  MA-DiferencRov_13.gif, x R.

Integrální křivka

Nakresleme integrální  křivku předchozího partikulárního řešení

MA-DiferencRov_14.gif

MA-DiferencRov_15.gif

MA-DiferencRov_16.gif

MA-DiferencRov_17.gif

Aplikační příklad

Izotermní vnitřní difuze v porézním katalyzátoru je popsána diferenciální rovnicí MA-DiferencRov_18.gif MA-DiferencRov_19.gif= MA-DiferencRov_20.gif MA-DiferencRov_21.gif s okrajovými podmínkami y ’(0)=0 a y(1)=1,  kde
y(x) je koncentrace,
a charakterizuje tvar částice katalizátoru (a=0,1,2 pro desku, váleček a kuličku)
n je řád reakce
φ Thieleho modul.
Nalezněme partikulární řešení pro parametry
n=1, a=2, φ=1. Nakresleme integrální  křivku  partikulárního řešení

MA-DiferencRov_22.gif

Nejdříve  najdeme obecné řešení diferenciální rovnice.

MA-DiferencRov_23.gif

MA-DiferencRov_24.gif

Partikulární řešení diferenciální rovnice, splňující dané okrajové podmínky.

MA-DiferencRov_25.gif

MA-DiferencRov_26.gif

Definujeme partikulární řešení jako funkci

MA-DiferencRov_27.gif

MA-DiferencRov_28.gif

Partikulární řešení  je  MA-DiferencRov_29.gif, x R \ {0}.

Integrální křivka partikulárního řešení diferenciální rovnice. Podmínka y ’(0)=0 je splněna limitně, stejně tak hodnotu y(0) musíme vypočítat limitně.

MA-DiferencRov_30.gif

MA-DiferencRov_31.gif

MA-DiferencRov_32.gif

MA-DiferencRov_33.gif

MA-DiferencRov_34.gif

MA-DiferencRov_35.gif


Diferenciální rovnice -  numerické (přibližné) řešení

Partikulární řešení

Najděme partikukární  řešení diferenciální rovnice y=MA-DiferencRov_36.gif
splňující počáteční podmínku y(0) = 1

MA-DiferencRov_37.gif

MA-DiferencRov_38.gif

Nejdříve zkusíme najít partikulární řešení diferenciální rovnice analyticky.

MA-DiferencRov_39.gif

MA-DiferencRov_40.gif

Protože  jsme přesné řešení nenašly, najdeme numerické řešení diferenciální rovnice, protože víme, že řešení existuje.

MA-DiferencRov_41.gif

MA-DiferencRov_42.gif

Nedostaneme analytické řešení, ale můžeme si vypočítat přibližnou hodnotu v libovolném bodě z intervalu <0,1>.

MA-DiferencRov_43.gif

MA-DiferencRov_44.gif

MA-DiferencRov_45.gif

MA-DiferencRov_46.gif

Nebo si můžeme  nakreslit přibližnou integrální křivku.

MA-DiferencRov_47.gif

MA-DiferencRov_48.gif

Aplikační příklad

Izotermní vnitřní difuze v porézním katalyzátoru je popsána diferenciální rovnicí MA-DiferencRov_49.gif MA-DiferencRov_50.gif= MA-DiferencRov_51.gif MA-DiferencRov_52.gif s okrajovými podmínkami y ’(0)=0 a y(1)=1,  kde
y(x) je koncentrace,
a charakterizuje tvar částice katalizátoru (a=0,1,2 pro desku, váleček a kuličku)
n je řád reakce
φ Thieleho modul.
Nalezněme partikulární řešení pro parametry
n=2, a=2, φ=1. Nakresleme integrální  křivku  partikulárního řešení

MA-DiferencRov_53.gif

Numerické partikulární řešení diferenciální rovnice. Úloha je okrajová a nelze jednoduše numericky řešit. Musíme použít metodu střelby. Stačí nám výsledek s přesností na 5 desetinných míst.

MA-DiferencRov_54.gif

Můžeme si vypočítat přibližnou hodnotu v libovolném bodě z intervalu <0,1>.

MA-DiferencRov_55.gif

MA-DiferencRov_56.gif

MA-DiferencRov_57.gif

MA-DiferencRov_58.gif

MA-DiferencRov_59.gif

MA-DiferencRov_60.gif

Můžeme si nakreslit přibližnou integrální křivku.

MA-DiferencRov_61.gif

MA-DiferencRov_62.gif

MA-DiferencRov_63.gif

MA-DiferencRov_64.gif

Přesvědčíme se, že jsme dostali “korektní řešení”. Malinko změníme okrajovou podmínku.

MA-DiferencRov_65.gif

Můžeme si nakreslit přibližnou integrální křivku.

MA-DiferencRov_66.gif

MA-DiferencRov_67.gif

MA-DiferencRov_68.gif

Zobrazíme si rozdíl obou řešení.

MA-DiferencRov_69.gif

MA-DiferencRov_70.gif

MA-DiferencRov_71.gif

Vzhledem k tomu, že jsme požadovali řešení s přesností MA-DiferencRov_72.gif dá se předpokládat, že řeření je “korektní”. Při malé změně počáteční podmínky se řešení moc nezměnilo.

Spikey Created with Wolfram Mathematica 9.0