автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование процессов теплопроводности и гидродинамики численно-аналитическими методами на основе использования дополнительных граничных условий

005004298 На правах рукописи

Кудинов Игорь Васильевич

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ ТЕПЛОПРОВОДНОСТИ И ГИДРОДИНАМИКИ ЧИСЛЕННО-АНАЛИТИЧЕСКИМИ МЕТОДАМИ НА ОСНОВЕ ИСПОЛЬЗОВАНИЯ ДОПОЛНИТЕЛЬНЫХ ГРАНИЧНЫХ УСЛОВИЙ

Специальность: 05.13.18 Математическое моделирование, численные методы и комплексы программ

- 1 ДЕК 2011

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Самара 2011

005004298

Работа выполнена на кафедре «Прикладная математика и информатика» Федерального государственного бюджетного образовательного учреждения высшего про фессионального образования «Самарский государственный технический университет»

Научный руководитель - доктор физико - математических наук, профессор

Радченко Владимир Павлович

Официальные оппоненты: доктор физико - математических наук, профессор

Жданов Александр Иванович

доктор технических наук, профессор Клюев Николай Ильич

Ведущая организация ФГБОУ ВПО «Ульяновский государственный

технический университет»

Защита состоится «22» декабря 2011г. в 13:00 часов на заседании Диссертационного совета Д 212.217.03 при ФГБОУ ВПО «Самарский государственный технический университет» по адресу: 443100, г. Самара, ул. Галактионовская 141, корпус 6, ауд. 33.

С диссертацией можно ознакомиться в научной библиотеке Самарского государственного технического университета по адресу: г. Самара, ул. Первомайская, 18.

Отзывы на автореферат в двух экземплярах, заверенные печатью, направлять по адресу: 443100, г.Самара, ул.Молодогвардейская, 244, Главный корпус, учёному секретарю Диссертационного совета Д 212.217.03; факс (846) 278-44-00.

Автореферат разослан ноября 2011г.

Ученый секретарь диссертационного совета

Н.Г. Губанов

Общая характеристика работы

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

Существенным фактором, ограничивающим возможности математического моделирования в области тепломассопереноса, является практическое отсутствие точных аналитических методов исследования нелинейных краевых задач, в том числе гиперболического типа. Привлечение гиперболического дифференциального оператора необходимо при исследовании полей различных потенциалов при сверхмалых значениях времени (обработка материалов лазерным излучением, распределение температуры в начальной стадии теплового удара, нагрев при динамическом распространении трещины, движение космических аппаратов в атмосфере и др.).

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

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

В связи с этим возникает потребность в разработке методов, позволяющих без введения каких-либо существенных допущений получать решения в аналитическом виде, пусть даже ценой выполнения большого количества предварительных вычислительных процедур, связанных с использованием компьютерной техники. Эти методы должны сочетать в себе положительные стороны классических аналитических (аналитический вид решения) и численных (решения сложнейших краевых задач без введения упрощающих допущений) методов при возможности избегать наиболее существенные их недостатки. В этом направлении наиболее эффективными являются методы, в которых используется понятие глубины термического слоя (интегральные методы теплового баланса). Их разработке и развитию посвящены труды БиоМ., ГудменаТ., Швеца М.Е., ВейникаА.И., Беляева Н.М., Дородницына A.A.,

3

Глазунова Н.Т., Карташова Э.М. и других. Введение при использовании этих методов дополнительных граничных условий во многих случаях позволяет получать приближенные аналитические решения практически с заданной степенью точности. Развитие этих методов применительно к решению сложных линейных и нелинейных краевых задач является актуальной научной проблемой, имеющей большое теоретическое и прикладное значение.

Цель диссертационной работы состоит в математическом моделировании процессов теплопроводности и гидродинамики на основе эффективных аналитических и численно-аналитических методов решения краевых задач с использованием интегральных методов теплового баланса, ортогональных методов Канторовича Л.В. и Бубнова-Галеркина, новых алгоритмов и специального программного обеспечения.

Для достижения указанной цели решались следующие задачи.

1. Разработка методов получения точных и приближенных аналитических решений гиперболических уравнений теплопроводности и движения жидкостей, описывающих распределение соответствующих полей потенциалов с учетом конечной скорости распространения возмущении.

2. Построение аналитического решения задач динамического и теплового пограничных слоев при граничных услоьиях первого и третьего рода на стенке.

3. Разработка численно-аналитического метода решения краевой задачи Стефана с абляцией (при удалении расплавляемого вещества).

4. Разработка и исследование метода получения вихревых полей потенциалов (тепловых, гидравлических, диффузионных и прочих) путем организации движения их источников по круговым орбитам с учетом найденных закономерностей этого движения.

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

Научная новизна полученн -гс в диссертационной работе результатов заключается в следующем:

1. Разработана новая математическая модель процесса теплообмена с учётом конечной скорости распространения теплоты, позволяющая, в отличие от существующих, снять противоречия, связанные с наличием скачков температуры внутри тела и с появлением отрицательных температур в обратной тепловой волне.

2. На основе введения фронта температурного возмущения и дополнительных граничных условий получены новые численно-аналитические решения задач динамического и теплового пограничных слоев при граничных условиях первого и третьего рода на стенках.

3. На основе использования дополнительных граничных условий с учетом движения фронта температурного возмущения и фронта плавления впервые получено аналитическое решение задачи Стефана с абляцией (с удалением расплавляемого вещества).

4. Теоретически разработан метод получения вихревых полей потенциалов, возникающих при вращении по определенным закономерностям источников этих полей по круговым орбитам.

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

Практическая значимость работы. Разраббтанные в диссертации методы и полученные аналитические решения отличаются заметной простотой конструкции при точности, достаточной для прикладных задач. Такие решения полезны в тех случаях, когда решения задач теплопроводности являются промежуточными стадиями других исследований, например, решения обратных задач, задач термоупругости, автоматизированного управления и проектирования. В частности, полученные в диссертации аналитические решения задач теплообмена при теченйи жидкостей в трубах и плоских каналах были использованы при разработке математических моделей и программного комплекса для теплосети центрального теплоснабжения г. Саратова (от СарГРЭС и ТЭЦ-5). 1

Математические модели теплосети и циркуляционной системы Новокуйбышевской ТЭЦ-1 позволили определить давления, скорости, расходы теплоносителя и потери напора в любой точке моделируемых систем, найти наиболее оптимальные режимы текущей работы, выполнить предварительные проекты реконструкций, а также составить планы построения новых участков трубопроводных систем.

Применительно к расчетам температурного состояния труб барабанов котлов Новокуйбышевской ТЭЦ-1 внедрена методика расчета, основывающаяся на введении фронта температурного возмущения и дополнительных граничных условий. Для расчетов температурного состояния труб, представленных в виде многослойных конструкций (тепловая изоляция - металлическая стенка трубы - отложения накипи на внутренней поверхности труб), внедрена методика, основанная на теории обобщенных функций.

На защиту выносятся следующие основные результаты и положения.

1. Точные и приближённые численно-аналитические решения гиперболических уравнений теплопроводности, учитывающих конечную скорость распространения возмущений; исследование сходимости и погрешности численно-аналитических решений.

2. Численно-аналитические решения задач динамического и теплового пограничных слоев при граничных условиях первого и третьего рода на стенке, полученные на основе введения фронта температурного возмущения и дополнительных граничных условий; оценки сходимости и погрешности решений".

3. Численно-аналитическое решение задачи СГтефана с абляцией (с удалением расплавляемого вещества), полученное на основе дополнительных граничных условий и учёта движения фронтов температурного возмущения и плавления; оценки сходимости и погрешности решения.

4. Математическая модель вихревых полей потенциалов и численный метод их расчёта, основанный на конечно-разностном методе расщепления и методе прогонки.

5. Программный комплекс для исследования математических моделей сложных разветвленных многокольцевых трубопроводных систем и численный метод их расчёта, позволяющий рассчитывать давления, скорости, расходы и температуру движущихся сред.

6. Результаты исследований на основе программного комплекса математической модели теплосети г. Саратова (от Саратовских ТЭЦ-5 и ГРЭС), позволившие установить причины недостаточной эффективности работы теплосети, оптимизировать планы реконструкции и выполнения ремонтных работ.

Апробация работы. Основные результаты работы доложены и обсуждены на Всероссийской научной конференции молодых ученых «Наука, технологии, инновации» (г. Новосибирск, 2007 г.); Тринадцатой Международной научно-технической конференции студентов и аспирантов «Радиоэлектроника, электротехника и энергетика» (г. Москва, 2007 г.); Четырнадцатой Международной научно-практической конференции студентов, аспирантов и молодых ученых «Современные техника и технологии» (г. Томск, 2008 г.); Третьей молодежной Международной научной конференции «Тинчуринские чтения» (г. Казань, 2008 г.); Седьмой Международной конференции «Математическое моделирование физических, экономических, технических, социальных систем и процессов» (г. Ульяновск, 2009 г.); Шестой, Седьмой и Восьмой Всероссийских научных конференциях с международным участием «Математическое моделирование и краевые задачи» (г. Самара, 2009, 2010, 2011 гг.); в Школе-семинаре молодых ученых и специалистов под рук. акад. РАН А.И. Леонтьева «Проблемы газодинамики и теплообмена в аэрокосмических технологиях» (г. Москва, 2009, 2011 гг.); в Школе-семинаре молодых ученых и специалистов под рук. акад. РАН В.Е. Алемасова «Проблемы тепломассообмена и гидродинамики в энергомашиностроении» (г. Казань,

2010 г.); на научном семинаре «Прикладная математика и механика» Самарского государственного технического университета под рук. профессора Радченко В.П. в 2010 и

2011 гг.

Связь диссертационной работы с планами научных исследований. Представленная работа является обобщением теоретических и экспериментальных исследований, выполненных автором в Самарском государственном техническом университете. Исследования проводились по планам госбюджетных тематик Минвуза РФ № 1.18.09 «Получение аналитических решений задач теплопроводности на основе теории обобщенных функций и дополнительных граничных условий», № 1.21.11 «Разработка методов получения точных аналитических решений дифференциальных уравнений гиперболического типа», а также по Аналитической ведомственной целевой программе «Развитие научного потенциала высшей школы», тематический план НИР № 551/02 «Разработка нового направления получения аналитических решений задач математической физики на основе введения фронта температурного возмущения и дополнительных граничных условий».

Результаты работы использовались при выполнении энергетического аудита Самарского государственного технического университета, хоздоговорных работ с Волжской территориальной генерирующей кампанией, с Куйбышевским и Новокуйбышевским нефтеперерабатывающими заводами, ОАО «Самараоргсинтез».

Экономический эффект, подтвержденный соответствующим актом внедрения, приведенным в приложениях диссертации, составляет 3 миллиона рублей.

Обоснованность выносимых на защиту научных положений, выводов и рекомендаций, а также достоверность полученных автором результатов подтверждается соответствием математических моделей физическим процессам, протекающим в энер-

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

Публикации. Материалы диссертации опубликованы в 38 печатных работах. В автореферате приведено 26 основных научных работ, из которых 9 статей - в рецензируемых журналах из перечня ВАК, 5 статей опубликованы в центральных и академических изданиях. По результатам исследований опубликованы две монографии, одна из которых - в центральном издательстве «URSS».

Личный вклад автора. Работы [7, 10, 20, 26] выполнены самостоятельно. В основных работах [1 - 6, 8, 9, 11] диссертанту принадлежит постановка проблем исследований, непосредственное выполнение основной части работы, которая выполнена совместно с другими авторами. В остальных работах [12 - 19, 21 - 25], также опубликованных в соавторстве, диссертанту в равной степени с другими авторами принадлежат постановки задач, получение решений, содержание работы и анализ результатов.

Структура и объем диссертации. Диссертационная работа состоит из введения, шести глав, выводов, списка используемой литературы, приложений; изложена на 150 страницах основного машинописного текста и 31 странице приложений, содержит 74 рисунка, 6 таблиц. Список использованной литературы включает 93 наименования.

Краткое содержание работы

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

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

В обзоре указывается, что большой интерес представляют приближенные аналитические методы, позволяющие получать решения, хотя и приближенные, но в аналитическом виде, причем с точностью, во многих случаях достаточной для инженерных приложений. Разработке и развитию таких методов посвящены труды Л.В. Канторовича, С.Г. Михлина, П.В. Цоя, Н.М. Беляева, Ю.Т. Глазунова, Э.М. Карташова, Е.В. Стефанюк и др.

Следует однако отметить, что, несмотря на определенный прогресс, эти методы пока еще недостаточно развиты применительно к решению нелинейных задач, задач

7

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

1. Плохая обусловленность матриц коэффициентов систем алгебраических линейных уравнений, к решению которых обычно приводит использование приближенных аналитических методов (Ритца, Бубнова-Галеркина, коллокаций, Треффтца, наименьших квадратов и др.).

2. В ряде методов краевая задача сводится к решению различных интегральных уравнений (Вольтерра, Фредгольма, интегралу теплового баланса и др.), даже точные решения которых приводят к выполнению исходных дифференциальных уравнений лишь в среднем.

3. При получении приближенных аналитических решений для большого числа приближений относительно собственных чисел получаются характеристические уравнения высокой степени, решения которых, если они возможны, имеют низкую точность.

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

Математическая постановка задачи теплопроводности для бесконечной пластины при симметричных граничных условиях первого рода с учетом конечной скорости распространения теплового возмущения имеет вид:

(*»о; о^<1); (2.1) ®М)=1; (2.2)

а©(^,0)/5Ро = 0; (2.3)

50(О,Ро)/а^ = О; (2.4)

0(1,Ро)=О, (2.5)

где Рог = а тг/б2 - время релаксации; 0 = (/'-?ст)/(го -Гст) - относительная избыточная температура; !~=х/8~ безразмерное время; Ро = а т/62 - безразмерное время.

Полученное в диссертации точное аналитическое решение задачи (2.1)-(2.5) имеет вид (2.6).

©*&*>) = £ [С,; ^р(%Ро)+С2,ехр(%Ро)]со/г^1 (г = 2/-1; ¡ = к = 1,2,...), (2.6)

где

1 ± Vi-4Forv,)/(2For)(/ = 1,2; * = 1,2,...); Cu =-C2i^-; C2k=±-/

:\k rn

к Zl* i

(знак «плюс» берётся при г = 1, 5, 9,13,...; «минус»-при г = 3, 7,11,15,...).

Непосредственной подстановкой можно убедиться, что все исходные уравнения задачи (2.1) - (2.5) удовлетворяются точно. В качестве примера на рис. 2.1 представлено распределение температуры при Рог = 6,25 • 10~3.

о 02 0,4 0,6 0,8 % 1,0

Рис. 2.1. Распределение температуры в пластине с учетом конечной скорости распространения теплоты при For = 6,25-КГ3 по формуле (2.6)

Анализ результатов, представленных на рис. 2.1, позволяет сделать вывод о том, что при малых значениях времени Fo и пространственной координаты на фронте температурного возмущения наблюдается скачок температуры от величины температуры на фронте тепловой волны до начальной температуры. Область тела, находящаяся за пределами фронта тепловой волны, оказывается невозмущенной и температура здесь равна начальной температуре.

Температура на фронте тепловой волны описывается формулой г(т) = t() - ехр[- т/(2тг)], полностью совпадающей с соотношением для температуры на фронте, полученным в работах A.B. Лыкова и К. Баумейстера, где решения были найдены для полупространства.

Отмечается линейная закономерность движения фронта тепловой волны по пространственной координате во времени: хф = 8-fVr. Этот факт подтверждается исследованиями К. Баумейстера.

После достижения на фронте волны температуры 0 = 1,0 полученные по формуле (2.16) результаты полностью совпадают с решением параболического уравнения теплопроводности.

Следует отметить, что сходимость ряда (2.6) существенно зависит от величины числа Фурье, для которого необходимо определять распределение температуры по толщине пластины. Например, для задачи с приведенными выше исходными данными в диапазоне числа Фурье 0,1<Ро<со для сходимости с погрешностью в 0,01 % от точного решения достаточно всего нескольких членов этого ряда. При 0,0001<Ро<0,1 сходимость имеет место при 10 + 1000 членах ряда. Для всех чисел Фурье, при которых происходит скачок в температурных кривых, число членов ряда (2.16), необходимых для его сходимости, существенно возрастает (от 10000 при Ро = 10~5 до 1000000 при Ро = 10"8). При дальнейшем уменьшении числа Фурье количество членов ряда (2.16) может достигать нескольких миллионов.

Обнаруженный из решения (2.6) скачок температурных кривых свидетельствует о возникновении изотерм внутри тела. Кроме того, исследования решения (2.6) позволяют сделать вывод о наличии отрицательных температур в обратной тепловой волне (см. рис. 2.2). Все эти факты свидетельствуют о том, что уравнение (2.1) неадекватно уравнению теплового баланса ввиду того, что не выполняется принцип максимума (по А.Н. Тихонову).

0,54 0,5

Го - 0 5477

1 3 2 2- § -0,01

.0,3

1 0,15

0.7 а 0

0,8 _ о 85 2

'о-рЛ / 1.1 / 1.0 / —

Рис. 2.2. Изменение температуры в пластине с учетом конечной скорости распространения теплоты при Рог = 0,3 по формуле (2.6)

Для преодоления указанных трудностей в диссертации была предложена формула для теплового потока

ВТ

д2Т

дц

ц--"к--Х.т, .

дх дхск дх

(2.7)

в которой учитывается релаксация во времени не только теплового потока, но и градиента температуры.

На основе формулы (2.7) с использованием уравнения теплового баланса выведено следующее гиперболическое уравнение теплопроводности:

эо^го) „ а20(5,Ро)_э2о(^,Ро), Гд дЩъРО)

! + Ро - = " ~V»»' + р0 " / (г = 2А-1; Л = 1,2,...). (2.8)

зро г аро2 г д^дРО к '

В диссертации приводится точное аналитическое решение уравнения (2.8) при симметричных граничных условиях первого рода, которое имеет вид

к=1

Р°)=Е1[С1*еХр(2иР°)+С2* ехр(г24р°)]С05[Г^

(2.9)

где :,к =

2Рог

(/ = 1,2; Л = 1,2,...); Сп = -С2к

(знак «плюс» берётся при г = 1, 5, 9,13,...; «минус» -

(к = 1,2,...); С2,=+ — /

I 21к)

при г = 3, 7,11,15,...).

Анализ решения (2.9) позволяет сделать вывод об отсутствии скачков и отрицательных значений температур внутри рассматриваемой области. 1,0

0,2 0,4 0,6 0,8 \ 1,0

Рис. 2.3. Изменение температуры в пластине с учетом конечной скорости распространения теплоты при Рог = 0,001 по формуле (2.8)

В третьей главе диссертации представлены результаты разработки метода решения задач теплообмена с учетом перемещения фронта плавления (задача Стефана с абляцией) на основе введения фронта температурного возмущения и дополнительных граничных условий.

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

В задачах Стефана с абляцией искомыми величинами являются распределение температуры в остающемся (нерасплавленном) слое вещества и закон перемещения

11

фронта фазового перехода во времени. Для их определения особенно эффективными оказываются методы, основанные на введении фронта температурного возмущения и дополнительных граничных условий. В этом случае кроме фронта фазового перехода необходимо определять еще закономерность перемещения во времени фронта температурного возмущения (глубину прогретого слоя).

Подобная задача с использованием интегрального метода теплового баланса решается в два этапа. На первом этапе рассматривается задача о прогреве тела от начальной температуры Т0 до температуры плавления Тт на его поверхности. На втором этапе задача решается уже с учетом перемещения фронта плавления.

Математическая постановка задачи определения температурного состояния для полубесконечного тела при граничных условиях первого рода до момента достижения температуры плавления на его поверхности имеет вид

^ = (х>0; 0*,«»); (3.1)

от дх

Т(х,0) = Т0; =

дх

где Г - тепловой поток, Вт/м2, Г0 - начальная температура, °С; а - коэффициент температуропроводности, м2/с.

В диссертационной работе задача (3.1) решена для четырнадцати приближений. Решение, например, в третьем приближении имеет вид:

Т(х,х) = Т0+ —

10 5+Ы гЫ + 51 в

(3.2)

где 6(т) = Г0 + 7гл/27ах/20 / X - глубина (толщина) прогретого слоя (фронт температурного возмущения). Из (3.2) получаем формулу для определения температуры на поверхности тела

Г(0,т) = Т0 +/>Ч/27ят/20/Х. (3.3)

Отличие решения (3.3) в третьем приближении от точного составляет 2,6 %. Решение в четырнадцатом приближении отличается от точного менее чем на 1 %. Для получения таких решений использовались дополнительные граничные условия, определяемые по следующим общим формулам, справедливым для любого числа приближений:

ЭТ(0,т)/сЬг"=0; а"чГ(6,т)/ах"4 =0; д"Т{Ь,х)/дх" =0 (и = 7,9,11,...).

После того как температура на поверхности х = 0 становится равной температуре плавления, наступает второй этап теплообмена с учетом перемещения фронта плавления. Отметим, что фронт температурного возмущения 5(т) при т = тпл, имеющий в момент времени начала плавления некоторую конкретную величину, определяется из решения задачи (3.1).

Математическая постановка задачи для второго этапа теплообмена имеет вид

= (х>0; 0 < х < со); (3.4)

от дх

г(*,0) = г.(*); (3.5)

П*,т) = Т„л; (3.6)

Г(б,т) = Г0;

ЭГ(5,т)

дх

= 0;

дх (к

(3.7) (3-8)

(3.9)

где s(z) - перемещение фронта плавления во времени; Га1 - температура плавления; <2 - теплота фазового перехода; р - плотность; к - коэффициент теплопроводности.

Начальное распределение температуры Тн(х) находится из решения задачи (3.1). Например, в первом приближении оно имеет вид

Тн (х) = Т0+ Р(х - 8(т„)У /(2Щтт)У.

(3.10)

В диссертации получено следующее аналитическое решение задачи (3.4) - (3.9):

0(*,т ) = к

(8-х)1

-{к-1)

(5-х)2

(¿ = 1,2,3,...). (3.11)

Результаты расчетов по формуле (3.11) приведены на рис. 3.1. Исходные данные для решения задачи были следующими:

X. = 0,19 Вт/(м к); <2 = 60000 Дж/кг; р = 1600 кг/м3; Т0 = 20 °С; Гш = 80 °С; F = 50000 Вт/м2; а = 1,35 • 10"6 м2/с.

Рис. 3.1. Изменение температуры в диапазоне ¿(т) < х < б(т): 1,2,3-соответственно первое, второе и четырнадцатое приближения

\

V -Юс ' т =30с

\

V 1/2из \ 1 2 3

0.02 0,03 0,04 X 0,05

Анализ результатов, связанных с задачей Стефана, позволяет сделать вывод, что температурное поле рассматриваемого тела заключено между двух фронтов - фронта плавления и фронта температурного возмущения, движущихся с различными скоростями во времени. Важным является тот факт, что разработанный в третьей главе метод позволяет выполнять оценку температурного состояния для сверхмалых значений временной и пространственной координат.

В четвертой главе диссертации приводятся результаты исследований, связанных с получением аналитических решений задач для ламинарного динамического (уравнения Прандтля) и теплового (уравнения Польгаузена) пограничных слоев, которые соответственно имеют вид:

дих 0иг д2их .... д!(х ,у) д((х,у) д21(х,у) ....

и,—- + о—- = у—т-; (4.1) иг +и„ =а—(4.2)

' дх у ду ду2 " дх у ду ду2

где и^, и,, - составляющие скорости по осям х и у; V - кинематическая вязкость;

/ - температура; а - коэффициент температуропроводности.

Граничные условия для уравнений (4.1) и (4.2) по форме идентичны и, например, для уравнения (4.1) они имеют вид

К = иЛ=о =0; = и=

8у

= 0;

= 0, (4.3)

о

где б(дг) - толщина динамического пограничного слоя.

Для получения решения уравнения (4.1), (4.2) путем их осреднения сводятся соответственно к интегральным уравнениям Кармана и Г.Н. Кружилина:

йу

у=0

(1х> * с1у

(4.4)

(4.5)

где Т(х,у) = 1(х,у)-^- 7ф=гср-гст.

Решение интегральных уравнений (4.4), (4.5) с соответствующими граничными условиями вида (4.3) принималось в виде ряда, неизвестные коэффициенты которого находились из основных (4.3) и дополнительных граничных условий, задаваемых на фронте динамического (теплового) возмущения или, что то же самое, на границах пограничных слоев.

Полученные аналитические решения позволили уточнить критериальные уравнения для коэффициентов теплоотдачи и касательных напряжений.

В диссертации также получено решение задачи для теплового пограничного слоя в случае, когда вместо граничного условия первого рода на стенке задано граничное условие третьего рода. Результаты расчетов для этой задачи в сравнении с известным точным решением Шлихтинга, полученным численным методом, отличаются менее, чем на 1 % и представлены на рис. 4.1.

о 1 2 3 4 ¡46

Рис. 4.1. Распределение температур @ = Т/Тср от координаты

r| = ^^N/u7vx . 1, 2, 3, 4- соответственно первое, второе, третье, четвертое приближения; 5 - точное решение: РН = а фх/и Д; Рг = 0,72 В пятой главе диссертации рассматриваются проблемы, связанные с разработкой способов получения вихревых полей потенциалов различной природы (температурных, диффузионных, фильтрационных, электрических и др.). В частности, показано, что для получения вихревых полей потенциалов в замкнутом цилиндрическом или

| сферическом пространстве необходимо организовать вращение источников этих полей по круговой орбите. При этом скорость вращения источников строго связана со скоростью распространения потенциалов указанных полей.

Для получения вихревого температурного поля было организовано вращение те-I плового источника по круговой орбите. Математическая постановка задачи в данном случае имеет вид:

1 дТ(х,у,х) = д2Т(х,у,х) д2Т(х,у,х)

I а дт дх2 ду2 ' (5Л)

Т(х,у,0) = Т0- (5.2)

Т(хт,уг,х) = Т0+Т^ 5ш(сот - у), (5.3)

где Т0 - начальная температура; Гм - амплитуда колебаний температуры стенки; хг'Ут ~ граничные значения координат цилиндрической камеры; © = 2т1У; V - частота колебаний; у - сдвиг фазы колебаний каждого последующего источника относитель-! но предыдущего.

| Декартова система координат в задаче (5.1)-(5.3) для цилиндрической камеры

была принята для удобства решения этой задачи численным методом расщепления и прогонки. Исходные данные для решения модельной задачи (5.1) - (5.3) были приняты следующими: а = 6-1(Г5м2/с; Л = 0,36м-радиус цилиндра; Го=50°С; 7„=250°С.

Имитацию вращения источника выполняли 146 источников теплоты, расположенных по окружности цилиндрической камеры. Изменение температуры в точках расположения источников соответствовало соотношению (5.3). Величина сдвига фазы колебаний каждого последующего (по часовой стрелке) источника по отношению к предыдущему была принята равной = 1 /146.

Распределение изотерм при т = 100 с (квазистационарное состояние) дано на рис. 5.1. Их анализ позволил установить, что изотермы имеют вихревой характер распределения. При этом происходит вращение изотерм в окружном направлении при одновременном их перемещении к центру цилиндра в радиальном направлении. Следует отметить, что движение тепловых потоков происходит по вихревым линиям, причем в противоположных направлениях в различных ветвях этих линий.

(т - 100 с, квазистационарное состояние): направление движения тепловых потоков

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

В диссертации приводится инструкция по использованию программы расчеты гидравлических систем, в которой дается общая характеристика программы, последовательность формирования элементов расчетной схемы (узлы сети, участки сети, насосы, запорные устройства и проч.), рекомендации по установке программы на компьютере, структура меню, действия, выполняемые программой.

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

В диссертации даны результаты построения компьютерной модели теплосети объединенного централизованного теплоснабжения г. Саратова (от Саратовской ТЭЦ-5 и Саратовской ГРЭС). Результаты расчетов выводятся как в цифровом, так и в графическом виде. Например, распределение давлений для участка теплосети г. Саратова приведено на рис. 6.1.

250 Р, м вод. ст.

160

130

100

70

40

10

1 В23

3 К-6 -и

......

. <3 ) НС-6

.......... .......... ..1

11

Ь,км 15

1 3 5 7 У Рис. 6.1. Эпюры распределения давлений в прямой и обратной магистралях первого вывода ТЭЦ-5 Расчеты, выполненные на основании предложенных моделей, показали крайне нерациональное распределение потребителей теплоты по конкретным ее источникам, неоптимальное «дерево» из тепловых выводов, построенное в Саратовских тепловых сетях, низкую эффективность вновь построенного участка трубопровода ТЭЦ-5 длиной 4,5 км и диаметром 1200 мм. 16

В заключении перечислены основные результаты, полученные в диссертационной работе.

1. Разработана новая математическая модель процесса теплообмена с учётом конечной скорости распространения теплоты, позволяющая, в отличие от существующих, снять противоречия, связанные с наличием скачков температуры внутри тела и с появлением отрицательных температур в обратной тепловой волне. На основе совместного использования методов Фурье и ортогонального метода Бубнова-Галеркина получены точные аналитические решения гиперболических уравнений теплопроводности для пластины при симметричных граничных условиях первого рода. Впервые показано, что процесс теплообмена с учетом конечной скорости распространения теплового возмущения (с учетом инерционности теплового потока) характеризуется чередованием прямых и обратных тепловых волн, на каждой из которых имеет место скачок температуры.

2. Используя интегральные методы теплового баланса на основе использования фронта температурного возмущения и дополнительных граничных условий разработан метод решения сложных линейных и нелинейных краевых задач тепломассопереноса, позволяющий получать решения для сверхмалых значений временной и пространственной координат.

3. На основе использования дополнительных граничных условий с учетом движения фронта температурного возмущения и фронта плавления впервые получено аналитическое решение задачи Стефана с абляцией (с удалением расплавляемой среды). Исследования показали практически линейный характер изменения фронта плавления во времени. При этом скорость перемещения фронта плавления возрастает с увеличением теплового потока на поверхности плавления и уменьшается с увеличением теплоты фазового перехода (теплоты плавления).

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

5. Разработаны теоретические основы метода получения вихревых полей потенциалов (тепловых, диффузионных, электрических и прочих), основанного на организации вращения источников этих полей по круговым орбитам. Выполненные в диссертации применительно к вихревому температурному полю исследования позволили установить, что распределение температур и тепловых потоков внутри окружности вращения источника имеет вихревой характер. Причем движение тепловых потоков, возникающих от выравнивания минимумов и максимумов температур происходит в противоположных направлениях.

6. Построена математическая модель теплосети централизованного теплоснабжения г. Саратова от Саратовской ТЭЦ-5 и ГРЭС и разработано программное обеспечение для её реализации. В основе построения модели лежат два закона Кирхгофа и теория графов. Модель позволяет определять давления, скорости, расходы и температуры теплоносителя в любой точке теплосети в текущий момент времени. Выполненные на модели исследования показали, что построенное в теплосетях «дерево» теплосети является неоптимальным, что приводит к значительным потерям напора, проявляющимся в уменьшении располагаемого перепада давлений между прямым и обратным трубопроводом.

Список монографий:

[1] Кудинов В.А., Кудииов И.В. Аналитические методы теплопроводности // Монография. Самара: Самарский государственный архитектурно-строительный университет. 2011. - 272 с.

[2] Кудинов В.А., Кудинов И.В. Методы решения параболических и гиперболических уравнений теплопроводности // Монография. - М.: Книжный дом «Либроком», 2011.-270 с.

Список основных публикаций в рецензируемых журналах из перечня ВАК:

[3] Аверин Б.В., Стефанюк Е.В., Кудинов КВ. Получение аналитического решения уравнений гидродинамического пограничного слоя на основе введения дополнительных граничных условий // Известия Самарского научного центра Российской Академии наук. Специальный выпуск. «Актуальные вопросы тепло- и массообмена, энергоэффективность, исследование вихревых закрученных потоков». - 2008. - С. 39 - 46.

[4] Стефанюк Е.В., Кудинов И.В., Ларгина Е.В. Построение приближенных аналитических решений нелинейных обыкновенных дифференциальных уравнений на основе использования дополнительных граничных условий // Вестник Самарского государственного технического университета. Серия Физико-математические науки. -№ 1 (18).-2009.-С. 122-132.

[5] Стефанюк Е.В., Кудинов И.В. Построение аналитических решений уравнений динамического и теплового пограничных слоев // Вестник Самарского государственного технического университета. Серия Технические науки.- №2(24).- 2009,-С. 205 -211.

[6] Радченко В.П., Кудинов И.В. Получение аналитических решений нелинейных задач теплопроводности на основе введения дополнительных граничных условий // Вестник Самарского государственного технического университета. Серия Физико-математические науки. - № 1 (20). - 2010. - С. 161 - 170.

[7] Кудинов И. В. Использование компьютерной модели для проектирования тепловых сетей // Вестник Самарского государственного технического университета. Серия Технические науки. - №4(27).-2010.-С. 174-181.

[8] Кудинов В.А., Кудинов И.В. Об одном методе получения точного аналитического решения гиперболического уравнения теплопроводности на основе использования ортогональных методов // Вестник Самарского государственного технического университета. Серия Физико-математические науки. - № 5 (21). - 2010. - С. 159 - 170.

[9] Стефанюк Е.В., Кудинов И.В. Получение аналитического решения уравнения теплового пограничного слоя на основе введения дополнительных граничных условий // Вестник Самарского государственного аэрокосмического университета имени академика С.П. Королева. - № 1 (21). - 2010. - С. 174 - 184.

[10] Кудинов И.В. Построение компьютерных моделей систем теплоснабжения больших городов // Вестник Самарского государственного технического университета. Серия Технические науки. - № 1 (29). - 2011. - С. 212 - 219.

[11] Кудинов В.А., Кудииов И.В. Теоретические основы получения вихревых полей потенциалов // Известия РАН. Энергетика. - № 1. - 2011. - С. 105 - 112.

В других изданиях:

[12] Стефанюк КВ., Кудинов И.В. Ан алитические решения уравнений динамического и теплового пограничного слоя при граничных условиях первого и третьего рода // Проблемы газодинамики и тепломассообмена в аэрокосмических технологиях: Труды XVII Школы-семинара молодых ученых и специалистов под руководством академика РАН Леонтьева А.И. Т.2. - Жуковский: ЦАГИ, 2009. - С. 139 -142.

[13] КудиновИ.В., Еремин A.B. Аналитическое решение уравнений теплового пограничного слоя с учетом диссипации энергии // Проблемы тепломассообмена и гидродинамики в энергомашиностроении. Материалы докладов VII Школы-семинара молодых ученых и специалистов под руководством академика РАН Алемасова В.Е.-Казань: КГТУ, 2010. - С. 184-187.

[14] Кудинов В.А., Кудинов И.В. Теплообмен при течении Куэтта с учетом диссипации энергии // Известия высших учебных заведений и энергетических объединений СНГ. Энергетика. - 2011. - № 2. - С. 43 - 51.

[15] Кудинов В.А., Кудинов И.В. Математическое моделирование процессов теплообмена и фазовых превращений с учетом абляции // Инженерно-физический журнал. Т. 84.-2011.-№5.-С. 1065-1074.

[16] Стефанюк Е.В., Кудинов И.В., ЛаргинаЕ.В. Обобщенные функции в задачах теплопроводности для многослойных конструкций // Вестник Самарского государственного технического университета. Серия: Математическая. - 2008. - № 2 (8). -С. 41-56.

[17] Аверин Б.В., Кудинов И.В., ЛаргинаЕ.В., Еремин A.B. Обобщенные функции в нелинейных задачах теплопроводности для многослойных конструкций // Повышение энергоэффективности зданий и сооружений. Межвузовский сборник научных трудов. Самара: СГАСУ. - 2009. - С. 45 - 52.

[18] КудиновИ.В., Ларгина Е.В., Еремин A.B., Будыльников М.Н. Теплообмен в плоском канале с учетом диссипации энергии // Вестник Самарского государственного технического университета. Серия: Математическая. -2009. - №2 (10). -С. 38-47.

[19] КудиновИ.В., ЛаргинаЕ.В., Котов В.В., Биктагирова Е.Ю. Дополнительные граничные условия в задачах теплопроводности с переменными физическими свойствами среды // Вестник Самарского государственного технического университета. Серия: Математическая. -2009. -№2 (10). - С. 48 - 55.

[20] Кудинов И.В. Аналитические решения уравнений теплового пограничного слоя с учетом диссипации энергии // Повышение энергоэффективности зданий и сооруже-

ний. Межвузовский сборник научных трудов. Самара: СГАСУ. - Вып. 5,- 2010. С. 146- 150.

[21] Стефанюк Е.В., КудиновИ.В., ЛаргинаЕ.В. Математическое моделировани теплопроводности в многослойных конструкциях на основе теории обобщенны функций // Математическое моделирование физических, экономических, технических социальных систем и процессов. Сборник трудов Седьмой Международной конферен ции. Ульяновск: УлГУ. -2009. - С. 254.

[22] Стефанюк Е.В., Кудинов И.В. Математические модели гидродинамического теплового пограничных слоев с учетом дополнительных граничных условий // Мате матическое моделирование физических, экономических, технических, социальнь систем и процессов. Сборник трудов Седьмой Международной конференции Ульяновск: УлГУ. - 2009. - С. 255.

[23] Стефанюк Е. В., Кудинов КВ., ЛаргинаЕ.В. Построение аналитических реше ний уравнений динамического и теплового пограничных слоев // Математическое мо делирование и краевые задачи. Сборник трудов Шестой Всероссийской научной кон ференции с международным участием. Часть 2,- Самара: СамГТУ. - 200 С. 187-191.

[24] Кудинов И.В., ЛаргинаЕ.В., Еремин A.B., Будильников М.Н. Теплообмен плоском канале с учетом диссипации энергии // Математическое моделирование краевые задачи. Труды седьмой Всероссийской научной конференции с международ ным участием. Часть 2. - Самара: СамГТУ. -2010. - С. 148 - 152.

[25] Кудинов И.В., ЛаргинаЕ.В., Котов В.В., Биктагирова Е.Ю. Математическо моделирование теплопроводности при переменных физических свойствах среды // Ма тематическое моделирование и краевые задачи. Труды седьмой Всероссийской науч ной конференции с международным участием. Часть 2.- Самара: СамГТУ, - 201 С. 163 -167.

[26] Кудинов И.В. Получение приближённого аналитического решения задачи Сте фана с удалением расплавляемой среды // Математическое моделирование и краевы задачи. Труды восьмой Всероссийской научной конференции с международным уча стием. Часть 2. - Самара: СамГТУ. - 2011. - С. 66 - 72.

Заказ № 1109. Формат 60 х 84 1/16. Уч. изд. л. 1,00. Тираж 100 экз.

Отпечатано в типографии Самарского государственного технического университета 443100, Самара, Молодогвардейская, 244, корпус 8.

Автореферат отпечатан с разрешения диссертационного совета Д 212.217.03 ФГБОУ ВПО «Самарский государственный технический университет» (протокол № 6 от 16 ноября 2011 г.)

Общая характеристика работы.

1. Аналитический обзор и постановка задачи.

2. Получение аналитических решений гиперболических дифференциальных уравнений.

2.1 Аналитические решения гиперболических уравнений теплопроводности (бесконечная пластина при симметричных граничных условиях первого рода).

2.2 Приближенный метод решения'гиперболических уравнений.

2.3 Получение точных аналитических решений гиперболических уравнений теплопроводности со смешанными производными.

2.4 Выводы по главе 2.

3. Исследование процесса теплопроводности на основе введения фронта температурного возмущения и дополнительных-граничных условий.

3.1 Задачи теплопроводности для1 бесконечной пластины при симметричных граничных условиях первого рода.

3.2 Математическое моделирование процессов теплообмена. и фазовых превращений с учетом абляции.

3.3* Выводы по главе 3.

4. Аналитические решения уравнению динамического и теплового пограничных слоев.

4.1 Динамический и тепловой пограничные слои.

4.2 Подобие динамического и теплового пограничных слоев гидродинамическая <теория теплообмена).

4.3 Получение аналитического решения краевой задачи.для динамического пограничного слоя.

4.4 Получение аналитического решения краевой задачи для теплового пограничного - слоя <.

4.5 Получение аналитического решения краевой задачи для теплового пограничного слоя при граничных условиях третьего рода на стенке.

4.6 Выводы по главе 4.

5. Теоретические основы получения*вихревых полей потенциалов.

5.1 Математическое моделирование вихревого температурного поля».

5.2 Исследование структуры вихревого движения при вращении источниковгполей потенциалов по круговым орбитам.

5.3 Формирование вихревых полей потенциалов при вращении точечных источников по круговым орбитам.

5.4 Выводы по главе 5.

6. Алгоритмическое и программное обеспечение для трубопроводных систем.

6.1 Теоретические основы расчета кольцевых разветвленных гидравлических сетей.

6.2 Основные принципы разработки и построения компьютерных моделей гидравлических сетей.

6.3 Применение компьютерных моделей для исследования систем теплоснабжения больших городов.

6.4 Выводы по главе 6.

Актуальность, работы; Разработка, исследование, обоснование математических моделей для процессов тепломассопереноса и гидродинамики, а также совершенствование методов, алгоритмов и программного обеспечения* для их реализации важны как для развития численных и аналитических методов решения дифференциальных уравнений, описывающих эти процессы, так и для создания эффективной информационной технологии получения, наполнения и прогнозирования новых знаний при исследовании широкого круга линейных и нелинейных явлений с учётом и без учёта конечной скорости распространения разного рода возмущений (тепловых, гидравлических, диффузионных и др.).

Существенным фактором, ограничивающим возможности математического моделирования в области тепломассопереноса, является практическое отсутствие- точных аналитических методов исследования нелинейных краевых задач, в том числе гиперболического типа. Привлечение гиперболического дифференциального оператора необходимо при исследовании полей различных потенциалов при сверхмалых значениях времени (обработка, материалов лазерным излучением, распределение температуры в начальной стадии теплового удара, нагрев при динамическом распространении трещины, движение космических аппаратов в атмосфере и ДР-)

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

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

В связи с этим возникает потребность в разработке методов, позволяющих без введения каких-либо существенных допущений получать решения в аналитическом виде, пусть даже ценой выполнения большого количества предварительных вычислительных процедур, связанных с использованием компьютерной техники. Эти методы должны сочетать в себе положительные стороны классических аналитических (аналитический вид решения) и численных (решения сложнейших краевых задач без введения упрощающих допущений) методов при возможности избегать наиболее существенные их недостатки. В этом направлении наиболее эффективными являются методы, в которых используется понятие глубины термического слоя (интегральные методы теплового баланса). Их разработке и развитию посвящены труды Био М., Гудмена Т., Швеца М.Е., Вейника А.И., Беляева Н.М.,. Дородницына A.A., Глазунова Н.Т., Карташова Э.М. и других. Введение при использовании этих методов* дополнительных граничных условий во многих случаях позволяет получать . приближенные аналитические решения практически с заданной степенью точности. Развитие этих методов применительно к решению сложных линейных и нелинейных краевых задач является-* актуальной научной проблемой, имеющей большое теоретическое и прикладное значение.

Цель диссертационной« работы состоит в математическом моделировании процессов теплопроводности и гидродинамики на основе эффективных аналитических и численно-аналитических методов* решения краевых задач с использованием интегральных методов теплового баланса, ортогональных методов. Канторовича JI.В: и Бубнова-Галеркина, новых алгоритмов и специального программного обеспечения.

Для достижения указанной цели решались следующие задачи.

1. Разработка методов получения точных и приближенных аналитических решений гиперболических уравнений теплопроводности и движения жидкостей, описывающих распределение соответствующих полей потенциалов учетом конечной скорости распространения возмущений.

2. Построение аналитического решения задач динамического и теплового пограничных слоев при граничных условиях первого и третьего рода на стенке.

3. Разработка численно-аналитического метода решения краевой задачи Стефана с абляцией (при удалении расплавляемого вещества).

4. Разработка и исследование метода получения вихревых полей потенциалов. (тепловых, гидравлических, диффузионных и прочих) путем организации движения их источников по круговым орбитам с учетом найденных закономерностей этого движения.

5. Создание ' математических моделей сложных разветвленных трубопроводных систем (водо- нефте- газо- проводов и проч.) и программного комплекса для исследования процессов, протекающих в реальных системах, с определением давлений, скоростей, расходов и температур движущейся среды.

Научная новизна полученных в диссертационной работе результатов заключается в следующем:

1. Разработана новая математическая модель процесса теплообмена с учетом конечной скорости распространения теплоты, позволяющая, в отличие от существующих, снять противоречия, связанные с наличием скачков температуры внутри тела и с появлением отрицательных температур в обратной тепловой волне.

2. На основе введения фронта температурного возмущения и дополнительных граничных условий получены новые численно-аналитические решения задач динамического и теплового пограничных слоев при граничных условиях первого и третьего рода на стенках.

3. На основе использования дополнительных граничных условий с учетом движения фронта температурного возмущения и фронта плавления впервые получено аналитическое решение задачи Стефана с абляцией (с удалением расплавляемого вещества).

4. Теоретически разработан метод получения вихревых полей потенциалов, возникающих при вращении по определенным закономерностям источников этих полей по круговым орбитам.

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

Практическая^ значимость работы. Разработанные в диссертации методы и полученные аналитические решения отличаются заметной простотой конструкции при точности, достаточной для прикладных задач. Такие решения полезны в тех случаях, когда решения задач теплопроводности являются, промежуточными стадиями других исследований, например,- решения обратных задач, задач термоупругости, автоматизированного- управления- и проектирования. В частности, полученные в диссертации аналитические решения задач теплообмена при течении жидкостей в трубах и плоских каналах были использованы при разработке математических моделей и программного комплекса для теплосети центрального теплоснабжения г. Саратова (от СарГРЭС и ТЭЦ-5).

Математические модели теплосети и циркуляционной системы Новокуйбышевской ТЭЦ-1 позволили определить давления, скорости, расходы теплоносителя' и потери напора в любой точке моделируемых систем, найти наиболее оптимальные режимы текущей работы, выполнить предварительные проекты реконструкций, а также составить планы построения новых участков трубопроводных систем:

Применительно к расчетам температурного состояния труб барабанов котлов Новокуйбышевской ТЭЦ-1 внедрена методика расчета, основывающаяся на введении фронта температурного возмущения и дополнительных граничных условий. Для расчетов температурного состояния труб, представленных в виде многослойных конструкций (тепловая изоляция — металлическая стенка трубы — отложения накипи на внутренней поверхности труб), внедрена методика, основанная на теории обобщенных функций.

На1 защиту выносятся следующие основные результаты и положения.

1. Точные и приближённые численно-аналитические решения гиперболических уравнений теплопроводности, учитывающих конечную скорость распространения возмущений; исследование сходимости и погрешности численно-аналитических решений.

2. Численно-аналитические решения задач динамического и теплового пограничных слоев при граничных условиях первого и третьего рода на стенке, полученные на основе введения фронта температурного возмущения и дополнительных граничных условий; оценки сходимости и погрешности решений. .

3. Численно-аналитическое решение задачи Стефана с абляцией (с удалением расплавляемого вещества), полученное на основе дополнительных граничных условий и учёта движениям фронтов температурного возмущения и плавления; оценки сходимости и погрешности решения.

4. Математическая модель вихревых полей потенциалов и численный метод их расчёта, основанный на конечно-разностном методе расщепления и методе прогонки.

5. Программный комплекс для исследования математических моделей сложных разветвленных многокольцевых трубопроводных систем и численный метод их расчёта, позволяющий рассчитывать давления, скорости, расходы и температуру движущихся сред.

6. Результаты исследований на основе программного комплекса математической модели теплосети г. Саратова (от Саратовских ТЭЦ-5 и ГРЭС), позволившие установить причины недостаточной эффективности работы теплосети, оптимизировать планы реконструкции и выполнения ремонтных работ.

Апробация работы. Основные результаты работы доложены и обсуждены на Всероссийской научной конференции молодых ученых «Наука, технологии, инновации» (г. Новосибирск, 2007 г.); Тринадцатой Международной научно-технической конференции студентов и аспирантов «Радиоэлектроника, электротехника и энергетика» (г. Москва, 2007 г.); Четырнадцатой Международной научно-практической конференции студентов, аспирантов и молодых ученых «Современные техника и технологии» (г. Томск, 2008 г.); Третьей молодежной Международной научной конференции «Тинчуринские чтения» (г. Казань, 2008 г.); Седьмой Международной конференции «Математическое моделирование физических, экономических, технических, социальных систем и процессов» (г. Ульяновск, 2009т.); Шестой, Седьмой и Восьмой Всероссийских научных конференциях с международным участием «Математическое моделирование и краевые задачи» (г.Самара, 2009, 2010, 2011 гг.); в Школе-семинаре молодых ученых и специалистов под рук. акад. РАН А.И. Леонтьева «Проблемы газодинамики и теплообмена в аэрокосмических технологиях» (г. Москва, 20092011 гг.); в Школе-семинаре молодых ученых и специалистов под рук. акад. РАН В:Е. Алемасова «Проблемы тепломассообмена и гидродинамики в энергомашиностроении» (г.Казань, 2010 г.); на научном семинаре «Прикладная математика и механика» Самарского государственного технического университета под рук. профессора Радченко В-П. в 2010 и 2011 гг.

Связь диссертационной работы с планами научных исследований; Представленная» работа является обобщением, теоретических и экспериментальных исследований, выполненных автором- в Самарском государственном техническом университете: Исследования проводились по планам госбюджетных тематик Минвуза РФ № 1.18.09 «Получение аналитических решений задач теплопроводности на основе теории обобщенных функций и дополнительных граничных условий», № 1.21.11 «Разработка методов получения точных аналитических решений дифференциальных уравнений гиперболического типа», а также по Аналитической? ведомственной« целевой программе «Развитие научного потенциала высшей школы», тематический план- НИР:№ 551702 «Разработка нового направления? получения? аналитических. решений задач математической- физики; на основе введения: фронта температурного возмущения и дополнительных граничных условий». ' :

Результаты работы, использовались при выполнении энергетического аудита Самарского государственного; технического университета, хоздоговорных работ с Волжской территориальной- генерирующей кампанией, с Куйбышевским и Новокуйбышевским нефте-. перерабатывающими заводами, ОАО «Самараоргсинтез».

Экономический эффект, подтвержденный^ соответствующим актом внедрения; приведенным в приложениях диссертации, составляет 3 миллиона рублей.

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

Публикации; Материалы диссертации опубликованы в 38 печатных работах. В* автореферате приведено 26 основных научных работ, из которых 9 статей опубликованы в рецензируемых журналах, из перечня ВАК, 5 статей - в центральных и академических изданиях. По результатам исследований опубликованы две монографии, одна из которых— в центральном издательстве «URSS».

Личный вклад автора. Работы [74, 77, 87, 93] выполнены самостоятельно. В основных работах [68 - 73, 75, 76, 78] диссертанту принадлежит постановка проблем исследований, непосредственное выполнение основной части работы, которая выполнена совместно с другими авторами. В остальных работах [79 - 86, 88 - 92], также опубликованных в соавторстве, диссертанту в равной степени с другими авторами принадлежат постановки задач, получение решений, содержание работы и анализ результатов.

Структура и объем диссертации. Диссертационная работа состоит из введения, шести глав, выводов, списка используемой литературы, приложений; изложена на 150 страницах основного машинописного текста и 31 странице приложений, содержит 74 рисунка, 6 таблиц. Список использованной литературы включает 93 наименования.

6.4 Выводы по главе 6

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

1. В сложных тепловых сетях, содержащих большое число кольцевых структур, необходимо строить оптимальное «дерево» Для этого перекрытие задвижек выполняется лишь на участках сети с минимальными скоростями течения теплоносителя либо в застойных зонах (зонах естественного разделения тепловыводов). Эффективным средством построения оптимального «дерева» теплосети является ее компьютерная модель, позволяющая определить необходимое число задвижек и участки, на которых их наиболее целесообразно установить.

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

ЗАКЛЮЧЕНИЕ

1. Разработана новая математическая модель процесса теплообмена с учетом конечной скорости распространения теплоты, позволяющая, в отличие от существующих, снять противоречия, связанные с наличием скачков температуры внутри тела и с появлением отрицательных температур в обратной тепловой волне. На основе совместного использования методов Фурье и ортогонального метода Бубнова-Галеркина получены точные аналитические решения гиперболических уравнений теплопроводности для пластины при симметричных граничных условиях первого рода. Впервые показано, что процесс теплообмена с учетом конечной скорости распространения теплового возмущения (с учетом инерционности теплового потока) характеризуется чередованием прямых и обратных тепловых волн, на каждой из которых имеет место скачок температуры.

2. Используя интегральные методы теплового баланса на основе использования фронта температурного возмущения и дополнительных граничных условий разработан метод решения сложных линейных и нелинейных краевых задач тепломассопереноса, позволяющий получать решения для сверхмалых значений временной и простраственной координат.

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

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

5. Разработаны теоретические основы метода получения вихревых полей потенциалов (тепловых, диффузионных, электрических и прочих), основанного на организации вращения источников этих полей по круговым орбитам. Выполненные в диссертации применительно к вихревому температурному полю исследования позволили установить, что распределение температур и тепловых потоков внутри окружности вращения источника имеет вихревой характер. Причем движение тепловых потоков, возникающих от выравнивания минимумов и максимумов температур происходит в противоположных направлениях.

6. Построена математическая модель теплосети централизованного теплоснабжения г. Саратова от Саратовской ТЭЦ-5 и ГРЭС и разработано программное обеспечение для ее реализации. В основе построения модели лежат два закона Кирхгофа и теория графов. Модель позволяет определять давления, скорости, расходы и температуры теплоносителя в любой точке теплосети в текущий момент времени. Выполненные на модели исследования показали, что построенное в теплосетях «дерево» теплосети является неоптимальным, что приводит к значительным потерям напора, проявляющимся в уменьшении располагаемого перепада давлений между прямым и обратным трубопроводом.

1. Баумейстер К., Хатипл Т. Гиперболическое уравнение теплопроводности. Решение задачи о полубесконечном теле // Теплопередача. 1969, № 4. С. 112 119.

2. Беляев Н.М., Рядно A.A. Методы нестационарной теплопроводности. М.: Высш. школа, 1978. 328 с.

3. Беликов В.К, Шаронова О.В., Бойко Г.П. Определение эффективного значения коэффициента температуропроводности плоской сложной системы // Теплообмен и гидродинамика: Сб. науч. тр. Красноярск, 1981. С. 35 43.

4. Бровкин Л.А., Гузов Л.А. Инженерный расчет нагрева многослойной пластины при граничных условиях 1-го рода //Изв. вузов СССР. Сер. Энергетика. 1985, № 9. С. 94 97.

5. Власов A.A. Статистические функции распределения. М.: Наука, 1966.

6. Гелъфанд И.М., Шилов Г.Е. Обобщенные функции и действия над ними. М.: Физматгиз, 1959. 470 с.

7. Гудмен Т. Применение интегральных методов в нелинейных задачах нестационарного теплообмена //Проблемы теплообмена. Сб. науч. тр. М.: Атомиздат, 1967. С. 41 96.

8. Жуковский B.C. Основы теплопередачи. М-Л.: Госэнергоиздат. 1960.

9. И. Канторович Л.В. Использование идеи метода Галеркина в методе приведения к обыкновенным дифференциальным уравнениям // Прикл. мат. и механ. Т. 6. № 1, 1942. С. 31 40.

10. Канторович Л.В., Крылов В.И. Приближённые методы высшего анализа. Л.: Физматгиз, 1962. 708 с.

11. Карташов Э.М. Аналитические методы в теории теплопроводности твёрдых тел. М.: Высш. школа, 1985. 480 с.

12. Карташов Э.М. Аналитические методы в теории теплопроводности твёрдых тел. М.: Высш. школа, 2001. 550 с.

13. Кеч В., Теодореску П. Введение в теорию обобщенных функций с приложениями в технике. М.: Мир, 1978.

14. Коздоба Л.А., ЧарныйН.Д. Эквивалентность тепловых режимов и эффективные характеристики однородных и многослойных оболочек // Промышленная теплоэнергетика. 1984. Т. 6. № 2. С. 14-21.

15. Коздоба Л.А. Принцип эквивалентности в теории теплопроводности // Тепломассообмен — VII. Т. 7. Минск: Ин-т тепло- и массообмена им. A.B. Лыкова АН БССР, 1984. С. 34 39.

16. КоляноЮ.М. Применение обобщенных функций в термомеханике кусочно-однородных тел. // В кн. Математические методы и физико-механические поля. Киев: Наукова думка. Вып. 7, 1978. С. 7 — 11.

17. Кудниов В.А., Аверин Б.В., СтефанюкЕ.В. Теплопроводность и термоупругость в многослойных конструкциях. Учебное пособие. М.: Высшая школа, 2008. 305 с.

18. Кудинов В.А., Карташов Э.М. Техническая термодинамика. Учебное пособие для вузов. 5-ое издание. М.: Высшая школа, 2007. 261 с.

19. Кудинов В.А., Аверин Б.В., СтефанюкЕ.В., Назаренко С. А. Анализ нелинейной теплопроводности на основе определения фронта температурного возмущения // Теплофизика высоких температур. Т. 44. № 5, 2006. С. 577-585.

20. Кудинов В.А., Стефанюк Е.В. Задачи теплопроводности на основе определения фронта температурного возмущения // Известия АН. Энергетика. № 4, 2008. С. 122 138.

21. Кудинов В. А., СтефанюкЕ.В. Аналитический метод решения задач теплопроводности на основе введения фронта температурного возмущения и дополнительных граничных условий // ИФЖ. Т. 82. № 3, 2009. С. 540 558.

22. Кудинов В.А. Метод координатных функций в нестационарных задачах теплопроводности // Изв. АН-Энергетика (обзор): № 3> 2004. С. 82 104.

23. Кудинов В.А., Карташов Э.М., Калашников В.В. Аналитические решения задач тепломассопереноса и термоупругости для многослойных конструкций. Учеб. пос. для втузов. М.: Высшая школа, 2005. 340 с.

24. Кудинов В.А., Карташов Э.М. Гидравлика. Учеб. пособ. для вузов. Третье издание. М.: Высшая школа,,2008. 200 с.

25. Кудинов-В.А., Кудинов A.A., Воробьев Б.В., Денисов А.Ю: Метод сведения задач теплопроводности для многослойных конструкций к однослойным // Изв. АН. Сер: Энергетика. 1993, № 3. С. 135 — 142.

26. Лыков A.B. Теория'теплопроводности. М.:„Высшая школа, 1967. 600 с.

27. Лыков A.B. Тепломассоперенос: Справочник. М.: Энергия, 1978. 480 с.

28. Лыков A.B., Михайлов Ю. А. Теория тепло- и массопереноса. М.: Госэнергоиздат, 1963. 535 с.

29. Лойцянский Л.Г. Механика жидкости и газа. Учеб. для вузов. 7-ое изд. М.: Дрофа, 2003. 840 с.

30. Михеев М.А., МихееваИ.М. Основы теплопередачи. М.: Энергия, 1977. 343 с.

31. ОбразцовИ.Ф., ОнановГ.Г. Строительная механика скошенных тонкостенных систем. М.: Машиностроение, 1973. 659 с.

32. Онанов Г.Г. Уравнения с сингулярными коэффициентами типа дельта-функция и ее производных //Докл. АН СССР. Т. 191. №5, 1970. С. 997-1000.

33. Пехович А.И., Жидких В.М. Расчеты теплового режима твердых тел. JL: Энергия, 1976.

34. Петухов B.C. Теплообмен и сопротивление при ламинарном течении жидкости в трубах. М.: Энергия, 1967. 412 с.

35. Подстригай Я. С., Ломакин В.А., Коляно Ю.М. Термоупругость тел неоднородной структуры. М.: Наука, 1984. 368 с.

36. Прибытков И. А., Левицкий И.А. Теоретические основы теплотехники. М.: Издательский центр «Академия», 2004. 464 с.

37. ПостольникЮ.С. Метод осреднения функциональных поправок в задачах теплопроводности // Тепло- и массоперенос: Сб. тр. Минск. Т. 8, 1972. С. 23-29.

38. Стефанюк Е.В. Модельные представления аналитических решений краевых задач теории теплообмена на основе введения дополнительных граничных условий // Дисс. докт. техн. наук. Москва. МАТИ. 2010.

39. Теория тепломассообмена. Учебник для . вузов. Под ред. А.И. Леонтьева. М.: Высшая школа, 1979. 495 с.

40. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. М.: Наука, Л974.43 . Форсайт Дж., Малькольм М., Моулер К. Машинные методы математических вычислений. М.: Мир, 1980. 279 с.

41. Цирельман Н.М. Прямые и обратные задачи тепломассопереноса. Энергоатомиздат, 2005. 392 с.

42. Цой П.В. Методы расчета отдельных задач тепломассопереноса. М.: Энергия; 1971. 382 с.

43. Шумаков Н.В. Метод последовательных интервалов в теплометрии нестационарных процессов. М.: Атомиздат, 1979. 212 с.

44. ШлихтингГ. Теория пограничного слоя. М'.: Наука, 1969. 472 с.

45. Швец М.Е. О приближенном решении некоторых задач гидродинамики пограничного слоя // Прикладная математики«и механика. Т. 13. № 3, 1949.

46. Шашков А.Г., Бубнов В.А., Яновский С.Ю. Волновые явления теплопроводности: системно-структурный- подход. Изд. 2-ое, доп. М.: Едиториал УРСС, 2004. 296 с.

47. Шиммел М.М., Бек, Доналдсон. Эффективный коэффициент теплопроводности многослойного композитного материала // Теплопередача. 1977, №3. С. 130-136.

48. Юдаев Б.Н. Теплопередача. М.: Высшая школа, 1981. 319 с.

49. Bui An Ton. An optimal control free boundary problem for the Navier-Stokes equations Nonlinear Analysis. Volume 63, Issues 5-7, 30 November 200515 December 2005, Pages 831-839.

50. BuriakJ. A pair of dual integral equation // Proc. Math. Soc. Edinburg, 1962. V. 13. №2. P. 179-189.

51. ErdelyiA., Sneddon I. Fractional integration and dual integral equations I I Canad. J. Math. 1962. V. 14. № 5. P. 685 698.

52. Guozhen Ln, Peiyong Wang. On the uniqueness of a solution of a two-phase free boundary problem. Journal of Functional Analysis. Article in Press.

53. Sneddon I., Srivastav R. Dual relation involving Fourier-Bassel siries// Proc. Roy. Soc. Edinburgh, 1964, V. 46; pt. 3, p. 150 160:

54. Tiegang Fang. A. note on the unsteady boundary layers over a. flat plate. International Journal of Non-Linear Mechanics. Volume 43, Issue 9, November 2008, Pages 1007-1011.

55. T. Wei, Y.S. LiAn. Inverse boundary problem for one-dimensional heat equation with a multilayer domain. Engineering Analysis with Boundary Elements. Volume 33, Issue 2, February 2009, Pages 225 232.

56. Ziqi Sun. Inverse boundary value problems for a class of semilinear elliptic equations. Advances in Applied Mathematics. Volume 32, Issue 4, May 2004, Pages 791 800.

57. Основные публикации автора,диссертационной работы^

58. Кудинов В.А., Кудинов И.В. Аналитические методы теплопроводности //Монография. Самара: Самарский . государственный архитектурно-строительный университет. 2011. — 212 с.

59. Кудинов В.А., Кудинов И.В. Методы решения параболических и гиперболических уравнений теплопроводности //Монография.- М.: Книжный дом «Либроком», 2011. 270 с.

60. Актуальные вопросы тепло- и массообмена, энергоэффективность, исследование вихревых закрученных потоков». — 2008. — С. 39 — 46.

61. Стефанюк Е.В., Кудинов И.В. Построение аналитических решений уравнений! динамического и теплового пограничных слоев // Вестник Самарского государственного технического университета. Серия Технические науки. № 2 (24). - 2009. - С. 205 - 211.

62. Кудинов И.В. Использование компьютерной модели для проектирования' тепловых сетей //Вестник Самарского государственного технического университета. Серия Технические науки. № 4 (27). - 2010. -С. 174-181.

63. Кудинов И.В. Построение компьютерных моделей систем теплоснабжения больших городов //Вестник Самарского государственного технического университета. Серия Технические науки. — № 1 (29). — 2011.— С. 212 — 219.

64. Кудинов В.А., Кудинов И.В. Теоретические основы получения вихревых полей потенциалов // Известия РАН. Энергетика. — № 1. — 2011. — С. 105 112.

65. Кудинов В.А., Кудинов И.В. Теплообмен при течении Куэтта с учетом диссипации энергии // Известия высших учебных заведений и энергетических объединений СНГ. Энергетика. — 2011. — № 2. — С. 43 — 51.

66. Кудинов В.А., Кудинов И.В. Математическое моделирование процессов теплообмена и фазовых превращений с учетом абляции // Инженерно-физический журнал. Т. 84. 2011. - №»5. - С. 1065 - 1074.

67. Стефанюк E.B!, Кудинов И.В., ЛаргинаЕ.В. Обобщенные функции в задачах теплопроводности для многослойных конструкций // Вестник Самарского государственного* технического университета. Серия: Математическая. 2008'. - № 2 (8). - С. 41 - 56.

68. Кудинов И.В., ЛаргинаЕ.В., Еремин A.B., Будыльников М.Н. Теплообмен в плоском канале с учётом диссипации энергии // Вестник Самарского государственного технического университета. Серия: Математическая. -2009. №2 (10). -С. 38 - 47.

69. Кудинов И.В. Аналитические решения уравнений теплового пограничного слоя с учетом диссипации энергии // Повышение энергоэффективности зданий и сооружений. Межвузовский сборник научных трудов. Самара: СГАСУ. Вып. 5. - 2010. - С. 146 - 150.