B.2.6 Ukázka možných způsobů řešení jednoduchých eliptických parciálních diferenciálních rovnic
Následující řešené příklady demonstrují některé základní způsoby počítání eliptických PDR:
Laplaceova rovnice
Laplaceova rovnice je v kartézských souřadnicích v nejjednodušší formě definovaná ve tvaru
V tomto příkladu budeme řešit Laplaceovu rovnici na obdélníkové oblasti s rozměry a,b, se smíšenými Dirichletovými a Neumannovými podmínkami (viz odstavec 3.2.1), v podobě
Separací proměnných u(x,y)=X(x)Y(y) dostáváme rovnici
kde konstanta λ může nabývat hodnot λ=0,λ>0,λ<0.
λ=0: Budeme předpokládat separované funkce X(x) a Y(y) ve tvaru polynomů, vzhledem k okrajovým podmínkám budou dostatečné polynomy 1. stupně, tedy X(x)=Ax+B,Y(y)=Cx+D. Z okrajové podmínky u(0,y)=0 vyplývá B=0∨C=D=0, pokud ovšem C=D=0, potom Y(y)=0 a tedy u(x,y)=0 všude. Pokračujeme-li s B=0, dostáváme Ax(Cy+D)=0 a tedy, zahrneme-li další okrajovou podmínku u(x,0)=0, musí být AxD=0. Případ A=0 dává X(x)=0, tedy u(x,y)=0 všude. Uvažujeme-li také D=0, potom u(x,y)=AxCy=0, z další okrajové podmínky ux(a,y)=0 vyplývá ACy=0, tedy A=0∨C=0, v obou případech ovšem u(x,y)=0. Případ λ=0 dává tedy pouze triviální řešení.
λ>0: Z rovnice B.163 dostáváme obecné řešení ve tvaru
u(x,y)=[Acosh(√λx)+Bsinh(√λx)][Ccos(√λy)+Dsin(√λy)].B.164Z okrajové podmínky u(0,y)=0 vyplývá A[Ccos(√λy)+Dsin(√λy)]=0, tedy C=D=0, potom ovšem Y(y)=0 a tedy u(x,y)=0 všude. Pokračujeme-li s A=0, dostáváme u(x,y)=Bsinh(√λx)[Ccos(√λy)+Dsin(√λy)]. Zahrneme-li další okrajovou podmínku u(x,0)=0, musí být BCsinh(√λx)=0. Případ B=0 dává X(x)=0, tedy u(x,y)=0 všude. Pokračujeme-li s C=0, potom u(x,y)=BDsinh(√λx)sin(√λy) a z další okrajové podmínky ux(a,y)=0 vyplývá √λBDcosh(√λa)sin(√λy)=0, tedy B=0∨D=0, v obou případech ovšem u(x,y)=0. Případ λ>0 dává tedy také pouze triviální řešení.
λ<0: Z rovnice B.163 dostáváme obecné řešení ve tvaru
u(x,y)=[Acos(√λx)+Bsin(√λx)][Ccosh(√λy)+Dsinh(√λy)].B.165Z okrajové podmínky u(0,y)=0 vyplývá A[Ccosh(√λy)+Dsinh(√λy)]=0, tedy C=D=0, potom Y(y)=0 a tedy u(x,y)=0 všude. Pokračujeme s A=0, dostáváme u(x,y)=Bsin(√λx)[Ccosh(√λy)+Dsinh(√λy)]. Zahrneme-li další okrajovou podmínku u(x,0)=0, musí být BCsin(√λx)=0. Případ B=0 dává X(x)=0, tedy u(x,y)=0 všude. Pokračujeme-li s C=0, potom u(x,y)=BDsin(√λx)sinh(√λy) a z další okrajové podmínky ux(a,y)=0 vyplývá √λBDcos(√λa)sinh(√λy)=0, tedy B=0∨D=0 (v obou případech ovšem u(x,y)=0∨cos(√λa)=0.
Poslední případ dává řešení
cos(√λa)=0, tedy √λ=(2k−1)π2a, a tedy
u(x,y)=∞∑k=1Kksin[(2k−1)π2ax]sinh[(2k−1)π2ay], kde Kk=BD.B.166
Uplatníme-li také čtvrtou okrajovou podmínku, dostáváme
Z argumentů funkce sinus vyplývá řešení pouze pro k=1, tedy K1=u0/sinh[πb/(2a)]. Výsledné řešení se zahrnutím všech okrajových podmínek bude
Další typickou eliptickou parciální diferenciální rovnicí může být například tzv. Poissonova rovnice typu Δu(x,y)=f(x,y), tedy nehomogenní eliptická rovnice, nejčastěji používaná veformě gravitační Poissonovy rovnice, ΔΦ=4πGρ, kde Φ je gravitační potenciál, ρ je hustota hmoty a G je gravitační konstanta, nebo Poissonovy rovnice elektrostatického potenciálu, ΔΦ=−ρ/ϵ, kde ρ je hustota elektrického náboje a ϵ je permitivita. Řešení vícerozměrné Poissonovy rovnice je analogické k řešení Laplaceovy rovnice a také například k řešení nehomogenní hyperbolické parciální diferenciální rovnice:
Poissonova rovnice s konstantní pravou stranou
Řešme jednoduchou rovnici, definovanou na oblasti x>y2, tj. na oblasti ohraničené parabolou x=y2 s vrcholem v bodě (0,0), jejíž osu tvoří kladná část osy x,
s Dirichletovou podmínkou na hranici oblasti, u(y2,y)=0. Předpokládejme řešení ve formě všech členů polynomu 2. stupně s neurčitými koeficienty,
kdy po jeho parciálním derivování ve smyslu rovnice B.169 snadno zjistíme: F=1−D. Po dosazení okrajové podmínky, tedy z rovnice
dostáváme nenulové hodnoty koeficientů pouze pro B=−1,F=1. Hledaná rovnice tedy bude
Poissonova rovnice s konstantní pravou stranou na kruhové oblasti, s nehomogenní okrajovou podmínkou
Uvnitř kruhové oblasti s poloměrem R platí následující rovnice,
kdy na hranici oblasti platí Dirichletova podmínka u(x,y1)=1, z níž vyplývá y1=±√R2−x2. Analogicky k parabolickým rovnicím s nehomogenními okrajovými podmínkami rozdělíme hledanou funkci u(x,y) na součet dvou funkcí, například U(x,y) a v(x,y), pro které bude platit:
Uxx+Uyy=4,
U(x,y1)=0,
vxx+vyy=0,
v(x,y1)=1.
Obdobně jako v předchozím případě předpokládáme pro každou funkci úplný polynom 2. stupně s neurčitými koeficienty,
což dává F=2−D, f=−d. Po dosazení okrajové podmínky můžeme rovnice B.175 a B.176 přepsat ve tvaru
jednotlivé nenulové koeficienty budou: A=−R2,D=1,F=1,a=1. Po sečtení rovnic B.177 a B.178 dostáváme hledanou výslednou funkci
Poissonova rovnice s obecnou pravou stranou, Newtonovy okrajové podmínky
Řešme obdobným způsobem Poissonovu rovnici ve tvaru
s okrajovými podmínkami u(0,y)=y2,ux(0,y)=0. Abychom po derivování dostali členy potřebného stupně, musíme nyní ovšem předpokládat řešení ve tvaru úplného polynomu 4. stupně s neurčitými koeficienty,
Příslušné druhé derivace tedy v tomto případě budou
Pro jednotlivé koeficienty dostáváme
Dosazením Dirichletovy okrajové podmínky dostáváme A=0,C=0,F=1,J=0,Q=0, z relací B.184 ihned vyplývá D=−4,H=0,K=0,M=0. Dosazením Neumannovy okrajové podmínky dostáváme B=0,E=0,I=0,N=0, z relací B.184 následně vyplývá G=0,L=1/6. Dosazením nenulových koeficientů do rovnice B.181 dostáváme hledanou výslednou funkci,
Podrobně je tato problematika popsána např. v učebnici Franců (2011).