Відмінності між версіями «Початково-крайова задача для рівняння теплопровідності. Метод Фур’є»

Матеріал з Вікі ЦДУ
Перейти до: навігація, пошук
(Створена сторінка: {{Механика сплошных сред}} '''Уравнение диффузии''' или '''уравнение теплопроводности''' пред…)
 
Рядок 1: Рядок 1:
{{Механика сплошных сред}}
+
Навідміну від попереднього пункта рівняння теплопровідності однорідне, але початкова умова неоднорідна.
'''Уравнение диффузии''' или '''уравнение теплопроводности''' представляет собой частный вид [[дифференциальное уравнение|дифференциального уравнения]] в частных производных. Бывает нестационарным и стационарным.
+
  
Математически уравнение диффузии и уравнение теплопроводности не различаются, и применение того или иного названия ограничено только конкретным приложением, причем второе представляется более частным, так как можно говорить, что в этом случае речь идет о диффузии тепловой энергии.
+
<math>{\begin{cases} {U}{t}-a^2{U}_{xx}=0, (1) & \\ U(0,t)=0,(2) & \\U(l,t)=0, (3)& \\U(x,0)=\varphi(x),(4)\end{cases}}</math>
 +
U-температура
  
В смысле интерпретации при решении ''уравнения диффузии'' речь идет о нахождении зависимости концентрации вещества (или иных объектов) от пространственных координат и времени, причем задан коэффициент (в общем случае также зависящий от пространственных координат и времени), характеризующий проницаемость среды для диффузии. При решении ''уравнения теплопроводности'' речь идет о нахождении зависимости температуры среды от пространственных координат и времени, причем задана теплоемкость и теплопроводность среды (также в общем случае неоднородной).
 
  
Физически в том и другом случае предполагается отсутствие или пренебрежимость макроскопических потоков вещества. Таковы физические рамки применимости этих уравнений. Также, представляя непрерывный предел указанных задач (то есть не более, чем некоторое приближение), уравнение диффузии и теплопроводности в общем не описывают статистических флуктуаций и процессов, близких по масштабу к длине и времени свободного пробега, также весьма сильно отклоняясь от предполагаемого точного решения задачи в том, что касается корреляций на расстояниях, сравнимых (и больших) с расстояниями, проходимыми звуком (или свободными от сопротивления среды частицами при их характерных скоростях) в данной среде за рассматриваемое время.
 
  
Это в подавляющей части случаев сразу же означает и то, что уравнения диффузии и теплопроводности по области применимости далеки от тех областей, где становятся существенными квантовые эффекты или конечность скорости света, то есть в подавляющей части случаев не только по своему выводу, но и принципиально, ограничиваются областью классической ньютоновской физики.
 
  
* Ближайшим формальным, а во многом и содержательным, аналогом уравнения диффузии является [[уравнение Шрёдингера]], отличающееся от уравнения диффузии множителем мнимая единица перед производной по времени. Многие теоремы о решении уравнения Шрёдингера и даже некоторые виды формальной записи его решений прямо аналогичны соответствующим теоремам об уравнении диффузии и его решениях, однако качественно их решения различаются очень сильно.
+
Задачу будемо розв'язувати методом Фур'є(методом відокремлюваних змінних)
  
== Нестационарное уравнение ==
+
<math>{U(x,t)=X(x)T(t), (5)}</math>
''Нестационарное'' уравнение диффузии классифицируется как [[Параболические уравнения|'''параболическое''' дифференциальное уравнение]]. Оно описывает распространение растворяемого вещества вследствие [[диффузия|диффузии]] или перераспределение [[температура|температуры]] тела в результате [[теплопроводность|теплопроводности]].
+
  
=== Одномерный случай ===
+
<math>{(5)\rightarrow \;(1)}</math>
В случае одномерного диффузионного процесса с коэффициентом диффузии (теплопроводности) <math>D</math> уравнение имеет вид:
+
: <math>\frac{\partial}{\partial t}c(x,\;t)=\frac{\partial}{\partial x}D\frac{\partial}{\partial x}{c(x,\;t)}+f(x,\;t).</math>
+
При постоянном <math>D</math> приобретает вид:
+
: <math>\frac{\partial}{\partial t}c(x,\;t)=D\frac{\partial^2}{\partial x^2}{c(x,\;t)}+f(x,\;t),</math>
+
где <math>c(x,\;t)</math> — концентрация диффундирующего вещества, a <math>f(x,\;t)</math> — функция, описывающая источники вещества (тепла).
+
  
=== Трёхмерный случай ===
+
<math>{{X}{T'}=a^2{X''}{T}}</math> <math>{\mid\frac{1}{X}{T}a^2}
В трёхмерном случае уравнение приобретает вид:
+
: <math>\frac{\partial}{\partial t} c(\vec{r},\;t)=(\nabla,\;D\nabla c(\vec{r},\;t))+f(\vec{r},\;t),</math>
+
где <math>\nabla=(\partial_x,\;\partial_y,\;\partial_z)</math> — [[оператор набла]], а <math>(\;,\;)</math> — скалярное произведение. Оно также может быть записано как
+
: <math>\partial_t c=\mathbf{div}\,(D\,\mathbf{grad}\,c)+f,</math>
+
а при постоянном <math>D</math> приобретает вид:
+
: <math>\frac{\partial}{\partial t} c(\vec{r},\;t)=D\Delta c(\vec{r},\;t)+f(\vec{r},\;t),</math>
+
где <math>\Delta=\nabla^2=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2}</math> — [[оператор Лапласа]].
+
 
+
=== ''n''-мерный случай ===
+
<math>n</math>-мерный случай — прямое обобщение приведенного выше, только под оператором набла, градиентом и дивергенцией, а также под оператором Лапласа надо понимать <math>n</math>-мерные версии соответствующих операторов:
+
: <math>\nabla=(\partial_1,\;\partial_2,\;\ldots,\;\partial_n),</math>
+
: <math>\Delta=\nabla^2=\partial_1^2+\partial_2^2+\ldots+\partial_n^2.</math>
+
Это касается и двумерного случая <math>n=2</math>.
+
 
+
=== Мотивировка ===
+
 
+
==== A. ====
+
Обычно уравнение диффузии возникает из эмпирического (или как-то теоретически полученного) уравнения, утверждающего пропорциональность потока вещества (или тепловой энергии) разности концентраций (температур) областей, разделённых тонким слоем вещества заданной проницаемости, характеризуемой коэффициентом диффузии (или теплопроводности):
+
: <math>\Phi=\varkappa\frac{\partial c}{\partial x}</math> (одномерный случай),
+
: <math>\mathbf j=\varkappa\nabla c</math> (для любой размерности),
+
в сочетании с уравнением непрерывности, выражающим сохранение вещества (или энергии):
+
: <math>\frac{\partial c}{\partial t}+\frac{\partial\Phi}{\partial x}=0</math> (одномерный случай),
+
: <math>\frac{\partial c}{\partial t}+\mathrm{div}\,\mathbf j=0</math> (для любой размерности),
+
с учетом в случае уравнения теплопроводности ещё теплоёмкости (температура = плотность энергия / удельная теплоемкость).
+
* Здесь источник вещества (энергии) в правой части опущен, но он, конечно же, может быть легко туда помещён, если в задаче есть приток (отток) вещества (энергии).
+
 
+
==== B. ====
+
Кроме того, оно естественно возникает как непрерывный предел аналогичного разностного уравнения, возникающего в свою очередь при рассмотрении задачи о случайном блуждании на дискретной решётке (одномерной или <math>n</math>-мерной). (Это простейшая модель; в более сложных моделях случайных блужданий уравнение диффузии также возникает в непрерывном пределе). Простейшей интерпретацией функции <math>c</math> в этом случае служит количество (или концентрация) частиц в данной точке (или вблизи неё), причём каждая частица движется независимо от остальных без памяти (инерции) своего прошлого (в несколько более сложном случае — с ограниченной по времени памятью).
+
 
+
=== Решение ===
+
'''В одномерном случае''' [[Функция Грина|фундаментальное решение]] однородного уравнения с постоянным — не зависящем от <math>x</math> и <math>t</math> — <math>D</math> (при начальном условии, выражаемом [[дельта-функция|дельта-функцией]] <math>c_f(x,\;0)=\delta(x)</math> и граничном условии <math>c_f(\infty,\;t)=0</math>) есть
+
: <math>c_f(x,\;t)=\sqrt{\frac{1}{4\pi Dt}}\exp\left(-\frac{x^2}{4Dt}\right).</math>
+
В этом случае <math>c_f(x,\;t)</math> можно интерпретировать как плотность вероятности того, что одна частица, находившаяся в начальный момент времени в исходном пункте, через время <math>t</math> перейдёт в пункт с координатой <math>x</math>. То же самое — с точностью до множителя, равного количеству диффундирующих частиц — относится к их концентрации, при условии отсутствия или пренебрежимости взаимодействия диффундирующих частиц между собой. Тогда (при таких начальных условиях) средний квадрат удаления диффундирующих частиц (или соответствующая характеристика распределения температуры) от начальной точки
+
: <math>\langle x^2\rangle=\int\limits_{-\infty}^{+\infty}x^2 c_f(x,\;t)\,dx=2Dt.</math>
+
 
+
<!--Такому параболическому закону подчиняется, например, рост окисных плёнок на поверхностях металлов при высоких температурах. Необходимым условием при этом является достаточное количество кислорода в окружающей среде.
+
(это я убрал, так как не выглядит прямо относящимся к делу - граничные условия другие, да и слишком частный пример к тому же. Сергей Сашов)-->
+
 
+
В случае произвольного начального распределения <math>c(x,\;0)</math> общее решение уравнения диффузии представляется в интегральном виде как [[свёртка]]:
+
: <math>c(x,\;t)=\int\limits_{-\infty}^{+\infty}c(x',\;0)c_f(x-x',\;t)\,dx'=\int\limits_{-\infty}^{+\infty}c(x',\;0)\frac{1}{\sqrt{4\pi Dt}}\exp\left(-\frac{(x-x')^2}{4Dt}\right)\,dx'.</math>
+
 
+
=== Физические замечания ===
+
Так как приближение, реализуемое уравнениями диффузии и теплопроводности, принципиально ограничивается областью низких скоростей и макроскопических масштабов (см. выше), то неудивительно, что их фундаментальное решение на больших расстояниях ведёт себя не слишком реалистично, формально допуская бесконечное распространение воздействия в пространстве за конечное время; надо при этом заметить, что величина этого воздействия так быстро убывает с расстоянием, что этот эффект как правило в принципе ненаблюдаем (например, речь идёт о концентрациях много меньше единицы).
+
 
+
Впрочем, если речь идёт о ситуациях, когда могут быть экспериментально измерены столь маленькие концентрации, и это для нас существенно, нужно пользоваться по меньшей мере не дифференциальным, а разностным уравнением диффузии, а лучше — и более подробными микроскопической физической и статистической моделями, чтобы получить более адекватное представление о реальности в этих случаях.
+
 
+
== Стационарное уравнение ==
+
В случае, когда ставится задача по нахождению установившегося распределения плотности или температуры (например, в случае, когда распределение источников не зависит от времени), из нестационарного уравнения выбрасывают члены уравнения, связанные с временем. Тогда получается '''стационарное уравнение теплопроводности''', относящееся к классу [[Эллиптические уравнения|эллиптических уравнений]]. Его общий вид:
+
: <math>-(\nabla,\;D\nabla c(\vec{r}))=f(\vec{r}).</math>
+
* При <math>D</math>, не зависящем от <math>\vec{r}</math>, стационарное уравнение диффузии становится [[уравнение Пуассона|уравнением Пуассона]] (неоднородное), или [[уравнение Лапласа|уравнением Лапласа]] (однородное, то есть при <math>f=0</math>):
+
: <math>\Delta c(\vec{r})=-\frac{f(\vec{r})}{D},</math>
+
: <math>\Delta c(\vec{r})=0.</math>
+
 
+
== Постановка краевых задач ==
+
* '''Задача с начальными условиями ([[задача Коши]]) о распределении температуры на бесконечной прямой'''
+
Если рассматривать процесс теплопроводности в очень длинном стержне, то в течение небольшого промежутка времени влияние температур на границах практически отсутствует, и температура на рассматриваемом участке зависит лишь от начального распределения температур.
+
 
+
''Найти решение уравнения теплопроводности в области <math>-\infty\leqslant x\leqslant +\infty</math> и <math>t\geqslant t_0</math>, удовлетворяющее условию <math>u(x,\;t_0)=\varphi(x)\quad(-\infty<x<+\infty)</math>, где <math>\varphi(x)</math> — заданная функция.''
+
* '''Первая краевая задача для полубесконечного стержня'''
+
Если интересующий нас участок стержня находится вблизи одного конца и значительно удалён от другого, то мы приходим к краевой задаче, в которой учитывается влияние лишь одного из краевых условий.
+
 
+
''Найти решение уравнения теплопроводности в области <math>-\infty\leqslant x\leqslant +\infty</math> и <math>t\geqslant t_0</math>, удовлетворяющее условиям
+
: <math>\left\{\begin{array}{l}
+
u(x,\;t_0)=\varphi(x),\quad(0<x<\infty) \\
+
u(0,\;t)=\mu(t),\quad(t\geqslant t_0)
+
\end{array}\right.</math>
+
где <math>\varphi(x)</math> и <math>\mu(t)</math> — заданнные функции.''
+
 
+
* '''Краевая задача без начальных условий'''
+
Если момент времени который нас интересует достаточно удалён от начального, то имеет смысл принебречь начальными условиями, поскольку их влияние на процесс с течением времени ослабевает. Таким образом, мы приходим к задаче, в которой заданы краевые условия и отсутствуют начальные.
+
 
+
''Найти решение уравнения теплопроводности в области <math>0\leqslant x\leqslant l</math> и <math>-\infty<t</math>, удовлетворяющее условиям
+
: <math>\left\{\begin{array}{l}
+
u(0,\;t)=\mu _1(t), \\
+
u(l,\;t)=\mu _2(t),
+
\end{array}\right.</math>
+
где <math>\mu_1(t)</math> и <math>\mu_2(t)</math> — заданнные функции.
+
 
+
* '''Краевые задачи для ограниченного стержня'''
+
Рассмотрим следующую краевую задачу:
+
: <math>u_t=a^2 u_{xx}+f(x,\;t),\quad 0<x<l,\;0<t\leqslant T</math> — уравнение теплопроводности.
+
Если <math>f(x,\;t)=0</math>, то такое уравнение называют ''однородным'', в противном случае — ''неоднородным''.
+
: <math>u(x,\;0)=\varphi(x),\quad 0\leqslant x\leqslant l</math> — начальное условие в момент времени <math>t=0</math>, температура в точке <math>x</math> задается функцией <math>\varphi(x)</math>.
+
: <math>\left.\begin{array}{l}
+
u(0,\;t)=\mu_1(t), \\
+
u(l,\;t)=\mu_2(t),
+
\end{array}\right\}\quad 0\leqslant t\leqslant T</math> — краевые условия. Функции <math>\mu_1(t)</math> и <math>\mu_2(t)</math> задают значение температуры в граничных точка 0 и <math>l</math> в любой момент времени <math>t</math>.
+
 
+
В зависимости от рода краевых условий, задачи для уравнения теплопроводности можно разбить на три типа. Рассмотрим общий случай (<math>\alpha_i^2+\beta_i^2\ne 0,\;(i=1,\;2)</math>).
+
: <math>\begin{array}{l}
+
\alpha_1 u_x(0,\;t)+\beta_1 u(0,\;t)=\mu_1(t), \\
+
\alpha_2 u_x(l,\;t)+\beta_2 u(l,\;t)=\mu_2(t).
+
\end{array}</math>
+
 
+
Если <math>\alpha_i=0,\;(i=1,\;2)</math>, то такое условие называют ''условием первого рода'', если <math>\beta_i=0,\;(i=1,\;2)</math> — ''второго рода'', а если <math>\alpha_i</math> и <math>\beta_i</math> отличны от нуля, то условием ''третьего рода''. Отсюда получаем задачи для уравнения теплопроводности — первую, вторую и третью краевую.
+
 
+
== Способы решения уравнений теплопроводности ==
+
 
+
=== Метод разделения переменных ===
+
==== Однородное уравнение теплопроводности с однородными граничными условиями ====
+
Рассмотрим следующую задачу
+
: <math>\begin{array}{l}
+
u_t=a^2 u_{xx},\quad 0<x<l,\;0<t\leqslant T \\
+
u(x,\;0)=\varphi(x);\quad 0\leqslant x\leqslant l \\
+
\left.\begin{array}{l}
+
u(0,\;t)=0, \\
+
u(l,\;t)=0. \\
+
\end{array}\right\}\quad 0\leqslant t\leqslant T \\
+
\end{array}</math>
+
Требуется найти функцию <math>u(x,\;t)</math> для <math>\forall(x,\;t):0\leqslant x\leqslant l,\;0\leqslant t\leqslant T</math>.
+
 
+
Представим искомую функцию в виде произведения
+
: <math>u(x,\;t)=X(x)T(t).</math>
+
Затем предполагаемую форму решения подставим в исходное уравнение, получим
+
: <math>X(x)T'(t)=a^2 X''(x)T(t).</math>
+
Разделим выражение на <math>a^2 X(x)T(t)</math>:
+
: <math>\frac{1}{a^2}\frac{T'(t)}{T(t)}=\frac{X''(x)}{X(x)}=-\lambda,\;\lambda=\mathrm{const}.</math>
+
Так как в левой части уравнения у нас находится функция зависящая только от <math>t</math>, а в правой — только от <math>x</math>, то, фиксируя любое значение <math>x</math> в правой части, получаем, что для любого <math>t</math> значение левой части уравнения постоянно. Таким же образом можно убедиться, что и правая часть постоянна, то есть равна некой константе <math>-\lambda</math> (минус взят для удобства). Таким образом, мы получаем два обыкновенных линейных дифференциальных уравнения:
+
: <math>\begin{array}{l}
+
X''(x)+\lambda X(x) = 0, \\
+
T'(t)+a^2\lambda T(t)=0.
+
\end{array}</math>
+
 
+
Обратим внимание на граничные условия исходной задачи и подставим в них предполагаемый вид уравнения, получим:
+
: <math>\begin{array}{l}
+
u(0,\;t)=X(0)T(t)=0, \\
+
u(l,\;t)=X(l)T(t)=0,
+
\end{array}</math>
+
откуда <math>X(0)=X(l)=0</math> (<math>T(t)\ne 0</math>, так как в противном случае мы имели бы решение <math>u(x,\;t)=0</math>, а мы ищем только нетривиальные решения).
+
 
+
С учетом полученных граничных условий мы получаем [[Задача Штурма — Лиувилля|задачу Штурма — Лиувилля]]:
+
: <math>\begin{array}{l}
+
X''(x)+\lambda X(x)=0; \\
+
X(0)=0, \\
+
X(l)=0. \\
+
\end{array}</math>
+
 
+
Её решение сводится к решению [[Линейное дифференциальное уравнение с постоянными коэффициентами|линейного дифференциального уравнения]] и рассмотрению трёх случаев:
+
# <math>\lambda<0.</math>
+
#: В этом случае общий вид решения будет следующим:
+
#: <math>X(x)=C_1 e^{\sqrt{-\lambda}x}+C_2 e^{-\sqrt{-\lambda}x}.</math>
+
#: Подставив граничные условия, мы убедимся, что решение будет <math>X(x)\equiv 0</math>, а мы ищем только нетривиальные решения, следовательно, этот случай не подходит.
+
# <math>\lambda=0.</math>
+
#: Общий вид решения
+
#: <math>X(x)=C_1 x+C_2.</math>
+
#: Несложно убедиться, что этот вариант нам также не подходит.
+
# <math>\lambda>0.</math>
+
#: Общий вид решения
+
#: <math>X(x)=C_1\cos(\sqrt\lambda x)+C_2\sin(\sqrt\lambda x).</math>
+
#: Подставим граничные условия:
+
#: <math>\begin{array}{l}
+
X(0)=C_1=0, \\
+
X(l)=C_2\sin(\sqrt\lambda l)=0.
+
\end{array}</math>
+
#: Так как мы ищем только нетривиальные решения, <math>C_2=0</math> нам не подходит, следовательно
+
#: <math>\begin{array}{l}
+
\sin(\sqrt\lambda l)=0, \\
+
\sqrt\lambda l=\pi n,\quad n=1,\;2,\;\ldots \\
+
\end{array}</math>
+
#: <math>\lambda_n=\left(\frac{\pi n}{l}\right)^2,\quad n=1,\;2,\;\ldots</math>
+
#: Отсюда
+
#: <math> X_n(x)=C_n\sin\left(\frac{\pi n}{l}x\right),\;\quad n=1,\;2,\;\ldots</math>
+
 
+
C учетом найденных <math>\lambda</math>, выведем общее решение линейного дифференциального уравнения
+
: <math>T'(t)+a^2\left(\frac{\pi n}{l}\right)^2 T(t)=0.</math>
+
Должен получиться ответ
+
: <math>T_n(t)=D_n\exp\left(-a^2\left(\frac{\pi n}{l}\right)^2 t\right),\quad D_n=\mathrm{const}.</math>
+
 
+
Теперь всё готово для того, чтобы записать решение исходной задачи:
+
: <math>u_n(x,\;t)=X_n(x)T_n(t)=C_n\sin\left(\frac{\pi n}{l}x\right)\exp\left(-a^2\left(\frac{\pi n}{l}\right)^2 t\right),\quad n=1,\;2,\;\ldots</math>
+
В результате у нас получилось бесконечное количество частных решений уравнения. Все эти частные решения линейно независимы, то есть линейная комбинация любого количества решений равна нулю, только если все коэффициенты при них равны нулю. Поэтому логично предположить, что суммируя все частные решения по <math>n</math> от единицы до бесконечности, мы получим общее решение исходной задачи.
+
: <math>u(x,\;t)=\sum\limits_{n=1}^\infty u_n(x)=\sum\limits_{n=1}^\infty C_n\sin\left(\frac{\pi n}{l}x\right)\exp\left(-a^2\left(\frac{\pi n}{l}\right)^2 t\right).</math>
+
Осталось определить значение константы <math>C</math> (зависящей от <math>n</math>) из начального условия
+
: <math>u(x,\;0)=\varphi(x).</math>
+
Для того, чтобы определить значение <math>C_n</math>, необходимо разложить функцию <math>\varphi(x)</math> в [[ряд Фурье]]:
+
: <math>\begin{array}{l}
+
\varphi(x)=\sum\limits_{n=1}^\infty A_n\sin\left(\dfrac{\pi n}{l}x\right), \\
+
A_n=\dfrac{2}{l}\displaystyle\int\limits_0^l \varphi(\xi)\sin\left(\dfrac{\pi n}{l}\xi\right)\,d\xi.
+
\end{array}</math>
+
Получаем:
+
: <math>\begin{array}{l}
+
u(x,\;0)=\sum\limits_{n=1}^\infty C_n\sin\left(\dfrac{\pi n}{l}x\right)=\sum\limits_{n=1}^\infty A_n\sin\left(\dfrac{\pi n}{l}x\right), \\
+
C_n=A_n=\dfrac{2}{l}\displaystyle\int\limits_0^l \varphi(\xi)\sin\left(\dfrac{\pi n}{l}\xi\right)\,d\xi.
+
\end{array}</math>
+
Откуда общее решение:
+
: <math>u(x,\;t)=\sum\limits_{n=1}^\infty \left(\dfrac{2}{l}\int\limits_0^l \varphi(\xi)\sin\left(\dfrac{\pi n}{l}\xi\right)\,d\xi\right)\sin\left(\dfrac{\pi n}{l}x\right)\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right).</math>
+
В курсе [[Математическая физика|математической физики]] доказывается, что полученный ряд удовлетворяет всем условиям данной задачи, то есть функция <math>u(x,\;t)</math> дифференцируема (и ряд сходится равномерно), удовлетворяет уравнению в области определения и непрерывна в точках границы этой области.
+
 
+
==== Неоднородное уравнение теплопроводности с однородными граничными условиями ====
+
Рассмотрим способ решения неоднородного уравнения:
+
: <math>\begin{array}{l}
+
u_t=a^2 u_{xx}+f(x,\;t),\quad 0<x<l,\;0<t\leqslant T \\
+
u(x,\;0)=0;\quad 0\leqslant x\leqslant l \\
+
\left.\begin{array}{l}
+
u(0,\;t)=0, \\
+
u(l,\;t)=0. \\
+
\end{array} \right\}\quad 0\leqslant t\leqslant T
+
\end{array}</math>
+
Пусть
+
: <math>
+
\begin{array}{l}
+
u_n(x,\;t)=X_n(x)T_n(t), \\
+
f_n(x,\;t)=X_n(x)F_n(t), \\
+
X_n(x)=\sin\left(\dfrac{\pi n}{l}x\right).
+
\end{array}</math>
+
Тогда, пользуясь очевидным соотношением <math>X''_n(x)=-\left(\frac{\pi n}{l}\right)^2 X_n(x)</math>, перепишем исходное уравнение как:
+
: <math>\begin{array}{l}
+
X_n(x)T'_n(t)=-a^2\left(\dfrac{\pi n}{l}\right)^2 X_n(x)T_n(t)+X_n(x)F_n(t), \\
+
T'_n(t)=-a^2\left(\dfrac{\pi n}{l}\right)^2 T_n(t)+F_n(t).
+
\end{array}</math>
+
Решим последнее линейное неоднородное уравнение методом [[Метод Лагранжа (дифференциальные уравнения)|вариации постоянной]]. Сначала найдём общее решение однородного линейного уравнения
+
: <math>\begin{array}{l}
+
T'_n(t)=-a^2\left(\dfrac{\pi n}{l}\right)^2 T_n(t), \\
+
T_n(t)=D\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right).
+
\end{array}</math>
+
В общем решении заменим постоянную <math>D</math> на переменную <math>D(t)</math> и подставим в исходное уравнение.
+
: <math>\begin{array}{l}
+
T_n(t)=D(t)\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right), \\
+
D'_n(t)\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)-a^2\left(\dfrac{\pi n}{l}\right)^2\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)D_n(t)=-a^2\left(\dfrac{\pi n}{l}\right)^2\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)D_n(t)+F_n(t), \\
+
D'_n(t)\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)=F_n(t), \\
+
D_n(t)=\displaystyle\int F_n(t)\exp\left(a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)\,dt+A_n, \\
+
T_n(t)=A_n\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)+\exp\left(-a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)\displaystyle\int F_n(t)\exp\left(a^2\left(\dfrac{\pi n}{l}\right)^2 t\right)\,dt.
+
\end{array}</math>
+
Из начального условия получаем:
+
: <math>\begin{array}{l}
+
u_n(x,\;0)=X_n(x)T_n(0)=0, \\
+
T_n(0)=0.
+
\end{array}</math>
+
С учетом условия для <math>T</math>, получаем
+
: <math>T_n(t)=\int\limits_0^t \exp\left(-a^2\left(\frac{\pi n}{l}\right)^2 (t-\tau)\right)F_n(\tau)\,d\tau.</math>
+
Так как
+
: <math>f_n(x,\;t)=X_n(x)F_n(t)=\sin\left(\frac{\pi n}{l}x\right)F_n(t),</math>
+
то <math>F_n(t)</math>, очевидно, является коэффициентом ряда Фурье, и равен
+
: <math>F_n(t)=\frac{2}{l}\int\limits_0^l f(\xi,\;t)\sin\left(\frac{\pi n}{l}\xi\right)\,d\xi.</math>
+
В результате, общая формула такова:
+
: <math>u(x,\;t)=\sum\limits_{n=1}^\infty X_n(x)T_n(t)=\sum\limits_{n=1}^\infty \left[\int\limits_0^t \exp\left(-a^2\left(\frac{\pi n}{l}\right)^2 (t-\tau)\right)\left\{\frac{2}{l}\int\limits_0^l f(\xi,\;\tau)\sin\left(\frac{\pi n}{l}\xi\right)\,d\xi\right\}\,d\tau\right]\sin\left(\frac{\pi n}{l}x\right).</math>
+
 
+
==== Общая первая краевая задача ====
+
Во многих случаях удаётся решить уравнение теплопроводности с неоднородными краевыми и начальным условиями
+
: <math>\begin{array}{l}
+
u_t=a^2 u_{xx}+f(x,\;t), \\
+
u(x,\;0)=\varphi(x), \\
+
u(0,\;t)=\mu_1(t), \\
+
u(l,\;t)=\mu_2(t)
+
\end{array}</math>
+
с помощью методов, описанных выше и следующего несложного приёма. Представим искомую функцию в виде суммы:
+
: <math>\begin{array}{l}
+
u(x,\;t)=\tilde u(x,\;t)+U(x,\;t), \\
+
\tilde u(x,\;0)=u(x,\;0)-U(x,\;0)=\varphi(x)-U(x,\;0), \\
+
\tilde u(0,\;t)=0, \\
+
\tilde u(l,\;t)=0.
+
\end{array}</math>
+
Найдём функцию <math>U(x,\;t)</math>:
+
: <math>\begin{array}{l}
+
U(x,\;t)=Ax+b, \\
+
U(0,\;t)=b=\mu_1(t), \\
+
U(l,\;t)=Al+\mu_1=\mu_2\Rightarrow A=\dfrac{\mu _2(t)-\mu_1(t)}{l}, \\
+
U(x,\;t)=\dfrac{\mu_2(t)-\mu_1(t)}{l}x+\mu_1(t).
+
\end{array}</math>
+
Таким образом, исходная задача свелась к следующей:
+
: <math>\begin{array}{l}
+
\tilde u_t=a^2\tilde u_{xx}+f(x,\;t)-\dfrac{\mu'_2(t)-\mu'_1(t)}{l}x-\mu'_1(t), \\
+
\tilde u(x,\;0)=\varphi(x)-\dfrac{\mu_2(0)-\mu_1(0)}{l}x-\mu_1(0), \\
+
\tilde u(0,\;t)=0, \\
+
\tilde u(l,\;t)=0.
+
\end{array}</math>
+
После того, как мы найдём функцию <math>\tilde u(x,\;t)</math>, искомую функцию найдём по формуле
+
: <math>u(x,\;t)=\tilde u(x,\;t)+\frac{\mu_2-\mu_1}{l}x+\mu_1.</math>
+
 
+
== [[Принцип максимума модуля|Принцип максимума]] ==
+
Пусть функция <math>u(x,\;t)</math> в пространстве <math>D\times[0,\;T],\;D\in\R^n</math>, удовлетворяет однородному уравнению теплопроводности <math>\frac{\partial u}{\partial t}-a^2\Delta u=0</math>, причем <math>D</math> — ограниченная область. Принцип максимума утверждает, что функция <math>u(x,\;t)</math> может принимать экстремальные значения либо в начальный момент времени, либо на границе области <math>D</math>.
+
 
+
[[Категория:Дифференциальные уравнения в частных производных|Диффузии уравнение]]
+
[[Категория:Физические законы и уравнения|Диффузии уравнение]]
+
 
+
[[cs:Rovnice vedení tepla]]
+
[[de:Wärmeleitungsgleichung]]
+
[[en:Diffusion equation]]
+
[[fr:Équation de la chaleur]]
+
[[he:משוואת החום]]
+
[[it:Equazione del calore]]
+
[[nl:Warmtevergelijking]]
+
[[pl:Równanie przewodnictwa cieplnego]]
+
[[pt:Equação do calor]]
+
[[sl:Difuzijska enačba]]
+
[[sv:Värmeledningsekvationen]]
+
[[zh:熱傳導方程式]]
+

Версія за 19:11, 19 травня 2010

Навідміну від попереднього пункта рівняння теплопровідності однорідне, але початкова умова неоднорідна.

Неможливо розібрати вираз (невідома помилка): {\begin{cases} {U}{t}-a^2{U}_{xx}=0, (1) & \\ U(0,t)=0,(2) & \\U(l,t)=0, (3)& \\U(x,0)=\varphi(x),(4)\end{cases}}

U-температура



Задачу будемо розв'язувати методом Фур'є(методом відокремлюваних змінних)

Неможливо розібрати вираз (невідома помилка): {U(x,t)=X(x)T(t), (5)}


Неможливо розібрати вираз (невідома помилка): {(5)\rightarrow \;(1)}


Неможливо розібрати вираз (невідома помилка): {{X}{T'}=a^2{X''}{T}}

 Неможливо розібрати вираз (невідома помилка): {\mid\frac{1}{X}{T}a^2}