Решение краевых задач для одномерного волнового уравнения методом разделения переменных (Фурье)




Эта задача интерпретировать малые поперечные колебания струны без воздействия внешних сил (свободные колебания), струна имеет конечную длину 0 ≤ x ≤ l, в каждой точке этой струны задано начальное ее отклонение равновесного положения и начальная скорость. Будем пока считать, что концы струны жестко закреплены, что соответствует краевому условию


Будем искать решение поставленной задачи
, тогда
и 
поставим эти уравнения в начальное уравнение.

для того, чтобы проинтегрировать раздел. переменные

поскольку последнее равенство выполнено для любых значений Х и Т, то найдется такое число К, что



из краевых условий
=> 
=> 
Следовательно, для нахождения функций Х(х) и Т(t) получаем две задачи

1-ый случай: К > 0, тогда

, где С1, С2, С3, С4 – произвольные постоянные.


Система однородная, определитель ≠ 0.
Полученная система на С3 и С4 будет иметь отличное от нуля решение только в том случае, если определитель = 0.
К = 0, что противоречит, что К = 0.
2-ой случай: К = 0
=> 




3-ий случай: К < 0
Положим 




– собственные значения поставленной задачи штурма – Муввиля, каждая из которых:



Пользуясь принципом суперпозиции решения однородных задач, получаем, что общее решение будет иметь вид:

Для определенных коэффициентов An и Bn воспользуемся начальными условиями
– ряд Фурье.
Функция f(x) по синусам, а An – коэффициент ряда Фурье по функциям f(x).


– это разложение функции g(x) в ряд Фурье.





Пример: Тугонатянутая гибкая струна. Книга Троицкой, стр. 23, задача 4
«Смешанная задача для волнового уравнения в прямоугольнике»



, где S – граница области, в которой ищется решение.
Пусть теперь мембрана имеет вид прямоугольника.

воздействие внешних сил не учитываем





Применим метод разделения переменных




Подставим эти выражения для производных в исходное уравнение:

разделим это равенство на
получим:

Причем последнее равенство выполнено для любых точек x, y, t => существует постоянное число λ, такое что



Причем, последнее равенство выполнено при любых x и y => найдутся числа λ1 и λ2 такие что:
λ = λ1 + λ2 и
, а 
Из граничных условий:
для любых точек x и t =>
Y(0) = 0
для любых точек y и t =>
X(l) = 0
для любой точки y, t =>
X(0) = 0
для любой точки x, t =>
Y(m) = 0
Следовательно, получаем две задачи Штурма – Люнвина для нахождения функций X(x) и Y(y).
Задачи:
1. Пусть λ1 > 0, тогда
C1 = C2, т.к.
λ1 > 0 быть не может
|
1. Аналогично левому столбику, можно показать, что λ2 > 0 быть не может
|
2. λ1 = 0 X(x) = ax + b
b = 0
a = 0
λ1 ≠ 0
решение тривиальное
| 2. Аналогично λ2 ≠ 0 |
3. λ1 < 0
Характеристическое уравнение
С1 = 0
– бесконечно много собственных значений задачи Штурмана – Люнвина
| λ2 < 0
Характеристическое уравнение
С3 = 0
– бесконечно много собственных значений
|
Каждое из найденных собственных значений порождает собственную функцию


Следовательно:

Найдем T(t)

Характеристическое уравнение:



– бесконечное число собственных функций задачи Штурмана – Люнвина, на нахождение функции T(t), тогда общее решение исходной задачи будем искать в виде:

где Ank и Bnk – определяются, исходя из начальных условий

Не трудно показать, что система функций
– ортогональная система функций на прямоуг
Bnk – коэффициент Фурье, разложения функции f(x,y) в ряд Фурье по этой системе =>

Воспользуемся заданной начальной скоростью для нахождения коэффициента Ank.

откуда получаем, аналогично предыдущему, что
– коэффициент Фурье разложения функции g(x,y) в ряд Фурье по системе функций


22 вариант (14)










|

U(t,x,y,z) – температура в точке
в момент времени t.
Воспользуемся законом Фурье для плотности потока тепла W в направлении нормали
в единицу времени.
, где
– производная функции температуры U(t,x,y,z) вдоль нормали
, k – коэффициент теплообмена.
k – может быть функцией температуры, точки, времени, т.е.
.
Рассмотрим часть тела V ограниченную поверхностью S.
Напишем уравнение баланса тепла в объеме D за малое время Δt.
– где М – точка объема D – функция плотности тепла от внешних источников, тогда Q1 – количество тепла от внешних источников.
– количество тепла, пришедшее в объем D за счет внешних источников за время Δt.
– расход тепла за счет выходящего из D потока.
– изменение количества тепла в области D за время Δt, где
С(x,y,z) – теплоемкость тела (вещества) в точке x,y,z;
ρ(x,y,z) – плотность вещества;
Ut(x,y,z,t) – изменение температуры.
Уравнение баланса тепла по закону сохранения энергии имеет вид:
Q3 = Q1 – Q2
или в интегральной форме

(*)
U(x,y,z) 

где 
Поток векторного поля через поверхность S:
– теорема Остроградского – Гаусса.
Применим теорему Остроградского – Гаусса к последнему интегралу (*).
,
тогда уравнение теплового баланса приобретает вид:


Последнее равенство выполнено для любой области D, целиком лежащей в объеме V; в виду произвольности области D, получаем уравнение теплопроводности:
уравнение теплопроводности, уравнение распространения тепла в объеме V.
ρ, C, k – функции от М(x,y,z) и времени t;
k – коэффициент теплообмена (Фурье).
Интересный случай, когда среда однородная, т.е. ρ, C, k – постоянные и не зависят ни от положения точки, ни от времени, тело – однородно и его характеристики с течением времени не меняются, тогда


, обозначим
– получим =>
уравнение теплопроводности для однородного тела.
– однородное уравнение теплопроводности.
– волновое уравнение.
Уравнения теплопроводности принадлежат к классу уравнений параболического типа.
Доказательство:



Совершенно аналогично выводится уравнение диффузии.
Воспользуемся законом Нернста для потока вещества W в направлении
.

U(M,t) – концентрация диффундированного вещества в точке М объема V в момент времени t.
, где
D – коэффициент диффузии.

f(M,t) – количество вещества, попадающего в точку М от внешних источников.
- диффунд. из тела, выходящий поток в-ва
где
d – коэффициент пористости среды, в котором происходит диффузия
Аналогично пред. можно показать, что из уравнения:

Уравнение диффузии, отлич. от уравнения теплопроводн. только тем, что коэффициенты D и d имеют другой физический смысл, если считать, что среда однородная, то уравнение приобретет вид
, обозначим 
Функция конечного вещества

Уравнение диффуз и распространенного тепла в некотором объеме одинаковы.
В 1-ом случае:
U(M,t) – конц. вещества в V
2-ом U(M,t) - темпер.
Уравнение теплопроводн. имеет множество решений, чтобы было единственным нужно задать начальные и краевые условия
1. Пусть λ1 > 0, тогда
C1 = C2, т.к.
λ1 > 0 быть не может
1. Аналогично левому столбику, можно показать, что λ2 > 0 быть не может
b = 0
a = 0
λ1 ≠ 0
решение тривиальное
– бесконечно много собственных значений задачи Штурмана – Люнвина
С3 = 0
– бесконечно много собственных значений