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

Алхасова, Джамиля Алибековна

Исследование и гидродинамические расчеты внутрискважинных теплообменников с продольными ребрами тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

Алхасова, Джамиля Алибековна АВТОР

кандидата технических наук УЧЕНАЯ СТЕПЕНЬ

Махачкала МЕСТО ЗАЩИТЫ

2009 ГОД ЗАЩИТЫ

01.04.14 КОД ВАК РФ

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

Автореферат диссертации на тему "Исследование и гидродинамические расчеты внутрискважинных теплообменников с продольными ребрами"

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

ООЗАЬЬ

Алхасова Джамиля Алибековна

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

Специальность: 01.04.14. - теплофизика и теоретическая теплотехника

АВТОРЕФЕРАТ

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

* «

Махачкала -2009

003466732

Диссертационная работа выполнена в Учреждении Российской академии наук Институте проблем геотермии Дагестанского НЦ РАН

Научный руководитель:

доктор технических наук, профессор Алишаев М.Г.

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

доктор технических наук, профессор Ильин А.К.

кандидат технических наук Джаватов Д. К.

Ведущая организация: Государственное образовательное учреждение высшего профессионального образования Московский государственный университет инженерной экологии

Защита состоится «__£» мая 2009 г. ъ/£££ч&соъ на заседании диссертационного совета ДМ 002.071.01 при Учреждении Российской академии наук Институте проблем геотермии Дагестанского НЦ РАН по адресу: 367030, г. Махачкала, пр.И.Шамиля, д. 39а.

С диссертацией можно ознакомиться в библиотеке Учреждения Российской академии наук Института проблем геотермии Дагестанского НЦ РАН.

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

Ученый секретарь диссертационного совета ДМ 002.071.01 доктор технических наук

¿А

БазаевА.Р.

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

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

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

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

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

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

В соответствии с целевым направлением работы задачами исследования явились:

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

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

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

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

5. Проведение оптимизационных расчетов по выбору конструктивных параметров оребренных скважинных теплообменников;

6. Представление рекомендаций по выбору конструктивных параметров скважинных теплообменников с продольным оребрением поверхности трубы.

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

1. Предложены новые методы расчета конструкций внутрисква-жинного теплообменника с продольным оребрением;

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

3. Предложен метод гидродинамического расчета перепада давления для турбулентного потока жидкости по кольцевому сечению с

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

Основные положения, выносимые на защиту.

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

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

3. Метод расчета перепада давления при движении потока жидкости по кольцевому сечению с ребрами. Зависимости перепада давления от высоты ребер и их числа для турбулентных течений. Случаи других сечений.

4. Рекомендации по выбору конструктивных параметров скважинных теплообменников для промысловых и лабораторных условий.

5. Расчет теплопередачи через ребро с применением метода Фурье. Влияние толщины и скорости потоков на передачу тепла в сква-жинном теплообменнике.

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

7. Рост теплопередачи от внутреннего оребренного кольца к внешнему кольцу через покоящуюся среду.

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

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

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

Региональной научно-технической конференции «Системы обеспечения тепловых режимов преобразователей энергии и системы транспортировки теплоты» (г.Махачкала, 2005), Международном симпозиуме молодых ученых, аспирантов и студентов «Инженерные и технологические исследования для устойчивого развития» (г.Москва, 2006), I и II Школах молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов» (г.Махачкала, 2006, 2008), Шестой Всероссийской научной молодежной школе «Возобновляемые источники энергии» (г.Москва, 2008), Всероссийской научно-практической конференции «Системы обеспечения тепловых режимов преобразователей энергии и системы транспортировки теплоты» (г. Махачкала, 2008).

Публикации. По теме диссертации опубликовано 8 работ, из них 1 патент Российской Федерации на изобретение.

Объем и структура диссертации. Диссертация состоит из введения, пяти глав, заключения и списка литературы, включающего 104 наименования. Общий объем работы 143 страницы машинописного текста. Работа содержит 53 рисунка и 13 таблиц.

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

В первой обзорной главе рассматриваются состояние и проблемы геотермальной энергетики в регионах Северного Кавказа, перспективы вовлечения низкопотенциальных геотермальных ресурсов в топливно-энергетический баланс Юга России.

В настоящее время потенциал геотермального тепла для большинства регионов Северного Кавказа остается недоиспользованным. Так, по некоторым данным количество эксплуатируемых скважин на месторождениях термальных вод составляет 172 для Республики Дагестан (РД) и 185 для Северного Кавказа (СК), добыча термальных вод 6,0 млн. м3/год для РД и 11,0 млн.м3/год для СК, производство тепловой энергии составляет 100 тыс. Гкал/год для РД и 500 тыс. Гкал/год для СК. В процентном отношении использование потенциала термальных вод составляет: для РД 12% от разведанных запасов и 3% от прогнозных запасов; для СК 12% и 2,4% соответственно. Такое положение дел наблюдается, несмотря на благоприятные предпосылки к интенсивному развитию, т.е. на значительную ресурсную "базу геотермальных вод, достаточно высокую экономическую эффективность, принадлежность геотермального тепла к возобновляемым источникам энергии. В этой связи, проблема развития геотермального производства для частичного замещения поставок природного газа в регионы при постоянном увеличении экспортного потока и роста тарифов представляется весьма актуальной проблемой, имеющей общегосударственное значение.

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

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

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

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

Осуществлена оценка теплопереноса в прямоточном скважинном теплообменнике для подогрева вышележащих низкопотенциальных артезианских вод (глубина залегания ~ 1000 м) путем отбора тепла термальных вод (глубина залегания ~ 3000 м) (рис. 1). Конструкция подобного теплообменника была применена на скв. № Зт Кизлярско-го месторождения термальных вод РД.

Рис. 1. Конструкция скважины с прямоточным теплообменником.

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

Г,(2) = ГФ +ГЛ-Г2 + ]ТС1.е"г; Т2(г) = Т + г{-+- Гг + ¿|

Се4

Л =

( 1 1 Л

1 | д] 1

<а Щг Р,

. Яг

/2±Л Д + ^Д + а

14-а/3

а—

2л

1, к 1, щ. 1. ад!

см — 1п—2—+—1п-2-+—1п— ^ ¿М Ъ К

2л

, 1 . # 1 . ЯГ 1 , +А-ЗД! 7 75 +7 »+Т ——

С, = т„ - т^ - гл -{р{тар - г„)+г+4тар -тю+гДк, - к,)-'

с2 =

Р^-Т^+Г + к,

-т +г

1 <7, 1 а д2а ¡3)

Л,1, Щ - внутренний и внешний радиусы металлической эксплуатационной колонны; ~ внутренний и внешний радиусы обсадной металлической трубы; Л(1) - условный радиус температурного влияния скважины; Н- глубина теплообменника; Г- геотермический градиент; Лв, Дм> Лг ~ коэффициент теплопроводности воды, металла, горной породы; ёп, ¿п - толщины температурного пристеночного слоя на внутренней поверхности внутренней трубы и на наружной поверхности внутренней трубы и внутренней поверхности обсадной трубы; св, сг - объемная теплоемкость воды и горной породы; Чи 42 - объемные расходы термальной воды в теплообменнике (по внутренней трубе) и артезианской воды в теплообменнике (по кольцевому зазору); Т« - температура горячей термальной воды на входе в теплообменник; Гш - пластовая температура горячей термальной воды; Тар - температура артезианской воды на входе в теплообмен-

ник; Тг{2) - невозмущенная температура горной породы вдали от скважины; Т;{2)=Тар - Гг\ Тг(г)= Тт-Гг.

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

Ш^-гЬ-г^+г^е»),

где

г=

2 тг

', Л, - дп Хи Л, Лг Л, >

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

7(2), °С

100-

1- 9;»/ООО.у3/г)Ж ф-500 лУега

2- ц^КЮО^Мт

3- ЦгШОм'/а-т; ц?*2000м?/ет

2,М

Рис. 2. Профили температур теплоносителей дня прямоточного скважинного теплообменника. Гш=110°С; Гар=45°С; Лв =0,63 Вт/м-°С; Ям=45 Вт/м-°С; Лг= 2 Вт/м-°С;

С/; =4,19 МДж/м3-°С; Л,''=28 мм; Я' =30 мм; Я'2 =64 мм; Щ =70 мм;

$г=0,125мм

При больших дебитах теплообмен термальной воды с окружающей горной породой до входа в теплообменник незначителен (потери составляют ~1 -2°С).

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

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

Рис. 3. Схема нагрева пресной воды в противоточном внутрискважинном

теплообменнике. 1 - внутренняя эксплуатационная колонна; 2 - среднее межтрубное пространство; 3 - наружное межтрубное пространство

Г,(г) - температура термальной воды по внутренней трубе; Т2{г) - температура закачиваемой пресной воды в среднем межтрубье; Т3(г) - температура добываемой пресной воды в наружном межтрубье.

Профили температур теплоносителей по высоте теплообменника при различных значениях расходов изображены на рис. 4. Температура термальной воды, поднимающейся с глубины 2000 м, равна 100 С, а начальная температура пресной закачиваемой воды равна 10 °С.

Т(г)°С 85 70 55 40

.л--------------------------------------------------------------------------------

О 50 100 150 гм

Рис. 4. Профили температур теплоносителей по высоте противоточного скважинного теплообменника.

=66 мм; =73 мм; Я'2 =113 мм; Яе2 =123 мм; ^ =174 мм; Д® =184 мм

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

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

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

внутренней трубы отличается от Г/, так как температурный пристеночный слой толщины &г оказывает тепловое сопротивление. Такой же толщины температурный пристеночный слой образуется и на внешней поверхности внутренней трубы, а также на поверхностях ребер. Процесс передачи тепла предполагается установившимся, т.е. не зависящим от времени.

Дифференциальное уравнение, описывающее изменение температуры вдоль прямоугольного ребра с постоянным по длине поперечным сечением получено в предположении малости толщины ребра 28 и температурного пограничного слоя толщины 5т

о>

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

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

Ае-кЦ-*) л ек1'~*) I Л 5 1 - *- к

2 2 Ае'а+еы УЛА ! + *■ 8

где Т. - температура на внешней поверхности внутренней трубы, которую считаем одинаковой, как между ребрами, так и их внутренними торцами.

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

42 8Т г' 8тк ' Ае~2"+1

Температура внешней поверхности трубы определена согласно условиям сохранения теплового потока, что приводит к формулам

8Т ~ Т^.+

Т. =•

Ч „ * '

—+ 2 —

V Л/)

21 + А

г.-Г, =

Т}~Тг 2(1 + ^)

ьл/ у

Поток тепла через межреберную поверхность трубы определяем по приближенной формуле, учитывая, что 8«Я и &г«Я

лЯ - п8

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

,3 =2п6^[т-Тг]=2п^Лв ТГ\ !+1 „ дт 8Т 2^1 + к )Ае + е

Общий поток тепла через трубу с продольными ребрами представляет собой сумму д= д!+д2+Чз

-Т2)

т 8 — + —

в ^м

' ЯД'! п (1-е~и)(1 + Ле'") | п8 А +1 яЯ) лЯк Ае'2" +1 + яй Ае"" + е"

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

д п8 | п (\-е'")(1 + Ае'ы) [ п8 А +1

<7о

лК пКк

Ае'2»* 1

яК Ле'" +е"

(2)

В расчетах принимались следующие данные: Л,=27 мм, 23=2 мм, Хн=0,63, 2.^=45 Вт/м-°С, и=16; 24; 32. Эффективность продольного оребрения при скоростях движения воды и= 1; 0,5 и 0,2 м/с оценена по формуле (2). Результаты приведены в табл. 1. Толщины пристеночного температурного слоя при данных скоростях составили ¿•/=0,125; 0,25 и 0,5 мм соответственно.

Таблица 1

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

Высота ребра /, мм Скорость потока и, м/с

1 0,5 0,2

п=16 п=24 п=32 п=16 п=24 п=32 п=16 п=24 п=32

3 1,304 1,456 1,609 1,401 1,602 1,803 1,471 1,707 1,943

6 1,365 1,547 1,730 1,553 1,829 2,106 1,742 2,112 2,483

9 1,374 1,560 1,747 1,595 1,892 2,189 1,862 2,293 2,724

12 1,375 1,562 1,750 1,605 1,908 2,210 1,910 2,365 2,820

15 1,375 1,562 1,750 1,608 1,911 2,215 1,928 2,392 2,856

18 1,375 1,563 1,750 1,608 1,912 2,217 1,935 2,402 2,870

21 1,375 1,563 1,750 1,608 1,913 2,217 1,937 2,406 2,874

Помимо стальной трубы, были проведены расчеты для случаев алюминиевой и медной конструкций скважинного теплообменника. Медные ребра толщиной 2 мм, высотой 21 мм в количестве 24 штук способствуют значительному увеличению теплового потока, до 5 раз при небольшой скорости потоков о=0,2 м/с.

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

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

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

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

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

лЯ2 лйг И 4я-2Д5 где X - так называемый коэффициент гидравлического сопротивления по трубе.

Экспериментально полученные профили скорости в трубе при развитом турбулентном движении можно приблизить степенным выражением, за исключением малой пристеночной области

и(г) = и0(1-г"/.г)1 л = 20* 50, (3)

где и0 - скорость потока на оси круглой трубы. Для ламинарного движения п=2.

Значение касательного напряжения т0 на стенке трубы получено дифференцированием (3) с учетом г=Я.

пи„ (п + 2Ъ V . Я к Я А п+2

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

Формула (4) указывает лишь на возможность введения понятия вязкого пристеночного слоя А и дает его ориентировочную толщину.

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

2яЛ£г0 =АрлИг; = Х^-—лК1 \ Д = (5)

А 2 2Й Яри

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

Для кольцевого сечения с ребрами вторая из формул (5) примет вид

(2яй, + 2л/ + 2лй2 )1 г0 = (я{Д\ - Я?) - пЩЬр,

где Я\- внешний радиус внутренней оребренной трубы; Я2 - внутренний радиус наружной трубы. Для динамического перепада давления в этом случае получено выражение

2 п(Кх+Кг)+п1

Др ------— ---Ц!-^-ри Ь

- Я, ) - п1д 4 л{Я\ -¡{*)-п18

Средняя скорость потока вдоль кольцевого сечения

и перепад давления при движении потока по кольцевому сечению с ребрами

Я Л{К+К)+Ш

Р А(я<Л1-%)-п15)г

Для кольцевого сечения без ребер при /=О оно переходит в известное в технической литературе выражение

АР, =

л рО>1 4 л-^-ЯЛЛг+Л.)2

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

Ар АР,

1 + -

п1

я(Я2+Я1))

1 —

п18

~ Щ))

(6)

Оребрение поверхности приводит к значительному росту перепада давления (до 5 раз и более) в зависимости от количества ребер и их высоты (рис. 5).

25 I, ММ

Рис. 5. Прирост перепада давления в зависимости от высоты ребер

и их числа. Л/=27 мм; /?2=б0 мм; 28=2 мм

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

На рис. 6 приводятся эти подсчитанные отношения для случая стальной конструкции.

Рис. 6. К выбору критерия оптимальности для стальной конструкции

теплообменника.

Для меди, когда Хм =300 Вт/м-°С, графики эти имеют более пологий вид, оптимум размазывается, и высоты рёбер в 2-3 раза выше по сравнению со случаем стали.

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

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

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

Далее найдено решение для элемента симметрии, принимаемого за Г-образную фигуру. Она включает в себя рассеченную посередине половину ребра и примыкающую к ней половинку куска трубы. Эта Г-образная фигура разбивается на три прямоугольные части (первая часть - половина ребра толщины <5; вторая часть - кусок трубы на стыке с ребром; третья часть - кусок трубы между ребрами высотой К).

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

\Ч4 Т

\яз

2 ф 1

Рис. 7. Элемент симметрии для оребренной трубы.

1-ребро; 2-кусок трубы на стыке с ребром; 3-кусок трубы между ребрами

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

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

Распределение температуры внутри ребра (область 1) определяется рядом

(( (

я=0 I А V

Ч\ \ « ^

гвтс^)

1+5тфРп)

г-^СОв^у); ^ЯГЛ + ОГ,

Значения а„ находим в виде корней уравнения у = (л-л + - А

Распределение температуры внутри куска трубы на стыке с ребром (область 2)

2 „=« М„Щ2-Мп) 2

Распределение температуры внутри межреберной части (область 3)

Х,~—-

/ хх 2(1-(1+2Л)ги з1{хпх)Мх„у) . _ , Пх,у)-1-{1-Ту)2+ я X,- —

Картина распределения температуры в трех областях контактной задачи представлена на рис. 8. С целью достижения сходимости рядов Фурье число их членов берется равным 1000 .

у, ми

0.« ал

0.401-

0.М

о:в т

о'},! Г (_¡_л.

-М.

г 1

1.Т Т

0.5 1 и а о а » е « х, мм

Рис. 8. Распределение температуры в областях контактной задачи. Лл*=45 Вт/м-°С; 8 =1 мм; А =5 мм

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

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

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

Я,

<5 т

7*0

7=0

Рис. 9. Разностная сетка для ребра.

I 8

Ах =—, Ау = —, *,=Дх-1, / = 0,1 у.=8-Ъу-], у = 0,1,..., и

пх пу

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

Матричная запись разностных уравнений ребра имеет вид

ль

(т \ ' ¡-1,0 \ '/+1,0

+ ЛС ^1+1,1

Т V '-'■»г > Л. , Т КМ*»)

= 0,где А1 = АР =

а + — О О 2

0 6 0

О О Ь

О 0 0 О 0 0

О О О о о о

ь о о ь

АС =

~(2а + Ь + с + \) 1 0

1 -2(1 + 6) 1 О 1 -2(1 + Ь)

О О

О 0

О О О

-2(1 +Ь) 1 2 -2(1 + 6)

для I = 1,2,...,пх -1; у = в...пу

Матричная запись уравнений на торце ребра при { = пх несколько изменится

= 0,где АЬп =

(т Л ^т Л

АЬп + /4Си т„л

Т Т V. «.V у

АСп=

1+Е.

О О

0 0 .. 0 0

0 2 Ь 0 .. 0 0

0 0 2Ь .. 0 0

0 0 0 .. 2 Ь 0

0 0 0 .. 0 Ъ

... 0

... -Ам+я+г^]

^ л

1 4+6+у)

яАЛу . . _ ЛУ

Дх2

2ЯмЛх

Ах

ЛмАх

Л„Ау , 1 + ,

; с= 8 * ; --Ь; е = с

ЛМ5Т \ + еу

ЛВ8Т Л/4у

1 , АУ + 2<5Г Ах 1 + ^

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

Далее приводится применение метода поперечно-продольной прогонки для элемента симметрии ребристой теплопередающей поверхности.

Элемент симметрии теплопередающей поверхности представляет собой дугу окружности толщины 23 и высоты ЪсШЪп с реберным продолжением толщины 8 и длины /. При большом числе ребер п (12 и более) можно пренебречь кривизной трубы и принять кусок трубы за прямоугольную область. На рис. 10 показан элемент симметрии и его замена на область из двух состыкованных прямоугольников.

щ ь» т„.

<5 1

¿я

*П

■23-

Нх х -I

Рис. 10. Сеточная область для элемента симметрии.

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

На границах у= 0 и у=к ставятся условия отсутствия оттока тепла - условия симметричного продолжения решения. На внутренней и внешней поверхностях ставятся условия теплообмена через пристеночный температурный слой толщины 8п и 6т2-

Граничные условия аппроксимировались со вторым порядком точности, так же как и оператор Лапласа. Однако вместо уравнения Лапласа для определения установившегося поля температур применялось двумерное уравнение теплопроводности,' решение которого отыскивалось для больших значений времен. Переход на стационарный режим происходит за конечное время, ибо граничные условия задачи жестко заданы и область решения имеет малые размеры. Уравнение теплопроводности

1 дТ д2Т дгТ 2 Л .. _ .,.. -г—= —г + —г> а =— ><>0 а2 81 дх ду2 рс

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

грП +1 грц

1Ы-ж =

1-1,> г-1ц + Л+и '¡,]-\ + 1 ¡.¡*\

— - — — + ------—----

Ах2 Ду2

Преимущество неявной схемы в том, что временной шаг Д1 можно брать любым, счет по такой схеме устойчив. Регулированием шага А1 можно уменьшить машинное время расчетов достижения установившегося режима теплообмена. Счет прекращается, когда поправка к решению для следующего слоя в относительных величинах становится менее 10"11.

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

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

Для проведения расчетов была составлена программа на языке Турбо-паскаль. Счет прекращался, когда максимальная разность температур в узлах сетки становилась меньше 10".

На рис. 11 изображено температурное поле элемента симметрии стальной трубы.

■ 0,4-0.45 0 0,35-0,4

■ 0,3-0,35 Ш 0,25-0,3

■ 0,2-0,25

□ 0,15-0,2

□ 0.1-0.15

■ 0,05-0,1 И 0-0,05

9 ю4

поперечная координата, мм

Рис. 11. Температурное поле элемента симметрии стальной трубы.

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

Рис. 12. Теплопередача через неподвижную среду.

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

* 87 0

" п' 8<р~

й * Е.т о &

Эр

Рис. 13. Разностная сетка для области задачи.

А<р = -^—, Аг = ———; 1 = п,Аг пп9 пг

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

Сеточное уравнение для внутренних точек области задачи + (т;.,-, - (2 + а1, + ар, % + Г, у+1)+ ар,Г|+и = 0,

Сеточное уравнение для точек на верхней границе

+ [2Т1Н -(2 + а/, + ар,щГми = О Сеточное уравнение для точек ребра на нижней границе

й _ К8гМ

+ И + Чк + )+ =0. ' " ЛвАг2

Сеточное уравнение для точек торцевой поверхности ребра на нижней границе

_1 г А(рг

¿Л-и + (- (1 + + аг, )Г.} + аг,Тм. = 0 . аг> 2Г< Дг*

Сеточное уравнение для точек на нижней границе вне ребра "М-и + (" (2 + а1< + аР, К + 2 )+ ар,ТМ1 = 0 Матричная запись разностных уравнений имеет вид

АЬ,

(т ^ '/-1,0

(т

+ АС,

л-АР,

(т ^ 11+1,0

^/+1,1

= 0

II, 0 0 0 а/,. О О 0 а1,

ООО ООО

О О

а], О

О а1.

,где //,. = с1, для г < п1; И, = а/,, для г > и/

АР,=

'АС, =

рр< 0 0 0 0

0 аР, 0 0 0

0 0 ар, 0 0

0 0 0 ар, 0

о 0 0 0 ар,

для I < п1; рр! = аг, для I

2 0 . . 0 0'

1 1 . . 0 0

0 Л . о 0

0 0 0 .. ■ 2,. 1

0 0 0 . . 2 г,-. ДЛЯ

-(2 + а/,. + ар;)

сс1 1' о. .'. . 0 0"

1 1 .. . 0 0

0 1 . 0 0 для

0 0 0 .. • ч 1

0 0 0 .. . 2

I > п1

АС: =

сс- = -(1 + 2(1,) для {< п1; сс, = -(1 + с11 + аг1) для г = п1.

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

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

скважинного теплообменника до нескольких раз (до 3-х раз при наличии 32 стальных ребер высотой 21 мм при Я;=30 мм, Л2=60 мм).

В ЗАКЛЮЧЕНИИ формулируются основные результаты и выводы:

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

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

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

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

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

Методика предусматривает введение вязкого пристеночного слоя турбулентных течений, толщина которого определяется гидравлическим сопротивлением в зависимости от числа Рейнольдса. Приведены зависимости приведённого перепада давления от высоты ребер для различного их числа. Оребрение поверхности приводит к значительному росту перепада давления при большом количестве ребер.

6. Сформулирован критерий оптимальности выбора параметров скважинного теплообменника с ребрами. За такой критерий принят максимум отношения приведенного теплового потока к приведенному дебиту скважины. Эти отношения подсчитаны для стальной и медной конструкций скважинного теплообменника.

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

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

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

1. Мамаев НИ., Алхасова Д.А. Повышение эффективности тепло-обменных аппаратов путем интенсификации процесса теплообмена.// Материалы Международной конференции «Возобновляемая энергетика: проблемы и перспективы». Т. 1. Махачкала, 2005. С. 326-330.

2. Мамаев Н.И., Алхасова Д.А. Расчет гидравлических сопротивлений и тепловых потоков в теплообменниках с продольными ребрами.// Материалы региональной научно-технической конференции «Системы обеспечения тепловых режимов преобразователей энергии и системы транспортировки теплоты». Махачкала, 2006. С. 20-26.

3. Алхасова Д.А., Мамаев Н.И. Исследование процесса теплообмена на модели внутрискважинного теплообменника.//Материалы I школы молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала, 2006. С. 172-177.

А. Мамаев Н.И., Алхасова Д.А. Способ подогрева воды для отопления и установка для его осуществления.//Патент RU 2336466 С2. 2008. Бюл. № 29.

5. Алхасова Д. А, Увеличение теплообмена между потоками ореб-рением перегородки.//Современные проблемы математики и смежные вопросы. Материалы Международной конференции «Мухтаровские чтения». Махачкала, 2008. С. 31-34.

6. Алхасова ДА. Расчет прямоточного внутрискважинного тепло-обменника.//Материалы шестой Всероссийской научной молодежной школы «Возобновляемые источники энергии». Москва, 2008. С. 13-19.

7. Алхасова ДА. Расчет влияния толщины на передачу тепла от внутреннего потока скважинного теплообменника к внешнему потоку через оребренную поверхность .//Материалы II школы молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала, 2008. С. 229-236.

8. Алхасова ДА. Оценка теплопереноса в геотермальном сква-жинном теплообменнике.//Материалы Всероссийской научно-практической конференции «Системы обеспечения тепловых режимов преобразователей энергии и системы транспортировки теплоты». Махачкала, 2009. С. 51-55.

Заказ №73. Тираж 100.

Типография Дагестанского научного центра РАН 367015, г. Махачкала, 5-й жилгородок, корпус 10

Содержание диссертации автор исследовательской работы: кандидата технических наук, Алхасова, Джамиля Алибековна

ВВЕДЕНИЕ.

ГЛАВА 1. СОСТОЯНИЕ И ПРОБЛЕМЫ ГЕОТЕРМАЛЬНОЙ ЭНЕРГЕТИКИ В РЕГИОНАХ СЕВЕРНОГО КАВКАЗА.

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

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

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

ГЛАВА 2. ТЕХНОЛОГИЧЕСКИЕ СХЕМЫ СЪЕМА ТЕПЛА ТЕРМАЛЬНЫХ ВОД.

2.1. Сравнение расчетных и конструктивных характеристик различных типов теплообменных аппаратов.

2.2. Расчет прямоточного внутрискважинного теплообменника.

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

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

ГЛАВА 3. ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ПРОЦЕССА ТЕПЛООБМЕНА ПРОДОЛЬНЫМ ОРЕБРЕНИЕМ ТЕПЛОПЕРЕДАЮЩЕЙ ПОВЕРХНОСТИ.

3.1. Теплопередача через ребристую стенку. Вывод уравнения теплопроводности для ребра постоянной толщины в установившемся режиме

3.1.1. Осреднение температуры ребра по поперечному сечению.

3.1.2. Понятие температурного пристеночного слоя.

3.2. Схема сосредоточенной емкости.

3.3. Расчет гидравлических сопротивлений оребренных поверхностей.

ГЛАВА 4. РЕШЕНИЕ ЗАДАЧ ТЕПЛОПЕРЕДАЧИ ОТ ВНУТРЕННЕГО ПОТОКА К ВНЕШНЕМУ ПОТОКУ С УЧЕТОМ ТОЛЩИН ТРУБЫ И РЕБРА.

4.1. Ребро неограниченной высоты.

4.2. Ребро конечной высоты.

4.3. Контактная задача.

4.3.1. Часть 1 - решение для ребра.

4.3.2. Часть 2 — решение для куска трубы.

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

4.4.1. Часть 3 - решение для куска трубы между ребрами.

ГЛАВА 5. КОНЕЧНО-РАЗНОСТНЫЕ МЕТОДЫ РЕШЕНИЯ.

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

5.2. Конечно-разностная схема для прямоугольного ребра. Применение метода матричной прогонки.

5.3. Применение метода поперечно-продольной прогонки для элемента симметрии.

5.4. Разностная задача для определения эффективности оребрения внутренней поверхности для теплопередачи через неподвижную кольцевую среду.

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

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

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

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

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

В соответствии с целевым направлением работы задачами исследования явились:

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

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

Диссертация состоит из введения, пяти глав, заключения и списка литературы, включающего 104 наименования. Общий объем работы 143 страницы машинописного текста. Работа содержит 53 рисунка и 13 таблиц.

Заключение диссертации по теме "Теплофизика и теоретическая теплотехника"

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

6. Сформулирован критерий оптимальности выбора параметров скважин-ного теплообменника с ребрами. За такой критерий принят максимум отношения приведенного теплового потока к приведенному дебиту скважины. Эти отношения подсчитаны для стальной и медной конструкций скважинного теплообменника.

7. Рассмотрено применение метода Фурье для расчета теплопередачи через оребренную поверхность. Учитывается как толщина ребра и изменение температуры поперек ребра, так и толщина трубы и изменение температуры, как по радиусу, так и по окружности трубы. То есть рассматривается общая двумерная задача распределения температуры. Сравнение с ранее полученными результатами показало, что для миллиметровых толщин влияние изменения температуры поперек трубы и ребра несущественно.

Список источников диссертации и автореферата по физике, кандидата технических наук, Алхасова, Джамиля Алибековна, Махачкала

1. Алишаев М.Г. О притоке к скважинам с щелевым фильтром // Известия АН СССР, МЖГ. 1967. №5.

2. Алхасов А.Б. Геотермальная энергетика: проблемы, ресурсы, технологии. М.: Физматлит, 2008. 376с.

3. Алхасов А.Б. Перспективы освоения геотермальных ресурсов Северного Кавказа // Перспективы энергетики. 2003. Т.7. С. 367-375.

4. Алхасов А.Б. Способ одновременно-раздельной эксплуатации двух термоводоносных пластов/ Патент RU 2105351 С1 РФ// Открытия. Изобретения. 1998. № 5.

5. Алхасов А.Б. Технологии освоения геотермальных ресурсов осадочных бассейнов // Материалы Международной конференции ИПГ ДНЦ РАН «Возобновляемая энергетика: проблемы и перспективы». Махачкала, 2005.Т. 1. С.173-185

6. Алхасов А.Б., Алишаев М.Г. Извлечение тепла грунта скважинным теплообменником в сезонном режиме работы // Известия РАН. Энергетика. 2007. №2. С. 137-144.

7. Алхасов А.Б., Алишаев М.Г. Скважинный теплообменник для съема тепла с горной породы // Материалы международной конференции ИПГ ДНЦ РАН «Возобновляемая энергетика: проблемы и перспективы», Махачкала, 2005. т.1, с.263-275.

8. Алхасов А.Б., Гайдаров Г.М., Магомедбеков Х.Г. Паротурбинная установка для геотермальной электростанции / Патент RU 2035588 С1 РФ / Открытия. Изобретения. 1995. №14.

9. Алхасов А.Б., Исрапилов М.И. Использование геотермальной энергии для подогрева подпиточной воды // Водоснабжение и санитарная техника. 1996. №4. С. 25-26.

10. Алхасов А.Б., Раджабов Р.И. Теплообмен при совместно-раздельной добыче термальных вод из двух горизонтов// Вестник ДНЦ РАН. 2000. №8. С. 28-33.

11. Алхасов А.Б., Рамазанов М.М., Абасов Г.М. Использование геотермальной энергии в горячем водоснабжении // Водоснабжение и санитарная техника. 1998. № 3. С. 24-25

12. Алхасова Д.А. Увеличение теплообмена между потоками оребрением перегородки // Материалы Международной конференции «Мухтаровские чтения». Современные проблемы математики и смежные вопросы. Махачкала. 2008. С. 31—34.

13. Алхасова Д.А. Расчет прямоточного внутрискважинного теплообменника // Материалы шестой Всероссийской научной молодежной школы «Возобновляемые источники энергии». Москва, 2008. С. 13-19.

14. Алхасова Д.А., Мамаев Н.И. Исследование процесса теплообмена на модели внутрискважинного теплообменника // Материалы I школы молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала, 2006. С. 172-177.

15. Антуфьев В.М. Эффективность различных форм конвективных поверхностей нагрева. M-JL: Энергия, 1966. 182 с.

16. Альтшуль А.Д. Гидравлические сопротивления. М.: Недра, 1982.218с.

17. Андреев В.А. Теплообменные аппараты для вязких жидкостей. JL: Энергия, 1971. 151 с.

18. Барановский Н.В., Коваленко JI.M., Ястребенецкий А.Р. Пластинчатые и спиральные теплообменники. М.: Машиностроение, 1973. 286 с.

19. Барило А.А., Маслова Н.А. Химико-аналитические аспекты использования геотермальных вод // Вщновлювана енергетика. 2006, №4. С.73-77.

20. Белан С.И., Панин JI.H., Гусейнов Н.М. Оценка перспектив использования геотермальных вод Юга России // Материалы Международной конференции ИПГ ДНЦ РАН «Возобновляемая энергетика: проблемы и перспективы», Махачкала, 2005. т.1. С.166-170.

21. Берман Э. Геотермальная энергия. М.: Мир, 1978. 416 с.

22. Бицадзе А.В., Калиниченко Д.Ф. Сборник задач по уравнениям математической физики. М.: Наука, 1977. 379с.

23. Богуславский Э.И., Певзнер JI.A., Хахаев Б.Н. Перспективы развития геотермальной технологии // Разведка и охрана недр, №7-8, 2000. С.43-48.

24. Болгарский А.В., Мухачев Г.А., Щукин В.К. Термодинамика и теплопередача. М.: Высшая школа, 1964. 458 с.

25. Васильев В.Я. Экспериментальное исследование рациональной интенсификации конвективного теплообмена в прямоугольных каналах трубчато-пластинчатых теплообменных поверхностей с поперечными выступами и ка-навкми // Теплоэнергетика. 2007, № 8. С.50-62.

26. Васильев В.Я. Экспериментальное исследование рациональной интенсификации конвективного теплообмена в прямоугольных рассеченных каналах пластинчато-ребристых теплообменных поверхностей // Теплоэнергетика. 2006, № 12. С.58-67.

27. Васильев Г.П., Шилкин Н.В. Использование низкопотенциальной тепловой энергии земли в теплонасосных системах. АВОК, №2, 2003. С.52-60.

28. Гортышов Ю.Ф., Олимпиев В.В., Попов И.А. Эффективность про-мышленно эффективных интенсификаторов теплопередачи (Обзор. Анализ. Рекомендации) // Изв. РАН. Энергетика. 2002, №3. С. 102-118.

29. Гортышов Ю.Ф., Попов И.А. Научные основы расчета высокоэффективных компактных теплообменных аппаратов с рациональными интенсифика-торами теплоотдачи //Теплоэнергетика. 2006, №4. С.2-13.

30. Готовский М.А. Интенсификация конвективного теплообмена и самоорганизация вихревых структур. // Теплоэнергетика. 1995 №3. С.55-59.

31. Доброхотов В.И., Поваров О.А. Использование геотермальных ресурсов в энергетике России // Теплоэнергетика. 2003. № 1. С. 2-11.

32. Дрейцер Г.А. Проблемы создания высокоэффективных трубчатых теплообменных аппаратов // Теплоэнергетика. 2006, №4. С.31-38.

33. Дрейцер Г.А. Проблемы создания компактных трубчатых теплообменных аппаратов // Теплоэнергетика. 1995, №3. С. 11-18.

34. Дрейцер Г.А., Мякочин А.С. Влияние геометрической формы турбу-лизаторов на эффективность интенсификации конвективного теплообмена в трубах// Теплоэнергетика. 2002, № 6. С.57-59.

35. Дрейцер Г.А., Исаев С.А., Лобанов И.Е. Расчет конвективного теплообмена в трубах с периодическими выступами // Проблемы гидродинамики и теплообмена в энергетических установках. М.: Изд-во МЭИ. 2003. Т.1. С. 57-60.

36. Дрейцер Г.А. Эффективность использования закрутки потока для интенсификации теплообмена в трубчатых теплообменных аппаратах // Теплоэнергетика. 1997, № 11.С.61-65.

37. Дядькин Ю.Д., Гендлер С.Г. Процессы тепломассопереноса при извлечении геотермальной энергии. JL: изд. ЛГИ, 1985. 93 с.

38. Идельчик И.Е. Гидравлические сопротивления М.: Госэнергоиздат, 1964. 316с.

39. Извлечение тепловой энергии термальных вод / Г.М.Гайдаров, А.Б.Алхасов, Ш.А.Гаджиев, З.Д.Абуев // Разведка и охрана нед. 1990. № 8. С. 41-43.

40. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача. M.-JL: Энергия, 1965. 423 с.

41. Исследование режимов совместной работы теплового насоса с вертикальным грунтовым теплообменником / Федянин В.Я., Утемесов М.А., Федин Л.Н., Горбунов Д.Л.// Теплоэнергетика, 1997, №4. С.21-23.

42. Калинин М.И., Баранов А.В. Метод расчета глубинных теплообменников для односкважинной технологии геотермального теплоснабжения // Охрана и разведка недр. 2003. №6. С. 53-60.

43. Калинин М.И., Хахаев Б.Н., Баранов А.В. Геотермальное теплоснабжение центральных регионов России с использованием мелких и глубоких скважин. Электрика, №4, 2004. С.8-13.

44. Калинин Э.К., Ярхо С.А. Влияние чисел Рейнольдса и Прандтля на эффективность интенсификации теплообмена в трубах // Инженерно-физический журнал. Т. XX, 1971, №4. С.592-599.

45. Калинин Э.К., Дрейцер Г.А., Ярхо С.А. Интенсификация теплообмена в каналах. М.: Машиностроение, 1990. 211с.

46. Кейс В.М., Лондон А.Л. Компактные теплообменники. М.: Госэнерго-издат, 1962.158с.

47. Клепанда А.С., Филиппов Э.Б., Пашко П.В. Методика расчета на ЭВМ парокомпрессионного теплового насоса//Холодильная техника, №7, 1990, с. 1013.

48. Кошкин В.К., Калинин Э.К. Теилообменные аппараты и теплоносители. М.: Машиностроение, 1971. 198с.

49. Курбанов М.К. Геотермальные и гидроминеральные ресурсы Восточного Кавказа и Предкавказья. М.: Наука, 2001. 206с.

50. Кутателадзе С.С. Основы теории теплообмена. Новосибирск, Наука, 1970. 660с.

51. КухлингХ. Справочник по физике. М.: Мир, 1985. 520 с.

52. Лабунцов Д.А. Вопросы теплообмена при пузырьковом кипении жидкости // Теплоэнергетика. 1972. № 9. С. 14 19.

53. Лабунцов Д.А. Обобщенные зависимости для теплоотдачи при пузырьковом кипении жидкостей // Теплоэнергетика 1960 . № 5. С. 76 81.

54. Лебедев П.Д., Щукин А.А. Теплоиспользующие установки промышленных предприятий. М.: Энергия, 1970. 408с.

55. Леонтьев А.И., Гортышов Ю.Ф., Олимпиев В.В., Попов И.А. Эффективные интенсификаторы теплоотдачи для ламинарных (турбулентных) потоков в каналах энергоустановок //Изв. РАН. Энергетика. 2005, №1. С.75-91.

56. Лыков А.В. Теория теплопроводности. М.: Гостехиздат, 1952.392 с.

57. Люстерник Л.А. Замечания к численному решению краевых задач уравнения Лапласа и вычислению собственных значений методом сеток.// Тр. Математического института им. В.А. Стеклова. Том 20. М.: АН СССР, 1947. С. 49-64.

58. Люстерник Л.А. О разностных аппроксимациях оператора Лапласа. //Успехи математических наук, 1954, т. 9, выпуск 2. С. 3-66.

59. Магомедов К.М. Теоретические основы геотермии. М.: Наука. 2001,277с.

60. Макаров A.M., Лунева Л.А., Залогина Н.И. Нестационарный процесс теплообмена в трубе с продольным оребрением // Инженерно-физический журнал. 1994. Том 66. № 6. С.673-680.

61. Мамаев В.А., Одишария Г.Э., Клапчук О.В. Гидродинамика газожидкостных систем. М.: Недра, 1981. 356с.

62. Мамаев Н.И., Алхасова Д.А. Повышение эффективности теплообменных аппаратов путем интенсификации процесса теплообмена // Материалы Международной конференции «Возобновляемая энергетика: проблемы и перспективы». Т. 1. Махачкала. 2005. С. 326—330.

63. Мамаев Н.И., Алхасова Д.А. Способ подогрева воды для отопления и установка для его осуществления //Патент RU 2336466 С2. 2008. Бюл. № 29.

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

65. Назмеев Ю.Г., Конахина И.А. Интенсификация теплообмена при течении вязкой жидкости в трубах с винтовой накаткой // Теплоэнергетика. 1993, №11. С.59-62.

66. Пермяков В.А., Пермяков К.В., Якименко А.Н., Нейбургер А.Н. К вопросу выбора типа водо-водяных подогревателей для систем теплоснабжения // Промышленная энергетика. 2000, №4. С.37-44.

67. Плановский А.Н., Николаев П.И. Процессы и аппараты химической и нефтехимической технологии. М.: Химия, 1972. 393с.

68. Поваров О.А., Томаров Г.В. Всемирный геотермальный конгресс WGC-2005 // Теплоэнергетика. 2006. № 3. С. 78-80.

69. Поваров О.А., Томаров Г.В., Семенов В.Н. Исследование и опыт внедрения от коррозии и эрозии металла паровых турбоустановок // Теплоэнергетика. 2002, №12. С.22-28.

70. Проблемы геотермальной энергетики Дагестана / Под ред. Амирхано-ва Х.И. и Ятрова С.Н. М.: Недра, 1980. 208с.

71. Промышленная теплоэнергетика и теплотехника: Справочник. Под Общ. ред. Клименко А.В. и Зорина В.М. М.: МЭИ, 2004. 632 е.

72. Ройзен Л.И., Дулькин И.Н. Тепловой расчет оребренных поверхностей. М.: Энергия, 1977. 114с.

73. Рихтмайер З.В. Разностные методы решения краевых задач. ИИЛ, М.: 1960. 429с.

74. Самарский А.А., Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 1978. 592 с.

75. Семенов В.Н. Исследование поведения примесей многокомпонентного теплоносителя в тракте геотермальной электростанции // Теплоэнергетика. 2006, №3. С. 17-22.

76. Системы извлечения тепла земной коры и методы их расчета / Щер-бань А.Н., Цирульников А.С., Мерзляков Э.И., Рыженко И.А. Киев, Наук думка, 1986, с.248.

77. Телегин А.С., Швыдкий B.C., Ярошенко Ю.Г. Тепломассоперенос. М.,: ИКЦ Академкнига., 2002 .455с.

78. Теоретические основы теплотехники. Теплотехнический эксперимент: под общ. ред. Клименко А.В. и Зорина В.М. М.: МЭИ, 2001. 564с.

79. Теплообмен и гидравлическое сопротивление при ламинарном течении вязкой жидкости в трубах с искусственной шероховатостью / Ю.Г. Назме-ев, A.M. Конахин, Б.А. Кумиров и др.// Теплоэнергетика. 1993. №4. С.66-69.

80. Теплотехнический справочник / Под ред. В.Н. Юренева, П.Д. Лебедева, Т. 1. М.: Энергия, 1975. 744 с.

81. Теплотехнический справочник / Под ред. В.Н. Юренева, П.Д. Лебедева, Т. 2. М.: Энергия, 1976. 896 с.

82. Тихонов А.Н., Самарский А.А. Об однородных разностных схемах. ЖВМиМФ, 1961, т. 1, №1.С. 278-291.

83. Тихонов А.Н., Самарский А.А. Однородные разностные схемы высокого порядка точности на неравномерных сетках. ЖВМ и МФ, 1961, т.1, №3. С. 425-440.

84. Тихонов А.Н., Самарский А.А. Однородные разностные схемы на неравномерных сетках. ЖВМ и МФ, 1962, т. 2, №5. С. 812-832.

85. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: Наука, 1977.736 с. '

86. Чекалюк Э.Б. Термодинамика нефтяного пласта. М.: Недра, 1965.238с.

87. Черменский Г.А. Геотермия. М., Недра, 1972. 265с.

88. Швец М.Ю. Исследование эффективности использования различных типов теплообменного оборудования для геотермальной энергетики// Вщновлювана енергетика.2007, №2. С.80-85.

89. Шлихтинг Г. Теория пограничного слоя. М.: ИИЛ., 1956. 528с.

90. Шпильрайн Э.Э. Проблемы и перспективы возобновляемой энергии в России / Перспективы энергетики, 2003, т.7. С.393-403.

91. Шпильрайн Э.Э. Экологические аспекты применения возобновляемых источников энергии для децентрализованного энергоснабжения / Перспективы энергетики, 2002, т.6. С.299-306.

92. Шрадер И.Л., Дашчян А.А., Готовский М.А. Интенсифицированные трубчатые воздухоподогреватели // Теплоэнергетика. 1999. №9. С.54-56.

93. Эффективные интенсификаторы теплоотдачи для ламинарных (турбулентных) потоков в каналах энергоустановок / А.И. Леонтьев, Ю.Ф. Гортышов, В.В. Олимпиев., И.А. Попов // Изв. РАН. Энергетика. 2005. №1. С.75-91.

94. Юдаев Б.Н. Техническая термодинамика. Теплопередача. М.: Высшая школа, 1988. 479 с.

95. Ягов В.В. Научное наследие Д.А. Лабунцова и современные представления о пузырьковом кипении // Теплоэнергетика. 1995. №3. С.2-9.

96. Curtis R., Lund J., Sanner В., Rybach L., Hellstrom G. Ground Source Heat Pumps Geothermal Energy for Anyone, Anywhere: Current Worldwide Activity//Proceedings World Geothermal Congress 2005, Antalya, Turkey, 24-29 April 2005, 9p.

97. Hellstrom G., Sanner B. PC-programs and modeling for borehole heat exchanger design // International Summer School on Direct Application of Geothermal Energy. Bad Urach, 2001, p.35-44.

98. Sanner В. Examples of GSHP and UTES System in Germany// Proceedings World Geothermal Congress 2005, Antalya, Turkey, 24 April 2005, 14p.

99. Sanner В., Kohlsch O. Examples of Ground Source Heat Pumps (GSHP) from Germany/Anternational Summer School on Direct Application of Geothermal Energy. Bad Urach, 2001, p.81-94.