Расчет труб отопления: Подбор диаметра труб для отопления по мощности контура

Уравнение теплопередачи

Уравнение теплопередачи

Уравнение теплопередачи

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

Уравнение теплового баланса

Уравнение энергии для несжимаемой жидкости, текущей в трубе (ссылка 24):

(3-1)

, где ρ — плотность жидкости (единица СИ: кг/м3), A — площадь поперечного сечения трубы (единица СИ: м2), доступная для потока, Cp (единица СИ: Дж/(кг·К)) — теплоемкость при постоянном давлении T (единица СИ: K) — это температура. u — поле скоростей. Для получения информации о тангенциальной скорости потока в трубе см. раздел Теория интерфейса потока в трубе. Кроме того, k (единица СИ: Вт/(м·K)) представляет собой теплопроводность. Второй член в правой части соответствует теплоте трения, рассеиваемой за счет вязкого сдвига. Q (единица СИ: Вт/м) представляет собой общий источник тепла, а Qwall (единица СИ: Вт/м) представляет внешний теплообмен через стенку трубы.

Обратите внимание, что термин Qwall подробно описан ниже.

В правую часть уравнения можно добавить дополнительный член Qp, установив флажок Работа под давлением:

(3-2)

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

Теплообмен через стену

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

(3-3) (Вт/м)

В уравнении 3-3 (hZ)eff – эффективное значение коэффициента теплопередачи h (единица СИ: Вт/(м2·К)) умноженное на периметр стенки Z (единица СИ: м) трубы. Текст (единица СИ: K) внешняя температура снаружи трубы. См. Рисунок 3-5. Qwall появляется как исходный член в уравнении теплопередачи трубы, уравнение 3-1.

Для функции теплопередачи стены требуется внешняя температура и, по крайней мере, добавленный к ней подузел сопротивления внутренней пленки. Отдельные вклады коэффициентов теплопередачи могут быть добавлены подузлами к функции «Теплопередача стены». Подузлы:

Сопротивление внутренней пленки

Настенный слой

Внешнее пленочное сопротивление

Текст в уравнении 3-3 может быть константой, параметром, выражением или заданным полем температуры, вычисленным другим физическим интерфейсом, обычно интерфейсом 3D Heat Transfer. h автоматически рассчитывается через сопротивления пленки и слои стен, которые добавляются как подузлы; см. уравнение 3-17 и далее.

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

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

Рис. 3-5: Распределение температуры по стенке трубы.

rn (единица СИ: м) — внешний радиус стены n

w = r − r0 (единица СИ: м) координата стены, начиная с внутреннего радиуса r0

Δwn = rn − rn−1 (единица СИ: м) толщина стенки n

Zn (единица СИ: м) — внешний периметр стены n

.

hint и hext — коэффициенты пленочного теплообмена внутри и снаружи трубы соответственно (единица СИ: Вт/(м2·K)).

кн – теплопроводность (единица СИ: Вт/(м·К)) стены n

Раковина весов

На рис. 3-5 рассмотрим короткий отрезок ΔL трубы, перпендикулярный плоскости рисунка. Тепло, уходящее из внутренней жидкости этого сегмента в стенку, равно

.

(3-4)(Ш)

Здесь AQ =  ΔL2πr0 (единица СИ: м2) — площадь, доступная для теплового потока в стену. Для стационарных условий такое же количество тепла должно пройти через любую цилиндрическую оболочку радиусом r в стене 1 (или любой стене).

(3-5)

Переставить и интегрировать от r0 до r1.

(3-6)

Выполнить интеграцию

(3-7)

и переставить

(3-8).

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

(3-9)

Замена

(3-10),

и составление линейной комбинации уравнений Уравнение 3-9 дает

(3-11)

, где (hAQ)eff — эффективная проводимость:

(3-12)

Для общего случая с N слоями стенки это выглядит как

.

(3-13)

Теперь пусть

(3-14),

, где Z (единица СИ: м) — средний периметр (окружность) трубы, взятый по толщине стенок трубы. Объедините уравнение 3-10 и уравнение 3-14 так, чтобы

и вставить в уравнение 3-13:

(3-15)

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

.

одинаковая температура по окружности

теплопередача через стенку квазистатическая. Последнее означает, что предполагается, что стена сразу принимает равновесное распределение температуры, соответствующее T и Text. Если бы это предположение не было сделано, потребовалось бы вспомогательное УЧП по координате стены.

Для труб квадратной и прямоугольной формы средняя проводимость может быть аппроксимирована более простой суммой сопротивлений на плоской стенке, которую можно найти, например, (ссылка 13):

(3-16)

Сопротивление пленки можно рассчитать по

(3-17)

, где k — коэффициент теплопроводности материала, а Nu — число Нуссельта. dh — гидравлический диаметр, определяемый как

(3-18).

Коэффициенты внутренней и внешней пленки оцениваются как (T + T0)/2 и (TN + Text)/2 соответственно.

Коэффициент теплопроводности kn может зависеть от температуры и оценивается как (Tn + Tn−1)/2.

Чтобы вычислить dh в уравнении 3-18, локальный периметр рассчитывается как Z = f(w) и площадь поперечного сечения как A = f(w). Автоматические расчеты для круглых труб выполняются физическим интерфейсом, как Z = 2πr и A = πr2. Для прямоугольных труб это Z = 2 (ширина + высота), а A = ширина · высота. Для определяемых пользователем форм труб пользователь может вводить произвольные выражения.

Локальные температуры в каждом радиальном положении стенки трубы (см. Рисунок 3-5) рассчитываются с учетом того факта, что уравнение 3-3 также может применяться для каждого отдельного слоя стенки:

(3-19)

Комбинация уравнения 3-3, уравнения 3-19 и уравнения 3-15 или уравнения 3-16 (в зависимости от формы трубы) для каждого слоя стенки дает явное значение Tn.

Сопротивление внутренней пленки

Для внутренней ламинарной вынужденной конвекции при полностью развитом потоке в трубе число Нуссельта является константой, зависящей от поперечного сечения трубы. Значения перечислены в таблице ниже (ссылка 1). Интерфейс Pipe Flow выполняет интерполяцию, чтобы найти значения отношения ширины/высоты, не указанные в списке. Настройки по умолчанию для расчета коэффициента пленки — «Автоматически», что означает, что ламинарная и турбулентная корреляции применяются в соответствии с числом Re.

Таблица 3-1: Числа Нуссельта для внутреннего ламинарного течения в трубе

Сечение

ширина/высота

Ню

циркуляр

3,66

квадрат

1

2,98

прямоугольный

1,43

3,08

прямоугольный

2

3,39

прямоугольный

3

3,96

прямоугольный

4

4,44

прямоугольный

8

5,60

параллельные пластины

бесконечность

7,54

Для поперечных сечений, определяемых пользователем, по умолчанию предлагается значение Nu 3,66.

Для внутренней турбулентной вынужденной конвекции (3000 < Re <  6·106, 0,5 < Pr <  2000) применяется уравнение Гнелинского (ссылка 18):

(3-20)

Где Pr — число Прандтля:

(3-21)

Сопротивление пленки внутреннему потоку можно рассчитать, используя свойства материала, определенные в функции «Теплопередача», и рассчитанный коэффициент трения. Свойства материала оцениваются при средней внутренней температуре пленки (T + T0)/2 (см. Выражения для коэффициента трения Дарси).

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

Внешнее пленочное сопротивление

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

Для внешней вынужденной конвекции вокруг трубы, справедливой для всех Re и для Pr > 0,2, Черчилля и Бернштейна (ссылка 19) используется соотношение:

(3-22).

Для внешней естественной конвекции вокруг трубы используется корреляция Черчилля и Чу (ссылка 20), которая рекомендуется для Ra < 1012:

(3-23)

, где число Рэлея задается как:

(3-24)

и номер Грасгофа:

(3-25)

Выше d — внешний диаметр трубы, а β — коэффициент объемного теплового расширения жидкости:

(3-26)

Свойства материала оцениваются как (TN + Text)/2.

Стабилизация уравнения теплопередачи

Уравнение переноса в интерфейсе «Теплообмен в трубах» численно стабилизировано.

Численная стабилизация в справочном руководстве COMSOL Multiphysics

Калькулятор производительности тепловых труб — как использовать наш бесплатный онлайн-инструмент

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

Вот ссылка на онлайн-калькулятор.

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

ВХОДЫ

Длина тепловой трубы – полная длина тепловой трубы, если испаритель находится на одном конце.

Длина испарителя – длина источника тепла.

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

Тепловая трубка, тип — выбор фитиля «Стандартный» или «Эффективный». Стандартные фитили позволят вам сплющить тепловую трубку, прежде чем повлиять на максимальную мощность (Qmax). Обратите внимание, что мы можем изменить пористость и толщину фитиля, чтобы точно соответствовать требованиям вашего применения.

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

РЕЗУЛЬТАТЫ

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

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

Вот несколько советов по дизайну, которые вам пригодятся. Во-первых, тепловая трубка Qmax является аддитивной при условии, что каждая тепловая трубка (трубки) находится над источником тепла. В нашем примере это означает, что одна 8-миллиметровая тепловая трубка имеет Qmax 62 градуса Цельсия в горизонтальном положении, а две 8-миллиметровые трубки имеют Qmax 124 градуса. C. Во-вторых, предусмотрите коэффициент безопасности, чтобы избежать эксплуатации трубы на максимальной мощности. Снижение мощности тепловой трубки Qmax на 20 % является хорошим отраслевым стандартом — в этом случае мощность одной трубки будет чуть менее 50 Вт. Случайные короткие скачки мощности выше этого допустимы, если они все еще ниже номинального значения Qmax.

График зависимости мощности от дельты-T и таблицу (показанные ниже) необходимо использовать вместе с приведенным выше графиком. Допустим, мы выбрали 8-миллиметровую тепловую трубку, о которой говорили ранее: в горизонтальном положении она будет безопасно нести чуть менее 50 Вт (после снижения номинальных характеристик). Если мы подадим 40 Вт на один конец, другой конец будет иметь дельта-Т 4,3 ° C (чем ниже, тем лучше). Предположим, мы хотим использовать две тепловые трубки диаметром 8 мм.

Добавить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *