Новый подход к
моделированию полей
Языковые версии сайта:
Language no-Pyccku Global English Deutsch Espanol Francais Italiano Danmark Ceske Chinese

>> >> >>

Естественная конвекция с поверхности вертикальной пластины

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

Тип задачи:
Плоская задача стационарной теплопередачи.

Геометрия:

Высота пластины l = 1000 мм, толщина пластины d = 2 мм, площадь поверхности пластины 1 м2.

Дано:
температура воздуха слева от пластины T1 = 20°C,
температура воздуха справа от пластины T2 = -10°C,
теплопроводность пластины λпласт = 40 Вт/(м·K),
теплопроводность наружного воздуха λвозд = 0.027 Вт/(м·K).

Задание:
Определить температуры левой и правой поверхностей пластины и тепловой поток через них.

Решение:

Тепловой поток идет от более тёплой области к более холодной (в нашей модели слева направо через стенку). Тепловой поток внутри стенки передаётся за счёт теплопроводности. Тепло передается от стенки к воздуху за счёт конвекции. Коэффициенты конвекции зависят от многих параметров, включая температуру среды, форму и ориентацию поверхностей.
Для вертикально ориентированной пластины средний коэффициент конвекции можно найти, используя теорию подобия. Температуры левой и правой поверхностей пластины для расчета коэффициентов конвекции будем оценивать равными средней температуре пластины, которую можно принять за Tпл = (T1+T2)/2 = 5°C. Модель представляет собой задачу теплопроводности в материале пластины с граничными условиями свободной конвекции на поверхностях пластины. Решение задачи проводится в следующем порядке:

  1. Число Прандтля (Pr), число Грасгофа (Gr), число Релея (Ra) и число Нуссельта (Nu) определяем по известным формулам*:
         Pr = μ·C / λвозд, где μ - динамическая вязкость воздуха, C - теплоемкость воздуха;
         Gr = L3ρ2g·ΔT·β / μ2, где L - характерный размер (в данном случае равен l), ρ - плотность воздуха, g - ускорение свободного падения, ΔT - разность температур, β - коэффициент теплового расширения воздуха;
         Ra = Gr·Pr;
         Nu = [0.825 + 0.387·Ra1/6 / (1 + (0.492/Pr)9/16)8/27]2.

    Параметры воздуха ρ, C, μ, β рассчитываются для средних температур:
    левая поверхность: Tср1 = (T1 + Tпл)/2 = (20+5)/2 = 12.5 °C;
    C1 = 1000 Дж/(кг·K), μ1 = 1.87·10-5 Н·с/м2, ρ1 = 1.25 кг/м3, β1 = 0.003501 1/K;
    правая поверхность: Tср2 = (T2 + Tпл)/2 = (-10+5)/2 = -2.5 °C**
    C2 = 1000 Дж/(кг·K), μ2 = 1.87·10-5 Н·с/м2, ρ2 = 1.32 кг/м3, β2 = 0.003695 1/K.

    Численные значения параметров:
    ΔT1 = 20 - 5 = 15 °C, ΔT2 = 5 - (-10) = 15 °C;
    Pr1 = 0.69; Pr2=0.69;
    Gr1 = 2.30·109, Gr2 = 2.71·109;
    Ra1 = 1.59·108, Ra2 = 1.88·108;
    Nu1 = 141, Nu2 = 149.

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

  2. Определяем средний коэффициент конвекции на поверхностях пластины
    α = Nu·λвозд /L
    Подставив значения, получаем средние коэффициенты конвекции с поверхностей пластины
    α1 = 3.82 Вт/(м2·K); α2 = 4.02 Вт/(м2·K).
  3. Решаем задачу в ELCUT и находим температуры поверхностей пластины и поток тепла через эти поверхности.

Результат:
Температура левой поверхности пластины: T1 = 277.77 K (4.62 °C).
Температура правой поверхности пластины: T2 = 277.76 K (4.61 °C).
Поток тепла: F = 58.756 Вт (площадь пластины 1 м2).

Загрузить файлы задачи.

*CIBSE Guide C: Reference Data (2010). Butterworth-Heinemann, ISBN: 0750653604. Таблица 3.5.

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

Моделирование высоковольтных систем в ELCUT


Карта сайта