Тепломассообмен (физико-математические основы): Учебное пособие, страница 41

При    p n = - a h n2/ s2     производная     f¢(p n)  =  [i h n sh ( ih n)] /2 =                          

= - (h n  sin h n)/2 . 

Если  cos h n = 0,  то  sin h n = ± 1 = (-1)n+1,  f¢(p n) = - h n (-1)n+1/2.  Полином

Ф(рn) = (t н – t п) ch (i h n x/s) = (t н – t п) cos (h nX),   X = x/s.

Таким  образом,  правая  часть  уравнения  температурного  поля  в  оригинале  будет  иметь  вид

(t н – t п ) –  (t н – t п)

Левая  же  часть  после  обратного  преобразования  получится  в виде

L-1[ (t н /p) – T(x,p)] = t н  – t (x, t).

Таким  образом,  в  оригинале  уравнение  температурного  поля  будет  иметь  вид

Как  видно  из  этого  примера,  для  решения  задач,  кроме  таблиц,  нужно  знать  некоторые  свойства  преобразований.  Начинающим   лучше  всего  познакомиться  с  этими  свойствами  и  поучиться  пользоваться  ими  на  примерах  в  [2] ,  где  приведены  решения  уравнения  теплопроводности  при  различных  граничных  условиях  как  методами  разделения  переменных,  так  и  методами  интегральных  преобразований.

Кроме  точных  методов  решения  дифференциальных  уравнений  в  частных  производных  имеется  ряд  приближенных  аналитических  методов,  на  которых  мы  останавливаться  не  будем   в  связи   с  ограниченным  временем,  тем  более,  что  изложение  этих  методов  даже  более  громоздко,  чем  объяснение  точных,  а  распространения  широкого  они  не  получили.

3.8.   Численные   методы

С  появлением   быстродействующей  вычислительной  техники,  обладающей  большим  объёмом  памяти,  для  исследования  и  расчётов  процессов  тепло-  и  массопереноса  всё  более  широко  применяются  численные  методы:  метод  сеток (конечных  разностей),  конечных  элементов,  зональные.  Их  названия  говорят  о  том,  что  они  имеют  одну  и  ту  же  основу – малые  изменения  параметров  процессов  во  времени  и  в  пространстве.  Отличия  их  определяются  областями  применения  и  особенностями  таких  процессов  как  аэродинамика,  теплообмен  излучением,  упругие  деформации,  особенно  в  телах  сложной  формы.

Наиболее  универсальным  и  разработанным  к  настоящему  времени  является  метод  конечных  разностей,  позволяющий  решать   как  линейные  так  и  нелинейные  дифференциальные  и  интегродифференциальные   уравнения  переноса  массы  и  энергии.  Его  использование   позволяет  отказаться  от  упрощений  математической  модели  процессов,  причём  решение  может  быть  получено  с  любой  наперёд  заданной  точностью.  Большим  достоинством   этого  метода  является  также  осуществление  расчётов  путём  многократного  повторения  одинаковых  математических  операций,  что  упрощает  разработку  машинных  программ.

Метод  конечных  элементов  предпочтителен  при  расчётах  процессов  в  системах  сложной  геометрической  формы.

Зональный  метод  близок  по  форме  к  методу  конечных  элементов,  но  применяется  для  расчётов  лучистого  теплообмена   между  телами,  в  том  числе  и  при  участии  полупрозрачных  сред.

3.8.1.   Метод  конечных  разностей

Сущность  метода  конечных  разностей  заключается  в  том,  что  в  исходном  дифференциальном  уравнении,  или  в  системе уравнений,  а  также  в  условиях  однозначности,  бесконечно  малые (неопределенные)  величины  заменяются  величинами  малыми,  но  конечными (имеющими  определенную  величину).  В  общем  виде  эта  замена  не  представляет  труда. Например,  уравнение  теплопроводности  в  дифференциальной  форме     (¶ t/¶t) = a (¶ 2t/¶x2)  изменяется  только  за  счёт  одного  символа  ¶,  который  заменяется  символом  D, то  есть уравнение  приобретает  вид (Dt/Dt) = a(D2t/Dx2).  Но  теперь  это  уравнение  стало  алгебраическим,  так  как  входящие  в  него  величины  должны  иметь  определенные  числовые  значения, и позволяющим получить приближенный  результат. Замена   дифференциальных уравнений  конечно – разностными  уравнениями  может  быть  осуществлена  несколькими  способами.