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

Ключевые слова:
дом из дерева, экологическая чистота, нестационарная теплопередача, периодический тепловой режим, энергетическая эффективность
Текст
Текст произведения (PDF): Читать Скачать

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

Одним из вариантов пассивного управления тепловыми потоками является использование стены Тромба [4], предполагающей применение каменных и бетонных массивов для обеспечения тепловой инерции конструкции и выравнивания внутренних температурных колебаний в течение суток. Более простым способом получения теплого дома является утепление стены с деревянным каркасом [5]. Существует широкий набор легкодоступных теплоизолирующих материалов [6], которые просто устанавливаются в полость между внутренней облицовкой и внешней облицовкой. Их типичным недостатком является недолговечность и необходимость замены через 10–15 лет эксплуатации. В связи с этим, особый интерес среди экологически чистых решений в рамках концепции энергетически эффективного жилища вызывают дома из деревянных конструкций [7–9]. Использование древесины, как и любого другого материала при проектировании ограждающих конструкций, требует проведения теплофизических расчетов для определения условий обеспечения комфортного проживания.

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

При теплообмене в ограниченном пространстве явления нагревания и охлаждения газа протекают вблизи друг от друга, и разделить их не всегда возможно. Вследствие ограниченности пространства и наличия восходящих и нисходящих потоков, условия движения усложняются, поскольку они зависят не только от температурного напора, но и от формы и геометрических размеров обтекаемых тел [11, 12].

Механизм конвекции существенно зависит от формы и размеров ограждения рассматриваемой ячейки. Такая конвекция в ограниченном объеме помещения, вызванная только естественными силами (без вынужденной вентиляции), называется естественной. Одним из стандартных объектов исследований теплообмена в помещениях является модельная ячейка в форме параллелепипеда с различным расположением тепловых источников, в том числе периодических во времени [13]. Соответствующие теоретические и экспериментальные исследования показывают необходимость учета конкретных условий теплообмена на границах и наибольшую роль числа Рэлея в качестве определяющего параметра [14-17].

 

 

Рис. 1. Тепловые и радиационные потоки в помещении зимой [1]

 

 

При проведении конкретных расчетов требуются данные о тепловых характеристиках воздуха [18], которые описываются в интервале температур 300–350 K эмпирическими соотношениями [16] (в СИ): динамическая вязкость, плотность, теплопроводность, теплоемкость, кинематическая вязкость. Для экспериментального изучения конвективного теплообмена используются тестовые помещения также чаще всего в форме параллелепипеда с линейными размерами в несколько метров [14]. Выбор размеров определяется их близостью к размерам реальных помещений. При размерах помещения меньше ~1 м не будет проявляться турбулентная стадия конвективного теплообмена, в то время как существенно большие модели не только труднее изготовить, но в них не будет наблюдаться эффект стеснения конвективного потока в ограниченных объемах, характеризующий отличие естественной конвекции от свободной конвекции. В реальных помещениях идеальная картина свободного конвективного теплообмена нарушается влиянием ограничения объема и наличием нескольких холодных и нагретых поверхностей теплообмена. Естественная конвекция в помещении вблизи нагретой стенки характеризуется наличием трех выраженных зон [19]: зоны ламинарного движения, зоны турбулентного движении и зоны торможения вблизи потолка. С учетом растущего спроса на дома из дерева, представляет значительный практический интерес развитие эффективных инженерных методов расчета их теплового режима. Важнейшим при этом является правильный учет одновременно инсоляции и конвективного теплообмена.

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

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

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

Материалы и методы. Для проведения конкретных расчетов теплообмена помещения в рамках сформулированной модели необходимо задать ряд геометрических и физических параметров, представленных в таблице 1. Следует учитывать, что материал дерева является сильно анизотропным, и теплопроводность древесины вдоль волокон выше на 40 %–90 %, чем поперек волокон, в зависимости от породы дерева. Надежные данные по теплофизическим свойствам древесины содержатся в справочнике [20]. Рассмотрим режим поддержания в помещении постоянной температуры 22 °С. В этом случае основной расчетной величиной является расход энергии на поддержание заданной температуры непосредственно в данном помещении. Остальные смежные помещения будем считать находящимися при такой же температуре за счет собственного кондиционирования воздуха. Основные параметры, необходимые для теплофизического расчета, приведены в таблице 1.

 

Таблица 1

Основные параметры расчетной модели

Наименование параметра

Числовое значение параметра

Теплоемкость воздуха при постоянном давлении cp

1,005 кДж/(кг·K)

Удельная плотность воздуха ρ

1,2041 кг/м3

Объем помещения V

81,9 м3

Коэффициент теплообмена внутренней поверхности стены αw

3(стена противоположная стене с остеклением)

Коэффициент теплообмена внешней поверхности стены αe

4.2 Вт/(м2·К)

Коэффициент теплообмена поверхности пола αf

3 Вт/(м2·К)

Коэффициент теплообмена поверхности светового пятна αs

4,2 Вт/(м2·К)

Коэффициент теплообмена поверхности потолка αp

2,8 Вт/(м2·К)

Размеры остекленной стены (высота* ширину):

(2 м х 1,7 м)х2

Тепловое сопротивление единицы площади поверхности остекления R

0,55 (м2·К/Вт)

Коэффициент воздухопроницаемости остекления β

5 кг/(м2·ч)

Максимальное значения светового потока для светового пятна q0  (Вт/м2)

980 Вт/ м2

Коэффициент поглощения света деревом ε

0,6

Размеры стен, смежных с остеклением (высота* ширина):

(2,8 м х 4 м)х2

Размеры светового пятна на полу:

0,9 м2

Время освещения светового пятна t в ч.

4,5 ч

Толщина стен L

0,2 м

Толщина пола Lf

0,2 м

Толщина потолка Lp

0,2 м

Теплоемкость материала стены cw

1,56 кДж/(кг·К)

Удельная плотность материала стены ρw

540 кг/ м3

Теплопроводность материала стены λw

0,15 Вт/(м·К)

Теплоемкость материала пола cf

1,56 кДж/(кг·К)

Удельная плотность материала пола ρf

540 кг/м3

Теплопроводность материала пола λf

0,15 Вт/(м·К)

Теплоемкость материала потолка cp

1,56 кДж/(кг·К)

Удельная плотность материала потолка ρp

540 кг/м3

Теплопроводность материала потолка λp

0,15 Вт/(м·К)

 

 

В качестве примера приведем расчет теплового режима в один из жарких дней июля 2021 года. По данным на 14.07.2021 погодного портала [21] максимальная температура воздуха составляла 33°С в 17:00, а минимальная температура 16°С в 5:00. Соответствующий суточный график температурной зависимости для воздуха от времени представлен на рис. 2. Средняя температура составляет T0=24,5 °С , амплитуда колебаний T=8,5 °С . Для определения потока тепла, поступающего в помещение, необходимо учесть доминирующий теплообмен через наружную стену и окна. В качестве исходного приближения будем считать, что поступающая через окно лучистая энергия поглощается фасадом с коэффициентом ε . Коэффициент пропускания света трехслойным стеклопакетом примем в соответствии с ГОСТ 11-90 равным 73,4%. Отметим, что надежность определения коэффициентов теплообмена поверхностей зависит, в том числе, от скорости ветра, а сами их величины могут значительно отличаться в различных источниках [22]. 

 

 

Рис. 2. Зависимость температуры воздуха по шкале Цельсия в Воронеже 14.07.2021 от времени

 

 

На рис.3 представлено распределение температур в помещении, полученное с использованием инфракрасного тепловизора. Оно подтверждает сделанное нами предположение о наличии больших плоских поверхностей с почти постоянной температурой и переходных областей с шириной порядка толщины пола. Так, температура холодной части пола и окна отличается примерно на градус, в то время, как неоднородность температуры в световом пятне составляет около 0,1 °С. Это позволяет далее придерживаться при численном расчете принятой нами составной одномерной модели, в которой теплообмен наружного ограждения помещения складывается из теплообмена через стену и теплообмена через окна с учетом солнечной радиации. Величина светового потока взята по результатам экспериментальных измерений [23] равной 980 Вт/м2. Восход солнца происходил в 4:22, а заход солнца – в 20:36. Астрономический полдень соответствовал t0= 12:29.

Теплофизическая характеристика стеклопакетов взята из работы [24] для двухкамерного стеклопакета в виде полного термического сопротивления R=0,87  м2·К/Вт с учетом конвективного и радиационного теплообмена. Полный тепловой поток через наружное ограждение находится как сумма тепловых потоков через стену и остекленное пространство. При этом тепловая инерция стеклопакета мала, и мы ее не учитываем, рассчитывая плотность теплового потока через остекление по формуле (19). Тепловая инерция стены учитывается с помощью наложения тепловых волн, как это описано выше для нестационарной теплопередачи на основе нестационарного уравнения теплопроводности.

Большая часть практических теплофизических расчетов зданий проводится в рамках одномерных [25, 26] или двумерных стационарных и нестационарных моделей [27-32]. Для определения точности расчетов в рамках одномерной модели мы рассмотрели решение стационарной двумерной задачи о распространении тепла в однородной плоской панели толщиной L . Температура на одной поверхности панели с координатой y=0  поддерживается постоянной и равной T0 , а на другой поверхности с координатой y=L  она отлична от T0  только на конечном интервале 0<x<A , где равна T . Трехмерная реализация такой двумерной постановки задачи означает, что область повышенной температуры на одной стороне панели имеет вид длинной полосы с постоянной шириной A.

 

 

 

Рис. 3. Тепловой портрет стеклопакета и светового пятна на деревянном полу (изображение с поворотом на 90)

 

 

Замены переменных y/L , x/L , u(x,y)=(T(x,y)-T0)/T , приводят исходное уравнение для распространения пятна нагрева к задаче Дирихле, решение которой имеет вид [33]:

ux,y=12-1πarctge-x-cosysiny           (1)

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

Распространение тепла в плоском слое теплопроводящей среды описывается одномерным уравнением теплопроводности [34]:

T∂t=χ2Tx2 , χ=λ                 (2)

Здесь λ  – теплопроводность среды, c  – удельная теплоемкость, ρ – объемная плотность. Уравнение (2) решается с учетом начальных и граничных условий.  Они зависят от условий теплообмена с окружающей средой [35]. Учтем, что на границе конструкции поглощается дополнительный тепловой поток с плотностью qs(t) , вызванной солнечным излучением. Тогда теплообмен на поверхности ограждения учитывается с помощью граничного условия:

αTextt-T-0,t=λ∂T∂xx=-0+qs(t)      (3)

Здесь Textt  может быть температурой наружного воздуха для внешней границы конструкции или температурой внутри помещения, λ  – коэффициент теплопроводности материала ограждения.  Поскольку коэффициент теплообмена  α  сам является нелинейной функцией температуры, то в последовательной постановке задача расчета теплообмена является нелинейной. Мы линеаризуем ее, считая известным диапазон изменения температур и проводя усреднение по времени с учетом гармонического изменения температуры в заданных пределах.

В качестве примера приведем расчет среднего значения эффективного коэффициента теплообмена  αeff  вертикальной стены высотой h=3 м из условия равенства тепловых потоков за полупериод для переменного (зависящего от температуры) и постоянного эффективного коэффициентов теплообмена:

0π/2αφsinφ=0π/2αeffsinφ.    (4)

Здесь фаза колебания φ=ωt . Температуру воздуха примем равной T0=22°С . Возьмем теплопроводность воздуха равную 𝜆=0.0268 Вт/(м·К), а его кинематическая вязкость ν=16.05·10-6 м2/с. Основной характеристикой самой конвективной теплопередачи является число Рэлея Ra=gѲL3Tνχ  , где g – ускорение свободного падения, L  – характерный размер задачи, Ѳ  - температурный перепад, T – температура воздуха, ν  – его кинематическая вязкость, χ  – температуропроводность.  Формула для расчета коэффициента теплопередачи имеет вид [36]:

α=0,825+0,3245Ra1/62 λL .        (5)

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

αeff=0π/4αφsinφ.          (6)

На рис. 2 показана зависимость отношения СѲ=αeff/α(Ѳ)  от амплитуды колебаний температуры Ѳ . Тем самым коэффициент СѲ  показывает насколько эффективный коэффициент конвективного теплообмена меньше его максимального значения при гармоническом законе изменения температуры. Для дальнейших расчетов мы примем C=0,607  и будем считать граничные условия линейными с заменой коэффициента теплообмена его средним значением.

Следствием линейности уравнения теплопроводности и принятых нами приближенных линейных граничных условий является принцип суперпозиции температурных возмущений [35]. Используя его, представим задачу в виде суммы двух более простых задач. Первая задача является стационарной и учитывает все стационарные составляющие в граничных условиях. Вторая задача описывает установившийся колебательный режим с учетом меняющейся со временем температуры наружного воздуха и потока солнечного излучения. Квазистационарный режим в слое толщиной L  устанавливается спустя некоторое время переходного процесса t , определяемое коэффициентом Фурье Fo=λtL2~1 , за которое начальное температурное распределение переходит в стационарное распределение.

Стационарная часть расчета производится с учетом перепада температур внутри и снаружи панели и условия постоянства потока тепла на границе. Решение ищется в виде линейной функции  Tx=ax+b . Подставляя это решение в граничные условия, получим систему линейных уравнений для определения постоянных коэффициентов a,b . Для решения нестационарной части задачи при Textt=Tcos(ωt)  перейдем к комплексному представлению [48] и ищем решение в виде комплексной функции u(x) :

θx,t=uxeiωt,                    (7)

где u(x)  – пространственная часть тепловой волны, ω  – частота теплового возмущения. Пространственная часть решения внутри однородной панели имеет вид

ux=Aexp-γx+Bexpγx.        (8)

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

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

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

qSt=kq0sinπt-t0  τ,               (9)

где k  – коэффициент пропускания остекления, q0 – амплитуда плотности светового потока, t0  – начальный момент освещения, τ  – длительность освещения. Формула (9) описывает временную зависимость светового потока на интервале полупериода (t0, t0+τ) , за пределами которого поток считается равным нулю. Ее параметры будут далее считаться заданными, как исходные для расчетной модели [39-41]. Кроме того, мы будем считать заданными углы падения на поверхности пола и стен, а также размеры световых пятен. Эти геометрические параметры потребуются для определения величины конвективного теплообмена нагретых солнцем участков.

Суточный ход температуры воздуха [42] хорошо воспроизводится гармонической зависимостью:

 Textt=T0+Tcosωt.         (10)

Что касается функции (9), то она задана на конечном интервале и не является гармонической, поэтому прямое применение изложенного выше метода гармонических волн невозможно. Для его использования необходимо сначала разложить функцию (9) в ряд Фурье по гармоническим функциям с периодом в t=2π/ω  (сутки):

qSt=-∞cnexpiωnt.        (11)

где

cn=1t0tqStexp-iωntdt, n=0,±1,±2 . (12)

Для временной зависимости в виде (9) получим:

cn=kq0tt0 t0+τ sinπt-t0 τ expiωntdt .   (13)

Вычисление этого интеграла приводит к выражению:

cn=kq0t2-β2exp-t0(exp-iβτ+1) ,  (14)

где введено обозначение β=ωn .

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

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

Таким образом, нагрев воздуха внутри помещения мы будем учитывать с помощью уравнения теплового баланса, выражающего закон сохранения энергии. Сумма тепловых потоков равна скорости изменения тепловой энергии массы M  воздуха в помещении, так что [43]:

cpMdTdt=jQj ,                      (15)

Здесь cp  – теплоемкость воздуха, T  – температура воздуха в помещении, Qj  – тепловые потоки, идущие внутрь и наружу помещения через ограждающие поверхности, а также внутренние тепловые источники. Тепловые потоки через ограждающую поверхность внутри помещения можно записать в виде:

Qj=Sjqj+sjpj ,                 (16)

где qj  – плотность теплового потока на границе освещенной солнцем участка поверхности ограждения с номером j  и площадью Sj , pj  – плотность теплового потока на границе неосвещенной солнцем участке поверхности ограждения с номером j  и площадью sj . Нагретые солнцем и теневые участки ограждения имеют разную температуру поверхности, которую мы обозначим как TjS(0)  и Tjs(0)  соответственно. Плотность тепловых потоков внутри помещения определяются разностью температур поверхности стены и температуры воздуха в помещении:

qj=αjS(TjS0-T) ,                (17)

pj=αjs(Tjs0-T) .                (18)

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

qG=Ѳ(1+β)/R ,            (19)

где β  – коэффициент, учитывающий инфильтрацию воздуха.

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

Результаты соответствующего расчета теплового потока в течение суток представлены на рис.4.

 

 

 

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

 

 

Из полученной зависимости видно, что наиболее резко поток энергии нарастает в период прямого солнечного освещения помещения через оконный проем. Полная компенсируемая за счет кондиционирования воздуха тепловая энергия за сутки составляет 3,57 кВтч. С учетом холодильного коэффициента охлаждения 3, затрачиваемая электрическая мощность составляет около 1,2 кВтч. С учетом тарифов для садового некоммерческого товарищества 2,55 руб. за 1 кВтч, суточная стоимость составит около 3 руб. для одной комнаты. Для сравнения нами проведен расчет для случая, когда солнечное излучение не падает на стеклопакет, а отклоняется защитным экраном. Результаты расчетов приведены на рис. 5.

 

 

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

 

 

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

На основании проведенных расчетов можно сделать вывод, что правильные мероприятия по затенению оконных проемов с солнечной стороны здания позволяет существенно уменьшить расходы на кондиционирование воздуха в помещениях. Как видно из результатов расчетов, эффективный коэффициент теплообмена всегда ниже максимального значения, и их отношение остается почти постоянным (рис. 6). Получено значение среднего за период температурных колебаний коэффициента конвекции 0,607, отличается от имевшегося ранее значения 0,64 [44], с учетом его нелинейной зависимости от разности температур.

 

 

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

 

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

  1. Показано, что для описания тепловых потоков в домах из дерева можно c контролируемой точностью использовать составные одномерные модели теплообмена и теплопередачи.
  2. Уравнения теплообмена содержат нелинейные граничные условия, которые можно приближенно заменить линейными средними соотношениями.
  3. Предложено установившиеся решения с периодическими граничными условиями находить методом наложения тепловых волн с использованием комплексных переменных. Наличие периодического нагрева помещения солнечной радиацией можно учесть с помощью разложения временной зависимости в ряд Фурье.
  4. Сформулированная математическая модель теплообмена помещения учитывает эффекты запаздывания распространения температурного поля в стенах при изменении температуры окружающего воздуха и динамическое изменение теплового потока от световых потоков снаружи и внутри помещения.
  5. Использование затенения световых проемов с солнечной стороны позволяет получить снижение затрат на одну треть для пиковых режимов кондиционирования.
Список литературы

1. Donn M.R., Thomas G., Joyce R., Hawkins T., Brown L. Designing comfortable homes: Guidelines on the use of glass, mass and insulation for energy efficiency. Cement & Concrete Association of New Zealand. Wellington, 2001. 78 p.

2. Fitton R., Swan W., Hughes T., Benjaber M. The thermal performance of window coverings in a whole house test facility with single-glazed sash windows // Energy Efficiency. 2017. Vol. 10. Pp.1419-1431. Doihttps://doi.org/10.1007/s12053-017-9529-0.

3. Борискина И.В. Здания и сооружения со светопрозрачными фасадами и кровлями. Теоретические основы проектирования светопрозрачных конструкций. С-Петербург: Инженерно-информационный Центр Оконных Систем, 2012. 400 с.

4. Dabaieh M., Elbably A. Ventilated Trombe wall as a passive solar heating and cooling retrofitting approach; a low-tech design for off-grid settlements in semi-arid climates // Solar Energy. 2015. Vol. 122. Pp. 820-833.

5. Canadian wood-frame house construction. Ottawa - Ontario: Canada Mortgage and Housing Corporation, 2013. 335 p.

6. Bisegna F., Mattoni B., Gori P., Asdrubali F., Guattari C., Evangelisti L., Sambuco S., Bianchi, F. Influence of Insulating Materials on Green Building Rating System Results // Energies. 2016. № 9. Pp. 712-728. Doihttps://doi.org/10.3390/en9090712.

7. Mitterpach J., Štefko J. An Environmental Impact of a Wooden and Brick House by the LCA Method // KEM. 2016. Vol. 688. Pp. 204-209.

8. Dokka T.H., Hemstad K. Energy efficient residential buildings for the future - a handbook for designing passive houses and low energy residential buildings, IEA SHC Task 28/ECBCS Annex 38 Sustainable Solar Housing (in Norwegian). 2006. 105 p.

9. Simonson C.J., Salonvaara M., Ojanen T. Improving indoor climate and comfort with wooden structures. Technical Research Center of Finland, ESPOO. Vuorimiehentie, 2001. 298 p.

10. Ma L., Shao N., Zhang, J., Zhao T. The Influence of Doors and Windows on the Indoor Temperature in Rural House // Procedia Engineering, 2015. Vol. 121. Pp. 621-627.

11. Mohammed A.A., Khudheyer A.F. Natural convection heat transfer inside enclosure with various geometries and enhancement methods a review // International Journal of Engineering & Technology, 2018. Vol. 7 (4). Pp. 4892-4898.

12. Kamiyo O.M., Angeli D., Barozzi G. S., Collins M. W., Olunloyo V., Talabi S. A Comprehensive Review of Natural Convection in Triangular Enclosures // Applied Mechanics Reviews. 2010. Vol. 63(6). 13р.

13. Rahimi A., Dehghan Saee A., Kasaeipoor A., Hasani Malekshah E. A comprehensive review on natural convection flow and heat transfer: The most practical geometries for engineering applications // International Journal of Numerical Methods for Heat & Fluid Flow, 2019. Vol. 29 (3). Pp. 834-877. Doihttps://doi.org/10.1108/HFF-06-2018-0272.

14. Awbi H.B., Hatton A. Natural convection from heated room surfaces // Energy and Buildings. 1999. Vol. 30(3). Pp. 233-244.

15. Ramesh N., Venkateshan S.P. Effect of Surface Radiation on Natural Convection in a Square Enclosure // Journal of Thermophysics and Heat Transfer. 1999. Vol. 13 (3). Pp. 299-301.

16. Dimassi N., Dehmani L. Natural Convection of a Solar Wall in a Test Room under Tunisian Climate // International Scholarly Research Network ISRN Renewable Energy. 2012. Vol. 2012. 9 р.

17. Zohuri B., Fathi N. Natural or Free Convection // Thermal-Hydraulic Analysis of Nuclear Reactors, 2015. Pp. 287-309. doi:https://doi.org/10.1007/s12209-019-00208-9.

18. Incropera F.P., DeWitt D.P. Fundamentals of Heat and Mass Transfer (4th edition). NY: John Wiley & Sons, 1996.

19. Богословский В.Н. Строительная теплофизика. М.: Высшая школа, 1982. 414 с.

20. Wood handbook: wood as an engineering material. General technical report FPL; GTR-113. Madison, WI: U.S. Department of Agriculture, Forest Service, Forest Products Laboratory, 1999. 463 p. doi:https://doi.org/10.2737/FPL-GTR-113.

21. Архив погоды на 21 июля в г. Воронеж [Электронный ресурс]. Систем. требования: Internet Explorer. URL: https://world-weather.ru/pogoda/russia/voronezh/14-july/ (дата обращения: 01.03.2022)

22. Jayamaha S.E.G., Wijeysundera N.E., Chou S.K. Measurement of the heat transfer coefficient for walls // Building and Environment. 1996. Vol. 31(5). Pp. 399-407. Doi:https://doi.org/10.1016/0360-1323(96)00014-5.

23. Мерщиев А.А., Шепс Р.А., Лобанов Д.В., Шашин А.В. Определение величины потока прямого солнечного излучения, направленного на горизонтальную поверхность // Региональная архитектура и строительство. 2020. №4. С. 137-143.

24. Басок Б.И., Новиков В.Г., Давыденко Б.В., Беляева Т.Г., Новицкая М.П., Сороковой Р.Я. Радиационно-конвективный теплообмен здания с окружающей средой при воздействии солнечного излучения // ИФЖ. Том 93. № 1. 2020. С. 48-57.

25. Lienhard J.H (IV)., Lienhard J.H. (V). A heat transfer textbook. Phlogiston Press, Cambridge Mass., 2006. 749 p.

26. Полежаев В.И., Верезуб Н.А., Простомолотов А.И. Математическое моделирование конвективного тепломассообмена на основе уравнений Навье-Стокса. М.: Наука. 1987. 272 с.

27. Скляров К.А., Мелькумов В.Н., Климентов А.В. Математическая модель вентиляционных процессов в помещениях сложной конфигурации // Вестник ВГТУ. Сер. Энергетика. 2006. Т.2 (6). С. 53-55.

28. Мелькумов В.Н., Кузнецов С.Н., Павлюков С.П., Кузнецов Р.Н. Формирование конвективных воздушных потоков при действии в помещении источника тепла // Вестник ВолгГАСУ. Сер.: Стр-во и архит. 2008. Вып. 12 (31). С. 76-80.

29. Кузнецов Г.В., Максимов В.И., Шеремет М.А. Естественная конвекция в замкнутом параллелепипеде при наличии локального источника энергии // Прикладная механика и техническая физика. 2013. Т. 54. №4. С. 86-95.

30. Мартюшев С.Г., Шеремет М.А. Математическое моделирование ламинарного режима сопряженного конвективного тепломассопереноса в замкнутой полости с источником энергии в условиях поверхностного излучения // Инженерно-физический журнал. 2013. Т.86. №1. С. 107-115.

31. Унаспеков Б.А, Сабденов К.О., Ерзада М., Ауелбеков С.Ш., Таубалдиева А.С., Сейтказиев О.Д. Исследование свободного конвективного движения воздуха в отапливаемом помещении // ИФЖ. 2020. Том 93. № 1. С. 170-177.

32. Zhou L., Liu J., Huang Q. Analysis of Combined Natural Convection and Radiation Heat Transfer in a Partitioned Rectangular Enclosure with Semitransparent Walls // Trans. Tianjin Univ. 2019. Vol. 25. Pp. 472-487.

33. Пикулин В.П., Похожаев С.И. Практический курс по уравнениям математической физики. М: МЦНМО, 2004. 208 с.

34. Тихонов А.Н., Самарский А.А. Уравнения математической физики. М.: МГУ, 1999. 798 с.

35. Минко В.А., Подпоринов Б.Ф., Семиненко А.С. Комплексное проектирование установок центрального водяного отопления зданий жилищно-гражданского назначения // Белгород: Изд-во БГТУ, 2013. 203 с.

36. Böckh P., Wetzel T. Heat Transfer. Basics and Practice. Springer, Karsruhe, 2012. 276 p.

37. Харкевич А.А. Спектры и анализ. М.: Изд-во ЛКИ, 2007. 240 с.

38. Gustavsen A., Jelle B.P., Arasteh D., Kohler C. State of the Art Highly Insulating Window Frames // Research and Market Review. United States: N.p., 2007.

39. Харкнесс Е., Мехта М. Регулирование солнечной радиации в зданиях. М.: Стройиздат, 1984. 176 с.

40. Методическое пособие. Расчеты теплопоступлений в здание от проникающей солнечной радиации за отопительный период. Москва, 2017. 111 с.

41. СП 345.1325800.2017 «Здания жилые и общественные. Правила проектирования тепловой защиты». М.: Mинстрой России, 2017. 51 c.

42. Погосян Х.П. Воздушная оболочка Земли. Ленинград: Гидрометеорологическое издательство, 1962. 294 с.

43. Mendes N., Gustavo H.C., de Araújo O. H. Building thermal performance analysis by using matlab/simulink // Seventh International IBPSA Conference Rio de Janeiro, Brazil. 2001. Pp. 473-480.

44. Самарин О.Д. О коэффициенте конвективного теплообмена на поверхностях ограждений помещений при периодическом тепловом режиме // Энергоснабжение и энергоподготовка. 2010. №1 (63). С. 48-49.


Войти или Создать
* Забыли пароль?