Тепловой расчет теплообменника
Выполнение теплового расчёта теплообменника — инженерный расчёт, позволяющий перейти от технического задания (например, «требуется ГВС 200 кВт») к конкретному аппарату с определённой площадью теплопередачи, количеством пластин, гидравлическим сопротивлением и резервом по мощности. Это не «подбор готовой модели по каталогу»: речь идёт о цепочке из 8 расчётных шагов, базирующихся на уравнении теплового баланса Q = G · c · ΔT, логарифмическая температурная разность LMTD, величина коэффициента теплопередачи U, фактор загрязнения Rf а также проверочный гидравлический расчёт. В материале собраны методика, формулы, справочные значения, пример вычисления для ГВС 200 кВт, обзор программных продуктов (SWEP SSP, Alfa Laval CAS, Kelvion KOPT, Sondex Pro) и 6 частых ошибок, из-за которых подобранный аппарат не обеспечивает паспортную мощность на площадке.
Ключевые величины: λ, α, U и фактор загрязнения Rf
Прежде чем определять площадь поверхности и количество пластин, специалисту нужно чётко различать четыре основополагающих параметра. Даже в технической документации их нередко смешивают: вместо «коэффициента теплопередачи» указывают «коэффициент теплоотдачи», а фактор загрязнения считают как % от площади, хотя это термическое сопротивление в м²·K/Вт. Рассмотрим это подробно.
Способность стенки пропускать тепловой поток. Единица измерения Вт/(м·K). У AISI 304 этот показатель равен 15 Вт/м·K, у AISI 316 — 16, у титана Gr.1 — 22, у меди — 380, у графита — 120. С ростом λ термическое сопротивление стенки δ/λ падает. У тонких пластин (0,5–0,6 мм) доля стенки в суммарном сопротивлении ≤ 3%.
Насколько активно тепло переходит «жидкость → стенка». Единица измерения Вт/(м²·K). Определяется скоростью потока, вязкостью, критериями Re и Pr. В гофрированных каналах для воды α = 8 000–25 000, при конденсации пара — 12 000–18 000, для масла — 800–2 500. Находят по критериальным уравнениям такого вида Nu = C·Rem·Prn.
Интегральный параметр «жидкость 1 → стенка → жидкость 2». Определяется по формуле 1/U = 1/α₁ + δ/λ + 1/α₂ + Rf1 + Rf2. Единица измерения — Вт/(м²·K). В главную формулу подставляется именно U F = Q/(U·LMTD·Φ). Для пластинчатых теплообменников обычные значения: вода-вода 3 500–6 500, вода-пар 8 000–12 000, вода-гликоль 2 500–4 800.
Сопротивление теплопереносу от отложений накипи, ила и биоплёнки. Единица измерения м²·K/Вт. Учитывается заранее, чтобы через 5–7 лет работы мощность аппарата не просела. Для чистой деаэрированной воды Rf = 0,0001, для городской подпитки ГВС — 0,0003–0,0005, для скважинной воды с железом — до 0,001. Или же закладывают запас по площади 15–30%.
Эти четыре параметра — фундамент всей методики. Когда в ТЗ уже указано, что α со стороны воды составит 14 000, а α со стороны гликоля — 6 500, инженер быстро оценит: U будет примерно 3 200–3 800 Вт/(м²·K), а поверхность попадёт в диапазон 4–6 м² на каждые 100 кВт мощности. Об конструкции пластинчатого теплообменника подробнее — в материале «как работает теплообменник».
Базовая формула: Q = G · c · ΔT
Любой тепловой расчёт пластинчатого или кожухотрубного расчёт аппарата стартует с уравнения теплового баланса. Это не просто «одна из формул», а основа. Остальные зависимости (LMTD, F, ε-NTU, критерий Re) нужны позже — чтобы от вопроса «сколько энергии переносим» перейти к «какая нужна площадь и при каком сопротивлении».
Формула применяется дважды — отдельно по горячему и холодному контуру; обе мощности должны сойтись с учётом теплопотерь (как правило, ≤ 1%). Так проверяют согласованность ТЗ. Если греющая сторона даёт 200 кВт, а нагреваемая 235 кВт, значит, у заказчика ошибка в температурах или расходах.
Вспомогательные формулы: расход и теплоёмкость
В практических ТЗ расход чаще указывают в м³/ч (объёмный), тогда как формула требует кг/с. Перевод: G [кг/с] = GV [м³/ч] · ρ [кг/м³] / 3600. У воды при 70 °C ρ = 977,8 кг/м³, при 95 °C — 961,9. У пропиленгликоля 30% при 60 °C ρ ≈ 1 022. Ошибка в плотности на 3% (если взять воду 1 000 кг/м³ вместо фактической) сдвинет мощность на те же 3%, а следом и площадь.
Особый случай — пар. Когда греющий агент конденсируется, работу совершает не «температурный перепад», а удельная теплота парообразования r (для пара 130 °C r = 2 174 кДж/кг). Тогда формула выглядит так Q = Gпар · r, а пар расходуется на 2–3 порядка меньше, чем вода на ту же мощность. Это принципиально важно учитывать в паровых ГВС-схемах.
Температурный напор по логарифму LMTD и эффективность ε
Когда мощность Q найдена, необходимо выяснить, с каким температурным напором аппарат передаёт её. Простое ΔTср = (Tгорср − Tхолср)/2 — слишком грубая оценка. Верная формула — логарифмический температурный напор LMTD (англ. Log Mean Temperature Difference), так как по длине аппарата ΔT изменяется нелинейно.
Если ΔT₁ ≈ ΔT₂ (отличаются меньше чем в 2 раза), LMTD можно с погрешностью < 4% заменить на среднеарифметическое (ΔT₁+ΔT₂)/2. Если же отношение больше 2 — обязательно логарифм. В контурах отопления с близкими входными температурами LMTD обычно составляет 7–12 K, для ГВС с холодной водопроводной водой — 18–28 K, для пар-вода — 60–90 K.
Коэффициент F для перекрёстного потока
В реальном пластинчатом теплообменнике идеальный противоток на 100% невозможен — потоки в каналах частично пересекаются. Поэтому в формулу площади добавляют поправочный F-фактор (коэффициент эффективности схемы движения), уменьшающий расчётный LMTD: LMTDэфф = F · LMTDтеор. У одноходовых пластинчатых аппаратов F = 0,95–0,99, у многоходовых (2-2, 1-2) F = 0,82–0,92, у перекрёстного потока в кожухотрубных F = 0,75–0,90 (принимается по номограммам HEI/TEMA).
Метод ε-NTU для поверочных расчётов
Другой подход — метод эффективности и числа единиц переноса теплоты (ε-NTU): ε = Qфакт / Qмакс, где Qмакс = Cmin·(Tгор.вх − Tхол.вх). Число NTU = U·F / Cmin. Он удобен для проверочного расчёта: «известны U, F, расходы — определить температуры на выходе». Применяется в SWEP SSP и Alfa Laval CAS для быстрых итераций при выборе модели.
Площадь теплопередающей поверхности F и количество пластин
Когда есть Q, LMTD, F-фактор и ориентировочный U (из таблиц или итераций), переходим к ключевому результату — площади теплопередающей поверхности, измеряемой в м². Формула:
К примеру, при Q = 200 000 Вт, U = 4 200 Вт/(м²·K), LMTD = 21 K, F = 0,97 получаем Fрасч = 200 000 / (4 200 · 21 · 0,97) = 2,34 м². Это площадь без запаса. С учётом 18% на загрязнение — 2,76 м². Это и будет необходимая поверхность нагрева.
Расчёт числа пластин N
Для разборного пластинчатого теплообменника: N = Fрасч / Fпластина + 1, где Fпластина — паспортная площадь теплопередачи одной пластины. Обычные значения по моделям:
Площадь одной пластины составляет от 0,04 м² (малые ЭТРА и Funke FP 04) до 2,0 м² (Alfa Laval M30 и WideGap). Чем больше пластина, тем меньше деталей при той же суммарной площади, но крупнее рама. Для нашего примера (2,76 м²) на РИДАН НН 14 понадобится N ≈ 2,76/0,14 + 1 = 21 пластина, для Alfa Laval M6-MFG — N ≈ 2,76/0,15 + 1 = 20 пластин. О выборе модели подробнее — в статье «пластинчатые разборные теплообменники».
Кожухотрубный: количество и длина труб
Для кожухотрубных аппаратов формула преобразуется: F = Nтр · π · dнаруж · L, где Nтр — количество трубок (как правило, 50–800), dнаруж — наружный диаметр (16, 19, 25 мм), L — длина (1,5–6 м). Стандарты TEMA AES, BEM, NEN. Толщина 1,5–2,5 мм. Для нашего примера ТР (2,76 м²) при трубке Ø16×1,5×3000 Nтр = 2,76 / (π·0,016·3) = ≈ 19 трубок — выполняется в наименьшем ТР-У с диаметром кожуха 159 мм.
Гидравлический расчёт Δp: потери напора в каналах
Тепловой расчёт определяет «требуемую площадь», а гидравлический отвечает на вопрос: «сможет ли заданный расход пройти через аппарат без избыточного перепада давления?». Если ΔP выше 80 кПа на горячей или холодной стороне, насос не справится — придётся брать модель крупнее или делать параллельные ходы. Если ΔP ниже 20 кПа, аппарат «недозагружен»: U окажется меньше расчётного, и мощность не достигнет паспортной.
У гофрированных пластин (chevron) ζ определяется углом гофры (β = 30°, 45°, 60°) и числом Рейнольдса. При Re выше 400 поток турбулентный, ζ ≈ (15…35)/Re0,25. В канале между пластинами скорость, как правило, составляет 0,3–0,8 м/с; эквивалентный диаметр равен 4·b (b — глубина гофры 2–6 мм). Каталоги производителей (РИДАН, Alfa Laval, ЭТРА) содержат графики ΔP(G) по каждой модели — их можно интерполировать.
Как поступать, когда ΔP выходит за пределы допустимого. Варианта три: (1) взять модель с более крупной пластиной (меньше скорость) — обычно увеличение габарита на 1 типоразмер снижает ΔP в 1,8–2,3×; (2) перейти на двухходовую схему 2-2 (но F-фактор уменьшится с 0,98 до 0,88); (3) поставить аппараты параллельно — дробление расхода в 2 раза снижает ΔP примерно в 4×. Подробно о выборе схемы — материал «как подобрать теплообменник».
Образец комплексного расчёта: ГВС 200 кВт, 8 этапов
Соединим все этапы на практическом примере. Исходные данные: ИТП жилого дома, ГВС 200 кВт. Контур нагрева (тепловая сеть): T1вх = 95 °C, T1вых = 70 °C. Контур, который нагревается (ГВС): T2вх = 5 °C (зимний городской водопровод), T2вых = 60 °C (норма СН РК 4.01-02-2013). С обеих сторон — вода. Подпитка из городской сети — Rf = 0,0003 м²·K/Вт.
Этап 1. Сбор и проверка ТЗ
Полный набор — 8 параметров: Q=200 кВт; T1вх=95 °C; T1вых=70 °C; T2вх=5 °C; T2вых=60 °C; горячая среда=вода; холодная среда=вода; рабочее давление 6 бар; Rf=0,0003. Проверяем: ΔTгор=25 K, ΔTхол=55 K. Поскольку Qгор=Qхол, расход воды по горячей стороне должен составлять 55/25 = 2,2 раза выше, чем на холодной.
Этап 2. Расход по двум контурам (Q=GcΔT)
Для горячей стороны: G₁ = Q/(c·ΔT) = 200 000 / (4 200·25) = 1,905 кг/с ≈ 6,86 м³/ч. Холодная сторона: G₂ = 200 000 / (4 195·55) = 0,867 кг/с ≈ 3,12 м³/ч. Эти величины используются для гидравлики и подбора циркуляционного насоса.
Шаг 3. LMTD и F-фактор
ΔT₁ = 95 − 60 = 35 K; ΔT₂ = 70 − 5 = 65 K; LMTD = (35 − 65)/ln(35/65) = (−30)/(−0,619) = 48,5 K. Для пластинчатого аппарата 1-1 принимаем F-фактор = 0,98 → LMTDэфф = 47,5 K.
Шаг 4. Оценка U
Для системы вода-вода в пластинчатом теплообменнике с гофрой 60°/30° и средней скоростью 0,5 м/с U₀ ≈ 4 500 Вт/(м²·K). С учётом Rf = 0,0003: 1/U = 1/4 500 + 0,0003 = 0,000522, отсюда U ≈ 3 830 Вт/(м²·K).
Этап 5. Площадь теплообмена и запас
Fтеор = 200 000 / (3 830 · 47,5) = 1,099 м². С запасом 18% на загрязнение и износ прокладок: Fрасч = 1,099 · 1,18 ≈ 1,30 м².
Шаг 6. Подбор модели
Этап 7. Проверочный гидравлический расчёт
У всех 4 моделей ΔP находится в пределах норматива 20–80 кПа с обеих сторон. Наименьший перепад у КС-25 (28/16 кПа), наибольший у Funke FP 14 (45/27). При слабом напоре циркуляционного насоса у заказчика — однозначно берём КС-25 или РИДАН НН 14.
Этап 8. Итоговая проверка по ε-NTU и оформление отчёта
Проверяем: NTU = U·F/Cmin = 3 830 · 1,54 / (0,867 · 4 195) = 1,62. По номограмме для противотока при Cr=0,455 — ε ≈ 0,712. Qфакт = ε · Cmin · (T1вх−T2вх) = 0,712 · 3 637 · 90 = 233 кВт в чистом виде; с учётом Rf уменьшается до паспортных 200 кВт — расчёт совпал. Результат: рекомендуется модель РИДАН НН 14 / 11 пластин AISI 316 / EPDM, стоимость примерно 145 000–185 000 ₸ в зависимости от партии. О ценах подробнее — в статье «стоимость пластинчатого теплообменника».
Расчётные программы: SWEP SSP, Alfa Laval CAS, KOPT, Sondex Pro
На практике вручную выполняют лишь «прикидочный» расчёт — основной же делают в специализированном ПО. Оно обеспечивает точность ±3–5% по площади и автоматически подбирают модель из библиотеки изготовителя.
Выбор программы зависит от бренда. Если объект проектируется под Alfa Laval — расчёт ведут в CAS. Под РИДАН — в Selector. Когда нужно сравнить 3–5 брендов в одной задаче, удобнее обратиться к инженеру нашей компании: мы выполняем параллельный расчёт во всех профильных программах и формируем сводное КП. О составе исходных данных подробнее — в статье «что требуется для расчёта теплообменника» и «расчёт оборудования для теплообмена».
Запас на загрязнение: Rf для разных сред
Аппарат, рассчитанный «впритык» под идеальную чистую воду, теряет 8–15% мощности уже за первый год эксплуатации и до 30% — за 5 лет. Чтобы аппарат всю свою службу выдавал паспортную Q, в расчёт вводят фактор загрязнения Rf (англ. fouling resistance). Стандарт — TEMA RGP-G-1.2, а также стандарты ГОСТ и СТ РК и СТ РК.
Как закладывать запас. Есть два равнозначных подхода: (1) включить R в формулу 1/Uf1 + Rf2 до вычисления площади — это «инженерный» путь, методологически более чистый; (2) определить F для идеальных условий и умножить на коэффициент 1,15–1,30 — «практический» путь, который чаще используют в КП. Оба варианта дают одинаковую итоговую площадь с точностью до 2%. Об устройстве и эксплуатации подробнее — в статьях «пластинчатые теплообменники» и «аппараты из нержавеющей стали».
6 распространённых ошибок в тепловом расчёте
Подобные ошибки постоянно встречаются в расчётах начинающих инженеров и в «расчётах» на скорую руку в Excel. Каждая объясняет, почему теплообменник «не выходит на режим» или почему спустя 2 года работы мощность снижается на 20–30%.
Принимают c=4 190 Дж/(кг·K) для воды «в целом», хотя для греющего теплоносителя при 95 °C это 4 213, а при 80 °C — 4 196. Погрешность 0,5%, но вместе с остальными — до 5% по площади. Особенно важно для гликоля, где c определяется концентрацией.
Берут LMTD = 47,5 K и используют без поправки. Если не учесть F=0,88 для многоходовой схемы, теплообменник получается на 12% меньше требуемого. Особенно рискованно для двухходовых пластинчатых и кожухотрубных схем TEMA 1-2.
Новый аппарат даёт расчётную мощность лишь в первые недели. Без Rf = 0,0003 для городского ГВС через 3 года площадь «съедается» на 22–28%, и инженер недоумевает, почему ёмкостный нагреватель «не справляется».
Q в кВт или Вт (×1000), c в кДж/(кг·K) или Дж/(кг·K) (×1000), расход в м³/ч или кг/с (×ρ/3600). Чаще всего получается «теплообменник в 1000 раз больше нужного» либо «в 1000 раз меньше». Выход — всё в СИ и аккуратно расписано в таблице.
У пара тепловая способность определяется не c·ΔT, а r (теплотой парообразования, ≈ 2 174 кДж/кг при 130 °C). Если рассчитать пар как воду по Q=GcΔT, расход пара окажется завышенным в 100–500 раз, и КП «выстрелит» под несуществующий котёл.
Площадь совпала, модель выбрана, однако ΔP = 110 кПа на горячей стороне — насос не «проталкивает». Аппарат либо недогружен по расходу (значит, U ниже расчётного), либо работает на грани кавитации. Гидравлический расчёт обязателен для каждого ТЗ.
Чтобы избежать этого — три простых шага: (1) применять аттестованные программы (CAS, SSP, KOPT), а не самодельный Excel; (2) обязательно делать перекрёстную проверку по ε-NTU; (3) согласовывать расчёт с производителем до запуска в производство. Если расчёт заказывают у нас — все три уровня контроля действуют по умолчанию.
FAQ — 10 ответов для инженеров и проектировщиков
1. Сколько времени занимает правильный тепловой расчёт?
Опытный инженер при наличии исходных данных и доступа к CAS/SSP/KOPT выполняет расчёт 30–60 минут: сбор ТЗ занимает 10 мин, основной расчёт и подбор — 15, гидравлика, проверка ε-NTU и оформление отчёта — 20. При пробелах в ТЗ (не указана точная концентрация гликоля, отсутствует Rf), срок увеличивается вдвое на уточнения.
2. Какой минимум из 8 параметров требуется для расчёта?
(1) Q — мощность тепловая; (2) T1вх — горячий контур на входе; (3) T1вых — выход горячего; (4) T2вх — вход холодного; (5) T2вых — холодный контур на выходе; (6) среда (вода/гликоль/пар/масло); (7) давление рабочее; (8) фактор загрязнения Rf либо регион и тип воды. Детальный опросный лист приведён в статье «что нужно для расчёта».
3. Допустимо ли подобрать аппарат «по аналогу», не выполняя полный расчёт?
Допустимо — при условии, что ТТХ аналога совпадают 1:1 (мощность, температуры, среды, давления и материалы те же), и он действительно эксплуатируется. Однако в 70% случаев «аналог» выбран не лучшим образом: или запас 200%, или недостаток 15%. Надёжнее выполнить расчёт заново — за 30–60 минут получите технически корректный подбор, а не аргумент «раньше же работало».
4. В чём отличие теплового расчёта ГВС от отопительного?
В ГВС холодная вода поступает из городского водопровода (T2вх=5–15 °C), ΔT большой (50–60 K), а R обязателенf=0,0003–0,0005. Для отопления в отоплении — вода закрытого контура с небольшим ΔT (15–25 K), Rf=0,00005. При одинаковой Q аппарат ГВС всегда на 30–60% крупнее.
5. В каких случаях выбирают паяный, а в каких — разборный пластинчатый?
Паяный вариант — для чистых сред, давления до 30 бар, мощности до 800 кВт, частного дома или небольшой котельной. Разборный — для городского ГВС, мощностей более 200 кВт, когда нужна чистка или переборка пакета. Сопоставление по 10 параметрам — в материале «паяные vs разборные».
6. Какой аппарат предпочтительнее: пластинчатый или кожухотрубный?
Пластинчатый в 3–5× компактнее и в 2–3× эффективнее по U, однако не терпит загрязнённых сред и резких ΔP. Кожухотрубный занимает много места и имеет большую массу, но пригоден для пара, масла, конденсата и грязной воды. Подробное сравнение — в статье «кожухотрубный либо пластинчатый».
7. Откуда брать надёжные значения U для предварительной оценки?
Доступные источники: справочники Cao «Heat Transfer in Process Engineering», Kakaç-Liu «Heat Exchangers», ASHRAE Handbook (HVAC Systems), HEI Standards. Для условий Казахстана — РТМ 26-01-72-82, стандарты ГОСТ и СТ РК и СТ РК. В программах CAS/SSP/KOPT таблицы U уже включены.
8. Требуется ли расчёт по ASME для объектов в Казахстане?
На объектах в Казахстане обязательны ТР ТС 010/2011 «О безопасности машин и оборудования» и ТР ТС 032/2013 для аппаратов с давлением более 0,5 МПа. Пластинчатые охватывает стандарты ГОСТ и СТ РК и СТ РК. ASME Section VIII Div.1 требуется для экспорта в США и для крупных нефтехимических установок.
9. Как поступить, если часть из 8 параметров в ТЗ отсутствует?
Уточнить у заказчика. Когда объект уже построен — поднять проектную документацию либо выполнить замеры манометрами и термометрами на работающем контуре. Наиболее «обходной» вариант — применить типовые значения для похожего объекта (к примеру, для жилого МКД в 9 этажей мощность ГВС принимают 25–35 Вт/чел или 1,5–2,5 кВт/квартиру). См. «нагрузки ИТП жилого дома».
10. На кого ложится ответственность при ошибочном расчёте?
Когда расчёт выполняет поставщик, по договору он отвечает за выход аппарата на паспортную Q в условиях из опросного листа. Если расчёт делала проектная организация с СРО — ответственность на ней. Если «инженер заказчика без подписи» — никто, а рисков нет лишь у согласовавших. Поэтому важно, чтобы расчёт был оформлен документом с подписью инженера и печатью поставщика или ПИ.
Связанные материалы
Тепловой расчёт за 30–60 минут от инженеров ТПЛ-Сервис
Направьте опросный лист либо назовите 8 параметров — инженер за 30–60 минут выполнит расчёт в SWEP SSP, Alfa Laval CAS и Kelvion KOPT, сопоставит 3–5 моделей по площади, ΔP и стоимости, подготовит сводное КП. Расчёт ведётся по ТР ТС 010/2011, ТР ТС 032/2013, стандарты ГОСТ и СТ РК и СТ РК. По запросу — ASME Section VIII.
- Комментарии
