- •М.Е. Семенов, н.Н. Некрасова математическое моделирование физических процессов
- •Оглавление
- •Предисловие
- •Введение
- •Общие принципы математического моделирования
- •1.1. Математическое моделирование на основе фундаментальных законов природы
- •1.2. Вариационные принципы и математические модели
- •Контрольные вопросы и задания
- •Классификация уравнений математической физики как моделей атмосферных процессов. Уравнение волновых движений
- •2.1. Понятие об общем решении уравнений в частных производных
- •2.2. Классификация уравнений
- •2.3. Волновые уравнения
- •Продольные колебания стержня
- •Контрольные вопросы и задания
- •Методы решения волновых уравнений
- •3.1. Метод Даламбера
- •3.2. Собственные колебания
- •Контрольные вопросы и задания
- •3.3. Метод Фурье
- •Контрольные вопросы и задания
- •Уравнение теплопроводности
- •4.1. Вывод уравнения теплопроводности. Первая краевая задача
- •4.2. Решение уравнения теплопроводности методом Фурье
- •4.3. Охлаждение бесконечного стержня
- •Контрольные вопросы и задания
- •5. Хаос в дискретных моделях
- •5.1. Одномерные отображения
- •5.2. Теорема Шарковского
- •5.3. Двумерные отображения, сохраняющие площадь
- •Контрольные вопросы и задания
- •Система лоренца
- •Задача о конвекции в подогреваемом снизу слое
- •6.2. Вывод уравнений Лоренца
- •Контрольные вопросы и задания
- •Динамика системы лоренца
- •Результаты численного моделирования уравнений Лоренца
- •7.2. Ограниченность и диссипативность системы Лоренца
- •Контрольные вопросы и задания
- •Неподвижные точки. Устойчивость. Бифуркация
- •. Неподвижные точки
- •Устойчивость неподвижных точек
- •Бифуркации в системе Лоренца
- •Контрольные вопросы и задания
- •Обобщенные размерности
- •9.1. Информационные размерности
- •Величину d1 называют информационной размерностью.
- •9.2. Корреляционная размерность
- •Контрольные вопросы и задания
- •Обработка реализаций. Характеристики хаотической динамики
- •10.1. Реконструкция фазового пространства. Оценка корреляционной размерности по наблюдаемой
- •10.2. Вычисление ляпуновских показателей
- •Контрольные вопросы и задания
- •Заключение
- •Библиографический список
- •3 94006 Воронеж, ул. 20-летия Октября, 84
Контрольные вопросы и задания
3.2. Определить тип и найти общее решение уравнений:
а) ; б) ; в) .
3.3. Найти общее решение уравнения: .
3.4. Используя замены , упростить следующие уравнения (привести к каноническому виду):
а) ;
б) ;
в) .
3.5. Найти общее решение уравнения: .
Указание. Перейти к характеристическим переменным .
3.6. Используя формулу Даламбера, найти решение уравнений:
а) ;
б) ;
в) ;
г) ;
д) .
3.7. Найти в области решение уравнения ,
удовлетворяющее однородным краевым условиям:
и следующим начальным условиям:
а) ; б) ;
в) ;
г) .
3.8. Найти решение уравнения в области , удовлетворяющее краевому условию и начальным условиям:
Указание. Использовать преобразование Фурье (косинус или синус).
Уравнение теплопроводности
4.1. Вывод уравнения теплопроводности. Первая краевая задача
Рассмотрим однородный теплоизолированный с боков стержень конечной длины , имеющий постоянную по длине толщину, и настолько тонкий, чтобы в любой момент времени температуру тела во всех точках поперечного сечения можно было бы считать одинаковой.
Выберем ось (направив ее по оси стержня) так, чтобы стержень совпадал с отрезком [0; ] оси (рис. 4.1).
Рис.4.1. Однородный теплоизолированный с боков стержень
Обозначим температуру стержня в сечении x в момент t через . Тогда функция задает закон распределения температуры в стержне. Выведем дифференциальное уравнение для этой функции.
Выделим элемент стержня [x, x+Δx] и составим для него уравнение теплового баланса, согласно которому скорость изменения количества тепла в рассматриваемом объеме (изменение количества тепла в единицу времени), обусловленная теплоемкостью материала, равна количеству тепла, поступившему в этот объем в единицу времени вследствие теплопроводности.
Скорость изменения тепла в выделенном элементе стержня равна
где с – теплоемкость материала стержня, ρ – плотность, s – площадь поперечного сечения стержня.
По теореме о среднем
, где 0 < θ1 < 1.
Найдем количество тепла, поступившее в выделенный элемент стержня за единицу времени. Так как стержень теплоизолирован с боков, то тепло может поступать только через сечения, ограничивающие выделенный элемент стержня. Известно, что количество тепла, протекающее через сечение с абсциссой x за единицу времени, равно , где k – коэффициент теплопроводности, s – площадь сечения.
Поэтому искомое количество тепла равно
где 0 < θ2 < 1. (Согласно формуле конечных приращений Лагранжа к функции
).
Тогда уравнение теплового баланса будет иметь вид
.
Разделим обе части этого уравнения на sΔx (объем выделенного элемента стержня) и устремим Δx к нулю (стягивая выделенный элемент стержня к сечению). Получим уравнение теплопроводности для однородного стержня
, (4.1)
где – коэффициент температуропроводности.
Искомая функция должна удовлетворять уравнению (4.1), начальному условию , где – заданная функция от координат (это условие выражает закон распределения температуры по длине стержня в начальный момент времени) и граничным условиям
,
где и – заданные функции от времени. Они определяют температуру, поддерживаемую на конце стержня. Данное уравнение не учитывает тепловой обмен между поверхностью стержня и окружающим пространством.
Рассмотрим первую краевую задачу для уравнения теплопроводности. Пусть Ω – конечная область трехмерного пространства, S – граница области Ω, Q –цилиндр с основанием Ω в пространстве переменных (x, y, z, t), образующие которого параллельны оси 0t. QТ – часть этого цилиндра, ограниченная плоскостями t = 0 и t = Т (Т > 0), часть границы цилиндра QТ , состоящую из его нижнего основания и боковой поверхности, обозначим через Г.
Найдем в цилиндре QТ решение уравнения теплопроводности
, (4.2)
удовлетворяющее начальному условию
(4.3)
и граничному условию
, (4.4)
где Р – точка поверхности S. Функции и непрерывны, причем значения при t = 0 совпадают со значениями на границе S.
Задача нахождения решения уравнения (4.2) при условиях (4.3), (4.4) называется первой краевой задачей для уравнения теплопроводности.
Теорема. Функция , удовлетворяющая однородному уравнению теплопроводности (4.2) внутри цилиндра QТ и непрерывная вплоть до его границы, принимает наибольшее и наименьшее значение на Г, т. е. или при t = 0, или на боковой поверхности цилиндра QТ.
Из этой теоремы непосредственно вытекает, что:
решение первой краевой задачи единственно;
решение первой краевой задачи непрерывно зависит от правых частей начального и граничных условий;
если функция Т, непрерывная на замыкании QТ и удовлетворяющая первой краевой задаче, равна нулю на Г, то она тождественно равна нулю в замыкании QТ.