Нестационарное одномерное температурное поле в полуограниченной среде с заданной постоянной температурой на поверхности.
1.Изменение температуры в неограниченной среде. Пусть в неограниченной среде задано начальное распределение температуры T = f(x). Для того, чтобы определить, как это распределение будет изменяться со временем, можно каждую точку среды считать источником с количеством тепла, равным: q = rcTdx' = rcf(x')dx'. Тогда температуру в любой точке среды в любой момент времени можно записать в виде интеграла от фундаментального решения по x', полагая в этой формуле t' = 0:
Формула (1) дает решение поставленной задачи, если задан конкретный вид функции f(x).
Вычислить этот интеграл удобно по отдельности для областей x > 0 и x < 0. В области x > 0 сделаем замену: (x-x')2/(4at) = a2. Тогда
= В области x < 0 обозначим: x = -| x | и сделаем замену: ( -| x | -x')2/(4at) = ( | x |+ x')2/(4at) = a2. Тогда
= Итак, вычисление интеграла (2) приводит к следующему результату:
На рисунке изображен вид кривых T(x) в различные моменты времени t1 и t2 > t1 (кривые 2 и 3). Тепло постепенно перетекает из более нагретой области в более холодную. В пределе при t ® ¥, как видно из формул (3) и (4), во всей среде установится одинаковая температура T = T1 /2 (линия 4), как и должно быть по закону сохранения энергии. В точке x = 0 (в плоскости соприкосновения тел) температура равна T1 /2 в любой момент времени, как и должно быть из соображений симметрии. 2. Нестационарное одномерное температурное поле в полуограниченной среде с заданной постоянной температурой на поверхности.
Рассмотрим теперь полуограниченную среду, на поверхности которой поддерживается постоянная температура T1. Решим сначала задачу для частного случая T1 = 0, а затем обобщим полученный результат на ненулевую температуру. Пусть среда занимает область x > 0 и имеет некоторое произвольное начальное распределение температуры T = f(x). Воспользуемся предыдущим результатом, для чего продолжим (мысленно) среду в область x < 0, причем будем считать, что в этой области f (-| x' |) = -f(x') (5) (см. рисунок). Тогда температурное влияние каждой точки из области x > 0 на границу x = 0 будет компенсироваться влиянием симметричной точки из области x < 0, и граничное условие
T(0,t) = T1 = 0 будет, очевидно, удовлетворено автоматически, а распределение температуры в любой в точке x > 0 в любой момент времени можно записать в виде:
Учтем во втором интеграле условие (15.5) и одновременно поменяем пределы этого интеграла:
Эта формула определяет температуру в любой точке x полуограниченной среды в любой момент времени t, если задана начальная температура среды f(x'), а на поверхности x = 0 поддерживается нулевая температура. Если f(x') = const = T0, то эту величину можно вынести из-под знака интеграла. Тогда, выполняя вычисления аналогично тому, как это было сделано при выводе формул (3) и (4), находим:
= Формула (6) определяет температуру в любой точке x полуограниченной среды в любой момент времени t, если начальная температура среды была T0 = const, а на поверхности x = 0 поддерживается нулевая температура T1 = 0. Обобщить полученный результат на ненулевую температуру на поверхности T1 ¹ 0 проще всего переходом к новой переменной T' = T - T1. Тогда T'0 = (T0 - T1), и
отсюда
Формула (7) определяет температуру в любой точке x полуограниченной среды в любой момент времени t > 0, если начальная температура среды была T0 = const, а на поверхности x = 0 начиная с момента времени t = 0 поддерживается постоянная температура T1. Если T0 = 0, то
|