Сегодня, 8 октября
- 10:00-18:00 Санкт-ПетербургПетербургский международный газовый форум – 2024
- 10:00-18:00 Санкт-ПетербургXXVII Международная специализированная выставка газовой промышленности и технических средств для газового хозяйства «РОС-ГАЗ-ЭКСПО-2024»
- 10:00-18:00 Санкт-ПетербургМеждународная выставка по теплоэнергетике «КОТЛЫ И ГОРЕЛКИ-2024»
- 10:00-18:00 Санкт-ПетербургМеждународная выставка «Энергосбережение и энергоэффективность. Инновационные технологии и оборудование»
Ближайшие мероприятия
-
08 октября - 11 октябряСанкт-ПетербургПетербургский международный газовый форум – 2024
-
08 октября - 11 октября
-
08 октября - 11 октября
-
08 октября - 11 октября
-
09 октября - 10 октября
-
21 октября - 25 октября
Об эффективности воздушно-тепловых завес
Ю. Н. Марр, советник генерального директора АО «НПО «Тепломаш»
Нагретые воздушные струи завес, защищающих открытые двери и ворота, предотвращают понижение температуры внутреннего воздуха при контакте его с холодным наружным. Поскольку этот эффект обозначается как снижение теплопотерь через открытые проемы, то определение эффективности завес, на первый взгляд, не представляет трудности. Однако, как только вопрос переходит в руки специалистов, кажущаяся простота сменяется неоднозначностью и противоречивостью в предлагаемых оценках эффективности. Именно по этой причине при разработке стандарта ГОСТ 32512-2013 «Воздушные завесы. Общие технические условия» сообществом специалистов было принято решение не включать в стандарт определение эффективности завес и требования к ней.
Попытки оценить эффективность завес предпринимались еще в [1] по минимуму приведенных затрат при соотношении цен на материалы, трудозатрат и энергоносителей середины прошлого столетия в СССР. Некоторые аспекты эффективности завес (скорее, в плане постановки вопроса) были рассмотрены в [2]. В [3] использована энергетическая характеристика завесы как машиностроительного изделия, и аэродинамический КПД завесы как инженерного сооружения. Авторы [4–6] также коснулись этой темы. Однако ни в одной из упомянутых работ, за исключением [4], вопрос не нашел однозначного и окончательного решения. Лишь в [4] коэффициент тепловой эффективности наиболее точно соответствует существу темы данной работы.
1. Общая постановка вопроса об эффективности завес допускает рассмотрение данной темы в следующих аспектах:
— предмет оценки: завеса — как машиностроительное изделие (устройство для формирования нагретой струи с той или иной величиной потока импульса) или как инженерное сооружение (проектное решение по защите проема на заданные атмосферные и архитектурно-строительные условия, состоящее из комплекта завес — машиностроительных изделий);
— способ защиты проема: защита смесительного типа (без противодействия струями втекающему снаружи потоку) или защита шиберующего типа (активное струйное противодействие затеканию наружного воздуха);
— способ тепловой компенсации теплопотерь: при шиберующей защите воздух струи нагрет непосредственно в завесе или организована раздельная аэродинамическая защита холодными струями и тепловая компенсация подачей нагретых отдельно масс воздуха;
— подход к оценке эффективности: для расчетных зимних условий по расчетной тепловой мощности или для всего отопительного периода с переменной наружной температурой по суммарным затратам энергии;
— организация защиты: предусмотрено регулирование работы завес в зависимости от наружной температуры, в частности, изменением частоты вращения вентиляторов (изменением расхода воздуха), или регулирование отсутствует.
Здесь будут рассмотрены завесы как инженерные сооружения для защиты шиберующего и смесительного типов.
2. Трактовка эффективности завес как инженерного сооружения имеет смысл относительного уменьшения энергетических затрат (включая остаточные теплопотери и дополнительные компенсационные затраты) при защите проема завесой в сравнении с незащищенным проемом [2]. Она принята в справочнике ASHRAE [4] и имеет вид
Э = [Q1 – (Q2 + Qз)]/Q1 = 1 – (Q2 + Qз)/Q1, (1)
где Q1 — тепловая мощность мгновенных теплопотерь через открытый проем при неработающей завесе;
Q2 — тепловая мощность мгновенных теплопотерь через открытый проем при защите завесой;Qз — тепловая мощность завесы.
Область применения выражения (1) ограничена защитой шиберующего типа, поскольку при защите смесительного типа Q1 = Q2 и величина Э принципиально не может быть положительной.
Приведем без выкладок составляющие выражения (1), опираясь на общие балансовые соотношения для завес любой установки (верхних и боковых) [2]. Определим величину Q1 как тепловую мощность, необходимую для синхронного подогрева затекающего в проем наружного воздуха от расчетной зимней температуры t1 до внутренней температуры t2. Расход втекающего воздуха определяется по разности давлений в проеме ΔР. Величину Q2 определим как тепловую мощность синхронного подогрева втекающей в проем смеси (при защите завесой) от температуры смеси tсм до температуры t2 . При этом второе слагаемое в (1) примет вид
‹Q› = (Q2 + Qз)/Q1 = β{(1 – θсм)[1/q + 0,5(λ – 1)] + (1/θ2 – 1)}. (2)
Здесь параметр β = Gз/G1 = (1/μF‾)(Т1/Т2)(2 ΔР/ρ2vз2)-0,5;
vз— скорость струи воздуха на выходе из завесы;
ΔР — разность давлений в проеме;
μ — коэффициент расхода через незащищенный проем;
F‾ = Fпр/ Fз, θсм = (tсм – t1)/( t2 – t1) — безразмерная температура смеси ;
λ — коэффициент эжекции струи завесы;
θ2 = (t2 – t1)/( tз – t1) — температурный параметр;
tз — температура воздуха на выходе из завесы.
Температура смеси из теплового баланса равна
θсм = [1/θ2 + 0,5(λ – 1) – Q̅пот/ θ2][1/q + 0,5(λ – 1)] -1, (3)
где Q̅пот = Qпот / СрGз(tз – t1) — безразмерные потери тепла с уходящими наружу частями эжектированного снаружи воздуха.
Выражение (3) записано с учетом составляющей, привносящей в баланс теплоту внутренней эжекции с температурой t2 . Так требуется для верхних завес [2]. Для боковых завес принимается, что эжектируемые изнутри массы имеют температуру tсм [1]. Поэтому в выражении (3) из числителя и знаменателя удаляются слагаемые с коэффициентом эжекции, после чего остается
θсм = [1 – Q̅пот]q/ θ2. (3-1)
Подставляя (3) в (2), получим составляющую (1) в виде
‹Q› = β[1/q – 1 + Q̅пот/ θ2]. (4)
Выражение (4) совместно с (3) позволяет анализировать эффективность способов тепловой компенсации теплопотерь. Рассмотрим случай, когда вся тепловая мощность для полной компенсации теплопотерь подводится через завесу: tсм = t2 , θсм = 1, откуда по (3) имеем
θ2 = (1 – Q̅пот) q (5)
и далее по (4)
‹Q›(зав) = β[1/ (1 – Q̅пот) q – 1]. (6)
Для численной оценки по (6) примем q = 1, Q̅пот = 0,5. Это дает ‹Q›(зав) = β. Величина β лежит в диапазоне β = 0,12 – 0,2. Отсюда получается, что при полной защите проема завесой с нагревом наибольшая эффективность Э = 0,80 – 0,88.
Другой случай — раздельная аэродинамическая и тепловая защита — подразумевает защиту холодной завесой (θ2 = 1) совместно с тепловой компенсацией вне зоны действия завесы. При этом из (4) имеем
‹Q›(разд) = β[1/q – 1 + Q̅ пот]. (7)
Разделив (6) на (7) и полагая равенство показателей защиты, после преобразования получим
‹Q›(зав) / ‹Q›(разд) = {1 – Q̅пот2[1 – (1 – Q̅пот) q]-1}-1. (8)
В режиме полной защиты (q = 1) из (8) имеем
‹Q›(зав) / ‹Q›(разд) = 1 / (1 – Q̅пот), (9)
откуда следует, в соответствии с [2], что при Q̅пот ~ 0,5 раздельная защита может достигать двукратного преимущества перед полной защитой через завесу. В предельном режиме завес из (8) имеем (q = q*, Q̅пот = 0)
‹Q›(зав) / ‹Q›(разд) = 1 (10)
и различие между способами защиты исчезает.
Выражение (1) допускает не только оценку текущей тепловой мощности при фиксированной температуре наружного воздуха (например, расчетной зимней), но и интегральную энергетическую оценку за весь отопительный период. В этом случае под составляющими (1) следует понимать не тепловые мощности, а затраты теплоты (энергии) соответствующего вида за отопительный период. Текущая и интегральная за отопительный период оценки могут существенно различаться. С повышением наружной температуры режим работы завесы уходит от расчетного и, если ничего не менять, струи завесы начнут вылетать из проема, унося на улицу тепловую мощность завесы. Разность (Q1 – Q2) начнет уменьшаться, и вместе с нею эффективность Э. Как минимум для сохранения эффективности в [2, 8] разработаны рекомендации по изменению частоты вращения вентиляторов (расхода воздуха через завесу) при повышении наружной температуры.
3. Оценка эффективности по выражению (1) имеет специфический характер: сопоставляется ситуация, когда открытый проем защищен завесой, с ситуацией, когда проем не защищен. В этом проявляется, до некоторой степени, ее абсурдность, как если бы эффективность отопления помещения оценивалась по ситуации полного отсутствия его теплозащиты. Иными словами, сравнивались бы объекты, имеющие стены и потолок и не имеющие их. Информация, которую несет выражение (1), принципиально отличается от принятой в системе энергетических оценок технических объектов и опирающейся на привычное всем представление о коэффициенте полезного действия (КПД). Оценка типа (1) — результат применения объекта в некоей системе, при этом качество самого объекта вторично. Смысл КПД состоит в выведении на первый план собственно объекта путем сопоставления идеализированной (теоретической) его специфики с реальностью. Оставаясь в энергетических представлениях, КПД завесы должен быть отношением минимально возможной тепловой мощности на защиту проема завесой к мощности, реально затрачиваемой. В такой постановке и аэродинамический КПД по [3] не является коэффициентом полезного действия, а, скорее, носит характер оценки типа (1).
Как следует из предыдущего, минимальные тепловые затраты на защиту проема обеспечивает раздельная аэродинамическая и тепловая защита. Тепловой КПД защиты можно было бы определить, в этом случае, как
ηq = Qразд / Qз, (11)
где под Qз понимается тепловая мощность оцениваемой завесы с проектной величиной температуры смеси tсм/з, а под Qразд — тепловая мощность раздельной защиты, обеспечивающей конечную температуру потока воздуха в ареале проема, равную tсм/з. Все исходные условия защиты проема для числителя и знаменателя (11) следует принимать одинаковыми. Иными словами, числитель описывает аэродинамическую защиту проема холодной завесой на тех же условиях, что и знаменатель.
Парадокс состоит в том, что в предельном режиме q = q*, т. е. при слабой аэродинамической защите, согласно (10), ηq = 1. В то же время при полной защите q = 1 КПД имеет наименьшую величину. Разрешение парадокса простое. Потери тепла, связанные с затеканием в проем наружного воздуха при q < 1, одинаковы в числителе и знаменателе (11), они максимальны при q = q * и равны нулю при q = 1. Потери с уходящими массами, напротив, равны нулю при q = q* и максимальны при q = 1, причем эти потери больше в знаменателе. В связи с этим выражение (11) не может служить оценкой эффективности.
4. На практике качество работы завес, чаще всего, оценивается по тому, как подогрет затекающий в проем воздух. Недаром в [6] прямо сказано, что «критерием эффективности завес любого типа является температура смеси», а в [4] в качестве эффективности защиты охлаждаемых помещений введена безразмерная температура смеси θ в виде, приспособленном для этого. Исходя из теплового баланса, температура смеси (3) есть сумма всех тепловых поступлений и потерь. Безразмерная температура смеси по определению
θсм = (tсм – t1)/( t2 – t1) (12)
говорит о том, что это есть отношение мер реальной тепловой мощности к такой идеальной, которая может быть интерпретирована, как целевая установка приближения к внутренней температуре при защите проема. Таким образом, выражение (12) следует считать термическим КПД любого типа шиберующей защиты
ηt = θсм. (13)
При tсм → t1 ηt = 0, при tсм → t2 ηt = 1. В своде правил СП 60.13330.2012 «Отопление, вентиляция и кондиционирование воздуха» указаны нормированные минимальные температуры смеси. Их величины не зависят от наружной температуры. Можно считать, что этим задаются нормированные минимальные термические КПД на любые наружные условия. Максимальные значения КПД не ограничены. Однако следует принять во внимание, что стремление к повышению КПД во всех видах техники всегда сопровождается удорожанием объекта. Поэтому необходим разумный баланс между КПД и затратами на объект.
Перепишем для удобства выражение (3) в виде
ηt = θсм = [(1– Q̅пот)/θ2 + 0,5(λ – 1)][1/q + 0,5(λ – 1)] -1. (14)
В (14) для верхней завесы имеем по [2]
Q̅пот = η{0,25(2 – θ2)[1 – exp( -2εL)] + 0,5 θ2 εL}. (15)
εL = 0,065(L/bз)0,556 L = H/cosα
η = 2[0,5(λ + 1) – 1/q]/(λ – 1). (16)
Для боковых завес используется выражение (3-1) с графиками Q̅пот от F̅ и от q по [9].
Если зафиксировать размеры проема, ширину сопел завесы, разность давлений в проеме, температуру струи завесы, то вариация расхода завесы приведет к изменению коэффициента эжекции λ и показателя зашиты q. С наращиванием расхода будет увеличиваться показатель защиты (уменьшение проникновения наружного воздуха) и относительные потери тепла. В целом термический КПД будет расти вплоть до единицы. Холодные завесы (θ2 = 1) имеют наименьший КПД.
Проиллюстрируем численно оценку по выражению (14) на примере стандартного проема 4х4 м с температурами воздуха –26 0С, +18 0С, скоростью ветра 4 м/с, водой 95/70 0С. Этим условиям соответствует разность давлений в проеме 6,06 Па. Верхняя завеса состоит из двух модулей КЭВ-230П7020W. Угол струи к плоскости проема 300. Ширина сопла 150 мм, скорость струи в сопле 12,5 м/с. Температура струи на выходе 43 0С. Неизменные параметры:
ΔР = 6,06 Па, F̅̅ = 26,7 , Q̃пот = 0,338. Температуру струи при варьировании будем удерживать постоянной изменением выходных параметров воды, поэтому полагаем θ2 = 0,638 = const. Результат, соответствующий принятым условиям, приведен во втором столбце табл. 1. В остальных столбцах показано влияние вариации расхода воздуха и скорости струи в этой завесе. В диапазоне режимов работы от предельного (q* = 0,49) до полной защиты (q = 1) термический КПД возрастает от 0,846 до 1 одновременно с ростом потока импульса. Этим подтверждается непротиворечивый характер введенного термического КПД по (14). Не вызывает удивления, что максимальная эффективность при заданных фиксированных условиях достигается в режиме полной защиты на максимальном расходе воздуха и при максимальной тепловой мощности. Равенство КПД единице в первых двух столбцах означает защиту высокой температурной эффективности tсм = t2. Если, к примеру, разность давлений в проеме увеличится до 10,0 Па, то режим вплотную приблизится к предельному и КПД станет равным ηt = 0,85 (tсм = 11,4 0С). В отсутствие нагрева КПД опустится до ηt = 0,66, что в условиях нашего примера соответствует tсм = 3,0 0С.
Таблица 1. Термический КПД при неизменной ширине сопла bз = const
Параметр | Величина параметра | ||||
vз, м/с | 13,0 | 12,5 | 11,6 | 10,6 | 9,5 |
Vз, м3/час | 28 080 | 27 000 | 25 060 | 22 900 | 20 500 |
Iз , Н (один модуль) | 60,8 | 56,3 | 48,4 | 40,5 | 32,4 |
q | 1,0 | 0,91 | 0,77 | 0,62 | 0,49* |
η | 1,0 | 0,91 | 0,73 | 0,43 | 0 |
Q̅пот = η Q̃пот | 0,339 | 0,31 | 0,25 | 0,146 | 0 |
ηt = θсм | 1,0 | 0,994 | 0,949 | 0,898 | 0,846 |
Qз, кВт | 235 | 226 | 210 | 192 | 172 |
Оценим влияние уменьшения ширины сопла с ростом скорости при неизменном потоке импульса выбранного единичного модуля Iз = 1,2 х 0,15 х 2,0 х 12,52 = 56,3 Н = const. Температура струи, по-прежнему, остается неизменной θ2 = 0,638 = const. Расчет для верхней завесы по [2] приведен в табл. 2.
Таблица 2. Термический КПД при условии I = 56,3 Н = const
Параметр | Величина параметра | ||||
bз, м | 0,15 | 0,12 | 0,10 | 0,08 | 0,06 |
vз, м/с | 12,5 | 14,0 | 15,3 | 17,1 | 19,8 |
Vз, м3/час | 27 000 | 24 190 | 22 030 | 19 700 | 17 100 |
F̅ | 26,7 | 33,3 | 40,0 | 50 | 66,7 |
q | 0,91 | 0,86 | 0,81 | 0,76 | 0,70 |
η | 0,91 | 0,88 | 0,84 | 0,82 | 0,79 |
εL | 0,44 | 0,49 | 0,55 | 0,62 | 0,73 |
Q̃пот | 0,34 | 0,37 | 0,40 | 0,44 | 0,49 |
Q̅пот = η Q̃пот | 0,31 | 0,32 | 0,34 | 0,36 | 0,39 |
ηt = θсм | 0,994 | 0,96 | 0,927 | 0,90 | 0,865 |
Из табл. 2 видно, что увеличение скорости приводит к очевидному падению показателя защиты и вместе с ним термического КПД. Это объясняется нарастанием масштабного фактора F‾ (гидравлической длины струи) при неизменном потоке импульса.
Во всех оценках табл. 1 и 2 коэффициент качества струи был принят равным ξ = 0,8. Напомним, что коэффициент качества характеризует внешнее влияние на начальное состояние струи, приводящее в дальнейшем к ее размыванию, которое выражается в степени снижения максимальной скорости на оси струи по отношению к скорости стандартной невозмущенной турбулентной струи [2] ξ = vm/‹vm›. Внешнее влияние на струи завес может выражаться в сильной турбулизации потока тангенциальными вентиляторами, внесении неравномерности скорости по сечению сопла препятствиями и поворотами, разделением струи на части при истечении из разнесенных сопел. Обычно в завесах ξ = 0,8–0,9. Однако в случаях плохо организованных сопел коэффициент качества опускается до 0,6. На примере модуля КЭВ-230П7020W для принятых выше условий покажем, как поведет себя термический КПД c понижением качества струи (см. табл. 3).
Таблица 3. Влияние коэффициента качества на термический КПД
Параметр | Величина параметра | |||
ξ | 0,9 | 0,8 | 0,7 | 0,6 |
λ | 2,98 | 3,26 | 3,58 | 3,95 |
q | 0,96 | 0,91 | 0,87 | 0,82 |
η | 0,96 | 0,91 | 0,88 | 0,85 |
Q̅пот = η Q̃пот | 0,33 | 0,31 | 0,30 | 0,29 |
ηt = θсм | 1,0 | 0,99 | 0,98 | 0,96 |
Из табл. 3 видно, что коэффициент качества значительно сильнее влияет на аэродинамику, чем на тепловые характеристики — термический КПД лишь слабо убывает. Это объясняется тем, что при почти неизменяющихся потерях Q̅пот увеличение поступления холодного наружного воздуха с понижением q компенсируется приростом внутренней эжекции (рост λ).
Таким образом, термический КПД является полноценным инструментом для оценки действительной эффективности завесы как инженерного сооружения и сравнения между собой различных вариантов.
5. Наиболее непростой вопрос — это оценка эффективности завес смесительного типа. Выражение (1) для этих целей непригодно по существу и дает отрицательное значение Э. Выражение (12) также непригодно, поскольку, в отличие от шиберующих завес, в смесительных происходит лишь перемешивание втекающего холодного и подаваемого нагретого воздуха и температура смеси определяется из обычного теплового баланса. Как было сформулировано еще в [7], энергетическая эффективность защиты смесительного типа по отношению к альтернативным способам компенсации теплопотерь через открытые двери проявляется как результат интенсивного струйного перемешивания воздушных масс для минимизации тепловой мощности на подготовку воздуха, подаваемого в помещение. К сожалению, этому очевидному качественному соображению невозможно придать форму количественной оценки. Не отвергая высказанного, представим защиту смесительного типа как часть раздельной аэродинамической и тепловой защиты, в которой аэродинамическая часть или полностью отсутствует, или заменена комплексом архитектурно-строительных препятствий втеканию наружного потока в помещение. В этом случае становится возможным ввести понятие теплового КПД защиты смесительного типа, в числителе которого будет условно минимизированная тепловая мощность при наличии идеализированной аэродинамической защиты, а в знаменателе реальная тепловая мощность оцениваемой ситуации
ηсмес = Qид / Qсмес. (17)
В качестве идеализированного случая примем защиту герметичного помещения. Согласно [7], его можно интерпретировать, как ситуацию с наддувом условного тамбура. Это минимизирует втекание наружного воздуха и вытекание внутреннего, соответственно, теплопотери (Qид) также станут минимальными. Расход наружного воздуха для этого случая определим как
Gнар/герм = μнизК2Fпр/герм(2ρ1ΔPпр/герм)0,5, (18)
где ΔPпр/герм = 9,8(ρ1 – ρ2)(Нпр/4) — разность давлений в нижней части проема герметичного помещения (в данном случае условно герметизированного наддувом тамбура);
μниз — коэффициент расхода втекающего через нижнюю половину дверей наружного воздуха; К2 — коэффициент, учитывающий количество людей, проходящих через дверь в течение часа по [9];
Fпр/герм = В(Нпр/2).
При этом идеализированные теплопотери равны
Qид = Cp Gнар/герм (tсм – t1), (19)
где tсм — заданная температура смеси втекающего наружного воздуха и нагретого, подаваемого завесой. Знаменатель (17) определим через расчетную разность давлений в проеме
Qсмес = Cp Gнар (tсм – t1)
Gнар = μпрК2Fпр(2ρ1ΔPпр)0,5, (20)
где Fпр= ВНпр, а ΔPпр определяется по [9] через расчетную высоту hрасч, зависящую от этажности здания. Подставляя (18) – (20) в (17) получим после упрощений
ηсмес = (μниз /μпр)( Fпр/герм/ Fпр)[(Нпр/4) /hрасч]0,5.(21)
Для иллюстрации, принимая Нпр = 2 м, высоту этажа 3 м и вычисляя hрасч по [9], а также принимая ориентировочно μниз = 0,6, μпр = 0,7, найдем по (21) величины КПД для продуваемых зданий:
ηсмес (1 этаж) = 0,21, ηсмес (3 этажа) = 0,11, ηсмес (5 этажей) = 0,094.
Понятно, что с ростом этажности здания тепловой КПД убывает. При защите помещения герметичного типа hрасч = Нпр/4 и ηсмес = 1. В случаях «мягкого» шиберования втекающего потока, описанных в [10], коэффициент расхода μпр снижается против принятого табличного значения, что повысит КПД такой защиты.
Выводы
- Общепринятая эффективность завес по выражению (1) основана на сопоставлении полного отсутствия защиты проема с эффектом действия завесы. При таком подходе невозможно оценить инженерное сооружение в представлении о коэффициенте полезного действия.
- Оценка по (1) непригодна для защиты смесительного типа.
- Предложенная в [4] тепловая эффективность является полноценным термическим коэффициентом полезного действия защиты шиберующего типа (14).
- Для защиты смесительного типа предложен тепловой КПД по (21).
Литература
- Эльтерман В. М. Воздушные завесы. Изд. 2-е. М.: Машиностроение. 1966.
- Марр Ю. Н. Воздушно-тепловые завесы. Расчет и проектирование завес для защиты проемов промышленных и общественных зданий. — СПб.: АО «НПО «Тепломаш», 2017, — 160 с.
- Караджи В. Г., Московко Ю. Г. Вентиляционное оборудование. Технические рекомендации для проектировщиков и монтажников. — М.: АВОК-ПРЕСС, 2010 — 432 с.
- Стронгин А. С. Расчетные параметры и эффективность применения воздушных завес у ворот охлаждаемых помещений // Инженерные системы. АВОК-Северо-Запад. № 2. 2019.
- Гримитлин А. М., Стронгин А. С. Воздушные завесы для зданий и технологических установок. Учебное пособие. 2018. — 136 с.
- Гримитлин А. М., Дацюк Т. А., Крупкин Д. Я., Стронгин А. С., Шилькрот Е. О. Отопление и вентиляция производственных помещений. Издательство «АВОК Северо-Запад». СПб. 2007. — 400 с.
- Марр Ю.Н. О завесах смесительного типа//Инженерные системы. АВОК-Северо-Запад. № 1. 2012.
- Марр Ю. Н. Проблемы регулирования шиберующих завес// Инженерные системы. АВОК-Северо-Запад. № 2. 2018.
- Справочник проектировщика. Внутренние санитарно-технические устройства. Часть 3. Вентиляция и кондиционирование воздуха. Книга 1. / В. Н. Богословский, А. И. Перумов, В. Н. Посохин и др. 4-е издание. М.: Стройиздат. 1992.
- Марр Ю. Н. Экономичная и комфортная защита входных дверей общественных зданий завесами с «мягким» шиберованием // Инженерные системы. АВОК-Северо-Запад. № 2. 2016.
Скачать PDF версию статьи «Об эффективности воздушно-тепловых завес»