Fabritseeritud lahendite meetod

Fabritseeritud lahendite meetod (ingl Method of Manufactured Solutions) on metodoloogia osatuletistega differentsiaalvõrrandite süsteemi analüütiliste lahendite leidmiseks. Fabritseeritud lahendite meetodit kasutatakse differentsiaalvõrrandite süsteeme lahendavate numbriliste algoritmide verifitseeriminel.

Differentsiaalvõrrandite numbriliste lahendajate verifitseerimine

Osatuletistega differentsiaalvõrrandisüsteemide numbriline lahendamine arvutiprogrammide abil on tänapäevaste teadusarvutuste oluline osa. Sealjuures on oluline kontrollida, kas kasutatav numbriline algoritm on korrektselt implementeeritud. Viimane tegevus nõuab omakorda teadaolevalt korrektset, analüütilisel kujul lahendit. Selle leidmine praktikas esilekerkivate keerulisemate võrrandite, määramispiirkondade ning ääretingimuste korral ei pruugi võrrandit "toore jõuga" lahendades võimalikuks osutuda. Sümmeetrilistest erijuhtudest tuletatud analüütilised lahendid võivad sealjuures jätta tuvastamata olulisi koodis esinevaid vigu.

Fabritseeritud lahendite meetod lahendab "lahendi leidmise probleemi", fikseerides lahendi ning leides seejärel viimase jaoks sobiva, kuid algse ülesandega piisavalt sarnase differentsiaalvõrrandi. Meetodi kasutamine reaalse arvutiprogrammi puhul võib seega nõuda väiksemaid muudatusi programmi koodis, tõstes seevastu oluliselt hinnangut koodi (matemaatilise) korrektsuse osas.

Meetodi krirjeldus

Fabritseeritud lahendite meetodi idee võib esitada järgmiselt. Olgu meil differentsiaalvõrrand kujul

$L[u(x,y,z,t]=0$

millele otsitakse numbriliselt korrektset, kuid mitte ilmtingimata sisulise (füüsikalise) tähendusega lahendit $u=M(x,y,z,t)$. Selleks valime differentsiaalvõrrandi operaatori $L$ asemele viimasega sarnase operaatori $L'$ nii, et eelnevalt valitud lahend oleks saadud differentsiaalvõrrandi lahendiks:

$L'[u(x,y,x,t]=0$

Lihtsaim viis $L'$ leidmiseks on teisendus:

$L'=L-Q$

kus $Q=L[M]$

Ääre- ja algtingimused määratakse fabritseeritud lahendi $M$ pidevust arvestades.

Näide lahendi fabritseerimisest diffusioonivõrrandile

Näitena vaatleme meetodi rakendamist ajast sõltumatu kahemõõtmelise diffusioonivõrrandi $u_{xx}+u_{yy}=0$ lahendamisele Jacobi iteratsioonimeetodiga. Valime fabritseeritud lahendiks $M(x,y)=2x^{4}-y^{3}$. Differentsiaalvõrrandi operaatoriks on

$L=\partial _{xx}+\partial _{yy}$

seega

$L\left[M(x,y)\right]=L\left[2x^{4}-y^{3}\right]=24x^{2}-6y$

Seega on funktsioon $M(x,y)$ järgmise differentsiaalvõrrandi lahendiks.

$u_{xx}+u_{yy}-\left(24x^{2}-6y\right)=0$

Kui ruudukujulistest elementidest küljepikkusega $h$ koosneval võrel Jacobi meetodiga algset differentsiaalvõrrandit lahendav iteratsioonireegel on järgmine:

$u_{i,j}^{(k+1)}={\frac {1}{4}}\left(u_{i+1,j}^{(k)}+u_{i-1,j}^{(k)}+u_{i,j+1}^{(k)}+u_{i,j-1}^{(k)}\right)$

siis fabritseeritud lahendite meetodi kasutamiseks tuleb iteratsioonireeglisse lisada jääkliige $q(x,y)$

$u_{i,j}^{(k+1)}={\frac {1}{4}}\left(u_{i+1,j}^{(k)}+u_{i-1,j}^{(k)}+u_{i,j+1}^{(k)}+u_{i,j-1}^{(k)}-h^{2}q(ih,jh)\right)$

Fabritseeritud lahendiga testimisel seatakse jääkliikme väärtuseks $q(x,y)=24x^{2}-6y$, algset ülesannet lahendades seevastu $q(x,y)=0$.