Obtaining exact analytical decisions of tasks heat conductions with variables in time boundary conditions

Full Text

В аналитической теории теплопроводности известны методы, основанные на определении глубины термического (прогретого) слоя, которые объединены под общим названием «интегральные методы теплового баланса» [1-12]. Их использование оказалось эффективным при решении ряда краевых задач, не допускающих получение точных аналитических решений классическими методами. Однако основной проблемой является низкая точность. Согласно этим методам процесс теплопроводности разделяется на две стадии по времени, первая из которых характеризуется постепенным продвижением фронта температурного возмущения от поверхности тела к его центру. Дополнительной искомой функцией здесь является зависящая от времени глубина прогретого слоя (величина фронта температурного возмущения). Во второй стадии изменение температуры происходит по всей толщине тела, и в рассмотрение вводится дополнительная искомая функция, характеризующая изменение температуры во времени в центре пластины. Особенностью данной группы методов является противоречие, связанное с допущением о конечной скорости распространения теплоты, несмотря на то, что решению подлежит параболическое уравнение теплопроводности, характеризующее бесконечную ее скорость. Однако результаты многочисленных исследований, приведенные в работах [1, 10-12], показали, что с увеличением числа приближений () первой стадии процесса время () достижения фронтом температурного возмущения центра пластины уменьшается и в пределе при . Следовательно, с увеличением числа приближений получаемое решение приближается к описанию процесса, протекающего с бесконечной скоростью распространения теплоты, и таким образом снимается указанное выше противоречие. Следует отметить, что с увеличением числа приближений диапазон времени первой стадии процесса уменьшается, а второй - возрастает. Следовательно, роль второй стадии в определении температурного состояния возрастает. В связи с этим в настоящей работе рассматривается метод получения точного аналитического решения, позволяющий избежать использования первой стадии процесса в интегральном методе теплового баланса. Основные его положения рассмотрим на примере решения нестационарной задачи теплопроводности для бесконечной пластины с переменными во времени граничными условиями первого рода (температура стенки - линейная функция времени) в следующей математической постановке: ; (; ); (1) ; (2) ; (3) , (4) где - температура, ; - координата, ; - время, ; - начальная температура, ; - коэффициент температуропроводности, ; - скорость нагревания стенки, ; - половина толщины пластины, . Введем следующие безразмерные переменные: ; ; ; , (5) где - безразмерная температура; - число Фурье; - безразмерная координата; - безразмерный комплекс. Задача (1)-(4) с учетом (5) примет вид ; (; ); (6) ; (7) ; (8) , (9) где - безразмерная скорость нагревания стенки. С целью упрощения процесса получения аналитического решения введем дополнительную искомую функцию , (10) представляющую изменение температуры во времени в центре пластины. Ввиду бесконечной скорости распространения теплоты, описываемой уравнением (6), температура в центре пластины начинает изменяться сразу после приложения граничного условия первого рода на ее поверхности. В связи с этим диапазон временнóго изменения функции охватывает весь диапазон времени нестационарного процесса (). Решение задачи (6)-(9) принимается в виде , (11) где , - неизвестные коэффициенты; , - координатные функции. Очевидно, что благодаря принятой системе координатных функций соотношение (11) в любом приближении точно удовлетворяет граничным условиям (8), (9). Неизвестные коэффициенты будем находить из условия (10) и некоторых дополнительных граничных условий, определяемых так, чтобы их выполнение соотношением (11) было эквивалентно выполнению уравнения (7) в граничных точках и . Отметим, что методы решения краевых задач путем выполнения уравнения на границах области рассматривались в работах [13-15]. Для получения дополнительных граничных условий, выполняемых в точках и , продифференцируем условия (8), (9) по переменной : ; (12) . (13) Соотношения (12) и (13) с учетом уравнения (6) приводятся к следующим дополнительным граничным условиям: ; (14) . (15) Продифференцируем соотношения (14), (15) по переменной : ; (16) . (17) Подставляя правую часть уравнения (6) в соотношения (16), (17), получаем дополнительные граничные условия вида ; (18) . (19) Аналогично получаются и последующие дополнительные граничные условия: ; (20) . (21) Общие формулы для них будут ; (22) . (23) Соотношение (15) не входит в общую формулу (23), так как с учетом оно, по сути, представляет математическую запись уравнения (6), которое в данной точке будет выполнено полученным решением в любом приближении. Решение (11) благодаря принятой системе координатных функций в любом приближении точно удовлетворяет дополнительным граничным условиям (22), (23). Отметим, что для краевых задач с более сложными дифференциальными операторами (нелинейных, с переменными коэффициентами и др.) дополнительные граничные условия могут иметь более сложный вид, и в случае, когда они не будут заранее выполняться искомым решением, их необходимо выполнять за счет определения неизвестных коэффициентов , . В данном случае неизвестные коэффициенты , решения (11) будем находить из соотношения (10) и получаемых на его основе дополнительных граничных условий. Для нахождения первого из них продифференцируем соотношение (10) по переменной : . (24) Сравнивая соотношение (24) с уравнением (6), получаем дополнительное граничное условие вида . (25) Дифференцируя (25) по переменной , находим . (26) Соотношение (26) с учетом уравнения (6) приводится к следующему дополнительному граничному условию: . (27) С учетом (25), (27) можно записать общую формулу для этих условий: . (28) Таким образом, из основных (8), (9) и дополнительных граничных условий (10), (22), (23), (28) невыполненными решением (11) остаются лишь условия (10) и (28), которые и будут использованы далее при определении неизвестных коэффициентов , () решения (11). Найдем решение задачи (6)-(9) в первом приближении. Подставляя (11) (ограничиваясь одним членом ряда) в соотношение (10), для определения неизвестного коэффициента будем иметь алгебраическое линейное уравнение, из решения которого находим . Соотношение (11) с учетом найденного значения принимает вид . (29) Потребуем, чтобы соотношение (29) удовлетворяло осредненному в пределах толщины пластины уравнению (6) - интегралу теплового баланса: . (30) Подставляя (29) в (30), после определения интегралов относительно функции будем иметь обыкновенное дифференциальное уравнение: . (31) Интегрируя уравнение (31), находим , (32) где - постоянная интегрирования. Подставляя (32) в (31), получаем . (33) Для определения постоянной интегрирования составим невязку начального условия (7) и потребуем ортогональности невязки к координатной функции : . (34) Подставляя (33) в (34), относительно постоянной интегрирования будем иметь алгебраическое линейное уравнение, из решения которого находим . Соотношение (33) с учетом найденного значения постоянной интегрирования принимает вид . (35) Соотношение (35) представляет решение задачи (6)-(9) в первом приближении. Оно точно удовлетворяет уравнению (6) и граничным условиям (8), (9) и приближенно - начальному условию (7). Результаты расчетов по формуле (35) в сравнении с точным решением [16] приведены на рисунке. Их анализ позволяет заключить, что в диапазоне расхождение не превышает 0,08 %. Распределение температуры в пластине: - первое приближение (по формуле (35)); - второе приближение (по формуле (42)); - третье приближение (по формуле (48)); - точное решение Для получения решения во втором приближении подставим (11) (ограничиваясь двумя членами ряда) в условия (10) и (28) (при ). Отсюда для определения неизвестных коэффициентов , () получаем систему двух алгебраических линейных уравнений, из решения которой находим ; , (36) где ; ; . Соотношение (11) с учетом соотношений (36) принимает вид , (37) где ; . Подставляя (37) в (30), относительно неизвестной функции будем иметь следующее обыкновенное дифференциальное уравнение: , (38) где . Интегрируя уравнение (38), находим (39) где , - константы интегрирования. Подставляя (39) в (37) и составляя невязку начального условия (7), получаем . (40) Требуя ортогональности невязки (40) к координатным функциям и , относительно постоянных и будем иметь систему двух алгебраических линейных уравнений: . (41) Ввиду ортогональности косинусов неизвестные и в системе уравнений (41) разделяются (в каждое уравнение входит лишь одно неизвестное). Из решения этих уравнений находим ; . С учетом (39) и найденных значений и соотношение (37) принимает вид . (42) Соотношение (42) представляет решение задачи (6)-(9) во втором приближении. Оно точно удовлетворяет уравнению (6) и граничным условиям (8), (9) и приближенно - начальному условию (7). Результаты расчетов по формуле (42) в сравнении с точным решением [16] приведены на рисунке. Их анализ позволяет заключить, что в диапазоне отклонение от точного решения уменьшается с 0,08 % (в первом приближении) до 0,03 % - во втором. Для получения решения в третьем приближении подставим (11) (ограничиваясь тремя членами ряда) в (10), (28) (при ). Относительно неизвестных коэффициентов , будем иметь систему трех алгебраических линейных уравнений, из решения которой находим: ; ; (43) , где ; ; ; ; ; ; ; . Подставляя (11) (с учетом найденных значений , ) в интеграл теплового баланса (30), после определения интегралов относительно неизвестной функции будем иметь следующее обыкновенное дифференциальное уравнение: (44) где . Интегрируя уравнение (44), находим , (45) где , , - постоянные интегрирования, определяемые из начального условия (7). Составляя его невязку, получаем (46) Требуя ортогональности невязки (46) к координатным функциям , , относительно постоянных , , будем иметь систему трех алгебраических линейных уравнений: (47) Ввиду ортогональности косинусов неизвестные , , в системе уравнений (47) разделяются и находятся в виде ; ; . С учетом (43), (45) и найденных значений постоянных интегрирования , , соотношение (11) в третьем приближении принимает вид . (48) Соотношение (48) точно удовлетворяет уравнению (6) и граничным условиям (8), (9). Начальное условие (7) в данном случае удовлетворяется приближенно (в третьем приближении). Анализ результатов расчетов, приведенных на рисунке, позволяет заключить, что температуры, найденные по формуле (48), в диапазоне практически совпадают с точными их значениями [16]. Из анализа формул (35), (42), (48) можно заключить, что появление нового слагаемого в каждом последующем приближении подчиняется определенной закономерности. Поэтому можно записать общую формулу решения для любого числа приближений, которая имеет вид (49) где . Исследование решений для последующих приближений показывает, что все они описываются соотношением (49), которое при полностью совпадает с классическим точным аналитическим решением [16]. Выводы 1. Путем использования дополнительной искомой функции и дополнительных граничных условий в интегральном методе теплового баланса получено точное аналитическое решение нестационарной задачи теплопроводности при переменных во времени граничных условиях первого рода. Введение дополнительной искомой функции, характеризующей изменение во времени температуры в центре пластины, основывается на свойстве параболического уравнения теплопроводности, связанном с бесконечной скоростью распространения теплоты. 2. Дополнительные граничные условия, определяемые с использованием исходного дифференциального уравнения и заданных граничных условий, находятся в таком виде, чтобы их выполнение искомым решением было эквивалентно выполнению в граничных точках дифференциального уравнения краевой задачи. Показано, что выполнение уравнения в граничных точках приводит к его выполнению и внутри рассматриваемой области. 3. Отсутствие необходимости непосредственного интегрирования исходного уравнения в частных производных по пространственной переменной, при замене его выполнением интеграла теплового баланса, позволяет свести решение краевой задачи к интегрированию обыкновенного дифференциального уравнения относительно дополнительной искомой функции. В связи с этим данный метод может быть использован для решения нелинейных краевых задач, задач теплопроводности с переменными физическими свойствами среды, с переменными начальными условиями и других задач.

About the authors

Igor Vasil'evich Kudinov

Samara State Technical University

Email: totig@yandex.ru
(Ph.D. (Techn.)), Associate Professor 244, Molodogvardeyskaya str., Samara, 443100, Russian Federation

References