Ecuația căldurii
Versiunea actuală a paginii nu a fost încă revizuită de colaboratori experimentați și poate diferi semnificativ de
versiunea revizuită la 13 ianuarie 2019; verificările necesită
5 modificări .
Ecuația căldurii este o ecuație diferențială parțială de ordinul doi care descrie distribuția temperaturii într-o anumită regiune a spațiului și schimbarea acesteia în timp.
Tip de ecuație
În spațiul cu un sistem de coordonate arbitrar, ecuația căldurii are forma
unde este o constantă pozitivă (numărul este difuzivitatea termică ),
este operatorul Laplace și este o funcție a surselor de căldură [1] . Funcția dorită setează temperatura în punctul cu coordonatele la momentul respectiv .
Această ecuație poate fi explicată după cum urmează. Rata schimbării temperaturii în timp este proporțională cu curbura distribuției temperaturii în spațiu (derivata a doua). Cu alte cuvinte, cu cât curbura „cocoașelor” de temperatură în corp este mai mare, cu atât mai rapid are loc egalizarea temperaturii în aceste locuri.
În spațiul cu coordonate carteziene , ecuația căldurii ia forma
Ecuația de conducere a căldurii se numește omogenă dacă , i.e. nu există surse și „chivoare” de căldură în interiorul sistemului.
Problema Cauchy pentru ecuația căldurii
Ecuație omogenă
Luați în considerare problema Cauchy pentru ecuația căldurii
omogene :
unde este funcția inițială , continuă și mărginită pe întreg spațiul, iar funcția dorită este continuă și mărginită pentru și toate valorile argumentului .
Următoarele proprietăți sunt valabile pentru problema Cauchy omogenă [2] :
- Principiul maximului (teorema maximului și minim): Rezolvarea problemei omogene Cauchy satisface inegalitățile pentru toate și . [3]
- Teorema existenței și unicității: Pentru orice soluție a problemei omogene Cauchy există, este unică și depinde continuu de funcția inițială din bandă . Cu alte cuvinte, această problemă Cauchy este bine pusă [4] .
- Nucleul ecuației căldurii este soluția problemei Cauchy pentru ecuația omogenă a căldurii cu condiția inițială , unde este funcția delta Dirac . Arată ca:
unde este
pătratul scalar standard al vectorului . Uneori nucleul ecuației de căldură este numit și soluția sa fundamentală , deși cel mai adesea
soluția fundamentală este înțeleasă ca o funcție care se obține din nucleu prin înmulțirea cu
funcția Heaviside .
- Coincidența formulei pentru nucleul ecuației de căldură cu densitatea distribuției normale cu așteptare matematică zero și dispersie proporțională cu nu este întâmplătoare. Se explică prin faptul că transferul de căldură este asociat cu mișcarea browniană a particulelor, care este descrisă matematic folosind procesul aleator Wiener .
- Integrală Poisson: Într-un spațiu cu coordonate carteziene, soluția problemei Cauchy omogene este dată sub forma unei formule integrale numită integrală Poisson . Și anume, pentru toate există o convoluție față de variabila spațiu a nucleului cu funcția inițială:
- Integrala Poisson definește o soluție unică continuă și mărginită a problemei Cauchy dată (rețineți că există infinit de soluții nemărginite).
- Paradoxul fizic: din formula Poisson rezultă că dacă funcția inițială este egală cu zero peste tot, cu excepția unei zone limitate, de exemplu, dată de condiția , în care este pozitivă, atunci după o perioadă de timp arbitrar mică soluția va fi strict pozitiv în toate punctele din spațiu, cu valori arbitrar mari . Aceasta implică o afirmație paradoxală din punct de vedere fizic că căldura se propagă cu o viteză infinită. Explicația paradoxului este că ecuația căldurii nu descrie cu exactitate procesul fizic real de propagare a căldurii. Practica arată că în majoritatea cazurilor această ecuație oferă încă o aproximare destul de bună [2] .
Ecuație neomogenă
Luați în considerare problema Cauchy pentru ecuația căldurii
neomogene :
În acest caz, integrala Poisson are forma [5] :
Ecuația unidimensională a căldurii
Pentru cazul unei variabile spațiale x (problema încălzirii sau răcirii tijei), ecuația căldurii ia forma
Pentru această ecuație, puteți stabili și rezolva diverse probleme cu valori la limită , una dintre metodele de rezolvare care a fost propusă de matematicianul francez Fourier și îi poartă numele [6]
Metoda de separare a variabilelor (metoda Fourier)
Ecuația de căldură omogenă cu condiții la limită omogene
Luați în considerare următoarea problemă:
Trebuie să găsiți o funcție pentru .
Reprezentăm funcția dorită ca produs
Apoi substituim forma propusă a soluției în ecuația originală, obținem
Să împărțim expresia în :
Deoarece în partea stângă a ecuației avem o funcție care depinde numai de , iar în dreapta - numai de , apoi, fixând orice valoare în partea dreaptă, obținem că pentru orice valoare a părții stângi a ecuației este constantă . În același mod, vă puteți asigura că partea dreaptă este constantă, adică egală cu o anumită constantă (minusul este luat pentru comoditate). Astfel, obținem două ecuații diferențiale liniare obișnuite:
Să acordăm atenție condițiilor la limită ale problemei inițiale și să substituim forma propusă a ecuației în ele, obținem:
de unde ( , deoarece altfel am avea o soluție , și căutăm doar soluții nebanale).
Ținând cont de condițiile la limită obținute, obținem problema Sturm-Liouville :
Soluția sa se reduce la rezolvarea unei ecuații diferențiale liniare și luând în considerare trei cazuri:
-
În acest caz, forma generală a soluției va fi următoarea:
Prin înlocuirea condițiilor la limită, ne asigurăm că soluția va fi , și căutăm doar soluții netriviale, prin urmare, acest caz nu este potrivit.
-
Vedere generală a soluției
Este ușor de observat că nici această opțiune nu ni se potrivește.
-
Vedere generală a soluției
Inlocuim conditiile la limita:
Întrucât căutăm doar soluții non-triviale, nu ni se potrivește, așadar
De aici
Ținând cont de rezultatul găsit , derivăm soluția generală a ecuației diferențiale liniare
Ar trebui să primească un răspuns
Acum totul este gata pentru a scrie soluția problemei inițiale:
Ca rezultat, avem un număr infinit de soluții particulare ale ecuației. Toate aceste soluții particulare sunt liniar independente , adică o combinație liniară a oricărui număr de soluții este egală cu zero numai dacă toți coeficienții lor sunt egali cu zero. Prin urmare, este logic să presupunem că prin însumarea tuturor soluțiilor particulare de la unitate la infinit, vom obține o soluție generală a problemei inițiale.
Rămâne de determinat valoarea constantei (în funcție de ) din condiția inițială
Pentru a determina valoarea lui , este necesar să se extindă funcția într- o serie Fourier :
Primim:
De unde vine soluția generală:
În cursul fizicii matematice , se demonstrează că seria rezultată îndeplinește toate condițiile acestei probleme, adică funcția este diferențiabilă (și seria converge uniform ), satisface ecuația din domeniul definiției și este continuă la punctele limitei acestui domeniu.
Ecuație de căldură neomogenă cu condiții la limită omogene
Luați în considerare următoarea problemă pentru o ecuație neomogenă :
Lăsa
Apoi, folosind relația evidentă , rescriem ecuația originală ca:
Să rezolvăm ultima ecuație liniară neomogenă prin metoda variației constantei . În primul rând, găsim soluția generală a ecuației liniare omogene
În soluția generală, înlocuim constanta cu o variabilă și o înlocuim în ecuația originală.
Din condiția inițială obținem:
Ținând cont de condiția pentru , obținem
pentru că
atunci , evident, este coeficientul seriei Fourier și este egal cu
Ca urmare, formula generală este:
Problemă generală a primei valori la limită
În multe cazuri, este posibil să se rezolve ecuația de căldură neomogenă cu limită neomogene și condiții inițiale
folosind metodele descrise mai sus și următorul truc simplu. Reprezentăm funcția dorită ca o sumă:
Să găsim funcția :
Astfel, problema inițială se reduce la următoarele:
După ce găsim funcția , găsim funcția dorită prin formula
Literatură
În rusă
În engleză
- Crank, J.; Nicolson, P. & Hartree, D. R. (1947), A Practical Method for Numerical Evaluation of Solutions of Partial Differential Equations of the Heat-Conduction Type , Proceedings of the Cambridge Philosophical Society vol. 43: 50–67 , DOI 10.1017/S401700502031700502
- Einstein, Albert (1905), Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen , Ann. Fiz. Leipzig 17 Vol. 322 (8): 549–560 , DOI 10.1002/andp.19053220806
- Evans, L.C. (1998), Ecuații cu diferențe parțiale , Societatea Americană de Matematică, ISBN 0-8218-0772-2
- John, Fritz (1991), Ecuații cu diferențe parțiale (ed. a patra), Springer, ISBN 978-0-387-90609-6
- Wilmott, P.; Howison, S. & Dewynne, J. (1995), The Mathematics of Financial Derivatives: A Student Introduction , Cambridge University Press
- Carslaw, H.S. & Jaeger, J.C. (1959), Conducerea căldurii în solide (ed. a doua), Oxford University Press, ISBN 978-0-19-853368-9
- Thambynayagam, RKM (2011), Manualul de difuzare: soluții aplicate pentru ingineri , McGraw-Hill Professional, ISBN 978-0-07-175184-1
- Perona, P & Malik, J. (1990), Scale-Space and Edge Detection Using Anisotropic Diffusion, IEEE Transactions on Pattern Analysis and Machine Intelligence vol . 12 (7): 629–639
- Unsworth, J. & Duarte, FJ (1979), Difuzia căldurii într-o sferă solidă și Teoria Fourier , Am. J Phys. T. 47 (11): 891–893 , DOI 10.1119/1.11601
Link -uri
Note
- ↑ Tikhonov A.N., Samarsky A.A. Ecuații ale fizicii matematice. - cap. III, § 1. - Orice editie.
- ↑ 1 2 Petrovsky I. G. Prelegeri despre ecuații cu derivate parțiale. - cap. IV, § 40. - Orice editie.
- ↑ Dacă, împreună cu soluțiile mărginite, le considerăm nemărginite, principiul maximului nu este adevărat: mărginirea soluției nu rezultă din mărginirea datelor inițiale. Prin urmare, nu există o soluție unică. Vezi, de exemplu, A. Tychonoff, „Théorèmes d'unicité pour l'équation de la chaleur”, Mat. Sat., 42:2 (1935), 199–216
- ↑ Afirmațiile despre unicitatea și dependența continuă a soluției sunt o simplă consecință a principiului maximului.
- ↑ Erich Miersemann. Partielle Differenzialgleichungen, p. 156 . Preluat la 11 iunie 2015. Arhivat din original la 27 martie 2016. (nedefinit)
- ↑ Tikhonov A.N., Samarsky A.A. Ecuații ale fizicii matematice. - cap. III, § 2. - Orice editie.