Лабораторная работа задача коши для уравнения теплопроводности

Задача Коши для уравнения теплопроводности

По этой ссылке вы найдёте полный курс лекций по математике:

Рассмотрим однородное уравнение теплопроводности отвечающее случаю . отсутствию источников. Задача Коши ставится так: найти функцию t), удовлетворяющую уравнению и начальному условию Задача Коши для уравнения теплопроводности Физический смысл задачи состоит в определении температуры однородного бесконечного стержня в любой момент времени по известной его температуре в момент времени . Считается, что боковая поверхность стержня теплоизолирована, так что через нее тепло из стержня не уходит.

Предположи м, что достаточно гладкие функции, убываююте при х2 +t2 +00 настолько быстро, что сущ ествуют преобразования Фурье 2) законны операции дифференцирования Тогда, применяя преобразование Фурье к обеим частям уравнения (1) и условию (2), от задачи (1)-(2) перейдем к задаче Коши для обыкновенного дифференциального уравнения (величина £ играет роль параметра). Решение задачи (5)-(6) имеет вид Ранее мы установили, что где преобразование Фурье функции .

Отсюда, полагая t = получаем Таким образом, в правой части равенства (7) стоит произведение преобразований Фурье функций Пользуясь теоремой о свертке, в силу которой равенство (7) можно представить в виде Левая часть формулы (8) есть преобразование Фурье (по аргументу х) искомой функции и(х, t) , так что формулу (8) можно переписать так: откуда, пользуясь выражением для свертки функций 4>(х) ие Л, имеем Полученная формула дает решение исходной задачи (1)-(2) и называется интегралом Пуассона.

В самом деле, можно доказать, что для любой непрерывной и ограниченной функции ipt), определяемая формулой (9), имеет производные любого порядка по х и по t при t > 0 и удовлетворяет уравнению (1) при t > 0 и Vx. Покажем, что функция удовлетворяет начальному условию . Положим Тогда так что откуда при получим так как Сформулируем следующий важный результат. Теорема 1. В классе ограниченных функций решение задачи Кош и (1)-(2) единственно и непрерывно зависит от начальной функции. Пример.

Найти решение задачи Коши Задача Коши для уравнения теплопроводности А Пользуясь формулой Пуассона (9), получаем Прообразуем интеграл в правой чести.

Имеем Сделаем замену переменного Тогда интефал в правой части последнего равенства примет вид Из формулы (и) (Здесь мы воспольэов опись тем, что получаем, что / Таким образом, решение поставленной задачи о предел и тся формулой Лелю видеть, что построен ноя функция u(x,f) удовлетворяет начальному условию (2′). Непосредственной проверкой легко убедиться в том, что ата фуниция при удовлетворяет уравнение SautWtl. Из формулы Пуассона (9) следует, что тепло расоросграня ется вдоль стержня мгновенно.

Возможно вам будут полезны данные страницы:

Действительно, пусть начальная температура ) положительна для и равна нулю вне этого отрезка. Тогда для последующего распределения температур получаем откуда видяо, что при сколь угодно малых t > 0 и сколь угодно больших |х| имеем tt(x,t) > 0. Это обьяс кяется неточностью теоретических предпосылок при выводе уравнения теплопроводности , не учитыва юших инерциальн ость движе ния молекул. Тем не менее, уравнение тепло про водности дает хорошее количественное согласование с опытом. Более точное описание процессов переноса тепла дается так называемыми уравнениями переноса. 2.1.

Фундаментальное решение уравнения теплопроводности Функция входящая в формулу Пуассона (9), называется фундаментальным решением уравнения теплопроводности. Рассматриваемая как функция аргументов х, t, она удовлетворяет уравнению щ = а2ихх, в чем можно убедиться непосредственной проверкой. Фундаментальное решение имеет важный физический смысл, связанный с понятием теплового импульса. Допустим, что начальное распределение ip(x) температур таково:

Тогда в силу (9) распределение температур и, в стержне будет иметь вид По теореме о среднем где имеем Переходя в последнем равенстве к пределу при е -* 0, получим Это означает, что функция G(x, t\ хо) представляет распределение температур в стержне в момент t > 0, если в начальный момент t = 0 в точке х = Хо имелся бесконечный пик температур (при е -* 0 функция 4>е<х) +оо), а в остальных точках стержня температура была равна нулю.

может быть приближенно реализовано следующим образом: в момент t = 0 к точке х = Хо стержня на очень короткий промежуток времени подносится узкое пламя очень высокой температуры (тепловой импульс плотности ср). Это начальное распределение температур описы вается так называемой 6 — функцией Дирака, обозначаемой символом 6(х — Хо).

Не являясь функцией в обычном смысле, б-функция определяется формально при помощи соотношений на любом интервале (а, Р), содержащем точку хо Основным свойством, определяющим б-функцию, является следующее: для всякой непрерывной функции f(x) Таким образом, фундаментальное решение G(x, t\xq) является решением уравнения теплопроводности в бесконечном стержне при начальном распределении температуры График функции G(x>t;xa) при в различных значениях t > Оимеетвид (рис. 1).

Кривые 1, 2, 3 соответствуют моментам времени • Рисунок показывает, каквыравнива-ется температура в стержне после теплового импульса. Решение Задача Коши для уравнения теплопроводности задачи теплопроводности в бесконечном стержне при начальном условии можно рассматривать как результат суперпозиции температур, возникающих в точке х в момент времени t вследствие непрерывно распределенных по стержню тепловых импульсов интенсивности у>(Л) в точке Л, приложенных в момент t = 0.

Присылайте задания в любое время дня и ночи в ➔

Официальный сайт Брильёновой Натальи Валерьевны преподавателя кафедры информатики и электроники Екатеринбургского государственного института.

Все авторские права на размещённые материалы сохранены за правообладателями этих материалов. Любое коммерческое и/или иное использование кроме предварительного ознакомления материалов сайта natalibrilenova.ru запрещено. Публикация и распространение размещённых материалов не преследует за собой коммерческой и/или любой другой выгоды.

Сайт предназначен для облегчения образовательного путешествия студентам очникам и заочникам по вопросам обучения . Наталья Брильёнова не предлагает и не оказывает товары и услуги.

Метод сеток решение задачи Коши для уравнения теплопроводности

Автор: kpm12 • Май 26, 2019 • Лабораторная работа • 2,347 Слов (10 Страниц) • 284 Просмотры

Тема: Метод сеток решение задачи Коши для уравнения теплопроводности.

Цель: Научиться находить численное решение задачи Коши для уравнения теплопроводности по явной двухслойной разностной схеме.

Задание 1. Решить в прямоугольнике [pic 1] задачу Коши методом сеток с помощью явной двухслойной разностной схемы. Шаг h по оси ох взять [pic 2] , а шаг [pic 3] по оси [pic 4] выбрать из условия устойчивости.

В отчет по этому заданию включить:

1. явную разностную схему, аппроксимирующую данную задачу Коши;
2. обоснование выбора шага [pic 5] по оси [pic 6] ;
3. алгоритм выполнения задания;
4. программу по составленному алгоритму;
5. тестовую задачу для проверки работы программы, её точное решение и решение, полученное по программе для этой задачи;
6. таблицу значений решения задачи Коши в прямоугольнике D. Таблица должна состоять из 11 строк и 11 столбцов
7. график решения задачи в D;
8. выводы относительно полученных результатов

Решение . Явную двухслойную разностную схему, аппроксимирующую данную задачу, можно записать в виде:

Выбираем шаг [pic 12] . Шаг [pic 13] выбираем из условия устойчивости явной двухслойной разностной схемы, аппроксимирующей данную задачу. В данном случае это условие: [pic 14] . Можно взять [pic 15]

Равенство (1) теперь можно переписать в виде:

Вычисления можно организовать по алгоритму, состоящему из следующих шагов.

1. Вычисляем значения [pic 17] по формуле (2) для значений [pic 18] от [pic 19] до [pic 20] , где в нашем случае [pic 21] .

2. Полагая в (3) [pic 22] , вычисляем [pic 23] для значений [pic 24] от [pic 25] до [pic 26] , затем, полагая в (3) [pic 27] , вычисляем [pic 28] для значений [pic 29] от [pic 30] до [pic 31] и т. д., наконец полагая в (3) [pic 32] вычисляем [pic 33] для значений [pic 34] от 0 до [pic 35] . Таким образом, каждому [pic 36] [pic 37] соответствуют по формуле (3) значения [pic 38] для значений [pic 39] от [pic 40] до [pic 41] .

Найденные значения [pic 42] решения [pic 43] для узлов [pic 44] располагаем в виде таблицы: