Metoda linija za rješenje enodimenzionalne valovne jednačine podvrgnuta integralnom stanju

konzervacije
Autor povezuje otvoreni panel zaoverlap Fatemeh Shakeri Mehdi Dehghan
Pokazati više

https://doi.org/10.1016/j.camwa.2008.03.055Preuzmite prava i sadržaj

Pod licencom korisnika Elsevier
otvorena arhiva

Sažetak
Hiperbolične parcijalne diferencijalne jednačine sa integralnim uslovima služe kao modeli u mnogim
granama fizike i tehnologije. Nedavno je puno pažnje proučavano u proučavanju ovih jednačina i za njih
je postojala značajna matematička interesa. U ovom radu, rešenje jednodimenzionalne nelokalne
hiperboličke jednačine prikazano je metodom linija. Metoda linija (MOL) je opšti način gledanja
parcijalne diferencijalne jednačine kao sistema običnih diferencijalnih jednačina. Parcialni derivati u
odnosu na prostorne varijable diskretizirani su da bi dobili sistem ODE-ova u vremenskoj varijabli, a
zatim se može koristiti odgovarajući softver početne vrijednosti za rješavanje ovog ODE
sistema. Predlažemo dva oblika MOL-a za rešavanje opisanog problema.
 Prethodni članak u emisiji
 Sledeći članak u emisiji
Ključne riječi
Hiperbolična parcijalna diferencijalna jednačina
Ne klasična granična situacija
Metoda linija
Sistem običnih diferencijalnih jednačina
Integralno stanje konzervacije

1 . Uvod
Tokom poslednjih nekoliko godina postalo je sve očiglednije da se mnogi fizički fenomeni mogu opisati u
smislu hiperboličnih parcijalnih diferencijalnih jednačina sa integralnim uslovima koji zamenjuju klasični
granični uslov [1] . Ova vrsta jednačina proizlazi, na primjer, u proučavanju termoelastičnosti [2,3] ,
fizike plazme [4] , hemijske heterogenosti [5,6] itd. Sve više pažnje se posvećuje razvoju, analizi i
primjeni numeričkih metoda za rešavanje ovih problema. Nekoliko autora [7,8,1,9-12] su proučavali
hiperbolične početne granične vrijednosti u jednoj dimenziji koja uključuje nlokalne granične uslove.. Za
paraboličke jednačine koje podležu nelokalnim graničnim uvjetima zainteresovani čitač može
videti [13] i reference u njima. Takođe, nonlocalni problemi uključuju probleme sa integralnim izrazom u
inicijalnom stanju [14] . U ovom radu nećemo razgovarati o ovoj grupi.
U ovom istraživanju razmatramo sledeći problem ove porodice jednačina
(1.1)∂2v∂t2-∂2v∂x2=q(x,t),0≤x≤l,0<t≤T,

sa početnim uslovima
(1.2)v(x,0)=f1(x),0≤x≤l,

i
(1.3)vt(x,0)=f2(x),0≤x≤l,

i pogranično stanje Dirichleta

(1.4)v(0,t)=g1(t),0<t≤T,

i nelokalnog stanja
(1.5)∫0lv(x,t)dx=g2(t),0<t≤T,

gde q, f1, f2, g1 i g2poznate su funkcije. Pretpostavljamo da su funkcijeq, f1, f2, g1, i g2zadovoljavaju
uslove kako bi rešenje ove jednačine postojalo i jedinstveno. Postojanje i jedinstvenost rešenja ovog
problema govori se u [15] .
Autor [1] je predstavio nekoliko konačnih razlika šema za numeričko rešenje problema (1.1) - (1.5) . Te
tri nivoa tehnike zasnovane su na dva shema drugog reda (jedna eksplicitna i jedna ponderisana) i tehnika
četvrtog reda (ponderisana eksplicitna) [1] . Takođe u [16]izmenjena tehnika Legendre Tau razvijena je
za rešenje proučavanog modela. Pristup u svom radu sastoji se od smanjenja problema na skup
algebarskih jednačina proširivanjem približnog rješenja kao promjenjene funkcije Legendre sa
nepoznatim koeficijentima. Navedene su integralne i derivativne operativne matrice. Ove matrice zajedno
sa tau metodom se zatim koriste za procenu nepoznatih koeficijenata promjenjenih funkcija
Legendre. Autor [17] je razvio numeričku tehniku zasnovanu na integro-diferencijalnoj jednačini i
lokalnim interpolirajućim funkcijama za rešavanje jednodimenzionalne talasne jednačine koja je predmet
nelokalnog stanja konzervacije i odgovarajuće propisanih početnih graničnih
uslova. Autori [18]kombinovane konačne razlike i spektralne metode za rešavanje jednodimenzionalne
valovne jednačine sa integralnim uslovima. Vremenska varijabla se aproksimira pomoću konačne šeme
razlika. Ali, spektralna metoda se koristi za diskretizovanje prostorne varijable. Glavna ideja koja stoji iza
njihovog pristupa je da se dobiju rezultati visokog reda. Oni su testirali novu metodu za dva primera iz
literature [18] . Numerička šema za rešavanje talasne jednačine drugog reda sa datim početnim uslovima
i graničnim uslovima i integralnim uslovom umesto klasičnog graničnog stanja predstavljena je
u [19]. Za rešavanje modela koristi se kubna funkcija skaliranja B-splina zajedno sa funkcijom skaliranja
granica na intervalu [0, 1]. Dobijeni rezultati pokazuju da ovaj pristup efikasno rešava problem. Vrijedno
je naglasiti da matrice nisu retke, kao one koje se mogu dobiti kada se koriste metode konačne razlike, ali
ovdje [19] trebamo manje matrice da bi se neki rezultati uporedili s metodama konačnih razlika.
U ovom radu koristi se drugačiji pristup, rješenje gornje jednačine izračunava metodom linija. Metoda
linija je alternativni računski pristup koji podrazumijeva aproksimaciju prostornih derivata i
smanjivanjem problema na rješavanje sistema početnih vrijednosti običnih diferencijalnih jednačina i
zatim korištenjem vremenskog integratora za rješavanje ODE sistema. Jednostavno se može povećati
tačnost metoda korišćenjem visoko efikasnih pouzdanih inicijalnih ODE solvera, što znači da se
istovremena porudžbina tačnosti može postići u integraciji vremena bez upotrebe izuzetno malih koraka.
Ovaj rad je organizovan na sledeći način. U odeljku 2 uvodimo metodu linija na kratko i primenimo ga
na jednačine. (1.1) - (1.5) na dva različita načina. Neki numerički rezultati i upoređivanja sa metodom
konačne razlike prikazane u [1] dati su u odeljku 3, a na kraju se zaključuje zaključak u odeljku 4 .

2 . Metoda linija
Metoda linija je polu-diskretna metoda [20-25] koja uključuje smanjenje inicijalne granične vrijednosti u
sistem običnih diferencijalnih jednačina (ODE) u vremenu korištenjem diskretizacije u prostoru. Dobijeni
ODE sistem se može rešiti pomoću standardnog softvera za početnu vrijednost, koji može koristiti
promenljivi pristup vremenski stepen / varijabilni poredak s vremenskom lokalnom kontrolom
greške. Najvažnija prednost pristupa MOL-a je da ne sadrži samo jednostavnost eksplicitnih
metoda [26]ali i superiornost (prednost stabilnosti) implicitnih ako se ne koristi loši numerički metod za
rješavanje ODE-ova. Moguće je postići približne procjene višeg reda u diskretizaciji prostornih derivata
bez značajnog povećanja kompleksne računske vrijednosti. Ova tehnika ima široku primenljivost na
fizičke i hemijske sisteme modela PDE. Modeli koji uključuju rešenje mešanih sistema algebarskih
jednačina, ODE i PDE, rezolucija strmih pokretnih frontova, procena parametara i optimalna kontrola,
drugi problemi kao što su diferencijalne jednačine kašnjenja [27] , dvodimenzionalna sinusna-Gordonova
jednačina [28 ] , jednodimenzionalna proširena Boussinesqova jednačina Nwogu [29], problemi
parcijalne diferencijalne jednačine koji opisuju nelinearne talasne pojave, npr. potpuno nelinearna
jednačina Korteweg-de Vries (KdV) treće reda, Boussinesq jednačina četvrtog reda, Kaup-Kupershmidt
jednačina peta reda i proširena KdV5 jednačina [30] problem nelinearne inverzne toplotne
provodljivosti [31] , problem interfejsa [32] , multi-komponentni model razmnožavanja atmosferskih
zagađivača sa razmatranjem fazne transformacije zagađivača [33] , problem Bingama u cilindričnim
cevima [34] , matematički model za formiranje kapilara [35 ] , trodimenzionalna tranzijentna zračna
prenosna jednačina [36], eliptične parcijalne diferencijalne jednačine koje opisuju stabilnost mase i
transport energije u čvrstim materijama [37] i mnoge druge fizičke probleme.
Autor [38] je koristio MOL da transformiše početni problem granične vrijednosti koji je povezan sa
nelinearnom hiperboličnom Boussinesq jednačinom, u problem prvog reda nelinearne početne
vrijednosti. Numeričke metode se razvijaju zamenom matrično-eksponencijalnog izraza u relativnu
relaciju racionalnih aproksimenata. MOL je zaposlen u [25] za rešavanje jednačine Korteweg-de
Vries. Autori [39] istraživali su performanse različitih pojmova upojanja kako bi pružili neke smernice u
odabiru metoda upadanja u MOL rešenju jako konvektivnih parcijalnih diferencijalnih jednačina. MOL se
koristi u [40] da bi se dobilo numeričko rješenje matematičkog modela za formiranje kapilara u
angiogenezi tumora. Autori [41]primenili su pokretne i adaptivne mreže rešenja za numeričko rešenje
proširenog modela Korteweg-de Vriesa petog reda, motivisanog vodenim talasima u prisustvu površinske
napetosti. Takođe su istraživali dinamiku i interakciju ugrađenih solitona. Autori [42] primenili su MOL
za rešavanje paraboličkih parcijalnih diferencijalnih jednačina. Pokazali su (eksperimentalno) da
povećanje tačaka u formuli u prostornoj aproksimaciji ne može dovesti MOL i tačna rješenja u blisko
saglasnost (za paraboličke probleme, ali ne i nužno za hiperbolične ili eliptičke jednačine). Metoda linija i
metoda konačnih razlika su testirane u [43]sa stanovišta tačnosti rješenja i vremena centralne procesne
jedinice tako što ih primjenjuje na rješavanje vremenskih dvodimenzionalnih Navier-Stokesovih
jednačina za prolazni laminarni tok bez / sa iznenadnim ekspanzijom i upoređivanjem njihovih rezultata
sa numeričkim predviđanjima i merenjima u stalnom stanju prijavljeni u literaturi. Utvrđeno je da je MOL
superiorniji od metode konačne razlike u odnosu na CPU i vreme postavljanja i njegovu fleksibilnost za
inkorporiranje drugih konzervativnih jednačina [43] . Metoda koncepcija linija kombinovana je sa
uklanjanjem prostornih derivata na graničnim elementima kako bi se dobila metoda rešenja za parcijalne
diferencijalne jednačine koje su parabolične u vremenu [44]. Metoda graničnih elemenata ublažava
potrebu za prostornom diskretizacijom i baca problem u integralnom formatu. Stoga su greške povezane
sa numeričkom aproksimacijom prostornih derivata potpuno eliminisane. Primjene metode linija na
hiperbolične parcijalne diferencijalne jednačine sa krutim nelinearnim izvornim izrazima istraživane su
u [45] . Neke od opcija dostupnih za integraciju vremena prilikom upotrebe metoda pokretnih mreža kod
linija se ispituje u [46] . Nova tehnika je predložena u [47]za numeričku integraciju sistema običnih
diferencijalnih jednačina koja se javlja u metodu rešenja linija vremenski zavisnih parcijalnih
diferencijalnih jednačina. Ovaj sistem je obično krut, pa je poželjno da numerički metod reši da ima dobru
osobinu koja se odnosi na stabilnost.
Kako bismo koristili ovaj pristup za rešavanje (1.1) - (1.5) , diskretizujemo koordinatux sa M(M čak i)
ravnomerno razmaknute točke mreže xi=xi-1+h, x0=0, xM=l, i=1,2,...,M. Zapiši toh=l/M i možemo i
pisati xi=ih. U početku, u delu (I), koristimo aproksimaciju drugog reda za drugi derivatx u mrežnim
tačkama xi, i=1,2,...,M-1 iu delu (II), primjenimo aproksimaciju drugog derivata u četvrtom redu x u
mrežnim tačkama xi, i=2,3,...,M-2.
2.1 . MOL I
Neka vi(t) približni v(xi,t). Na slici 1 , redove na kojima se približavajuvi(t)su definisani,
prikazani. Korišćenje aproksimacije srednje razlike drugog reda za drugi derivat ux rezultira
(2.1)d2vidt2-vi-1-2vi+vi+1h2=q(xi,t),i=1,...,M-1,0<t≤T.
1. Preuzmite sliku u punoj veličini
Slika 1 . Obojana traka je domen na kojem je PDE definisan. Aproksimacijevi(t) su definisane duž
isprekidanih linija.

Takođe su uslovi (1.2) - (1.4) svedeni

(2.2)vi(0)=f1(xi),0≤x≤l,
(2.3)(vi)t(0)=f2(xi),0≤x≤l,

i
(2.4)v0(t)=g1(t),0<t≤T.

Primjenjujući Simpsonovo numeričko pravilo za integraciju, aproksimira nelokalno stanje (1.5) na

sljedeći način
(2.5)∫0lv(x,t)dx=hΣi=0Mcivi(t)=g2(t),0<t≤T,

gde c0=cM=13, c2k-1=43 for k=1,...,M2 i c2k=23, za k=1,...,M2-1. Sada, koristeći (2.4) i (2.5) , možemo
dobiti sledeću formulu zavM(t)
(2.6)vM(t)=-1hcM[hΣi=0M-1civi(t)-g2(t)]=-1hcM[hc0g1(t)+hΣi=1M-1civi(t)-g2(t)].

Pre nego što rešimo sistem ODE (2.1) sa početnim uslovima (2.2) i (2.3) , primenjujemo neke
transformacije da konvertujemo ovaj sistem u sistem jednačina prvog reda
u1(t)=dv1dt,
u2(t)=dv2dt,
⋮
uM-1(t)=dvM-1dt,

što rezultira sledećim sistemom 2M-2 jednačine

duidt-vi-1-2vi+vi+1h2=q(xi,t),i=1,...,M-2,0<t≤T,
duM-1dt-vM-2(1-cM-2cM)-vM-1(2+cM-1cM)-1cMΣi=1M-3civi(t)h2=q(xM-1,t)-1h3cM[hc0g1(t)-g2(t)],
dvidt=ui(t),i=1,...,M-1,0<t≤T,

sa početnim uslovima
vi(0)=f1(xi),ui(0)=f2(xi).

2.2 . MOL II
U ovom delu, za drugi derivat u x u mrežnim tačkama xi, i=2,3,...,M-2, koristimo sljedeću aproksimaciju
razlike u četvrtom redu
(vi)xx=-vi+2-16vi+1+30vi-16vi-1+vi-212h2,i=2,...,M-3,
(vM-2)xx=1hcM[hc0g1(t)-g2(t)]+vM-1(16+cM-1cM)-vM-2(30+cM-2cM)12h2+vM-3(16+cM-3cM)-vM-
4(1+cM-4cM)+1cMΣi=1M-5civi(t)12h2.

Koristeći transformacije uvedene u prethodnom dijelu, ODE sistem se dobija na sledeći način
du1dt-v0-2v1+v2h2=q(x1,t),0<t≤T,
duidt+vi+2-16vi+1+30vi-16vi-1+vi-212h2=q(xi,t),i=2,...,M-3,0<t≤T,
duM-2dt-vM-1(16-cM-1cM)-vM-2(30+cM-2cM)12h2-vM-3(16-cM-3cM)-vM-4(1+cM-
4cM)+1cMΣi=1M-5civi(t)12h2=q(xM-2,t)+112h3cM[hc0g1(t)-g2(t)],0<t≤T,
duM-1dt-vM-2(1-cM-2cM)-vM-1(2+cM-1cM)-1cMΣi=1M-3civi(t)h2=q(xM-1,t)-1h3cM[hc0g1(t)-
g2(t)],0<t≤T,
dvidt=ui(t),i=1,...,M-1,0<t≤T.

Sada možemo rešiti rezultujući ODE sistem pomoću ODE rešenja sa vremenskom lokalnom kontrolom
greške. Kao što je pomenuto u [23]osnovna ideja metode linija je zamena prostornih derivata u
parcijalnoj diferencijalnoj jednačini sa algebarskim aproksimacijama. Kada se ovo izvede, prostorni
derivati više nisu eksplicitno navedeni u smislu prostorno nezavisnih varijabli. Stoga, ustvari ostaje samo
varijabla početne vrijednosti, obično vrijeme u fizičkom problemu. Drugim rečima, sa samo jednom
preostalom nezavisnom promenljivom, imamo sistem običnih diferencijalnih jednačina koji približavaju
originalni PDE. Tada je izazov da se formuliše aproksimativni sistem ODE-ova. Kada se to izvede,
možemo primijeniti bilo koji integracioni algoritam za ODE početne vrijednosti za izračunavanje
približnog numeričkog rješenja PDE-u. Dakle, jedna od glavnih karakteristika MOL-a je korišćenje
postojećih, i generalno dobro uspostavljenih, [22] . U trenutnom radu, za rešavanje dobijenog ODE
sistema, koristimo ode45 solver u MATLAB-u koji se zasniva na eksplicitnoj formi Runge-Kutta (4, 5),
paru Dormand-Prince.

3 . Numerički primjeri
U ovom odeljku će se istražiti neki primeri kako bi se pokazalo pouzdanost i efikasnost predloženih šema
u ovom radu.
3.1 . Test 1
Kao prvi primjer, razmotrite Eqs. (1.1) - (1.5) sal=1, T=4,
q(x,t)=0,
f1(x)=0,f2(x)=πcos(πx),
g1(t)=greh(πt),g2(t)=0.

Tačno rešenje ove jednačine jeste

v(x,t)=cos(πx)greh(πt).

Upoređujemo rezultate dobijene procedurom u prethodnom odeljku metodom konačnih razlika uvedenim
u [1] u Tabeli 1 . Za svrhu upoređivanja postavili smoh=0.01 i izračunati apsolutne greške
za t=12 (konačno vreme Tu [1] je12).
Tabela 1 . Računarski rezultati za Primer 1

xi Točna vrijednost Apsolutna greška Apsolutna greška MOL MOL II apsolutna

prema [1] I greška
0.10 0.95105651629515 3,3 × 10 -5 1.5970339260263 × 10 -5 9.1695313453322 × 10 -8
0.20 0.80901699437495 3,0 × 10 -5 2.2782454249470 × 10 -5 8.0261219004285 × 10 -8
0.30 0.58778525229247 3,2 × 10 -5 2.0901699376741 × 10 -5 5.9441112476577 × 10 -8
0.40 0.30901699437495 3,1 × 10 -5 1.2281581424911 × 10 -5 3.3507949526168 × 10 -8
0.50 0.00000000000000 3,3 × 10 -5 1.612449505348880 × 10- 1.087626749329525 × 10-
17 16

0.60 - 3,4 × 10 -5 1.2281581424800 × 10 -5 3.3507949692702 × 10 -8

0.30901699437495
0.70 - 3,1 × 10 -5 2.0901699376519 × 10 -5 5.9441111255332 × 10 -8
0.58778525229247
0.80 - 3,2 × 10 -5 2.2782454249470 × 10 -5 8.0261216561794 × 10 -8
0.80901699437495
0.90 - 3,4 × 10 -5 1.5970339258931 × 10 -5 9.1695316228879 × 10 -8
0.95105651629515
1.00 - 3,2 × 10 -5 8,437694987151190 × 10- 2.664535259100376 × 10-
15 15
1.00000000000000
Takođe, maksimalne norme grešaka koje se razmatraju u nastavku su prikazane u tabelama 2 i 3 .
Tabela 2 . Norma‖vtačno-vMOL‖xi,∞ u Primeru 1

xi MOL I MOL II
0.10 1.9404299881232 × 10 -5 1.0078024292870 × 10 -7
0.20 2.9082139621384 × 10 -5 1.0395382088468 × 10 -7
0.30 2.7979999464800 × 10 -5 9.519132848634 × 10 -8
0.40 1.6949536840560 × 10 -5 6.254244772075 × 10 -8
0.50 1.109874786267721 × 10 -14 2.177963786387436 × 10 -14
xi MOL I MOL II
0.60 1.6949536832678 × 10 -5 6.254244200310 × 10 -8
0.70 2.7979999469574 × 10 -5 9.519131310975 × 10 -8
0.80 2.9082139622272 × 10 -5 1.0395382143980 × 10 -7
0.90 1.9404299884340 × 10 -5 1.0078024492710 × 10 -7
1.00 5.933031843596837 × 10 -13 2.680078381445128 × 10 -13
Tabela 3 . Norma‖vtačno-vMOL‖∞ u Primeru 1

MOL I MOL II
2.9082139622272 × 10 -5 1.0395382143980 × 10 -7
‖vtačno-vMOL‖xi,∞=max0<t≤T|v(xi,t)-vi(t)|.
‖vtačno-vMOL‖∞=max1≤i≤Mmax0<t≤T|v(xi,t)-vi(t)|.

Neki drugi numerički rezultati su prikazani na Sl. 2-4 . Vrijedno je napomenuti da se diskretizacija u
MOL-u primjenjuje samo na prostornu varijablu, što povećava konačno vrijemeT ne smanjuje tačnost
rješenja i možemo dobiti poboljšano rješenje korištenjem preciznijeg integratora za rješavanje ODE
sistema, ali u metodama konačnih razlika obično postoje ograničenja kada se konačno T povećava.

1. Preuzmite sliku u punoj veličini

Slika 2 . Prikaz tačnog rešenja u Primeru 1.
1. Preuzmite sliku u punoj veličini
Slika 3 . Parcelavi(t) for i=5k, k=1,...,10 u Primeru 1, grafikoni od vi(t) for i=5k, k=11,...,19, su
simetrije vi(t), i=5k, k=9,...,1, u odnosu na horizontalnu osu.
1. Preuzmite sliku u punoj veličini
Slika 4 . Parcelav(xi,t)-vi(t) za (a1) xi=14, (a2) xi=12, (a3) xi=34, (a4) xi=1 u Primeru 1.

Kao što se vidi iz Tabele 1 , rezultati dobijeni korišćenjem MOL I imaju generalno istu tačnost kao oni
dobijeni metodom konačne razlike u [1] , ali MOL II je tačniji od ove dve sheme.
3.2 . Test 2
Razmotrimo problem (1.1) - (1.5) sal=1, T=5,
q(x,t)=0,
f1(x)=cos(πx),f2(x)=0,
g1(t)=cos(πt),g2(t)=0.

v(x,t)=12(cos(π(x+t))+cos(π(x-t)))je tačno rešenje ovog sistema. Računarski rezultati koji su dobili MOL
I i MOL II sah=0.01upoređeni sa onima istraženim u [1] u Tabeli 4 . Zbog poređenja, apsolutne greške su
tabelirane nat=14. Iz ove tabele vidimo da je rješenje zasnovano na MOL II tačnije u odnosu na MOL I ili
metod konačne razlike razvijen u [1] . Norme grešaka su rezultirali MOL I i MOL II, grafikoni tačnog
rješenja, približno rješenjevi(t) za više vrijednosti ii greške koje je dobila MOL II data su u tabelama 5 i
6 i Fig. 5-7 respektivno.
Tabela 4 . Računarski rezultati za Primer 2

xi Točna vrijednost Apsolutna greška Apsolutna greška MOL MOL II apsolutna

prema [1] I greška
0.10 0.67249851196396 5,2 × 10 -5 1.2923596914405 × 10 -5 8.6007681532330 × 10 -8
0.20 0.57206140281768 5,1 × 10 -5 1.7467564661144 × 10 -5 8.9389747826019 × 10 -8
0.30 0.41562693777745 5,1 × 10 -5 1.3423564604209 × 10 -5 1.731384369208 × 10 -9
0.40 0.21850801222441 5,3 × 10 -5 7.057215508699 × 10 -6 9.28636417763 × 10 -10
0.50 0.00000000000000 5,0 × 10 -5 1.312806175591734 × 10- 6.396931098296924 × 10-
17 17

0.60 - 5,2 × 10 -5 7.057215509060 × 10 -6 9.28636584296 × 10 -10

0.21850801222441
0.70 - 5,4 × 10 -5 1.3423564604320 × 10 -5 1.731384868808 × 10 -9
0.41562693777745
0.80 - 5,3 × 10 -5 1.7467564658924 × 10 -5 8.9389752599978 × 10 -8
0.57206140281768
0.90 - 5,5 × 10 -5 1.2923596913073 × 10 -5 8.6007680755174 × 10 -8
0.67249851196396
1.00 - 5,4 × 10 -5 3.952393967665557 × 10- 2.198241588757810 × 10-
14 14
0.70710678118655
Tabela 5 . Norma‖vtačno-vMOL‖xi,∞ u Primeru 2
xi MOL I MOL II
0.10 3.204176544324699 × 10 -5 1.984972012314401 × 10 -7
0.20 4.558331840376351 × 10 -5 1.679581814739706 × 10 -7
0.30 4.182611406278181 × 10 -5 1.219577727695764 × 10 -7
0.40 2.458643023162122 × 10 -5 1.137402102918683 × 10 -7
0.50 1.242518910410138 × 10 -14 2.811945523278168 × 10 -14
0.60 2.458643024139118 × 10 -5 1.137402084530614 × 10 -7
0.70 4.182611405878500 × 10 -5 1.219577708821973 × 10 -7
0.80 4.558331840875951 × 10 -5 1.679581880242864 × 10 -7
0.90 3.204176544702175 × 10 -5 1.984971881308084 × 10 -7
1.00 3.066435994014682 × 10 -13 3.248512570053208 × 10 -13
Tabela 6 . Norma‖vtačno-vMOL‖∞ u Primeru 2

MOL I MOL II
4.558331840875951 × 10 -5 1.984972012314401 × 10 -7

1. Preuzmite sliku u punoj veličini

Slika 5 . Prikaz tačnog rešenja u Primeru 2.
1. Preuzmite sliku u punoj veličini
Slika 6 . Parcelavi(t) for i=10k, k=1,...,5 u Primeru 2, grafikoni od vi(t) for i=10k, k=6,...,9, su
simetrije vi(t), i=10k, k=4,...,1, u odnosu na horizontalnu osu.
1. Preuzmite sliku u punoj veličini
Slika 7 . Parcelav(xi,t)-vi(t) za (a1) xi=14, (a2) xi=12, (a3) xi=34, (a4) xi=1 u Primeru 2.

3.3 . Test 3
U ovom primjeru primjenjujemo predložene metode u ovom radu kako bi pronašli rješenje (1.1) -
(1.5) sal=1, T=4,
q(x,t)=-2(x-t)exp(-x-t),
f1(x)=0,f2(x)=xexp(-x),
g1(t)=0,g2(t)=-2texp(-t-1)+texp(-t).

v(x,t)=xtexp(-(x+t))je tačno rešenje ovog problema koji se može lako proveriti. Kao što su prethodni
primeri, rezultati su prikazani u tabelama 7 i 8 i Fig. 8-10 .
Tabela 7 . Norma‖vtačno-vMOL‖xi,∞ u Primeru 3

xi MOL I MOL II
0.10 8.790697957105 × 10 -7 1.3567209786181 × 10 -9
0.20 1.2140191431882 × 10 -6 1.2268780713587 × 10 -9
0.30 1.1133469118643 × 10 -6 9.826361113685 × 10 -10
0.40 7.043162150097 × 10 -7 6.474974688364 × 10 -10
0.50 4.225513185602 × 10 -7 1.327798704320 × 10 -10
0.60 7.709645756804 × 10 -7 7.196475637627 × 10 -10
0.70 1.0072074030434 × 10 -6 1.0519506238316 × 10 -9
0.80 1.2092548875610 × 10 -6 1.2563831364165 × 10 -9
0.90 1.1058150297394 × 10 -6 1.3559769695970 × 10 -9
1.00 1.3828775877189 × 10 -6 2.4905782158857 × 10 -9
Tabela 8 . Norma‖vtačno-vMOL‖∞ u Primeru 3

MOL I MOL II
1.3828775877189 × 10 -6 2.4905782158857 × 10 -9
1. Preuzmite sliku u punoj veličini
Slika 8 . Prikaz tačnog rešenja u Primeru 3.

1. Preuzmite sliku u punoj veličini

Slika 9 . Parcelavi(t) for i=10k, k=1,...,10 u Primeru 3.
1. Preuzmite sliku u punoj veličini
Slika 10 . Parcelav(xi,t)-vi(t) za (a1) xi=14, (a2) xi=12, (a3) xi=34, (a4) xi=1 u Primeru 3.

3.4 . Test 4
U četvrtom primeru rešavamo problem (1.1) - (1.5) sal=1, T=5,
q(x,t)=2x5+2x3-2x2-(20x3+6x-2)(t2-t),
f1(x)=0,f2(x)=-x5-x3+x2,
g1(t)=0,g2(t)=t(t-1)12.

Tačno rešenje ove jednačine jeste

v(x,t)=(x5+x3-x2)(t2-t).

Neki numerički rezultati su rezimirani u tabelama 9 i 10 i Fig. 11-14 .

Tabela 9 . Norma‖vtačno-vMOL‖xi,∞ u Primeru 4
xi MOL I MOL II
0.10 1.4589452381 × 10 -4 2.27205158021 × 10 -8
0.20 2.7172977818 × 10 -4 1.59301172808 × 10 -8
0.30 3.7488135133 × 10 -4 2.98653940467 × 10 -8
0.40 4.3322433196 × 10 -4 4.39071972114 × 10 -8
0.50 4.1659346493 × 10 -4 5.65441311551 × 10 -8
0.60 2.9994415140 × 10 -4 7.65543484160 × 10 -8
0.70 5.996250296 × 10 -5 9.44377208656 × 10 -8
0.80 3.3337581084 × 10 -4 1.134464433505 × 10 -7
0.90 8.9980978039 × 10 -4 1.326242955457 × 10 -7
1.00 1.66653416461 × 10 -3 2.4346786631213 × 10 -6
Tabela 10 . Norma‖vtačno-vMOL‖∞ u Primeru 4

MOL I MOL II
1.66653416461 × 10 -3 2.4346786631213 × 10 -6

1. Preuzmite sliku u punoj veličini

Slika 11 . Prikaz tačnog rešenja u Primeru 4.
1. Preuzmite sliku u punoj veličini
Slika 12 . Parcelavi(t) for i=5k, k=1,...,10 u Primeru 4.

1. Preuzmite sliku u punoj veličini

Slika 13 . Parcelavi(t) for i=5k, k=11,...,20 u Primeru 4.
1. Preuzmite sliku u punoj veličini
Slika 14 . Parcelav(xi,t)-vi(t) za (a1) xi=14, (a2) xi=12, (a3) xi=34, (a4) xi=1 u Primeru 4.

MOL pristup obično nam omogućava relativno lako i sa prihvatljivom efikasnošću da rešimo prilično
opšte i komplikovane parcijalne diferencijalne jednačine. Primjenjuje se na širok spektar problema u
mnogim područjima [47] . Čitač može videti Dodatak [22] za neke probleme u fizici, dinamici fluida,
reaktorskim modelima, automatskoj kontroli i još mnogo toga.

4 . Zaključak
Metoda linija je generalno prepoznata kao sveobuhvatan i moćan pristup numeričkom rješenju
diferencijalnih jednačina parcijalnih varijacija. Ova metoda se nastavlja u dva odvojena koraka: Prostorni
derivati se prvo zamenjuju konačnom razlikom, konačnim zapreminom, konačnim elementom ili drugim
algebarskim aproksimacijama, a zatim dobijeni sistem običnih diferencijalnih jednačina koji je obično
čvrst, integriran je u vremenu. Uspeh ove metode objašnjen je dostupnošću visokokvalitetnih numeričkih
algoritama za rješavanje krutih sistema ODE-ova. Ovaj rad istražio je MOL pristup za rešavanje
jednodimenzionalne hiperboličke jednačine s integralnim uslovom [1,48]. Takođe, novi algoritmi su
testirani na nekoliko problema. Rezultati računanja su potvrdili efikasnost, pouzdanost i tačnost ove
procedure. Vrijedno je naglasiti da rješenja za MOL 1 i MOL 2 (konačne razlike drugog i četvrtog reda)
jasno pokazuju poboljšanu tačnost metode višeg reda (MOL 2). Takođe, ove superiorne performanse
postižu se s vrlo malo povećanim računskim naporom (pošto aproksimacija MOL 2 konačne razlike u
četvrtom redu zahtjeva samo pet pojmova u ponderiranom sume u poređenju sa tri termina u MOL
1). Prema tome, preporučuje se upotreba konačnih razlika višeg reda u MOL-ovom rešenju jednačine
talasa, i generalno za druge parcijalne diferencijalne jednačine. [49] . Drugi zadatak je primijeniti tehniku
predloženu u tekućem istraživanju kako bi pronašla rješenja problema proučavanih u [50-53] .

Zahvale
Autori su veoma zahvalni recenzentu A za pažljivo čitanje rada i njegovih (njenih) konstruktivnih
komentara i sugestija koje su poboljšale rad.

Reference
[1]
M. DehghanO rešenju problema početne granice koji kombinuje Neumann i integralni uslov
za jednačinu talasa
Numer. Metode parcijalnih diferencijalnih jednačina , 21 ( 2005 ) , str. 24 - 40
CrossRefView Record in Scopus
[2]
DE CarlsonLinearna termoelastičnost
Enciklopedija fizike, vol. 2 , Springer , Berlin ( 1972 )
[3]
WA DayA smanjuje svojstvo rješenja paraboličke jednačine sa primjenama na
termoelastičnost i druge teorije
Quart. Appl. Matematika. , 4L ( 1983 ) , str. 468 - 475
CrossRefView Record in Scopus
[4]
AA SamarskiiNeki problemi moderne teorije diferencijalnih jednačina
Razlike. Uravn. , l6 ( 1980 ) , str. 122l - 1228
[5]
JH Cushman , TR GinnNeklokalni disperzioni neredovni mediji sa kontinuirano razvijanjem
vaga heterogenosti
Transport Porous Media , l3 ( 1993 ) , str. 123 - 138
CrossRefView Record in Scopus
[6]
JH Cushman , H. Xu , F. DengNonlocalni reaktivni transport sa fizičkom i hemijskom
heterogenosti: greška lokalizacije
Resurs resursa res. , 31 ( 1995 ) , str. 2219 - 2237
CrossRefView Record in Scopus
[7]
A. BouzianiMešani problem za određene nonclassičke jednačine koji sadrže mali parametar
Acad. Roj. Belg. Bik. CI, Sci. , 6 ( 5 ) ( 1994 ) , str. 389 - 400 hotela
View Record in Scopus
[8]
A. Bouziani Jasnorešenje za mešani problem sa nelokalnim uslovima za određene
pluriparaboličke jednačine
Hirošima matematika. J. , 27 ( 3 ) ( 1997 ) , str. 373 - 390
View Record in Scopus
[9]
IS Gordeziani , GA AvalishviliO konstrukciji rešenja nelokalnih inicijalnih graničnih
problema za jednodimenzionalne jednačine srednjih oscilacija
Matem. Modelirovanie , 12 ( 1 ) ( 2000 ) , str. 94 - 103
View Record in Scopus
[10]
NI Kavalloris , DS TzanetisPonašanje kritičnih rešenja nelokalnog hiperboličnog problema u
ohmicnom zagrevanju hrane
Appl. Matematika. E-Notes , 2 ( 2002 ) , str. 59 - 65
[11]
S. Mesloub , A. BouzianiNa klasu singularne hiperbolne jednačine sa ponderiranim
integralnim uslovima
Internat. J. Math. Matematika. Sci. , 22 ( 1999 ) , str. 511 - 520
CrossRefView Record in Scopus
[12]
LS PulkinaOsolvabilnostiuL2 nelokalnog problema sa integralnim uslovima za hiperbolične
jednačine
Differeti. Uravn. VN , 2 ( 2000 ) , str. 1 - 6
[13]
M. DehghanRačunska studija jednodimenzionalne paraboličke jednačine koja je predmet
neklasičnih graničnih specifikacija
Numer. Metode parcijalnih diferencijalnih jednačina , 22 ( 2006 ) , str. 220 - 257
CrossRefView Record in Scopus
[
1
4
]
M. DehghanImplicitna tehnika kolokacije za jednačinu toplote saneklasicniminicijalnim
uslovima
Int. J. Nelinearna Sci. Numer. Simul. , 7 ( 2006 ) , str. 447 - 450
View Record in Scopus
[
1
5
]

SA BeilinPostojanje rešenja za jednodimenzionalne talasne jednačine sa nelokalnim

uslovima
Electron. J. diferencijalne jednadžbe , 76 ( 2001 ) , str. 1 - 8
View Record in Scopus

A. Saadatmandi , M. DehghanNumeričko rešenjeenodimenzionalne valovnejednačine sa

integralnim uslovima
Numer. Metoda parcijalne diferencijalne jednačine , 23 ( 2007 ) , str. 282 - 292
CrossRefView Record in Scopus

WT Ang. Numerička metoda za jednačinu talasa pod uslovom ne-lokalnog stanja

konzervacije
Appl. Numer. Matematika. , 56 ( 2006 ) , str. 1054 - 1060
ČlanakDownload PDFView Record u Scopus-u

M. Ramezani , M. Dehghan , M. RazzaghiKombinovane konačne razlike i spektralne metode

za numeričko rešenje hiperbolične jednačine sa integralnim uslovima
Numer. Metode parcijalnih diferencijalnih jednačina , 24 ( 2008 ) , str. 1 - 8
CrossRefView Record in Scopus

M. Dehghan , M. LakestaniUpotreba kubnih funkcija za skaliranje B-splinaza rešavanje

jednodimenzionalne hiperbolične jednačine sa neklokalnim uslovima konzerviranja
Numer. Metode parcijalnih diferencijalnih jednačina , 23 ( 2007 ) , str. 1277 - 1289
CrossRefView Record in Scopus

G. Hall , JM Watt (Eds.) , Moderne numeričke metode za obične diferencijalne

jednačine , Clarendon Press , Oxford ( 1976 )

AM Loeb, WE Schiesser, krutost i tačnost u postupku integracije linija parcijalnih diferencijalnih

jednačina, u: Proc. Ljetne kompjuterske simulacije 1973, 2, 1973, str. 25-39

WE SchiesserNumerički metod linija

Academic Press , Njujork ( 1991 )

S. Hamdi, WE Schiesser, GW Griffiths, Metoda linija, iz naučnih

radova. http://www.scholarpedia.org/article/Method_of_Lines
S. Hamdi , WH Enright , WE Schiesser , JJ GottliebTačna rešenja i zakoni o konzervaciji za
spojene generalizovane Korteweg deVries i quintic regulirane jednačine dugog talasa
Nelinearni anal. , 63 ( 2005 ) , str. 1425 - 1434

WE SchiesserMetoda linijskog rješenja jednačine Korteweg-de Vries

Comput. Matematika. Appl. , 28 ( 1994 ) , str. 147 - 154
ČlanakDownload PDFView Record u Scopus-u

M. Dehghan Postupciraznih razlika za rešavanje problema koji nastaju pri modeliranju i

projektovanju određenih optoelektronskih uređaja
Matematika. Comput. Simulacija , 71 ( 2006 ) , str. 16 - 30
ČlanakDownload PDFView Record u Scopus-u

T. KotoMetoda približavanja linija diferencijalnih jednačina kašnjenja

Comput. Matematika. Appl. , 48 ( 2004 ) , str. 45 - 59
ČlanakDownload PDFView Record u Scopus-u

AG BratsosRešenjedvodimenzionalne sinusno-Gordonove jednačine metodom linija

J. Comput. Appl. Matematika. , 206 ( 2007 ) , str. 251 - 277
ČlanakDownload PDFView Record u Scopus-u

S. Hamdi , WH Enright , Y. Ouellet , WE SchiesserMetoda linijskih rešenja proširenih

Boussinesq jednačina
J. Comput. Appl. Matematika. , 183 ( 2005 ) , str. 327 - 342
ČlanakDownload PDFView Record u Scopus-u

P. Saucez , AV Wouwer , WE Schiesser , P. ZegelingMetoda proučavanja linija nelinearnih

disperzivnih talasa
J. Comput. Appl. Matematika. , 168 ( 2004 ) , str. 413 - 423
ČlanakDownload PDFView Record u Scopus-u

J. Taler , P. DudaRešenje nelinearnih inverznih problema toplotne provodljivosti koristeći

metodu linija
Prenos toplotne mase , 37 ( 2001 ) , str. 147 - 155
CrossRefView Record in Scopus

H. Han , Z. HuangDirektan metod linija za numeričko rešenje problema interfejsa

Comput. Metode Appl. Mech. Engrg. , 171 ( 1999 ) , str. 61 - 75
ČlanakDownload PDFView Record u Scopus-u
VD EgorovMetoda linija za rešavanje problema kinetike kondenzacije
J. Aerosol Sci. , 30 ( . Suppl I ) ( 1999 ) ., Str S247 - S248
ČlanakDownload PDFView Record u Scopus-u

G. Torres , C. TurnerMetoda pravih linija za problem Bingama u cilindričnim cevima

Appl. Numer. Matematika. , 47 ( 2003 ) , str. 543 - 558
ČlanakDownload PDFView Record u Scopus-u

A. Erdem , S. PamukMetoda linija za numeričko rešenje matematičkog modela za

formiranje kapilara: Uloga tumorskog angiogenog faktora u vanceličnoj matrici
Appl. Matematika. Comput. , 186 ( 2007 ) , str. 891 - 897
ČlanakDownload PDFView Record u Scopus-u

I. Ayranci , N. SelçukMOL rješenje DOM-a za prolazni zračni prenos u 3-D rasterskim

medijima
J. Kvantitativna spektroskopija prenosa zračenja , 84 ( 2004 ) ., Str 409 - 422
ČlanakDownload PDFView Record u Scopus-u

VR Subramanian , RE BelaPolu-analitička metoda linija za rešavanje eliptičnih parcijalnih

diferencijalnih jednačina
Chem. Eng. Sci. , 59 ( 2004 ) , str. 781 - 788
ČlanakDownload PDFView Record u Scopus-u

AG BratsosRešenje Boussinesq jednačine koristeći metodu linija

Comput. Metode Appl. Mech. Eng. , 157 ( 1998 ) , str. 33 - 44
ČlanakDownload PDFView Record u Scopus-u

P. Saucez , WE Schiesser , A. Vande WouwerIskrivljavanje u metodi linija

Matematika. Comput. Simulacija , 56 ( 2001 ) , str. 171 - 185
ČlanakDownload PDFView Record u Scopus-u

S. Pamuk , A. ErdenMetoda linija za numeričko rešenje matematičkog modela za formiranje

kapilara: Uloga endotelnih ćelija u kapilarnom
Appl. Matematika. Comput. , 186 ( 2007 ) , str. 831 - 835
ČlanakDownload PDFView Record u Scopus-u

P. Saucez , A. Vande Wouwer , PA ZegelingAdaptivni metod rešenja linija za proširenu

jednačinu Korteweg deVries
J. Comput. Appl. Matematika. , 183 ( 2005 ) , str. 343 - 357
ČlanakDownload PDFView Record u Scopus-u
Amr A. Sharaf , HO BakodahDobra prostorna diskretizacija u metodi linija
Appl. Matematika. Comput. , 171 ( 2005 ) , str. 1253 - 1263
ČlanakDownload PDFView Record u Scopus-u

N. Selcuk , T. Tarhan , S. TanrikuluUpoređivanje metode linija i konačnih razlika rješenja 2-

D Navier-Stokesove jednačine za prolazni laminarni protok cevi
Internat. J. Numer. Metode Eng. , 53 ( 2002 ) , str. 1615 - 1628
CrossRef

PA RamachandranMetoda linija sa graničnim elementima za prolazne probleme difuzije i

reakcije
Numer. Metode parcijalnih diferencijalnih jednačina , 22 ( 2006 ) , str. 831 - 846
CrossRefView Record in Scopus

I. Ahmad , M. BerzinsMOL solveri za hiperbolične PDEs sa izvornim uslovima

Matematika. Comput. Simulacija , 56 ( 2001 ) , str. 115 - 125
ČlanakDownload PDFView Record u Scopus-u

JR CashEfikasni integratori vremena u numeričkoj metodi linija

J. Comput. Appl. Matematika. , 183 ( 2005 ) , str. 259 - 274
ČlanakDownload PDFView Record u Scopus-u
H. Ramos , J. Vigo-AguiarSkoro L-stabilna BDF-tip metoda za numeričko rješenje krutih
ODE-a koje proizilaze iz metoda linija
Numer. Metode parcijalnih diferencijalnih jednačina , 23 ( 2007 ) , str. 1110 - 1121
CrossRefView Record in Scopus

M. Dehghan, A. Saadatmandi, Varijaciono iteracioni metod za rešavanje talasne jednačine pod

uslovom integralnog stanja očuvanja, Chaos Solitons Fraktali, 2008 (u štampi)

M. Dehghan Određivanjeparametara u parcijalnoj diferencijalnoj jednačini iz prekoračenih

podataka
Matematika. Comput. Modeliranje , 41 ( 2005 ) , str. 196 - 213
ČlanakDownload PDFView Record u Scopus-u

M. Dehghan , M. TatariPrepoznavanje nepoznate funkcije u paraboličkoj jednačini sa

prekomerno specificiranim podacima pomoću He's varijacijske iteracijske metode
Chaos Solitons Fraktali , 36 ( 2008 ) , str. 157 - 166
ČlanakDownload PDFView Record u Scopus-u

M. DehghanEfikasne tehnike za paraboličku jednačinu drugog reda podležu nelokalnim

specifikacijama
Appl. Numer. Matematika. , 52 ( 2005 ) , str. 39 - 62
ČlanakDownload PDFView Record u Scopus-u

M. DehghanIdentifikacija vremenski zavisnog koeficijenta u parcijalnoj diferencijalnoj

jednačini koja se podvrgava dodatnom merenju
Numer. Metode za parcijalne diferencijalne jednačine , 21 ( 2005 ) , str. 611 - 622
CrossRefView Record in Scopus

M. Dehghan Jednodimenzionalnatoplotna jednačina koja podleže graničnoj integralnoj

specifikaciji
Chaos Solitons Fraktali , 32 ( 2007 ) , str. 661 - 675
ČlanakDownload PDFView Record u Scopus-u