Εξίσωση θερμότητας

Η **εξίσωση θερμότητας** είναι παραβολική διαφορική εξίσωση μερικών παραγώγων, η οποία περιγράφει την κατανομή της θερμότητας (ή την μεταβολή της θερμοκρασίας) σε μια ορισμένη περιοχή καθώς μεταβάλλεται ο χρόνος.

## Διατύπωση της εξίσωσης
Για μια συνάρτηση *u*(*x*,*y*,*z*,*t*) τριών χωρικών μεταβλητών (*x*,*y*,*z*) (βλ. καρτεσιανές συντεταγμένες) και μιας μεταβλητής χρόνου *t*, η εξίσωση της θερμότητας είναι

:\frac{\partial u}{\partial t} -\alpha\left(\frac{\partial^2u}{\partial x^2}+\frac{\partial^2u}{\partial y^2}+\frac{\partial^2u}{\partial z^2}\right)=0
:
Πιο γενικά σε οποιοδήποτε σύστημα συντεταγμένων:

\frac{\partial u}{\partial t} - \alpha \nabla^2 u=0 

όπου *α* είναι μια θετική σταθερά, και το Δ ή ∇2 συμβολίζει τον τελεστή Λαπλάς. Στο φυσικό πρόβλημα της μεταβολής της θερμοκρασίας,η συνάρτηση *u*(*x*,*y*,*z*,*t*) δίνει την θερμοκρασία ενώ η *α* είναι η θερμική διαχυτικότητα.  Για τη μαθηματική μελέτη του προβλήματος είναι επαρκές να θεωρήσουμε την περίπτωση  *α*&nbsp;=&nbsp;1.

Παρατηρήστε ότι η εξίσωση της θερμότητας ,η οποία δίνεται από τον πρώτο θερμοδυναμικό νόμο (δηλαδή η αρχή διατήρησης της ενέργειας), μπορεί να γραφεί στην ακόλουθη μορφή (με την υπόθεση ότι δεν υπάρχει κάποια μεταφορά μάζας ή εκπομπή ακτινοβολίας). Αυτή η μορφή της εξίσωσης είναι πιο γενική και ιδιαίτερα χρήσιμη στην αναγνώριση της επιρροή της κάθε ιδιότητας στον κάθε όρο. π.χ *c<sub>p</sub>* ή *<math>\rho</math>*

\rho c_p \frac{\partial T}{\partial t} - \nabla \cdot \left( k \nabla T \right) = \dot q_V 

όπου \dot q_V  είναι η ογκομετρική θερμική ροή.

Η εξίσωση της θερμότητας έχει θεμελιώδη αξία σε διάφορα επιστημονικά πεδία. Στα μαθηματικά, είναι το αρχέτυπο των παραβολικών διαφορικών εξισώσεων με μερικές παραγώγους. Στη Θεωρία πιθανοτήτων, η εξίσωση της θερμότητας σχετίζεται με την μελέτη της κίνησης Brown μέσω της εξίσωσης  Fokker–Planck. Στα οικονομικά μαθηματικά χρησιμοποιείται για τη λύση της διαφορικής εξίσωσης μερικών παραγώγων Black–Scholes. Η εξίσωση διάχυσης, μια πιο γενική μορφή της εξίσωσης της θερμότητας, έρχεται σε σύνδεση με τη μελέτη της χημικής διάχυσης και άλλων σχετικών διεργασιών.

## Γενική Περιγραφή

Ας υποθέσουμε ότι έχουμε μια συνάρτηση *u* η οποία περιγράφει την θερμοκρασία σε ένα δοσμένο σημείο (*x*, *y*, *z*). Η τιμή της συνάρτησης *u* μεταβάλλεται με την πάροδο του χρόνου καθώς η θερμότητα διαδίδεται στον χώρο. Η εξίσωση θερμότητας χρησιμοποιείται για να προσδιορίσουμε την μεταβολή της  συνάρτησης  *u* με το χρόνο.  **Ο ρυθμός μεταβολή της ''u'' είναι ανάλογος με την "καμπυλότητα" της ''u***. 

Στο διπλανό σχήμα περιγράφεται ο τρόπος με τον οποίο μεταβάλλεται η θερμότητα σε μια μεταλλική ράβδο με την πάροδο του χρόνου. Μια από τις ενδιαφέρουσες ιδιότητες της εξίσωσης θερμότητας είναι η [[αρχή του μεγίστου]] η οποία λέει ότι η μέγιστη τιμή της *u* 

Άλλη μια ενδιαφέρουσα ιδιότητα είναι ότι ενώ η *u* έχει ασυνέχεια σε ένα αρχικό (χρονικό) σημείο  *t*&nbsp;=&nbsp;*t*<sub>0</sub>, η θερμοκρασία εξομαλύνεται καθώς  *t*&nbsp;>&nbsp;*t*<sub>0</sub>. Παραδείγματος χάρη ,αν μια μεταλλική ράβδος έχει θερμοκρασία 0 και μια άλλη θερμοκρασία 100 και τις φέρουμε σε επαφή άκρο με άκρο, τότε πολύ γρήγορα η θερμοκρασία στο σημείο επαφής θα φτάσει τους 50 και η γραφική παράσταση της θερμοκρασίας θα διαπεράσει ομαλά τις τιμές από το 0 έως το 100.

Η εξίσωση θερμότητας χρησιμοποιείται στην θεωρία πιθανοτήτων για να περιγράψει τον τυχαίο περίπατο. Για το λόγο αυτό έχει εφαρμογή στα οικονομικά μαθηματικά.

Είναι επίσης σημαντική στη [[γεωμετρία Riemann]] και συνακόλουθα στη [[τοπολογία]]: προσαρμόστηκε (η εξίσωση θερμότητας) από τον Richard Hamilton όταν όρισε τη [[Ricci flow]] η οποία χρησιμοποιήθηκε αργότερα από τον [[Γκριγκόρι Πέρελμαν]] για να αποδείξει την τοπολογική [[η εικασία του Πουανκαρέ|εικασία του Πουανκαρέ]].

== Το φυσικό πρόβλημα και η εξίσωση ==

===Παραγώγιση σε μια διάσταση===

Η εξίσωση θερμότητας προέρχεται από το νόμο του [[Ζοζέφ Φουριέ|Φουριέ]] και από την [[αρχή διατήρησης της ενέργειας]].  
Σύμφωνα με το νόμο του Φουριέ, ο ρυθμός ροής θερμικής ενέργειας ανά μονάδα επιφάνειας μέσω μιας επιφάνειας είναι ανάλογος με την αρνητική βαθμίδα θερμοκρασίας κατά μήκος της επιφάνειας ,

:<math>\mathbf{q} = - k \, \nabla u \ </math>

όπου *k* είναι η [[θερμική αγωγιμότητα]] και *u* είναι η θερμοκρασία.  Σε μία διάσταση, η βαθμίδα θερμοκρασίας είναι μια συνήθης χωρική παράγωγος, κι έτσι ο νόμος του Φουριέ θα είναι:

:<math>q = -k \frac{\partial u}{\partial x} \,</math>

Απουσία παραγωγής έργου, η μεταβολή της [[Εσωτερική ενέργεια|εσωτερικής ενέργειας]] ανά μονάδα όγκου στο υλικό , Δ*Q*, είναι ανάλογη με την μεταβολή της θερμοκρασίας,  Δ*u* (σε αυτό το τμήμα και μόνο, με Δ θα συμβολίζουμε τη διαφορά και όχι τον τελεστή Λαπλάς). Δηλαδή,

:<math>\Delta Q = c_p\rho \, \Delta u \,</math>

όπου *cp* είναι η [[ειδική θερμοχωρητικότητα|θερμοχωρητικότητα]] και *ρ* η πυκνότητα μάζας του υλικού. Θέτοντας την ενέργεια μηδέν στην απόλυτη μηδενική θερμοκρασία, η προηγούμενη εξίσωση μπορεί να γραφεί ως εξής 
:<math>Q = c_p\rho u. \,</math>

Η αύξηση εσωτερικής ενέργειας σε μια μικρή χωρική περιοχή του υλικού

:<math>x-\Delta x \le \xi \le x+\Delta x</math>

κατά τη διάρκεια της χρονικής περιόδου

:<math>t-\Delta t\le \tau \le t+\Delta t</math>

δίνεται από

:<math>c_p\rho \int_{x-\Delta x}^{x+\Delta x} [u(\xi,t+\Delta t)-u(\xi,t-\Delta t)]\, d\xi = c_p\rho\int_{t-\Delta t}^{t+\Delta t}\int_{x-\Delta x}^{x+\Delta x} \frac{\partial u}{\partial\tau}\,d\xi \, d\tau</math>

όπου χρησιμοποιήσαμε το [[θεμελιώδες θεώρημα του λογισμού]].  Απουσία παραγωγής έργου και αν δεν υπάρχουν πηγές θερμότητας ούτε καταβόθρες, η μεταβολή εσωτερικής ενέργειας στο διάστημα [*x*−Δ*x*, *x*+Δ*x*] αντιπρoσωπεύεται εξ ολοκλήρου από τη ροή θερμότητας κατά μήκος του συνόρου της επιφάνειας.  Από το νόμο του Φουριέ ,αυτό σημαίνει

:<math>k\int_{t-\Delta t}^{t+\Delta t}\left[\frac{\partial u}{\partial x}(x+\Delta x,\tau)-\frac{\partial u}{\partial x}(x-\Delta x,\tau)\right]\,d\tau = k\int_{t-\Delta t}^{t+\Delta t}\int_{x-\Delta x}^{x+\Delta x}\frac{\partial^2u}{\partial\xi^2}\,d\xi\, d\tau</math>

λόγω και πάλι του θεμελιώδους θεωρήματος του λογισμού.  Λόγω τώρα της [[αρχή διατήρησης της ενέργειας|αρχής διατήρησης της ενέργειας]] θα έχουμε,

:<math>\int_{t-\Delta t}^{t+\Delta t}\int_{x-\Delta x}^{x+\Delta x} [c_p\rho u_\tau - k u_{\xi\xi}]\, d\xi \, d\tau = 0.</math>

Το οποίο αληθεύει για οποιοδήποτε ορθογώνιο [*t*&nbsp;&minus;Δ*t*, *t*&nbsp;+&nbsp;Δ*t*]&nbsp;×&nbsp;[*x*&nbsp;&minus;&nbsp;Δ*x*, *x*&nbsp;+&nbsp;Δ*x*].  Λόγω του θεμελιώδους λήμματος του [[Λογισμός των μεταβολών|λογισμού των μεταβολών]] , οι ολοκληρωτέες θα πρέπει να μηδενιστούν ταυτοτικά:

:<math>c_p\rho u_t - k u_{xx} = 0. \,\!</math>

Η οποία μπορεί να ξαναγραφεί ως εξής:

:<math>u_t = \frac{k}{c_p\rho}u_{xx},</math>

ή:

:<math>\frac{\partial u}{\partial t} = \frac{k}{c_p\rho} \left(\frac{\partial^2u}{\partial x^2}\right) </math>

Η οποία είναι η εξίσωση θερμότητας, όπου ο συντελεστής (που συνήθως συμβολίζεται με α)

:<math>\alpha = \frac{k}{c_p\rho}\,\!</math>

ονομάζεται θερμικη διαχυτικοτητα.

Ένας πρόσθετος όρος που μπορεί να εισαχθεί στη εξίσωση για να λογαριάσουμε την απώλεια θερμότητας λόγω εκπομπής ακτινοβολίας , ο οποίος εξαρτάται από την υπερβολική θερμοκρασία *u* = *T* − *Ts* σε ένα δοσμένο σημείο συγκρινόμενη με την θερμοκρασία του περιβάλλοντος χώρου. Σε υπερβολικά χαμηλές θερμοκρασίες , η απώλεια θερμότητας είναι περίπου μ*u*, δίνοντας μια μονοδιάστατη εξίσωση μεταφοράς θερμότητας της μορφής
:<math>\frac{\partial u}{\partial t}= \frac{k}{c_p\rho} \left(\frac{\partial^2u}{\partial x^2}\right)-\mu u. </math>

Σε υπερβολικά υψηλές θερμοκρασίες, ωστόσο, ο [[νόμος Στέφαν-Μπόλτζμαν]] δίνει μια καθαρή απώλεια θερμότητας λόγω ακτινοβολίας ,η οποία είναι ανάλογη προς το  <math>T^4-T_s^4</math>, και η παραπάνω εξίσωση είναι ανακριβής. Για μεγάλες θερμοκρασίες έχουμε, <math>T^4 - T_s^4 \approx u^4</math>, το οποίο δίνει μια εξίσωση μεταφοράς θερμότητας υψηλής θερμοκρασίας  της μορφής

:<math>\frac{\partial u}{\partial t}=\alpha\left(\frac{\partial^2u}{\partial x^2}\right)-\ m u^4</math>

όπου <math>\ m= \epsilon\sigma p/\rho Ac_p</math>.  Εδώ, σ είναι η σταθερά Στέφαν-Μπόλτζμαν, ε είναι η χαρακτηριστική σταθερά του υλικού, *p* είναι η περίμετρος διατομής της ράβδου και *A* είναι το εμβαδόν της εγκάρσιας διατομής. Ωστόσο, χρησιμοποιώντας  *T* αντί για *u* παίρνουμε μια καλύτερη προσέγγιση στην περίπτωση αυτή.

===Το πρόβλημα στις τρεις διαστάσεις===
Στις ειδικές περιπτώσεις της κυματικής διάδοσης της θερμότητας σ' ένα ισότροπο και ομογενές μέσο σ' έναν 3-διάστατο χώρο , αυτή η εξίσωση είναι

:<math>{\partial u\over \partial t} =
\alpha \left({\partial^2 u\over \partial x^2 } +
{\partial^2 u\over \partial y^2 } +
{\partial^2 u\over \partial z^2 }\right)</math>&ensp;<math> = \alpha ( u_{xx} + u_{yy} + u_{zz} ) \quad </math>

όπου:

* *u* = *u*(*x*, *y*, *z*, *t*) είναι η θερμοκρασία ως συνάρτηση χώρου και χρόνου·
* <math>\frac{\partial u}{\partial t}</math> είναι ο ρυθμός μεταβολής της θερμοκρασίας σ' ένα σημείο με την πάροδο του χρόνου·
* *uxx*, *uyy*, and *uzz* είναι οι δεύτερες χωρικές [[παράγωγος|παράγωγοι]] (*thermal conductions*) της θερμοκρασίας στις κατευθύνσεις *x*, *y*, και *z*  αντίστοιχα·
* <math>\alpha = \frac{k}{c_p\rho}</math> είναι η θερμική διαχυτικότητα, μια χαρακτηριστική ιδιότητα του υλικού  που εξαρτάται από την *θερμική αγωγιμότητα* *k*, την *πυκνότητα μάζας* *ρ*, και την *ειδική θερμοχωρητικότητα* *cp*.

Η εξίσωση θερμότητα είναι συνέπεια του νόμου του Φουριέ για την αγωγιμότητα (βλέπε [[θερμική αγωγιμότητα]]).

Στην περίπτωση που το μέσο δεν είναι ολόκληρος ο χώρος, πρέπει να ορίσουμε  συνοριακές συνθήκες για *την u* ώστε να έχουμε μονοσήμαντη λύση της εξίσωσης θερμότητας. 

=== Εσωτερική παραγωγή θερμότητας ===
Η παραπάνω συνάρτηση *u* αναπαριστά την θερμοκρασία ενός σώματος. Εναλλακτικά, είναι ορισμένες φορές βολικό να αλλάξουμε τις μονάδες μέτρησης και να αναπαραστήσουμε την  *u* ως την θερμική ροή ενός μέσου. Και από την στιγμή που η θερμική ροή είναι ανάλογη με την θερμοκρασία σε ένα ομογενές μέσο, η εξίσωση θερμότητας έχει ισχύ και στις νέες μονάδες μέτρησης.

Ας υποθέσουμε ότι ένα σώμα υπακούει στην εξίσωση θερμότητας και, επιπλέον παράγει την δικιά του θερμότητα ανά μονάδα όγκου (π.χ., σε watts/litre - W/L) σε ένα ρυθμό που δίνεται από μια γνωστή συνάρτηση *q* που ποικίλει ανάλογα με τον χώρο και τον χρόνο. Τότε η θερμότητα ανά μονάδα όγκου *u* ικανοποιεί μια εξίσωση της μορφής

: <math>\frac{\partial u}{\partial t} = \alpha \left({\partial^2 u\over \partial x^2 } + {\partial^2 u\over \partial y^2 } + {\partial^2 u\over \partial z^2} \right) + \frac{1}{c_p\rho}q.</math>

Για παράδειγμα ,το νήμα πυράκτωσης ενός λαμπτήρα βολφραμίου παράγει θερμότητα, έτσι θα έχει μη μηδενική τιμή για την *q* όταν είναι αναμμένος. Ενώ όταν ο λαμπτήρας είναι  απενεργοποιημένος, η τιμή της *q* για τον νήμα πυράκτωσης θα είναι μηδέν.

== Επίλυση της εξίσωσης θερμότητας χρησιμοποιώντας σειρές Fourier ==


Η ακόλουθη τεχνική λύσης για την εξίσωση θερμότητας προτάθηκε από τον Ζοζέφ Φουριέ στην πραγματεία του *Théorie analytique de la chaleur (Αναλυτική θεωρία της θερμότητας)*, η οποία δημοσιεύτηκε το 1822. Ας θερωρήσουμε την εξίσωση θερμότητας για μια χωρική μεταβλητή. Η οποία μπορεί να προσομοιώσει την θερμική αγωγιμότητα σε μία ράβδο. Η εξίσωση είναι



όπου η  *u* = *u*(*x*, *t*) είναι μια συνάρτηση δύο μεταβλητών *x* και *t*. Εδώ 
* *x* είναι η χωρική μεταβλητή, έτσι *x* ∈ [0, *L*], όπου *L* είναι το μήκος της ράβδου.
* *t* είναι η μεταβλήτη χρόνου, έτσι *t* ≥ 0.

Υποθέτουμε την αρχική συνθήκη



όπου η συνάρτηση  *f*  είναι δοσμένη, και τις συνοριακές συνθήκες



Ας επιχειρήσουμε να βρούμε μια λύση της  που να μην είναι ταυτοτικά μηδέν , να ικανοποιεί της συνοριακές συνθήκες  και να έχει την παρακάτω  ιδιότητα:η *u* είναι ένα γινόμενο στο οποίο στο οποίο η εξάρτηση της *u* από τις *x*, *t* χωρίζεται, δηλαδή:



Αυτή η τεχνική λύσης ονομάζεται χωριζομένων μεταβλητών. Αντικαθιστώντας την *u* πίσω στην εξίσωση  έχουμε,

:<math>\frac{T'(t)}{\alpha T(t)} = \frac{X*(x)}{X(x)}.

Και από την στιγμή που το δεξί μέλος εξαρτάται μόνο από το  *x* και το αριστερό μόνο από το *t*, τα δύο μέλη θα είναι ίσα με μια σταθερή ποσότητα  −λ. Κι έτσι:

και

Θα δείξουμε τώρα ότι δεν υπάρχουν μη τετριμμένες λύσεις για την  για λ ≤ 0 δεν υπάρχουν:

Έτσι λύνεται η εξίσωση θερμότητας στην ειδική περίπτωση όπου η εξάρτηση τηf *u* έχει την ειδική μορφή .

Γενικά, το σύνολο των λύσεων της  που ικανοποιούν τις συνοριακές συνθήκες  ικανοποιεί επίσης την  και την . Μπορούμε να δείξουμε ότι η λύση της ,  και της  δίνεται από την

:u(x,t) = \sum_{n = 1}^{\infty} D_n \sin \left(\frac{n\pi x}{L}\right) e^{-\frac{n^2 \pi^2 \alpha t}{L^2}}

όπου

:D_n = \frac{2}{L} \int_0^L f(x) \sin \left(\frac{n\pi x}{L}\right ) \, dx.

Υπάρχουν και άλλες λύσεις κλειστής μορφής.

### Γενικεύοντας την τεχνική λύσης
Η τεχνική λύσης που χρησιμοποιήσαμε παραπάνω μπορεί κάλλιστα να επεκταθεί σε πολλές άλλου είδους εξισώσεις. Η ιδέα είναι ότι ο τελεστής *u<sub>xx</sub>* με τις μηδενικές συνοριακές συνθήκες μπορεί να παρασταθεί σε όρους των ιδιοδιανυσμάτων του . Αυτό οδηγεί με φυσικό τρόπο σε μια βασική ιδέα της φασματικής θεωρίας ,αυτή των γραμμικών αυτοσυζυγών ή αυτοπροσαρτημένων τελεστών.

Θεωρήστε τον γραμμικό τελεστή Δ*u* = *u<sub>xx</sub>*. Η άπειρη ακολουθία συναρτήσεων

: e_n(x) = \sqrt{\frac{2}{L}}\sin \left(\frac{n\pi x}{L}\right)

για *n* ≥ 1 είναι τα ιδιοδιανύσματα της Δ. Πράγματι

: \Delta e_n = -\frac{n^2 \pi^2}{L^2} e_n. 

Επιπλέον, οποιοδήποτε ιδιοδιάνυσμα  *f* της Δ με τις συνοριακές συνθήκες *f*(0) = *f*(*L*) = 0 είναι της μορφής *e**n* για κάποιο *n* ≥ 1. Οι συναρτήσεις *e**n* για *n* ≥ 1 σχηματίζουν μια ορθοκανονική ακολουθία ως προς ένα συγκεκριμένο εσωτερικό γινόμενο στον χώρο των πραγματικών συναρτήσεων στο διάστημα [0, *L*]. Αυτό σημαίνει

: \langle e_n, e_m \rangle = \int_0^L e_n(x) e^*_m(x) dx = \delta_{mn}

Τελικά, η ακολουθία {*e**n*}*n* ∈ **N** παράγει έναν πυκνό υπόχωρο του *L*2((0, *L*)). Αυτό δείχνει στην πράξη ότι έχουμε διαγωνοποιήσει τον τελεστή Δ.

## Θερμική αγωγή σε μη ομογενή ανισότροπα μέσα

Γενικά, η μελέτη της θερμικής αγωγιμότητας βασίζεται σε ορισμένες αρχές. Η θερμική ροή είναι μια μορφή ενεργειακής ροής, και συνεπώς έχει νόημα να μιλήσουμε για τoν χρονικό ρυθμό ροής θερμότητας σε μια περιοχή του χώρου.

* Ο χρονικός ρυθμός θερμικής ροής προς μια περιοχή *V* δίνεται από μια εξαρτώμενη απ' τον χρόνο ποσότητα *q**t*(*V*). Υποθέτουμε ότι η  *q* έχει μια πυκνότητα Q, έτσι ώστε

:: q_t(V) = \int_V Q(x,t)\,d x \quad 

* Η θερμική ροή είναι μια εξαρτώμενη απ' τον χρόνο διανυσματική συνάρτηση **H**(*x*) η οποία έχει τα ακόλουθα χαρακτηριστικά: ο χρονικός ρυθμός της θερμικής ροής μέσω ενός απειροστού τμήματος μιας επιφάνεια με εμβαδόν *dS* και με το κανονικό μοναδιαίο διάνυσμα  **n** είναι

:: \mathbf{H}(x) \cdot \mathbf{n}(x) \, dS 

:Συνεπώς ο ρυθμός ροής της θερμότητας προς την περιοχή *V* δίνεται επίσης απ' το επιεπιφάνειο ολοκλήρωμα

::  q_t(V)= - \int_{\partial V} \mathbf{H}(x) \cdot \mathbf{n}(x) \, dS 

:όπου  **n**(*x*)  είναι το μοναδιαίο διάνυσμα στο σημείο *x* στραμμένο  προς τα έξω.

* Ο νόμος του Φουριέ δηλώνει ότι η ροή θερμικής ενέργειας έχει την ακόλουθη γραμμική εξάρτηση από την βαθμίδα θερμοκρασίας

:: \mathbf{H}(x) = -\mathbf{A}(x) \cdot \nabla u (x) 

:όπου **A**(*x*) είναι ένας  3&nbsp;×&nbsp;3 πραγματικός πίνακας που είναι συμμετρικός θετικά ορισμένος.

*Από το θεώρημα απόκλισης (βλ. θεώρημα Gauss-Ostrogradsky), το προηγούμενο επιεπιφάνειο ολοκλήρωμα για την ροή θερμότητας προς μια περιοχή  *V* μπορεί  να μετασχηματιστεί στο ολοκλήρωμα όγκου

::\begin{align}
q_t(V) &= - \int_{\partial V} \mathbf{H}(x) \cdot \mathbf{n}(x) \, dS \\
&= \int_{\partial V} \mathbf{A}(x) \cdot \nabla u (x) \cdot \mathbf{n}(x) \, dS \\
&= \int_V \sum_{i, j} \partial_{x_i} \bigl( a_{i j}(x) \partial_{x_j} u (x,t) \bigr)\,dx 
\end{align}

* Ο χρονικός ρυθμός μεταβολής της θερμοκρασίας στο σημείο *x* είναι ανάλογος με την ροή θερμικής ενέργειας προς ένα τμήμα απειροστού όγκου, ;όπου η σταθερά αναλογίας εξαρτάται από μια σταθερά *κ*

:: \partial_t u(x,t) = \kappa(x) Q(x,t)\, 

Συγκεντρώνοντας όλες αυτές τις εξισώσεις, παίρνουμε την γενική εξίσωση θερμικής ροής:

:  \partial_t u(x,t) = \kappa(x) \sum_{i, j} \partial_{x_i} \bigl( a_{i j}(x) \partial_{x_j} u (x,t)\bigr) 

**Παρατηρήσεις**.

* Ο συντελεστής *κ*(*x*) είναι η αντίστροφη της θερμοχωρητικότητας της ουσίας στο σημείο *x* επί την πυκνότητα της ουσίας στο σημείο *x*.
* Στη περίπτωση που το μέσο είναι ισότροπο, ο πίνακας **A** είναι ένας βαθμωτός πίνακας ίσος προς την θερμική αγωγιμότητα.
* Στη περίπτωση που το μέσο είναι ανισότροπο όπου ο πίνακας **A** δεν είναι βαθμωτός (π.χ., εξαρτάται από το *x*), τότε σπάνια μπορούμε να έχουμε έναν αναλυτικό τύπο για την εξίσωση θερμότητας. Παρ ' όλα αυτά ,είναι συνήθως δυνατό να λάβουμε υπ' όψιν μας το σχετικό αφηρημένο πρόβλημα Cauchy και να δείξουμε ότι είναι ένα καλά ορισμένο πρόβλημα και/ή να δείξουμε κάποιες ποιοτικές ιδιότητες (όπως διατήρηση των αρχικών θετικών δεδομένων, άπειρη ταχύτητα διάδοσης, σύγκλιση γύρω από ένα σημείο ισορροπίας, ιδιότητες εξομάλυνσης). Αυτό συνήθως γίνεται με την θεωρία ημιομάδων μονής παραμέτρου: για παράδειγμα, αν ο *A* είναι ένας συμμετρικός πίνακας, τότε ο ελλειπτικός τελεστής που ορίζεται παρακάτω

::Au(x):=\sum_{i, j} \partial_{x_i} a_{i j}(x) \partial_{x_j} u (x)

:είναι αυτοσυζυγής  και σκεδαστικός, συνεπώς λόγω της φασματικής θεωρίας παράγει μια ημιομάδα μονής παραμέτρου.

## Θεμελιώδης λύση
Μία θεμελιώδης λύση, που επίσης λέγεται πυρήνας θερμότητας , είναι μία λύση της εξίσωσης της θερμότητας αντίστοιχη στην αρχική συνθήκη από ένα αρχικό σημείο- πηγή της θερμότητας σε μία γνωστή θέση.Αυτό μπορεί να χρησιμοποιηθεί για να βρούμε μία γενική λύση της εξίσωσης της θερμότητας πάνω στο πεδίο ορισμού΄δες για παράδειγμα ( Evans,1998) για εισαγωγική αναφορά.
Σε μία μεταβλητή,η συνάρτηση του Green είναι η λύση του προβλήματος αρχικής τιμής 
 
	
:\begin{cases}
 
	
u_t(x,t)-k u_{xx}(x,t)=0& (x, t)\in \mathbf{R} \times (0, \infty)\\
 
	
u(x,0)=\delta(x)&
 
	
\end{cases}
 
 		 	
	
όπου δ είναι η συνάρτηση Dirac.Η λύση του προβλήματος είναι η θεμελιώδης λύση
 
	
:\Phi(x,t)=\frac{1}{\sqrt{4\pi kt}}\exp\left(-\frac{x^2}{4kt}\right).  
Σε αρκετές χωρικές μεταβλητές, η θεμελιώδης λύση λύνει τα ανάλογα προβλήματα
 
	
:\begin{cases}
 
	
u_t(\mathbf{x},t) - k \sum_{i=1}^nu_{x_ix_i}(\mathbf{x},t) = 0 & (\mathbf{x}, t) \in \mathbf{R}^n \times (0, \infty)\\
 
	
u(\mathbf{x},0)=\delta(\mathbf{x})
 
	
\end{cases}
 
	
Η "n"- ιοστή θεμελιώδης λύση είναι το προϊόν των θεμελιωδών λύσεων για κάθε μεταβλητή'για παράδειγμα, 
 
	
:\Phi(\mathbf{x},t) = \Phi(x_1,t)\Phi(x_2,t)\dots\Phi(x_n,t)=\frac{1}{\sqrt{(4\pi k t)^n}}\exp \left (-\frac{\mathbf{x}\cdot\mathbf{x}}{4kt} \right ).

Η γενική λύση της εξίσωσης της θερμότητας για R *n* βρίσκεται με συνέλιξη,έτσι ώστε να λυθεί το πρόβλημα αρχικής τιμής με *u*(**x**, 0) = *g*(**x**), :u(\mathbf{x},t) = \int_{\mathbf{R}^n}\Phi(\mathbf{x}-\mathbf{y},t)g(\mathbf{y})d\mathbf{y}.

Το γενικό πρόβλημα σε ένα πεδίο ορισμού Ω στο **R***n* είναι 
 
	
: \begin{cases}
 
	
u_t(\mathbf{x},t) - k \sum_{i=1}^nu_{x_ix_i}(\mathbf{x},t) = 0& (\mathbf{x}, t) \in \Omega\times (0, \infty)\\
 
	
u(\mathbf{x},0)=g(\mathbf{x})&\mathbf{x}\in\Omega
 
	
\end{cases}
 
	
είτε με όριο δεδομένων τo πρόβλημα Dirichlet είτε το πρόβλημα του Neumann.Η συνάτηση του Green πάντοτε υπάρχει,αλλά αν το πεδίο ορισμού δεν μπορεί να αναλυθεί σε πρόβλημα μίας μεταβλητής(δες παρακάτω),πιθανόν να μην μπορεί γραφτεί ρητά.Άλλες μέθοδοι περιλαμβάνουν τη μέθοδο των εικόνων, το διαχωρισμό των μεταβλητών και το μετασχηματισμό Λαπλάς (Cole,2011).

## Μερικές λύσεις της συνάρτησης Green μονοδιάστατες

Μία ποικιλία από λύσεις της συνάρτησης Green σε μία διάσταση αναφέρονται εδώ άλλες υπάρχουν αλλού. Σε μερικές από αυτές το πεδίο ορισμού είναι (−∞,∞). Σε άλλες είναι το (0,∞)είτε με τις οριακές συνθήκες του Neumann ή του Ντίριχλετ. Μία άλλη ποικιλία είναι εκείνη η οποία λύνει την ανομοιογενή εξίσωση :u_{t}=ku_{xx}+f.
όπου η f είναι μία συνάρτηση του χ και t.
## Ομοιογενής εξίσωση θερμότητας

## Πρόβλημα στο (-∞,∞) αρχικών συνθηκών
 
	
:\begin{cases} 
 
	
u_{t}=ku_{xx}+f(x,t) & (x, t) \in \mathbf{R} \times (0, \infty) \\ 
 
	
u(x,0)=0 & IC 
 
	
\end{cases} 
 
	
:u(x,t)=\int_{0}^{t}\int_{-\infty}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)f(y,s)\,dy\,ds 
 
	
Σχόλιο:Αυτή η λύση είναι η συνέλιξη στο * 'R' *2,δηλαδή σε σχέση με τις μεταβλητές *χ* και *t* της θεμελιώδους λύσης
 
	
:\Phi(x,t):=\frac{1}{\sqrt{4\pi kt}}\exp\left(-\frac{x^2}{4 kt}\right)
 
	
και της συνάρτησης *f*(*x, t*),όπως ορίζεται σε όλο το  **R**2 και 0 για όλα τα *t* → 0. Επαληθεύει πως 
 
−	
:\left (\partial_t-k \partial_x^2 \right )(\Phi*f)=f,
 
	
που εκφράζονται στη γλώσσα των διανομών
 
	
:\left (\partial_t-k \partial_x^2 \right )\Phi=\delta,
 
	
όπου η κατανομή δ είναι η Dirac συνάρτηση που είναι η αξιολόγηση σε 0.

## Πρόβλημα στο (0,∞) με ομοιογενείς οριακές συνθήκες Dirichlet
	
	
:\begin{cases} u_{t}=ku_{xx}+f(x,t) & (x, t) \in [0, \infty) \times (0, \infty)  \\ u(x,0)=0 & IC \\ u(0,t)=0 & BC \end{cases} 
 
	
:u(x,t)=\int_{0}^{t}\int_{0}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \left(\exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)-\exp\left(-\frac{(x+y)^2}{4k(t-s)}\right)\right)
 
	
f(y,s)\,dy\,ds 
 
	
Σχόλιο:Αυτή η λύση λαμβάνεται από τον προηγούμενο τύπο όπως εφαρμόζεται στα δεδομένα *f*(*x*, *t*) κατάλληλα επεκτεταμένα στο **R** × [0,∞),ώστε να είναι μία περιττή συνάρτηση της μεταβλητής *χ*,δηλαδή επιτρέποντας *f*(−*x*, *t*) := −*f*(*x*, *t*) για όλα τα *χ* και *t*.Αντίστοιχα,η λύση του ανομοιογενούς προβλήματος στο (−∞,∞) είναι μία περιττή συνάρτηση σε σχέση με τη μεταβλητή *χ*για όλες τις τιμές του *t* και ειδικότερα πληροί τις ομοιογενείς οριακές συνθήκες Dirichlet *u*(0, *t*) = 0.

## Πρόβλημα στο (0,∞) με ομοιογενείς οριακές συνθήκες του Neumann
 
	
:\begin{cases} u_{t}=ku_{xx}+f(x,t) & (x, t) \in [0, \infty) \times (0, \infty)  \\ u(x,0)=0 & IC \\ u_x(0,t)=0 & BC \end{cases} 
 
	
:u(x,t)=\int_{0}^{t}\int_{0}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \left(\exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)+\exp\left(-\frac{(x+y)^2}{4k(t-s)}\right)\right)
 
	
f(y,s)\,dy\,ds 
 
	
*Σχόλιο*.Αυτή η λύση λαμβάνεται από τον πρώτο τύπο όπως εφαρμόζεται στα δεδομένα *f*(*x*, *t*) κατάλληλα επεκτεταμένο στο **R** × [0,∞),ώστε να είναι μία ακόμη συνάρτηση της μεταβλητής *χ*,δηλαδή αφήνοντας *f*(−*x*, *t*) := *f*(*x*, *t*) για όλα τα *χ* και *t*.Αντίστοιχα,η λύση του ανομοιογενούς προβλήματος στο (−∞,∞) είναι μία συνάρτηση σε σχέση με τη μεταβλητή *χ*για όλες τις τιμές του *t* και ειδικότερα είναι μία ομαλή συνάρτηση που πληροί τις ομοιογενείς οριακές συνθήκες του Neumann *u<sub>x</sub>*(0, *t*) = 0.

## Πρόβλημα στο (0,∞) ομοιογενών και μη ομοιγενών αρχικών συνθηκών του Dirichlet
:\begin{cases} u_{t}=ku_{xx} & (x, t) \in [0, \infty) \times (0, \infty) \\ u(x,0)=0 & IC \\ u(0,t)=h(t) & BC \end{cases} 
:u(x,t)=\int_{0}^{t} \frac{x}{\sqrt{4\pi k(t-s)^3}} \exp\left(-\frac{x^2}{4k(t-s)}\right)h(s)\,ds, \qquad\forall x>0

Σχόλιο.Αυτή η συνάρτηση είναι συνέλιξη με αναφορά στη μεταβλητή t 
:\psi(x,t):=-2k \partial_x \Phi(x,t)=\frac{x}{\sqrt{4\pi kt^3}}\exp\left(-\frac{x^2}{4kt}\right)
και στη συνάρτηση *h*(*t*).Από τότε που η Φ(*x*, *t*) είναι η θεμελιώδης λύση της 
:\partial_t-k\partial^2_x,
η συνάρτηση ψ(*x, t*) είναι επίσης μία λύση της ίδιας εξίσωσης θερμότητας, και τελικά είναι *u* := ψ ∗ *h*, 
με τις γενικές ιδιότητες της συνέλιξης.Περισσότερο, 
:\psi(x,t)=\frac{1}{x^2}\,\psi\left(1,\frac{t}{x^2}\right)
:\int_0^{\infty}\psi(x,t)\,dt=1,

έτσι ώστε από τη γενική αλήθεια σχετικά με προσέγγιση στην ταυτότητα
ψ(*x*, ⋅) ∗ *h* → *h* as *x* → 0 για τη συγκεκριμένη *h*.
## Ανομοιογενής εξίσωση της θερμότητας

## Το πρόβλημα στο (-∞,∞) με ομοιογενείς αρχικές συνθήκες

:\begin{cases} 
u_{t}=ku_{xx}+f(x,t) & (x, t) \in \mathbf{R} \times (0, \infty) \\ 
u(x,0)=0 & IC 
\end{cases} 
:u(x,t)=\int_{0}^{t}\int_{-\infty}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)f(y,s)\,dy\,ds 
΄΄΄Σχόλιο*.Αυτή η εξίσωση είναι η συνέλιξη *'R**<sup>2</sup>,σε σχέση με δύο μεταβλητές ''x'' και ''t'', της θεμελιώδης λύσης
:<math>\Phi(x,t):=\frac{1}{\sqrt{4\pi kt}}\exp\left(-\frac{x^2}{4 kt}\right)</math>
και της συνάρτησης ''f''(''x, t''),όπως ορίζονται σε όλο το **R**<sup>2</sup> και για όλα τα t που τείνουν στο 0 δίνουν την τιμή 0.
:<math>\left (\partial_t-k \partial_x^2 \right )(\Phi*f)=f,</math>
που στη γλώσσα των κατανομών μετατρέπεται σε 
:<math>\left (\partial_t-k \partial_x^2 \right )\Phi=\delta,</math>
όπου δ είναι η [[συνάρτηση Dirac]],η οποία εκτιμάται στο 0.

==Πρόβλημα σε (0,∞) με οριακές αρχικές συνθήκες Dirichlet και αρχικές συνθήκες==

:<math>\begin{cases} u_{t}=ku_{xx}+f(x,t) & (x, t) \in [0, \infty) \times (0, \infty)  \\ u(x,0)=0 & IC \\ u(0,t)=0 & BC \end{cases} </math>
:<math>u(x,t)=\int_{0}^{t}\int_{0}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \left(\exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)-\exp\left(-\frac{(x+y)^2}{4k(t-s)}\right)\right)
f(y,s)\,dy\,ds </math>

''Σχόλιο''.Αυτή η λύση παρατηρείται από τον τύπο όπως εφαρμόζεται για δεδομένα ''f''(''x'', ''t'') κατάλληλα επεκτεταμένα στο **R** × [0,∞),ώστε γίνεται περιττή συνάρτηση της μεταβλητής ''x'',θέτοντας ''f''(−''x'', ''t'') := −''f''(''x'', ''t'') για όλα τα ''x'' και ''t''.Η λύση του ανομοιογενούς προβλήματος στο (−∞,∞) είναι μία περιττή συνάρτηση για τη μεταβλητή  ''x'' και όλες τις τιμές του ''t'',και ικανοποιεί τις ομοιογενείς οριακές συνθήκες Dirichlet ''u''(0, ''t'') = 0.

==Πρόβλημα στο (0,∞) με ομοιογενείς οριακές συνθήκες Neumann και αρχικές συνθήκες==

:<math>\begin{cases} u_{t}=ku_{xx}+f(x,t) & (x, t) \in [0, \infty) \times (0, \infty)  \\ u(x,0)=0 & IC \\ u_x(0,t)=0 & BC \end{cases} </math>
:<math>u(x,t)=\int_{0}^{t}\int_{0}^{\infty} \frac{1}{\sqrt{4\pi k(t-s)}} \left(\exp\left(-\frac{(x-y)^2}{4k(t-s)}\right)+\exp\left(-\frac{(x+y)^2}{4k(t-s)}\right)\right)
f(y,s)\,dy\,ds </math>

''Σχόλιο''.Η λύση παρατηρείται από τον πρώτο τύπο όπως εφαρμόζεται σε δεδομένα ''f''(''x'', ''t'') καταλλήλως επεκτεταμένα σε **R** × [0,∞),ώστε να γίνει μία συνάρτηση του ''x'',θέτοντας ''f''(−''x'', ''t'') := ''f''(''x'', ''t'') για όλα τα ''x'' και ''t''.Η λύση από το ανομοιογενές πρόβλημα στο (−∞,∞) είναι μία συνάρτηση με τη μεταβλητή ''x'' για όλες τις τιμές του ''t'',και ικανοποιεί τις ομοιογενείς οριακές συνθήκες του Neumann ''u<sub>x</sub>''(0, ''t'') = 0.

== Παραδείγματα ==
 
	
Επειδή η εξίσωση της θερμότητας είναι γραμμική,λύσεις από άλλους συνδυασμούς οριακών συνθηκών,ανομοιογενούς όρου και αρχικών συνθηκών μπορούν να βρεθούν παίρνοντας έναν κατάλληλο [[γραμμικό συνδυασμό]]από τις παραπάνω λύσεις της συνάρτησης του Green.
 
	
Για παράδειγμα, για να λύσεις
 
	
 :<math>\begin{cases} u_{t}=ku_{xx}+f & (x, t) \in \mathbf{R} \times (0, \infty) \\ u(x,0)=g(x) & IC\end{cases} </math>

	
βάζοντας ''u'' = ''w'' + ''v'' όπου ''w'', ''v'',και ''r'' λύνουν τα προβλήματα
 
 :<math>\begin{cases} 
v_{t}=kv_{xx}+f, \, w_{t}=kw_{xx} \, & (x, t) \in \mathbf{R} \times (0, \infty) \\ 
v(x,0)=0,\, w(x,0)=g(x) \, & IC
\end{cases} </math>
	 	 	
==Ιδιότητα της μέσης τιμής για την εξίσωση της θερμότητας==
 
	
Λύσεις της εξίσωσης θερμότητας
 
	
:<math>(\partial_t -\Delta)u=0</math>
 
	
ικανοποιούν την ιδιότητα της μέσης τιμής ανάλογα με τις ιδιότητες των μέσων τιμών των [[αρμονικών συναρτήσεων]],λύσεις
 
	
 :<math>\Delta u=0</math>,
 
	
αν και λίγο πολύπλοκη.Ακριβώς,αν ''u'' λύνει
 
	
:<math>(\partial_t -\Delta)u=0</math>
 
	
και
 
	
:<math>(x,t)+E_\lambda\subset\mathrm{dom}(u)</math>
 
	
τότε
 
	
:<math>u(x,t)=\frac{\lambda}{4}\int_{E_\lambda}u(x-y,t-s)\frac{|y|^2}{s^2}ds\,dy,</math>
 
	
όπου ''E<sub>λ</sub>'' είναι μία ''θερμή μπάλα'',που αποτελεί υπερσύνολο της θεμελιώδους λύσης της εξίσωσης θερμότητας:
 
	
:<math>E_\lambda:=\{(y,s)\,:\,\Phi(y,s)>\lambda\},</math>
 
	
:<math>\Phi(x,t):=(4t\pi)^{-\frac{n}{2}}\exp\left(-\frac{|x|^2}{4t}\right).</math>
 
	
Σημείωσε πως 
 
	
:<math>\mathrm{diam}(E_\lambda)=o(1)</math>
 
	
καθώς το λ → ∞ ώστε ο πάνω τύπος να ισχύει για κάθε (''x, t'') στο ανοιχτό σύνολο dom(''u'') για λ αρκετά μεγάλο.Αντίθετα,κάθε συνάρτηση ικανοποιεί την παραπάνω ιδιότητα της μέσης τιμής σε μία ανοιχτή περιοχή του **R**<sup>''n''</sup> × **R*' είναι μία λύση της εξίσωσης θερμότητας.Αυτό μπορεί να δειχθεί με ένα επιχείρημα παρόμοιο με το ανάλογο ενός από [[των αρμονικών συναρτήσεων#Η ιδιότητα της μέσης τιμής|αρμονικές συναρτήσεις]].

==Σταθερή εξίσωση θερμότητας==
 
	
Η σταθερή εξίσωση θερμότητας δεν εξαρτάται από το χρόνο.Με άλλα λόγια υπάρχουν συνθήκες τέτοιες ώστε:
 
	
:<math>\frac{\partial u}{\partial t}=0</math> 	

Αυτή η συνθήκη εξαρτάται από το σταθερό χρόνο και από το χρόνο που πέρασε από τότε που επιβλήθηκαν οι οριακές συνθήκες.Έτσι,η προϋπόθεση αυτή πληρούται στις περιπτώσεις στις οποίες ο χρόνος είναι αρκετά γρήγορος ώστε η πιο περίπλοκη χρονοεξαρτώμενη εξίσωση της θερμότητας από τη σταθερή περίπτωση.Ισοδύναμα,η σταθερή κατάσταση υπάρχει για όλες τις περιπτώσεις για τις οποίες έχει περάσει αρκετός χρόνος ώστε το θερμικό πεδίο u πλέον δεν εξελίσσεται στο χρόνο.

Στην στατική περίπτωση , μια χωρική θερμική κλίση μπορεί ( ή δεν μπορεί ) να υπάρχει , αλλά αν υπάρχει , αυτό δεν αλλάζει στο χρόνο . Επομένως, η παρούσα εξίσωση περιγράφει το τελικό αποτέλεσμα σε όλα τα θερμικά προβλήματα στα οποία μία πηγή ενεργοποιείται (για παράδειγμα, ένας κινητήρας που ξεκίνησε σε ένα αυτοκίνητο ) , και αρκετός χρόνος έχει περάσει για όλες τις μόνιμες βαθμίδες θερμοκρασίας ώστε να εγκατασταθούν στο χώρο, μετά δε αυτά τα χωρικές βαθμίδες δεν αλλάζουν πλέον στο χρόνο ( όπως και πάλι , με ένα αυτοκίνητο στο οποίο ο κινητήρας έχει λειτουργήσει για αρκετά μεγάλο χρονικό διάστημα ) . Η άλλη ( τετριμμένη ) λύση είναι για όλες τους βαθμούς θερμοκρασίας για να εξαφανιστούν , καθώς, στην περίπτωση αυτή η θερμοκρασία γίνει ομοιόμορφh στο χώρο.
Η εξίσωση είναι πιο απλή και μπορεί να βοηθήσει να κατανοήσουμε καλύτερα τη φυσική των υλικών , χωρίς να εστιάζει στη δυναμική της διαδικασίας μεταφοράς θερμότητας . Χρησιμοποιείται ευρέως για απλά προβλήματα μηχανικής υποθέτοντας πως υπάρχει ισορροπία των πεδίων θερμοκρασίας και θερμότητα μεταφοράς , με το χρόνο.

Σταθερή συνθήκη:
:<math>\frac{\partial u}{\partial t}=0</math>

Η στατική εξίσωση θερμότητας για έναν όγκο που περιέχει μια πηγή θερμότητας ( την ανομοιογενή περίπτωση ) , είναι η [[εξίσωση του Poisson]]:
:<math>k \nabla^2 u=q</math>

όπου u είναι η [[θερμοκρασία]] , k είναι η [[θερμική αγωγιμότητα]] και q η θερμική πυκνότητα ροής της πηγής .

Στην [[ηλεκτροστατική]], αυτό είναι ισοδύναμο με την περίπτωση όπου ο χώρος υπό εξέταση περιέχει ένα [[ηλεκτρικό φορτίο]] .

Η στατική εξίσωση θερμότητας χωρίς μια πηγή θερμότητας εντός του όγκου (η ομοιογενής περίπτωση ) είναι η εξίσωση στην ηλεκτροστατική για όγκο ελεύθερου χώρου που δεν περιέχει ένα φορτίο . Περιγράφεται από την [[εξίσωση Λαπλάς]]:
:<math>\nabla^2 u=0</math>

==Θερμική διάχυση σε πολυμερή==
 
	
Μία άμεση εφαρμογή της εξίσωσης της θερμότητας,σε συνδυασμό με τη [[Fourier θεωρία]], σε σφαιρικές συντεταγμένες,είναι η πρόβλεψη των προφίλ θερμικής μεταφοράς και η μέτρηση της [[θερμικής διάχυσης]] σε [[πολυμερή]] (Unsworth και [[F. J. Duarte|Duarte]]). Αυτή η διπλή θεωρητική-πειραματική μέθοδος εφαρμόζεται σε καουτσούκ,διάφορα άλλα πολυμερή υλικά πρακτικού ενδιαφέροντος και σε υγρά.Οι συγγραφείς αυτοί προέρχονται από μία έκφραση της θερμοκρασίας στο κλεντρο της σφαίρας 
 
	
:<math>\frac{T_C - T_S}{T_0 - T_S} =2 \sum_{n = 1}^{\infty} (-1)^{n+1} \exp\left({-\frac{n^2 \pi^2 \alpha t}{L^2}}\right)</math>
 
	
όπου  είναι η αρχική θερμοκρασία της σφαίρας και  η θερμοκρασία στην επιφάνεια της σφαίρας,ακτίνας  .Αυτή η εξίσωση έχει επίσης βρει εφαρμογές στη μεταφορά πρωτεϊνης και στα θερμικά μοντέλα στη φυσική.

==Εφαρμογές==
==Διάχυση σωματιδίων==
		
*Κύριο άρθρο:Εξίσωση διάχυσης
*	
Κάποιος μπορεί να μοντελοποιήσει [[διάχυση]] με μία εξίσωση περιλαμβάνοντας:

*την ογκομετρική [[συγκέντρωση]] των σωματιδίων, συμβολίζοντας *c*,στην περίπτωση της[[συλλογικής διάχυσης]] ενόα μεγάλου αρθμού σωματιδίων, ή 
 
	
*τη [[συνάρτηση πυκνότητας πιθανότητας]] που σχετίζεται με τη θέση ενός μεμονωμένου σωματιδίου,που συμβολίζεται με *p*.
 
	
Σε κάθε περίπτωση,χρησιμοποιείς την εξίσωση της θερμότητας 
 
	
:<math>c_t = D \Delta c, \quad </math>
 
	
ή
 
	
:<math>P_t = D \Delta P. \quad </math>
 
	
Οι *c* και *P*είναι συναρτήσεις της θέσης και του χρόνου.*D* είναι ο συντελεστής διάχυσης που ελέγχει την ταχύτητα της διαδικασίας της διάχυσης,και τυπικά εκφράζεται σε μέτρο στο τετράγωνο ανά δευτερόλεπτο.Αν ο συντελεστής διάχυσης *D* δεν είναι σταθερός,αλλά εξαρτάται από τη συγκέντρωση *c*(ή *P*στη δεύτερη περίπτωση),τότε παίρνει την [[εξίσωση διάχυσης|μη γραμμική εξίσωση διάχυσης]].

===Κίνηση Brown===
 
	
Βάζοντας τη [[σταχαστική διαδικασία]] <math>X</math> να είναι η λύση της [[διαφορετικής στοχαστικής εξίσωσης]]
 
	
: <math>\left\{\begin{array}{ll}\text{d}X_t=\sqrt{2k}\text{d}B_t\\ X_0=0 \end{array}\right.</math>
 
	
όπου <math>B</math> είναι [[Wiener διαδικασία]] (καθορισμένη Brownian κίνηση). Τότε η [[συνάρτηση πυυκνότητας πιθανότητας]]  <math>X</math> δίνεται κάθε φορά <math>t</math> από
 
	
: <math>\frac{1}{\sqrt{4\pi kt}}\exp\left(-\frac{x^2}{4kt}\right)</math>
 
	
που είναι η λύση του προβλήματος αρχικής τιμής
 
	
: <math>\left\{\begin{array}{ll}u_t(x,t)-ku_{xx}(x,t)=0,\quad (x,t)\in\mathbb R\times(0,+\infty)\\ u(x,0)=\delta(x) \end{array}\right.</math>
 
	
όπου <math>\delta</math> είναι η [[Dirac delta συνάρτηση]].	

===Η εξίσωση Schrödinger για ένα ελεύθερο σωματίδιο===
 
	
Με μία απλή διαίρεση,η [[Schrödinger εξίσωση]] για ένα μονό σωματίδιο από [[μάζα]] m σε απουσία κάθε εφαρμοσμένου πεδίου δύναμης μπορεί να ξαναγραφτεί με τον ακόλουθο τρόπο:
 
	
:<math>\psi_t = \frac{i \hbar}{2m} \Delta \psi</math>,
 
		 
	
όπου *i* είναι η [[φανταστική μονάδα]], *ħ* είναι [[reduced Planck's σταθερά]], και *ψ* είναι η  [[κυματική συνάρτηση]] του σωματιδίου.
 
	
Αυτή η εξίσωση είναι τυπικά όμοια με την εξίσωση διάχυσης των σωματιδίων που λαμβάνεται μέσω του ακόλουθου μετασχηματισμού:
 
	
:<math>\begin{align}
 
	
c(\mathbf R,t) &\to \psi(\mathbf R,t) \\
 
	
D &\to \frac{i \hbar}{2m}
 
	
\end{align}</math>
 
	
Εφαρμόζοντας αυτό το μετασχηματισμό στις συναρτήσεις του Green ορισμένες στην περίπτωση των σωματιδίων διάχυσης δίνει την [[Schrödinger εξίσωση]],η οποία με τη σειρά της χρησιμοποιείται για την [[κυματική συνάρτηση]] μέσω ενός ολοκληρώματος της [[κυματικής συνάρτησης]] για 
 
	
*t'' = 0:
 
	
: \psi(\mathbf R, t) = \int \psi(\mathbf R^0,t=0) G(\mathbf R - \mathbf R^0,t) dR_x^0\,dR_y^0\,dR_z^0,
 
	
με
 
	
:G(\mathbf R,t) = \bigg( \frac{m}{2 \pi i \hbar t} \bigg)^{3/2} e^{-\frac {\mathbf R^2 m}{2 i \hbar t}}.
 
	
Παρατήρηση:αυτή η αναλογία μεταξύ της κβαντομηχανικής και της διάχυσης είναι καθαρά τυπική.Φυσικά,η εξέλιξη της κυματικής συνάρτησης ικανοποιώντας την Schrödinger εξίσωση θα μπορούσε να έχει προέλευση διάφορη της διάχυσης.

## Θερμική διάχυση σε πολυμερή

Μία άμεση εφαρμογή της εξίσωσης της θερμότητας,σε συνδυασμό με τη Fourier θεωρία, σε σφαιρικές συντεταγμένες,είναι η πρόβλεψη των προφίλ θερμικής μεταφοράς και η μέτρηση της θερμικής διάχυσης σε πολυμερή (Unsworth και Duarte). Αυτή η διπλή θεωρητική-πειραματική μέθοδος εφαρμόζεται σε καουτσούκ,διάφορα άλλα πολυμερή υλικά πρακτικού ενδιαφέροντος και σε υγρά.Οι συγγραφείς αυτοί προέρχονται από μία έκφραση της θερμοκρασίας στο κλεντρο της σφαίρας 
 
−	
:\frac{T_C - T_S}{T_0 - T_S} =2 \sum_{n = 1}^{\infty} (-1)^{n+1} \exp\left({-\frac{n^2 \pi^2 \alpha t}{L^2}}\right)
 
−	
όπου  είναι η αρχική θερμοκρασία της σφαίρας και  η θερμοκρασία στην επιφάνεια της σφαίρας,ακτίνας  .Αυτή η εξίσωση έχει επίσης βρει εφαρμογές στη μεταφορά πρωτεϊνης και στα θερμικά μοντέλα στη φυσική.

## Κι άλλες εφαρμογές
 
	
Η εξίσωση της θερμότητας προκύπτει στο modeling από ένα αριθμό φαινομένων και συχνά χρησιμοποιείται στα οικονομικά μαθηματικά στη μοντελοποίηση options.Η διάσημη Black–Scholes διαφορική εξίσωση μπορεί να μετατραπεί σε εξίσωση θερμότητας που επιτρέπει σχετικά εύκολα λύσεις από ένα οικείο όργανο μαθηματικών.Πολλές από τις επεκτάσεις των απλών μοντέλων επιλογής δεν έχουν λύσεις και ως εκ τούτου πρέπει να λυθούν αριθμητικά για να αποκτήσουν μία μοντελοποιημένη τιμή.Η εξίσωση που περιγράφει την πίεση της διάχυσης σε ένα πορώδες μέσο είναι πανομοιότυπη στη μορφή με την αξίσωση θερμότητας.Προβλήματα διάχυσης που ασχολούνται με Dirichlet, Neumann και Robin οριακές συνθήκες έχουν κλειστή μορφή αναλυτικών λύσεων .
 
	
Η εξίσωση της θερμότητας χρησιμοποιείται ευρέως στην ανάλυση εικόνας   και στη μηχανική μάθηση όπως τη θεωρία οδήγησης πίσω από τις scale-space ή graph Laplacians μεθόδους.Η εξίσωση της θερμότητας μπορεί να λυθεί αποτελεσματικά αριθμητικώς χρησιμοποιώντας τη  Crank–Nicolson method του .Αυτή η μέθοδος μπορεί να επεκταθεί σε πολλά μοντέλα χωρίς κλειστού τύπου λύσεις,βλέπε για παράδειγμα .
 
	
Μία αφηρημένη μορφή της εξίσωσης της θερμότητας manifolds παρέχει μία σημαντική προσέγγιση για το Atiyah–Singer index theorem,και έχει οδηγήσει σε περεταίρω εργασίες για την εξίσωση της θερμότητας στη Riemannian geometry.

## Παραπομπές


Κατηγορία:Διάχυση
Κατηγορία:Μεταφορά θερμότητας