Tepelná rovnice je parciální diferenciální rovnice druhého řádu , která popisuje rozložení teploty v dané oblasti prostoru a její změnu v čase.
V prostoru s libovolným souřadnicovým systémem má rovnice tepla tvar
|
kde je kladná konstanta (číslo je tepelná difuzivita ), je Laplaceův operátor a je funkcí zdrojů tepla [1] . Požadovaná funkce nastaví teplotu v bodě se souřadnicemi v okamžiku času .
Tuto rovnici lze vysvětlit následovně. Rychlost změny teploty v čase je úměrná zakřivení rozložení teploty v prostoru (druhá derivace). Jinými slovy, čím vyšší je zakřivení teplotních „hrbů“ v těle, tím rychleji dochází v těchto místech k vyrovnání teplot.
V prostoru s kartézskými souřadnicemi má rovnice tepla tvar
|
Rovnice vedení tepla se nazývá homogenní , pokud , tzn. uvnitř systému nejsou žádné zdroje a „propady“ tepla.
Zvažte Cauchyho problém pro rovnici homogenního tepla:
kde je počáteční funkce spojitá a ohraničená celým prostorem a požadovaná funkce je spojitá a omezená pro všechny hodnoty argumentu .
Pro homogenní Cauchyho problém platí následující vlastnosti [2] :
Zvažte Cauchyho problém pro rovnici nehomogenního tepla:
V tomto případě má Poissonův integrál tvar [5] :
Pro případ jedné prostorové proměnné x (problém ohřevu nebo chlazení tyče) má rovnice tepla tvar
Pro tuto rovnici můžete nastavit a řešit různé okrajové úlohy , což je jedna z metod řešení, kterou navrhl francouzský matematik Fourier a nese jeho jméno [6]
Zvažte následující problém:
Je třeba najít funkci pro .
Požadovanou funkci reprezentujeme jako produkt
Poté dosadíme navržený tvar řešení do původní rovnice, dostaneme
Rozdělme výraz na :
Protože na levé straně rovnice máme funkci, která závisí pouze na , a na pravé - pouze na , pak, když libovolnou hodnotu stanovíme na pravou stranu, dostaneme, že pro jakoukoli hodnotu levé strany rovnice je konstantní . Stejně tak se můžete ujistit, že pravá strana je konstantní, tedy rovna určité konstantě (mínus se bere pro pohodlí). Dostaneme tedy dvě obyčejné lineární diferenciální rovnice:
Dáme pozor na okrajové podmínky původní úlohy a dosadíme do nich navržený tvar rovnice, dostaneme:
odkud ( , protože jinak bychom měli řešení a hledáme pouze netriviální řešení).
Vezmeme-li v úvahu získané okrajové podmínky, získáme Sturmův-Liouvilleův problém :
Jeho řešení je redukováno na řešení lineární diferenciální rovnice a uvažování tří případů:
Vezmeme-li v úvahu zjištěné , odvodíme obecné řešení lineární diferenciální rovnice
Měl by dostat odpověď
Nyní je vše připraveno napsat řešení původního problému:
V důsledku toho máme nekonečný počet konkrétních řešení rovnice. Všechna tato konkrétní řešení jsou lineárně nezávislá , to znamená, že lineární kombinace libovolného počtu řešení je rovna nule pouze tehdy, jsou-li všechny jejich koeficienty rovny nule. Proto je logické předpokládat, že sečtením všech partikulárních řešení od jednoty do nekonečna získáme obecné řešení původního problému.
Zbývá určit hodnotu konstanty (v závislosti na ) z počáteční podmínky
Aby bylo možné určit hodnotu , je nutné rozšířit funkci do Fourierovy řady :
Dostaneme:
Odkud pochází obecné řešení:
V průběhu matematické fyziky je dokázáno, že výsledná řada splňuje všechny podmínky tohoto problému, to znamená, že funkce je diferencovatelná (a řada konverguje rovnoměrně ), splňuje rovnici v oboru definice a je spojitá v body hranice této domény.
Rovnice nehomogenního tepla s homogenními okrajovými podmínkamiZvažte následující problém pro nehomogenní rovnici :
Nechat
Poté pomocí zřejmého vztahu přepíšeme původní rovnici jako:
Vyřešme poslední lineární nehomogenní rovnici metodou variace konstanty . Nejprve najdeme obecné řešení homogenní lineární rovnice
V obecném řešení nahradíme konstantu proměnnou a dosadíme ji do původní rovnice.
Z počáteční podmínky dostaneme:
Vezmeme-li v úvahu podmínku pro , získáme
Protože
pak je samozřejmě koeficient Fourierovy řady a rovná se
V důsledku toho je obecný vzorec:
Obecný první okrajový problémV mnoha případech je možné řešit rovnici nehomogenního tepla s nehomogenními okrajovými a počátečními podmínkami
pomocí výše popsaných metod a následujícího jednoduchého triku. Požadovanou funkci reprezentujeme jako součet:
Pojďme najít funkci :
Původní problém je tedy redukován na následující:
Poté, co najdeme funkci , najdeme požadovanou funkci podle vzorce
Matematická fyzika | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
Typy rovnic | |||||||||||
Typy rovnic | |||||||||||
Okrajové podmínky | |||||||||||
Rovnice matematické fyziky |
| ||||||||||
Metody řešení |
| ||||||||||
Studium rovnic | |||||||||||
související témata |