Использование дифференциальных уравнений в частных производных для моделирования реальных процессов
я уравнения (6) могут быть различные.
Условия, которые соответствуют так называемой первой краевой задаче для
[pic], следующие:
u (x, 0) = ?(x),
(7)
u (0, t) = ?1(t),
(8)
u ([pic], t) = ?2(t).
(9)
Физическое условие (7) (начальное условие) соответствует тому, что при
[pic] в разных сечениях стержня задана температура, равная ?(x). Условия
(8) и (9) (граничные условия) соответствуют тому, что на концах стержня при
х = 0 и при х = [pic] поддерживается температура, равная ?1(t) и ?2(t)
соответственно.
Доказывается, что уравнение (6) имеет единственное решение в области
[pic], удовлетворяющее условиям (7) – (9).
2.1.2. Распространение тепла в пространстве.
Рассмотрим процесс распространения тепла в трехмерном пространстве.
Пусть u (x, y, z, t) – температура в точке с координатами (x, y, z) с
момент времени t. Опытным путем установлено, что скорость прохождения тепла
через площадку [pic]s, т. е. количество тепла, протекающего за единицу
времени, определяется формулой (аналогично формуле (1))
[pic] (10)
где k – коэффициент теплопроводности рассматриваемой среды, которую мы
считаем однородной и изотропной, n – единичный вектор, направленный по
нормали к площадке [pic]s в направлении движения тепла. Таким образом,
можем записать:
[pic]
где [pic] – направляющие косинусы вектора n, или
[pic]
Подставляя выражение [pic] в формулу (10), получаем:
[pic]Q = -k n grad u [pic]s.
Количество тепла, протекающего за время ?t через площадку ?s, будет
равно:
[pic]Q[pic]t = -k n grad u [pic]t [pic]s.
Вернемся к поставленной задаче. В рассматриваемой среде выделим малый
объем V, ограниченный поверхностью S. Количество тепла, протекающего через
поверхность S, будет равно:
[pic] (11)
где n – единичный вектор, направленный по внешней нормали к поверхности S.
Очевидно, что формула (11) дает количество тепла, поступающего в объем V
(или уходящего из объема V) за время [pic]t. Количество тепла, поступившего
в объем V, идет на повышение температуры вещества этого объема.
Рассмотрим элементарный объем [pic]v. Пусть за время [pic]t его
температура поднялась на [pic]u. Очевидно, что количество тепла,
затраченное на это повышение температуры элемента [pic]v, будет равно
[pic]
где с – теплоемкость вещества, ? – плотность. Общее количество тепла,
затраченное на повышение температуры в объеме V за время [pic]t, будет
[pic]
Но это есть тепло, поступающее в объем V за время [pic]t; оно
определено формулой (11) . Таким образом, имеет место равенство
[pic]
Сокращая на [pic]t, получаем:
[pic] (12)
Поверхностный интеграл, стоящий в левой части этого равенства,
преобразуем по формуле Остроградского (в векторной форме, где F –
дивергенция векторного поля, [pic] – замкнутая поверхность)
[pic]
полагая F = k grad u:
[pic]
Заменяя двойной интеграл, стоящий в левой части равенства (12), тройным
интегралом, получим:
[pic]
Применив теорему о среднем к тройному интегралу, стоящего слева,
получим :
[pic] (14)
где P (x, y, z) – некоторая точка объема V.
Так как мы можем выделить произвольный объем V в трехмерном
пространстве, где происходит распространение тепла, и так как мы
предполагаем, что подынтегральная функция в равенстве (13) непрерывна, то
равенство (14) будет выполняться в каждой точке пространства. Итак,
[pic] (15)
Но
[pic]
Подставляя в уравнение (15), получаем:
[pic] (16)
Если k – постоянное, то
[pic]
и уравнение (15) в этом случае дает:
[pic]
или, положив [pic]
[pic] (17)
Коротко уравнение (17) записывается так:
[pic]
где [pic]u – оператор Лапласа. Уравнение (17) и есть уравнение
теплопроводности в пространстве. Для того чтобы найти единственное решение,
отвечающее поставленной задаче, нужно задать краевые условия.
Пусть имеем тело [pic], поверхность которого [pic]. В этом теле
рассматривается процесс распространения тепла. В начальный момент
температура тела задана. Это соответствует тому, что известно значение
решения при t = 0 – начальное условие:
u (x, y, z, 0) = ? (x, y, z). (18)
Кроме того, должна быть известна температура в любой точке М
поверхности [pic] тела в любой момент времени t – граничное условие:
u (М, t) = ? (М, t). (19)
(Возможны и другие граничные условия.)
Если искомая функция u (x, y, z, t) не зависит от z, что соответствует
тому, что температура не зависит от z, то получаем уравнение:
[pic] (20)
- уравнение распространения тепла на плоскости. Если
рассматривается распространения тепла в плоской области D с
границей С, то граничные условия, аналогично (18) и (19),
формулируются так:
u (x, y, 0) = ? (x, y),
u (М, t) = ? (М, t),
где ? и ? – заданные функции, М – точка границы С.
Если же функция u не зависит ни от z, ни от y, то получаем уравнение
[pic]
- уравнение распространения тепла в стержне.
§2.2. Температурные волны.
Задача о распространении температурных волн в почве является одним из
первых примеров приложения математической теории теплопроводности, развитой
Фурье, к изучению явлений природы.
Температура на поверхности земли носит, как известно, ярко выраженную
суточную и годовую периодичность. Обратимся к задаче о распространении
периодических температурных колебаний в почве, которую будем рассматривать
как однородное полупространство [pic]. Эта задача является характерной
задачей без начальных условий, так как при многократном повторении
температурного хода на поверхности влияние начальной температуры будет
меньше влияния других факторов, которыми мы пренебрегаем (например,
неоднородность почвы). Таким образом, приходим к следующей задаче:
найти ограниченное решение уравнения теплопроводности
[pic] (1)
удовлетворяющее условию
u (0, t) = A cos [pic]t. (2)
Предполагается, что функции u (x, t) и ( (t) ограничены всюду, т.е.
[pic]
Запишем граничное условие в виде
[pic] (2’)
Из линейности уравнения теплопроводности следует, что действительная и
мнимая части некоторого комплексного решения уравнения теплопроводности
каждая в отдельности удовлетворяет тому же решению.
Если найдено решение уравнения теплопроводности, удовлетворяющее
условию (2’), то его действительная часть удовлетворяет условию (2), а
мнимая – условию
[pic]
Итак, рассмотрим задачу:
[pic] (3)
Ее решение будем искать в виде
[pic] (4)
где [pic] и [pic] - неопределенные пока постоянные.
Подставляя выражение (4) в уравнение (3) и граничное условие, находим:
[pic],
откуда
[pic]
Для u (x, t) имеем:
[pic] (5)
Действительная часть этого решения
[pic] (6)
удовлетворяет уравнению теплопроводности и граничному условию (2). Формула
(6) в зависимости от выбора знака определяет не одну, а две функции. Однако
только функция, соответствующая знаку минус, удовлетворяет требованию
ограниченности. Таким образом, решение поставленной задачи получаем в виде
[pic] (7)
На основании полученного решения можно дать следующую характеристику
процесса распространения температурной волны в почве. Если температура
поверхности длительное время периодически меняется, то в почве также
устанавливаются колебания температуры с тем же периодом, причем:
1.Амплитуда колебаний экспоненционально убывает с глубиной
| |  скачать работу |
Использование дифференциальных уравнений в частных производных для моделирования реальных процессов |
