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

КОСТЕНКО АНАТОЛИЙ ВАСИЛЬЕВИЧ

Численное моделирование и разработка комплекса

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

в регулярных продольнооребренных коридорных структурах

Специальность – 05.13 18 «Математическое моделирование, численные методы и комплексы программ.»

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

доктора технических наук

Комсомольск-на-Амуре - 2009

Работа выполнена в государственном образовательном учреждении высшего профессионального образования «Комсомольский-на-Амуре государственный технический университет» (ГОУВПО «КнАГТУ»)

Официальные оппоненты:

доктор физико-математических наук, профессор

Булгаков Виктор Кирсанович (г. Хабаровск)

доктор физико-математических наук, профессор

Карпов Александр Иванович (г. Ижевск)

доктор технических наук, профессор

Седельников Геннадий Дмитриевич (г.Комсомольск-на-Амуре)

Ведущая организация – Институт прикладной математики (ИПМ) ДВО РАН г.Владивосток

Защита состоится часов на заседании диссертационного совета Д 212. 092.03 при ГОУВПО, «КнАГТУ», по адресу: 681013, г. Комсомольск-на-Амуре, пр.Ленина, 27. Факс: (42172) 54 –08–87,E-mail: [email protected]

С диссертацией можно ознакомиться в библиотеке ГОУВПО «КнАГТУ».

Автореферат разослан «___» ___________ 2009 г.

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

Ученый секретарь

диссертационного совета,

ДМ 212.092.03

Кандидат физико-математических наук, профессор М.М. Зарубин

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Специалистам в области машиностроения, энергетики, теплотехники и теплоснабжения известно несоответствие уровней развития экспериментальных и численных методов решения задач, возникающих при проектировании теплообменных аппаратов (ТА) и их устройств. К сказанному можно привести факты. Во – первых, применяемые в настоящее время методы расчета ТА базируются, как правило, на результатах физических экспериментов; примером обобщения экспериментальных данных с помощью теории подобия и полуэмпирических методов, приводящих к простым критериальным зависимостям и соотношениям для коэффициентов сопротивления и теплоотдачи служат работы А.А. Жукаускаса, Р.В. Улинскаса, В.И. Катинаса, относящиеся к вопросам конвективного переноса, гидродинамики и вибрациям, особенностям теплообмена, сопротивления и теплоотдачи поперечно обтекаемых пучков труб в области низких и высоких чисел Рейнольдса; работы исследователей В.Ф. Юдина, Л.С. Тохтаровой, В.М. Антуфьева, Г.С. Белецкого посвящены аналогичным проблемам поперечно обтекаемым пучкам ребристых труб; позже работы И.Д. Лисейкина, В. А. Локшина, И.Д. Аронова, В.Н. Фоминой – мембранным конвективным поверхностям нагрева. Во – вторых, численное моделирование течений и теплообмена с помощью ЭВМ и компьютерных технологий, использование возможностей которых сулит повышение точности расчетов ТА, и, что более важно, переход к их автоматизированному проектированию, пока еще не нашло должного широкого применения при решении задач, имеющих практический интерес. Естественно, это объяснимо как сложностью самих рассматриваемых задач, так и тем, что большинство существующей литературы по численным методам рассчитано на вычислителей – математиков, и для специалистов в области проектирования техники вряд ли доступны без соответствующей подготовки. Исключением является крайне ограниченное число работ, которым следует отнести монографию П.Роуча «Вычислительная гидродинамика», хотя не в полной мере отвечает требованиям предъявляемым к расчету при проектировании ТА и их элементов (не все режимы течений рассматриваются, отсутствуют примеры конкретных реализаций и др.). В третьих, прогресс в развитии ЭВМ, компьютерных технологий и потребность доведения расчетных этапов проектирования ТА и их элементов до современного уровня стимулировала интерес к численному моделированию, но появление ряда статей по проблемам тепломассопереносу при вынужденной конвекции в журналах и обзорных материалах имеют недостатки, которые носят ущербный характер точности решения, игнорируя вопросы адекватности математической модели и эффективности расчетной процедуры и др.

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

Расчетные исследования ТА и их устройств, основанные на решениях уравнений Навье – Стокса и энергии, до сих пор носят единичный характер, причем первые из этих работ относятся к начальному периоду развития вычислительной гидродинамики и теплообмена (CFD), а обобщившая эти исследования энциклопедическая монография И.А. Белова и Н.А.Кудрявцева «Теплоотдача и сопротивление пакетов труб», опубликована почти двадцать лет назад. Появление и современное развитие ЭВМ и компьютерной техники может существенно повышать эффективность решения научных инженерных задач. Новые результаты достигаются при разработке программ, предназначенных для решения сложных задач,

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

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

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

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

пучков круговых продольно оребренных труб ТА поперечным потоком вязкой несжимаемой жидкости в ламинарном диапазоне изменения числа Рейнольдса говорят об актуальности темы исследования.

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

Целью работы является исследование с помощью численного моделирования неизотермического течения и теплообмена в регулярных коридорных прямолинейно оребренных (ассиметрично и симметрично) пучках труб кругового сечения поперечным потоком вязкой несжимаемой жидкостью в ламинарном диапазоен изменения числа Рейнольдса при постоянном значении числа Прандтля (Pr = 0,73) и разработка комплекса программ за счет :

- совершенствования методики реализации математических моделей течения и теплообмена;

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

- составления программных продуктов реализации задачи исследования;

- дискретизации расчетной области, ее разностной аппроксимации исходной системы уравнений и использовании вычислительного алгоритма расчета течения и теплообмена в коридорных регулярных гладкотрубных пучках;

- установления метода «среднего» Фурье эффективного получения параметров корреляционных зависимостей течения теплообмена по двум базовым;

- исследования сравнительной эффективности гладкотрубных и продольно оребренных пучков труб.

Для достижения поставленной цели работа содержит.

- Методы расчета течения теплообмена в регулярных структурах поперечно обтекаемых оребренных пучках труб кругового сечения и анализ результатов физического и численного экспериментов.

- Расчет течения теплообмена в регулярном коридорном пучке труб кругового сечения в ламинарном диапазоне изменения чисел Рейнольдса.

- Численное исследование стационарного обтекания регулярного коридорного пучка труб кругового сечения.

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

- Использование комплексного интегрального преобразования Фурье в исследовании течения и теплообмена и установление метода «среднего» Фурье в численном моделировании.

- Численное моделирование обтекания и теплообмена в продольно оребренных регулярных коридорных пучках труб поперечным потоком вязкой несжимаемой жидкости.

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

- Разработку комплекса программ численного исследования течения и теплообмена в регулярных коридорных структурах ТА и их элементов.

- Использование результатов численного исследования течения и теплообмена регулярных коридорных пучков труб кругового сечения и метода «среднего» Фурье в

практике исследований других наук.

- Приложение результатов исследований в практике.

Методы исследования. Основные результаты в работе получены с помощью реализации математических моделей течения и теплообмена (нелинейная система уравнений: неразрывности, Навье – Стокса и энергии) осуществлялось численными конечно-разностными итерационными с использованием метода контрольного объема (КО) и процедуры SIMPLE методами на ЭВМ. Ряд результатов получены с помощью анализа результатов численного моделирования и проведенных сравнительных данных, существующих в литературных источниках, других авторов, а также с использованием метода «среднего» Фурье. Экспериментальные данные подвергались статистической обработке с помощью программного продукта, разработанного автором «Корреляция Эйлер - Нуссельт», компьютерных программ Mathcad, Exel и др.

Достоверность результатов работы подтверждается :

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

- хорошим согласованием тестовых расчетов;

- использованием современных методов визуализации течений и теплообмена в виде картин течения функции тока и изотерм;

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

Научная новизна работы состоит в следующем:

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

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

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

- впервые показана перспектива использования продольного оребрения пучка труб для повышения тепловой эффективности ТА и их элементов;

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

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

- разработаны программы для ЭВМ «Коридор», «Трубные пучки», «Корреляция Эйлер - Нуссельт» и получены соответствующие регистрационные Свидетельства РФ;

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

- предложен метод оценки и получена формула расчета сравнительной эффективности трубных пучков кругового сечения;

- получен патент на изобретение «Ледяная переправа», который может быть использован в районах Севера и Дальнего Востока при сооружении мостовых переправ;

- предложено использование результатов численных расчетов задач фильтрации,

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

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

Научная и практическая ценность. Результаты, полученные в работе, имеют большое теоретическое и практическое значение. Можно констатировать, что решена крупная научно- практическая проблема неизотермического течения и теплообмена твердых тел поперечным потоком вязкой несжимаемой жидкостью в виде прямолинейно оребренных пучков труб кругового сечения. Они позволяют углубить и раскрыть представление о физической модели рассматриваемых течений и более четко описать ее определяющие механизмы. Разработанная методика численного моделирования течения и теплообмена позволяет определять тепловые и динамические характеристики ТА, ТЭУ и их элементов и других устройств и установить их оптимальные формы с целью достижения максимальной тепловой активности и минимального сопротивления при затратах одной и той же мощности на перемещение теплоносителя. Разработанные комплексы программ и алгоритмы позволяют решать также новые задачи вычислительного моделирования, связанные с разработкой и совершенствованием других структур, режимных течений, рабочих органов машиностроения, систем теплового управления потоками. Построенные комплексы программ позволяют использовать их при переходе к автоматизированному проектированию ТА, ТЭУ и другие. Алгоритмическая блок-сема решения задачи представляет собой самостоятельную научно-методическую ценность и внедрена в учебном процессе в КнАГТУ. Результаты нашли применение на «КнААПО» им Ю А Гагарина; ЗАО «Двльметаллургстрой»; ООО «Спецстройиндустрия»; МУП «Трамвайное управление»; ПО «Хлебная база № 59».

Автор защищает.

1.Оценку за счет выбора экстремально мелкой сетки качества полученной расчетной информации о течении и теплообмене в регулярном коридорном продольно оребренном пучке труб кругового сечения при использовании условий симметрии и периодичности в расчетной области.

2. Результаты расчета обтекания и теплообмена регулярного коридорного пучка труб кругового сечения поперечным потоком вязкой несжимаемой жидкостью.

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

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

5. Метод «среднего» Фурье получения корреляционных зависимостей течения и теплообмена.

6. Оценку сравнительной эффективности трубчатых поверхностей теплообмена кругового сечения.

7. Выводы о перспективности продольного оребрения на основе анализа сравнительной эффективности гладкотрубного и оребренного пучков.

8. Использование методов численного моделирования и их результатов в практической и научной деятельности.

Апробация результатов работы. Основные результаты диссертационной работы

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

1.На научно-технических конференциях Комсомольского на Амуре государственного технического университета (КнАГТУ) (г. Комсомольск на Амуре, 1988, 1989, 1990г.г.).

2. На научном семинаре по проблеме аэродинамике вихревых течений кафедры аэродинамики и динамики полета КНИЛ - 14 ОЛАГА (г. Ленинград,1988, 1989, 1990 г.г.).

3. На научном семинаре в лаборатории турбулентности ИТМО АН БССР (г. Минск, 1989, 1990г.г.).

4. На научном семинаре кафедры аэродинамики летательных аппаратов ЛИАП (г. Ленинград, 1990г.).

5. На научно-техническом семинаре кафедры аэродинамики КуАИ по вопросам механики жидкости и газов (г.Куйбышев, 1969, 1970 г.г.).

6. На математических конференциях в городах: Чебоксары (1967г.), Ярославле (1968г.), Горьком (1969г.).

7. На научном семинаре Института Кибернетики АН УССР (г. Киев, 1969г.).

8. На научно-технической конференции по проблемам механики сплошной среды (г. Комсомольск на Амуре, КнАГТУ, 1988г.).

9. На международном научно-техническом симпозиуме (г.Комсомольск на Амуре, КнАГТУ, 1995г.).

10. На международной научно-технической конференции ГОУВПО «КнАГТУ» (г. Комсомольск на Амуре, 2003, 2006 г. г.).

11. На международных научно-практических конференциях в области экологии и безопасности жизнедеятельности (г. Комсомольск на Амуре, 2007, 2008г.г.).

Ряд полученных расчетных результатов работы докладывались на профилирующих кафедрах КнАГТУ (1993 – 2008 г.г.) и семинарах Центра математического моделирования и информации КнАГТУ (1999 – 2008 г.г.).

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

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

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

Структура и объем работы. Диссертация состоит из оглавления, введения, 7 глав, заключения и выводов, списка 224 литературных источников и приложения. Диссертация изложена на 292 страницах, включая 75 рисунков и 19 таблиц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Г.С. Белецкий, В.Ф. Юдин, А.А. Жукаускас, Р,В. Улинскас, В.П. Исаченко, Я.Л. Полыновский, В.Г. Фастовский, Л.С. Тохтарова, А.И. Мицкевич, М.В. Кирпичев, Ю.В. Петровский, Д.Керн, А. Краус, В.А. Локшин, А.А. Гухман, В.Т. Буглаев, Л.А. Анисин, И.А. Белов, Н.А. Кудрявцев и др., которые сформулировали основные положения и разработали первые инженерные методики расчетов и установки исследований. Благодаря работам М.В. Кирпичева и его последователей М.В. Антуфьевым и В.Ф. Юдиным установлены критериальные зависимости коэффициентов сопротивления и теплопередачи для различных типов пучков (в том числе регулярных) и даны практические рекомендации для подсчета коэффициента теплоотдачи в условиях охлаждения и нагревания потоков газа. Экспериментальный материал обобщен и представлен в виде расчетных формул, графиков, номограмм, удобных для практического использования. Работы М.В. Кирпичева внесли также пионерские исследования по вропросам сравнительной эффективности теплообменных процессов пучках труб и развиты в работах В.М. Антуфьева, В,Ф. Юдина и др.

Существенно, что работы И.А. Белова, Н.А. Кудрявцева, С,А. Исаева дали новое, более перспективное направление численного моделирования процессов расчета течения и теплообмена поперечным потоком жидкости (газа) с помощью использования ЭВМ. Это направление сулит большую экономическую выгоду и переход к автоматизированным методам использования численных расчетов при проектировании ТА, ТЭУ и других их элементов.

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

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

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

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

Рассчитывается плоское ламинарное течение несжимаемой вязкой жидкости и конвективный теплообмен в регулярном коридорном продольно оребренном пучке круговых труб, расположенных с продольными и поперечными шагами l и h (в нашем случае l = h = S1 = S2). Диаметр цилиндра d принят за линейный масштаб. В коридорном пучке выделяется расчетная ячейка (модуль) ABCDEF, на проточных границах BC и ED которого ставим периодические граничные условия. На границах CD и BE (рис. 3.а) – условие симметрии; на омываемой стенке AF - условия прилипания. Решение тепловой задачи рассматривается отдельно от динамической с использованием ранее рассчитанных полей скоростей. Число Re = 40,100, 250 и 500, а число Pr = 0,73. Геометрические размеры принимаются равными l = 1,25; 1,50 и 2,0, а также другие (см. рис.3.а; 6.г.).

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

Характерной особенностью работы является адресность проведения обзора литературных источников в главах. Первая глава является основополагающей в критическом обзоре работ. Обзор работ по пучкам труб и их оребрению изначально восходит от работ Н. Гребера, Е. Шмидта, С.Богерта и Р. Майера, Е. Нейля, О.Е. Власова, Л.Н. Ильина, М.А. Стриковича,

Харпера и Броуна, Т. Боша, Н. Деча, М. Ланга, благодаря чему можно продвинуться в решении задач теплообмена вообще и, в частности, оребренных трубчатых поверхностей при наличии удовлетворительных экспериментальных данных по коэффициентам теплоотдачи ребер и несущих их поверхностей нагрева. Отмечаются работы по повышению коэффициентов теплоотдачи и методам интенсификации тепловой активности, к которым следует отнести работы А.Е. Берглиса, Р.Л. Веба и других зарубежных авторов. Большой вклад отечественных исследователей: А.М. Антуфьева, В.Ф. Юдина, М.В. Кирпичева, А.И. Мицкевича и других. В развитие экспериментальных методов исследования течения и теплообмена привело к подходу сравнительной оценки поверхностей нагрева, при этом можно ее упрастить, что предлагается в данной работе, имея корреляционные зависимости чисел Рейнольдса, Эйлера и Нуссельта. Численные исследования пучков труб позволили сделать вывод о том, что экспериментальные данные, полученные названными авторами достаточно хорошо согласуются с расчетными, полученными с помощью численного моделирования. Следовательно, актуальность проблемы численного моделирования течения и теплообмена продольно оребренных пучков труб можно решать в плоской постановке, пренебрегая концевыми эффектами. Из всего отмеченного следует, что выбранный метод исследования течения и теплообмена, который основан на численном решении полных уравнений Навье – Стокса и энергии, содержание которого предложено П. Роучем и И.А. Беловым, А.С. Гиневским, Н.А. Кудрявцевым, С. Патанкаром и С.А. Исаевым, Е.Спэрроу является эффективным с версией реализации периодического характера течения и теплообмена. Эта методология получила развитие в приложениях к проектированию трубчатых ТА, контейнерному транспотру, аналитическому приборостроению, терморегулированию и экологическим вопросам.

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

, (1)

где - вектор скорости (u, v); Гф - коэффициент диффузии; t – время; Ф – зависимая переменная; Sф – источниковый член. Уравнение (1) записано в нестационарном виде и включает в себя четыре члена: нестационарный, конвективный, диффузионный и источниковый.

Учитывая наложенные ограничения на характер течения, (1) запишем в виде:

, (2)

и включает три члена: конвективный, диффузионный и источниковый. Переменная Ф может принимать значения: 1; u или v, Т, где u и v – составляющие вектора скорости; Т – температура. Конкретный вид Гф и Sф зависит как от Ф, так и от выбранной системы координат

, (3)

где r и v - радиальная и тангенциальная координаты; u и v - радиальная и тангенциальная составляющие скорости соответственно.

Уравнение (2) используется для расчета течения и теплообмена как в ламинарном, так и в турбулентном режимах течения, используя дополнительные режимные условия. Следовательно, в случае ламинарного режима течения переменная Ф обозначает составляющие скорости (u, v) и температуру Т.

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

Ф= 1; ГФ = Sф = 0 (4)

- уравнение неразрывности;

Ф=u: Гu=1/Re; Su=- (5)

- уравнение изменения количества движения в радиальном направлении;

Ф = v: Гv = 1/Re; Sv= (6)

- уравнение изменения количества движения в тангенциальном направлении;

Ф = Т: ГТ = 1/ Re Pr ; SТ = 0 (7)

- уравнение энергии. Здесь Re = u d/ v - характерное число Рейнольдса; Pr = v/a - число Прандтля; u - среднемассовая скорость в минимальном проходном сечении пучка АА!

(см. рис.1.а); d – диаметр трубы; v - коэффициент кинематической вязкости; a - коэффициент температуропроводности. Система (3) - (7) может быть использована для расчета как стационарных, так и нестационарных течений в ламинарном диапазоне изменения числа Рейнольдса.

Обтекание пучка продольно оребренных труб кругового сечения поперечным потоком жидкости (газа) характеризуется неравномерностью распределения скорости, давления, температуры, локальных характеристик течения и теплообмена, сложностью геометрии расчетной области (см. рис.6.). При расчете таких течений С.А. Исаев, П.А. Баранов, Н.А. Кудрявцев, Т.А. Баранова в работе «Численное моделирование влияния чисел Рейнольдса» и

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

Отметим, что свойства метода КО и детализация решения описаны в работах С.Патанкара, П.Роуча и И.А. Белова. Важным достоинством которого является обеспечение точного

интегрального выполнения физических законов для любой группы КО, следовательно, и для всей расчетной области. Разностная схема является консервативной по массе, количеству движения и энергии. Важным является и принцип построения сетки – размещение узлов сетки в расчетной области. От этого зависит точность информации и ее физическая достоверность. Используем сетку, в которой тангенциальная и радиальная составляющие скорости хранятся в узлах, смещенных относительно узлов хранения давления и температуры на полшага (см. рис.7.а, б.). Значения переменной Фi можно определить, используя различные аппроксимационные схемы. В работе используется схема Леонарда против потока второго порядка точности. Численные эксперименты показывают, что схема Леонарда второго порядка точности близка по точности получения интегральных характеристик течения и теплообмена, а также по обеспечению устойчивости процесса численного решения. Рассматривая более общую задачу нестационарного обтекания и теплообмена, для способа аппроксимации нестационарного члена уравнений (3) – (7) можно на подобии уравнения (1) записать в следующем виде

, (8)

где F – стационарная часть, включающая в себя конвективные, диффузионные и источниковые члены. Практика расчетов показывает, что для обеспечения устойчивости и сходимости процесса численного моделирования выбирают неявные схемы, одну из которых рекомендует И.А. Белов – схему Перйе второго порядка точности по времени:

, (9)

которая подавляет нефизические пульсации параметров течения и достаточно проста в алгоритмизации. Для обеспечения записи (9) стационарной части F на (N + 1) - м временном слое применяется итерационная процедура.

В общем виде разностный аналог конвективно-диффузионного уравнения представляется следующим образом (см. рис.5 и 7)

, (10)

где Ф = u, v, Т; aN, aE, … - коэффициенты, определяемые способом нахождения конвективных и диффузионных потоков через грани КО; bФ - источниковая часть разностного уравнения. Учитывая факторы стабилизации окончательно получим, что

, (11)

где ; - соответственно диффузионные и конвективные коэффициенты; vol Ф - объем КО; - шаг по времени;

, (12)

Алгоритм решения конвективно-диффузионного уравнения (11) основан на LU - факторизации. Расчет давления основан на решении уравнения Пуассона для поправки давления, найдя которую можно определить и само давление, используя при этом процедуру SIMPLE.

Рисунок 7 иллюстрирует способ определения интерполяционных коэффициентов при заданных граничных условиях. Эти коэффициенты установлены автором в работе «Использование интерполяционных соотношений при нахождении параметров течения и теплообмена в расчетной области » можно также применять при разработке комплексов программ для построения многоблочного факторизованного алгоритма решения уравнений Навье – Стокса и энергии на базе структурированных пересекающихся сеток

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

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

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

Анализ расчетных результатов также позволяет сделать вывод о том, что с уменьшением плотности регулярных коридорных пучков труб ведет к снижению уровня силовых нагрузок на глубинную трубу, продольная составляющая которых вносит основной вклад в сопротивление. Приводятся картины течения и температурного поля для компоновки 1,50 х 1,50 и Re = 40; 500, а также графики эволюции с ростом числа Рейнольдса поверхностного распределения давления. Результаты численных экспериментов и полученные с помощью метода «среднего» Фурье для регулярных пучков 2,25 х 2,25 и 2,50 х 2,25, говорят о том, что уменьшение плотности пучков труб за пределы компоновки 2,0 х 2,0 ведет к автомодельному процессу течения во всем диапазоне изменения чисел Re = 40…500 и величина Eu стабилизируется и находится в пределах 0,2 – 0,1 (рис.10). Представленные кривые на графиках распределения местного числа Нуссельта, напряжение трения, давления – топологически подобны для ламинарного режима течения (рис.11, рис. 12).

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

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

коридорных пучках труб кругового сечения и их анализ, рассматривая мембраны как прямолинейные ребра – проставки.

Начиная с семидесятых годов ХХ столетия потребности государства в конструкторских и технологических разработках по промышленному освоению новых МКПН для парогенераторов поставили задачу исследований по расчету аэродинамического сопротивления и теплообмена. Большой вклад в экспериментальных исследованиях принадлежит Локшину В.А., Лисейкину И.Д., Аронову Д.И., Ленькову Ю.А., Назаренко В.С., Моргун А.В., Мигай В.К., Быстрову П.Г., Сотникову И.А.

Анализ результатов выявил расхождения при сравнении опытных данных с расчетными, что побудило автора работы к численным исследованиям для выработки рекомендаций. Следует обратить внимание на факт отсутствия сведений в литературных источниках по численному моделированию расчетов течения и теплообмена в МКПН, за исключением работы Костенко А.В. по использованию алгебраической блок-схемы при расчете аэродинамического сопротивления и теплоотдачи мембранного коридорного пакета труб. В данной главе все расчеты выполнены с использованием q – расчетной области (см. рис. 3). Проведен анализ расчетной информации, полученной по программе комплекса «Коридор» и методу «среднего» Фурье, разработанных автором. Приводятся картины течения и теплообмена, графики корреляционных зависимостей (см. рис.13 ) и др., а также использование результатов расчетов для практического применения их в практике проектирования ТА.

На рисунке 18 приведены изображения регулярного коридорного пучка труб, схема компоновки пучка и картина течения в виде линий постоянных значений функции тока при числе Re =40.

Среднее по периметру трубы пучка 2,0 х 2,0 значение Num при всех режимах течения сведены в таблицу 1. Расчет значений Eu – динамической характеристике этого регулярного коридорного пучка труб в таблице 2. Расчетные данные сравнительной эффективности гладкотрубных и оребренных пучков, как пример для рассмотренного случая пучка 2,0 х 2,0 и Pr = 0.73 (Re =500), приведены в таблице 3.

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

В пятой главе представлены результаты исследований по применению рядов и комплексного интегрального преобразования Фурье с конечными пределами, которое базируется на работах Г.Е. Пухова «Комплексное исчисление и его применение» и автора настоящей работы, по свойствам, относящимся к комплексному исчислению как к предмету по интегральным преобразованиям. Эти результаты использованы в практике исследования течения и теплообмена регулярных структур ТА и их устройств для получения корреляционных зависимостей по сопротивлению и теплоотдаче и другим интегральным и локальным характеристикам. Метод названный методом «среднего» Фурье, так как доказательство его основано на достаточном признаке сходимости ряда Фурье (что представляет собой обратное преобразование Фурье).

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

использование метода важный для практики проектирования ТА случай выбора корреляционной зависимости в переходном режиме течения, числа Нуссельта от числа Рейнольдса. Кривая 3 получен по методу «среднего» Фурье.

Для построения корреляционных зависимостей и построения их графиков автором создан вычислительный программный продукт «Корреляция Эйлер - Нуссельт».

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

Интерес представляют картины течения и теплообмена, графики распределения давления и трения, местного числа Нуссельта и другие.

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

Получены формулы

н = 2с - в, (13)

в = 2с - н. (14)

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

Глава содержит общие выводы и рекомендации для конструкторов-проектировщиков ТА и других устройств. Предложенные методы корреляционных зависимостей (13) и (14) могут пополнить по новым сведениям по теплообмену «Справочник по теплообменникам »: в 2 т. под редакцией Б.С. Петухова, В.К. Шикова. – М.: Энергоатомиздат, 1987.

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

= Num / Eu, (15)

где p и q – вещественные числа. Кроме (15) возможны и другие более сложные функциональные зависимости. Этому вопросу и его освещению посвящены результаты работ А.И. Мицкевича, В.М. Антуфьева, В.Ф. Юдина, М.В. Кирпичева. Здесь ограничимся рассмотрением (15), предложенном автором данной работы для проведения предварительного качественного анализа эффективности продольного оребрения. Числа p и q должны зависеть от ряда факторов. Учитывая, что все предыдущие расчеты проведены при постоянном значении числа Прандтля Pr = 0,73 (отсюда следует подобие теплового и динамического пограничных слоев, что отражено в трудах Л.И. Седова, Б.С. Петухова, Л.Г. Генина, Л.Г. Лойцянского) примем p=q=1 зависимость (15) примет вид

= Num / Eu, (16)

При решении проблемы сравнительной оценки поверхности нагрева М.В. Кирпичев предложил использовать для оценки так называемый энергетический коэффициент. В настоящей работе предлагается делать заключение об эффективности поверхности нагрева с помощью геометрического смысла, построив прежде корреляционные зависимости для сравниваемых поверхностей. Для заданного режима течения и данной компоновки труб по графикам Nu = f (Re) найти отношение площадей, которое даст ответ о наивыгоднейшей поверхности нагрева. Так имея зависимости Nu = S1 (Re) и Nu = S2 (Re), коэффициент эффективности трубного пучка определим соотношением :

= S1 (Re) / S2 (Re), (17)

где S1 и S2 – площади фигур, ограниченные зависимостями S1 (Re) и S2 (Re).

Таким образом, в общем случае

=d Re, (18)

где Re 0 и Re k - начальное и конечное значение числа Рейнольдса рассматриваемого режима течения.

Проведенные расчетные данные для коэффициента для гладкотрубного и продольно оребренных пучков труб собраны в таблицу. Из данных таблицы сравнительной эффективности пучка 2,0 х 2,0 при Pr = 0,73 видно, что во всем диапазоне расчетных чисел Рейнольдса можно подобрать вариант оребрения пучка труб, при котором он будет эффективен, чем гладкотрубный пучок. В частности, при Re = 40 … 250 симметрично оребренный пучок с четырьмя ребрами дает выигрыш в эффективности в сравнении с гладкотрубным в среднем на 30%. В тоже время при числе Re = 500 наиболее эффективном является вариант ассимметричного оребрения (выигрыш 14%); а продольно оребренный пучок с восьмью ребрами дает выигрыш 50%, а дальнейшее оребрение приводит к падению по эффективности (начиная с варианта десяти ребер).

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

На рис. 16.б – представлены распределения коэффициентов давления pw на развитой поверхности для двух вариантов оребрения. Для случая четырех ребер наблюдается резкое падение давления по достижении абсолютного максимума, который с ростом Re смещается к передней критической точке. В тоже время для второго варианта оребрения характерно резкое падение давления в местах расположения ребер на поверхности трубы, причем видно, что месторасположение абсолютного максимума давления pw зависит от числа Re. Зависимости коэффициента трения связаны с коэффициентом давления и геометрически.

На рис. 17.а – б - показано распределение числа Нуссельта для обоих вариантов оребрения. Отмечается, что наибольшее различие теплоотдачи оребренного и гладкотрубного пучков связано с наличием при оребрении большого числа локальных максимумов теплоотдачи (по числу ребер). Симметричное продольное оребрение приводит к более равномерной теплоотдаче поверхности трубы. В этой главе представлены также интегральные характеристики течения и теплообмена для продольного оребрения для обоих вариантов симметричного оребрения.

В итоге делаются выводы, из которых следует отметить, что:

- продольное оребрение является перспективным способом увеличения тепловой активности ТА при неизменной мощности на перемещение теплоносителя;

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

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

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

- следствием асимметрии распределения коэффициентов давления и трения на поверхности несущей трубы является возникновение помимо лобового сопротивления и поперечной нагрузки - подъемной силы, о чем отмечено в работе А.И. Олейникова и А.В. Костенко;

- даны рекомендации для практики проектирования ТА использовать асимметричное и отказ от симметричного оребрения;

- вычислительные комплексы программ «Коридор», «Трубные пучки» использовать для перехода к автоматизированному проектированию ТА и их элементов.

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

На рисунке 21 приведена процедура построения картин течения переменной, в виде этапов и систематизации точек изолинии = 0, отличающейся от предложенной П. Роучем в его книге «Вычислительная гидродинамика», введя понятие нормы сеточной ячейки.

ВЫВОДЫ

В результате исследований, проведенных в нестоящей работе:

1. Разработанный вычислительный алгоритм и комплекс программ «Коридор», «Трубные пучки», «Корреляция Эйлер-Нуссельт»для исследования и расчета течения и теплообмена в регулярных коридорных структурах пучков труб, позволяет эффективно получать достоверную и надежную информацию по интегральным и локальным характеристикам течения и теплообмена во всем расчетном диапазоне изменения числа Рейнольдса

(Re = 40…500).

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

3. Использование результатов течения и теплообмена в исследованных пучках позволило придти к методу «среднего» Фурье, позволяющем экономно и эффективно получать расчетную информацию и в других науках.

4. Исследование течения и теплообмена с помощью метода «среднего» Фурье в мембранных пучках труб показало хорошее согласование расчетных характеристик с соответствующими экспериментальными данными других авторов.

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

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

7. Выявлено, что следствием асимметричного оребрения является возникновение помимо лобового сопротивления и поперечной нагрузки – подъемной силы.

8. Предложены практические выводы о возможностях интенсификации теплообменных процессов в регулярных структурах ТА за счет асимметричного оребрения, позволяющего достичь значения теплосъема на 30% больше.

9. Предложен метод для установления коэффициента сравнительной эффективности пучков труб и получена простая в употреблении формула расчета.

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

11. Использовать результаты численных экспериментов при автоматизированном проектировании ТА и их устройств.

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

Рис.1. Схема расположения труб в пучке и их конфигурация: а- коридорное гладкотрубное; б – шахматное гладкотрубное; в – однорядный гладкотрубный пучок; е - мембранный коридорный пучок; ж - оребренный коридорный пучок; з – схема изображения гладкотрубного коридорного пучка; и – шахматное изображение пучка гладких труб; к- схема мембранного пучка; л – схема коридорного оребренного пучка.

Рис.2. Продольное оребрение: Рис. 3. Расчетные области:

а- симметричное; б- ассиметриченое; а-«q» расчетная область;

в- кольцевое. б – «Q» расчетная область

и полярная сетка на ней.

Рис. 4. Виды сечений и оребренные трубы: Рис.5. Размещение узлов: а, б – в декартовых

а – поперечное; б – поперечное и координатах; в – элемент KO, охватывающ

продольное сечения; в – г – виды оребре- чий узел P; г – узлы в полярной системе;

ния труб. ( ; х) – узлы хранения; – скалярных

переменных.

Рис.6. Изображение коридорного гладкотрубного пучка круговых труб и расчетной ячейки:

а – АВСД расчетная ячейка ( модуль) корудорного гладкотрубного кругового пучка;

б - расчетная область (ячейка) гладкотрубного кругового пучка;

в – расчетная область (ячейка или модуль) гладкотрубного кругового пучка в другой конфигурации; г – схема расчетной области (расчетного модуля); s1 и s2 – относительные поперечный и продольный шаги; d – диаметр трубы.

Рис. 7. Определения интерполяционных коэффициентов при заданных граничных условиях;

а – пристеночный КО для радиальной составляющей скорости; б – для всех остальных переменных; в – расчет интерполяционных коэффициентов.

Рис. 8. Картины течения с симметрично расположенным продольным оребрением регулярных коридорных пучков труб 2,0 х 2,0 в виде линий постоянных значений функции тока : а – вариант четырех ребер; б- вариант восьми ребер.

Рис. 9. Блок-схема решения нестационарной задачи в ламинарном диапазоне изменения числа Рейнольдса.

Рис. 10. Численные результаты расчета теплоотдачи трубы регулярных коридорных пучков при числе Pr = 0,73.

Рис. 11. Топологически подобные профили напряжения трения в регулярных коридорных пучках труб в ламинарном режиме изменения числа Рейнольдса. Развитие с ростом числа Re поверхностного трения : 1 – Re=40; 2 - Re=100; 3 - Re=250; 4 – Re=500.

Рис. 12. Эволюция с ростом числа Re поверхностного распределения давлениядля регулярных коридорных пучков труб различной плотности: 1,25 х 1,25; а – плотные пучки труб; б – свободные пучки труб

Рис. 13. Влияние числа Рейнольдса на число Эйлера для мембранных и ассиметрично оребренных регулярных коридорных пучков труб. Пучок : 1,25 х 1,25 – верхняя полоса; 1,50 х 1,50 – средняя полоса; 1.38 х 1,38 – верхняя полоса; 2,0 х 2,0 – полоса 5; 2,0 х 2,0 – нижняя полоса (4 и 8 ребер)

Рис. 14. Схема «среднего» Фурье: а – точки сходимости корреляционной зависимости

= F (Re) 1; 2 – базовые зависимости ; 3 – искомая кривая; б – число Нуссельта в зависимости от числа Рейнольдса в переходном режиме течения; - 3 – кривая, полученная по методу «среднего» Фурье

Рис. 15. Ассиметричное оребрение: а, б – картины течения ассиметрично и продольно оребренной трубы пучка 2,0 х 2,0 в виде линий постоянных значений функции тока и линий постоянных значений температуры при числе Re=500; в - распределение чмсла Nu / Num на поверхности ассиметрично и продрольно оребренной трубы пучка 2,0 х 2,0.

Рис. 16. Симметричное продольное оребрение : а, б, в, г – распределение коэффициентов трения и давления (случаи четырех и восьми ребер) соответственно для пучка труб 2,0 х 2,0

Рис. 17. Распределение среднего числа Нуссельта для симметрично оребренного регулярного коридорного пучка труб 1,88 х 1,88 : а- случай четырех ребер; : б- случай восьми ребер.

Рис 18. Регулярный коридорный мембранный пучок труб.

а – изображение пучка; б – схема компоновки пучка; в – картина течения в виде линий постоянных значений функции тока для пучка 2,0 х 2,0 и Re = 40

Рис. 19. Распределение местного числа Нуссельта на поверхности трубы мембранного регулярного коридорного пучка: а – эксперимент (Жукаускас, Re = 57); – настоящий расчет (Re = 40); –настоящий расчет ( Re= 500); - эксперимент (Улинскас, Re = 502); - эксперимент (Жукаускас, Улинскас Re = 500); б – тоже, что и в а, везде Pr =0,73

Рис. 20. Влияние числа Re на число Эйлера для регулярных коридорных гладкотрубных и продольнооребренных мембранных пучков труб. Обозначения : - настоящий расчет; – эксперимент (Локшин В.А., Лисейкин И.Д. ).

Рис. 21. Процедура построения картин течения для переменной.

Расчетные данные Num и Eu

Для пучка 2,0 х 2,0 и Pr = 0.73

Таблица 1

Re	40	100	250	500
Num	4,49	4,58	4,66	6,43

Таблица 2

Re	40	100	250	500
Eu	- 0,76	- 0,34	- 0,20	- 0,14

Расчетные данные сравнительной эффективности гладкотрубных и оребренных пучков труб

( 2,0 x 2,0, Pr =0,73)

Таблица 3

Re		2	4	8
500	46,0	52,5	50,2	23,8

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ ОПУБЛИКОВАНО В СТАТЬЯХ

1. Костенко А.В. Построение корреляционных зависимостей с использованием рядов и комплексного преобразования Фурье в исследовании течения и теплообмена в круговых пучках труб //Проблемы машиностроения и автоматизации. Международный журнал, № 2, - 2007. – С. 103-107.

2. Костенко А.В., Гусаков А.А. О теореме запаздывания в комплексном исчислении // В сб. Вопросы механики жидкости и газов. Тр.,вып.2. Куйбышев, 1969.- С. 57-61.

3. Костенко А.В. Некоторые замечания о порядке убывания комплексных изображений// В сб. «Механика», Куйбышев, 1969. –С. 62-67.

4. Костенко А.В. О некоторых теоремах комплексного исчисления // В сб. «Волжский математический сборник», вып. 7. Куйбышев, 1969. –С.57-61.

5. Костенко А.В. Нахождение комплексного изображения от произведения функции на факториальную // В сб. «Методы математического моделирования и теория электрических цепей». Вып.3,АН УССР. Киев,1969. –С.89-91.

6. Костенко А.В. Обобщенная теорема о комплексном изображении произведения двух функций // В сб. «Математика». Куйбышев,1970. –С.224-226.

7. Костенко А.В. Численное моделирование обтекания и теплообмена коридорного пакета труб с ассиметричным продольным оребрением. Деп. в ВИНИТИ, 15.01.92, 161-В92, 1992.-С.13.

8. Костенко А.В. Использование интерполяционных соотношений при нахождении параметров течения и теплообмена в расчетной области. Деп. в ВИНИТИ, 07.07.1988,№ 5476-В88.-М.:1988.-С.12.

9. Костенко А.В. Численное моделирование течения и теплообмена коридорных пакетов труб с симметричным продольным оребрением. Деп. в ВИНИТИ 24.01.90, № 485 – В 90. –М.: 1990. –С. 12.

10 Костенко А.В., Костенко В.А. Расчет обтекания и теплоотдачи коридорного пакета труб с симметрично расположенным продольным оребрением // Механика строительных конструкций из новых материалов и проблемы практического внедрения в производство: Материалы межд. научн. - техн. симпозиума, Комсомольск на Амуре: КнАГТУ, 1995.-С.85-88.

11. Костенко А.В., Костенко В.А. Расчет обтекания и теплоотдачи коридорного пакета труб с ассиметрично расположенным продольным оребрением // Механика строительных конструкций из новых материалов и проблемы практического внедрения в производство: Материалы межд. научн.- техн. симпозиума. Комсомольск на Амуре, КнАГТУ, 1995,-С.163-165.

12. Костенко А,В., Усольцев Ю.А., Черепюк И.Д. Численное моделирование течения вблизи поверхности обтекаемого тела симметричной формы // Проблемы механики сплошной среды. Ч.1: Материалы тр. научн.-техн. конф. Комсомольск на Амуре: КнАГТУ, 1988, -С. 163-165.

13. Костенко А.В. Численное моделирование обтекания и теплообмена в оребренных коридорных пучках труб поперечным потоком жидкости // Проблемы машиностроения и автоматизации. Международный журнал. №3, 2006. –С. 60-67.

14. Костенко А.в. Обнаружение подъемной силы течения в коридорном пакете труб с ассиметрично расположенным оребрением // Нелинейная динамика и прикладная

синергетика: Ч.1: Материалы тр. межд. научно-техн. конф. Комсомольск на Амуре: ГОУВПО «КнАГТУ», 2003.-С.88-81.

15. Костенко А.В. Сопротивление и теплоотдача коридорного пакета труб в ламинарном диапазоне числа Рейнольдса (результаты численного эксперимента). Автореферат дис. на соиск. ученой степ. к.т.н. АН Белоруссии. АНК ИТМО им. А.В. Лыкова. – Минск.1992,-С.21.

16. Костенко А.В. Об автомодельности процесса течения и теплообмена в ламинарном режиме изменения числа Рейнольдса // Повышение эффективности инвестиционной и инновационной деятельности в ДВ регионе и странах АТР:Ч.2: Материалы тр. межд. научн.-техн. конф. Комсомольск на Амуре :ГОУВПО «КнАГТУ», 2006. –С.143-145.

17. Костенко А.В., Костенко В,А., Черепюк И.Д, Расчет течения и теплообмена коридорного пакета труб в проблеме сооружений переправ через водные преграды // Проблемы механики сплошной среды: Ч.1: Материалы тр. межд. научн.-техн. конф. Комсомольск на Амуре: КнАГТУ, 1998. –С.163-164.

18. Патент изобретения № 2260648 РФ, «Ледяная переправа», зарегистр. в Гост.реестре изобр. РФ 20.09.2005; Заявлено 07.04.2003, заявка №2003109921. г. Москва.

19. Свидетельство об офиц. регистр. Программы для ЭВМ № 2003611264 «Коридор» / А.В. Костенко, заявка № 2003610709,03.04.2003. Зарегистр. в реестре программ для ЭВМ 28.05.2003, г. Москва.

20. Свидетельство об офиц. регистр. программы для ЭВМ № 2006613241 «Трубные пучки» / А.В. Костенко, заявка № 2006611329, 24.04.2006. Зарегистр. в реестре программ для ЭВМ 14.09.2006.г. Москва.

21. Свидетельство об офиц. регистр. программы для ЭВМ № 2007612330.

«Корреляция Эйлер-Нуссельт» / А.В.Костенко, заявка № 2007611491, 20.04.2007. Зарегистр. в реестре программ для ЭВМ 04.06.2007.г. Москва.

22. Костенко А.В. Расчет обтекания и теплоотдачи коридорного пакета труб с симметрично расположенным оребрением // Повышение эффективности инвестиционной и инновационной деятельности в ДВ регионе и странах АТР: Ч.2: Материалы межд. научн.-техн. конф. Комсомольск на Амуре: ГОУВПО «КнАГТУ». 2006. –С.135-142.

23. Костенко А.В., Воротников С.М. Элементы теории аналитических функций и операционного исчисления в задачах и примерах. Комсомольск на Амуре, 2003, -С. 256.

24. Костенко А.В., Воротников С М. Элементы теории аналитических функций и операционного исчисления в задачах и примерах. Комсомольск на Амуре, 2004 –С.274.

25. Костенко А.В. Использование алгоритмической блок-схемы при расчете аэродинамического сопротивления и теплоотдачи мембранного коридорного пакета труб. Деп. в ВИНИТИ 29.02.88, № 1628-В 88. –М.: 1988.-С.32

26. Костенко А.В. К набору экспериментальных данных в процессе обтекания и теплообмена труб с использованием критерия подобия // Повышение эффективности инвестиционной и инновационной деятельности в ДВ регионе и странах АТР: Ч.2: Материалы межд. научн.-практ. конф. Комсомольск на Амуре : ГОУВПО «КнАГТУ», 2006. –С.131-134.

27. Костенко А.В. Использование комплексов программ численного моделирования течения и теплообмена в регулярных коридорных пучках труб в задачах фильтрации // В сб. «Дальневосточная весна -2007»: Материалы межд. научн.-практ. конф. в области экологии и безопасности жизнедеятельности: Комсомольск на Амуре, 7-8 июня 2007.-С. 184-187.

28. Костенко А.В. Использование результатов численного моделирования течения и теплообмена в регулярных коридорных пучках труб в задачах фильтрации и теплопотерь органами дыхания // В сб. «Дальневосточная весна - 2008» : Материалы межд. научно-практической конф. в области экологии и безопасности жизнедеятельности: Комсомольск на Амуре, 2008 –С.203 - 205.

29. Олейников А.И., Костенко А.В. Обнаружение подъемной силы в задаче обтекания и теплообмена ассиметрично и продольно оребренных регулярных пучков труб поперечным потоком жидкости с помощью численного моделирования // Проблемы машиностроения и автоматизации. Международный журнал № 1, 2008. – С. 105 – 110.

30. Костенко А.В. К вопросу оценки сравнительной эффективности гладкотрубных и продольно оребренных коридорных пучков труб // Проблемы машиностроения и автоматизации. Международный журнал. №4, 2008. – С 103 - 104.

31. Костенко А.В. Алгоритмы и программы численного моделирования неизотермического течения в элементах теплообменников // Информатика и системы управления. - 2008. - № 1 (15). – С. 46 -55.

32. Костенко А.В. Использование программ численного моделирования мембранных коридорных теплообменных учстройств // Информатика и системы управления.- 2008. № 3 (17). – С. 40 – 46.

33. Костенко А.В. О методике коррекции градиента среднемассовой температуры теплообмена коридорных пучков труб в разработке программ компьютерных технологий //Информатика и системы управления. – 2008. - № 4 (18), - С.33 – 39.

34. Костенко А.В. Результаты численного моделирования течения и теплообмена в мембранных регулярных коридорных пучках труб // Вестник КнАГТУ. Вып. 12 Часть 1. - Комсомольск-на-Амуре. ГОУВПО «КнАГТУ», 2008. – 185 С.

35. Костенко А.В., Купченко А.В. Применение результатов численного моделирования течения и теплообмена мембранных пучков труб при построении устройств для использования газов при переработке нефти// В сб. «Дальневосточная весна - 2008». Материалы межд. научн.-практ. конф. в области экологии и безопасности жизнедеятельности: Комсомольскна Амуре. 2008. – С. 242-243.

36. Костенко А.В., Зигун К.В. Применение результатов численного моделирования к построению устройств фильтрации водных массивов// В сб. «Дальневосточная весна - 2008». Материалы межд. научн.-практ. конф. в области экологии и безопасности жизнедеятельности: Комсомольск на Амуре. 2008. – С. 243 – 244.

37. Костенко А.В. Численное моделирование течения и теплообмена в мембранных коридорных пучках труб поперечным потоком жидкости // Проблемы машиностроения и автоматизации. Международный журнал № 1, 2009. – С 95 – 99.