CC BY
55
УДК 693.5:536.24+532.51
ГНЫРЯ АЛЕКСЕЙ ИГНАТЬЕВИЧ, докт. техн. наук, профессор, gnyria@mail ru
КОРОБКОВ СЕРГЕЙ ВИКТОРОВИЧ, канд. техн. наук, докторант, доцент, korobkov@hotmail ru
КОШИН АНТОН АЛЕКСАНДРОВИЧ, аспирант, dawghood@mail. ru
МОКШИН ДМИТРИЙ ИЛЬИЧ, магистрант, mokshin@sibmail com
Томский государственный архитектурно-строительный университет, 634034, г. Томск, пл. Соляная, 2
ТЕРЕХОВ ВИКТОР ИВАНОВИЧ, докт. техн. наук, профессор,
terekhov@itp. nsc. ru
Институт теплофизики СО РАН,
630090, г. Новосибирск, пр. Академика Лаврентьева, д. 1
КОМПЛЕКСНЫЕ ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ АЭРОДИНАМИКИ И ТЕПЛООБМЕНА МОДЕЛЕЙ ЗДАНИЙ И СООРУЖЕНИЙ*
Статья посвящена комплексным экспериментальным исследованиям аэродинамики и теплообмена моделей высотных зданий. Рассмотрены: метод сажемасляной визуализации воздушных течений вдоль ряда из двух квадратных призм при увеличении калибра между ними, расположенных друг за другом, метод измерения полей статического давления с использованием многоканального дифференциального манометра, метод получения распределения коэффициента теплоотдачи по высоте модели здания.
Ключевые слова: аэродинамика, теплоотдача, поля статического давления, модели высотных зданий.
GNYRYA, ALEKSEY IGNATJEVICH, Dr. of tech. sc., prof., gnyria@mail. ru
KOROBKOV, SERGEY VIKTOROVICH, Cand. of tech. sc., assoc. prof., korobkov@hotmail. ru
KOSHIN, ANTON ALEKSANDROVICH, P.G. dawghood@mail. ru
MOKSHIN, DMITRIY ILJICH, graduate students, mokshin@sibmail com
Tomsk State University of Architecture and Building,
2 Solyanaya sq., Tomsk; 634003, Russia TEREKHOV, VICTOR IVANOVICH, Dr. of tech. sc., prof., terekhov@itp. nsc. ru Institute of Thermo-physics,
1 Academika Lavrentyev avenue, Novosibirsk, 630090, Russia
* Исследования проводятся при финансовой поддержке Российского фонда фундаментальных исследований (проект № 09-08-00523-а).
© А.И. Гныря, С.В. Коробков, А. А. Кошин, Д.И. Мокшин, В.И. Терехов, 2011
COMPREHENSIVE EXPERIMENTAL STUDY OF AERODYNAMICS AND HEAT TRANSFER OF MODELS OF BUILDINGS AND STRUCTURES
The article presents a brief report on the work devoted to a comprehensive experimental study of aerodynamics and heat transfer of models of high-rise buildings. The following methods are considered: the method of oil-black visualization of air flows along a row of two square prisms at caliber increase between them, arranged one after other; the method of measuring the static pressure fields using a multichannel differential manometer; the method of obtaining the distribution of the heat transfer coefficient on the height of the building model.
Keywords: aerodynamics, heat transfer, static pressure fields, models of high-rise buildings.
Введение
Аэродинамика зданий и сооружений детально исследовалась в большом числе работ, но в силу сложности структуры трехмерных отрывных потоков в данной проблеме остается много слабо изученных вопросов. В настоящее время в научной литературе имеется множество работ, посвященных изучению структуры турбулентных отрывных потоков, полей давления и коэффициентов теплоотдачи моделей высотных зданий. Эта проблема актуальна в связи с интенсивным строительством высотных сооружений, для которых ветровые нагрузки являются одним из определяющих факторов в процессе их эксплуатации. В то же время экспериментальные методы исследований трехмерных потоков достаточно сложны и дорогостоящи. Численные методы расчета турбулентных пространственных течений, которые активно развиваются в последнее десятилетие, также имеют ряд принципиальных недостатков. В частности, к таковым можно отнести отсутствие единой универсальной модели турбулентности, необходимость тщательного их тестирования и др. Поэтому большое значение приобретает разработка простых методов визуализации, дающих детальное и наглядное представление о структуре трехмерных потоков, которое одновременно с данными о полях давлений позволяют определить зоны повышенного и пониженного теплопереноса. Об этом свидетельствуют опытные данные авторов по изучению аэротермодинамики обтекания куба. Такие методы особенно важны при практическом решении вопросов оптимального расположения зданий городской застройки [1-9].
Для проведения комплексного исследования структуры течения воздушного потока были изготовлены экспериментальные модели квадратных призм для сажемасляной визуализации; измерения распределения давления; исследования теплообмена.
Сажемасляная визуализация
В настоящее время выполнен цикл исследований, посвященный визуа-лизационному изучению турбулентного обтекания квадратных призм (моделей дома) высотой H/h = 3-6 калибров. Одиночные призмы, а также их тандем располагались на плоскости при углах атаки 0 и 45°. Кроме того, были измерены коэффициенты давления по всем граням призмы, а также на плоскости
до призмы и в ее кормовой зоне. Для исследования аэродинамики при помощи сажемасляной визуализации были изготовлены экспериментальные модели квадратных призм, имитирующих высотные здания. Опыты проводились в аэродинамической трубе, установленной на кафедре технологии строительного производства ТГАСУ (рис. 1).
Рис. 1. Общий вид аэродинамической трубы испытательного стенда
Для визуализационных исследований структуры течения вдоль ряда из двух квадратных призм, расположенных друг за другом, было изготовлено шестнадцать испытуемых моделей в виде квадратных призм сечением 50x50 мм и высотой 150 и 300 мм по восемь на каждый из углов ф. Исследовалась картина обтекания потоком воздуха при углах атаки воздушного потока ф = 0° и ф = 45° и при калибрах между ними L/H = 0; 0,5; 1 и 2 (рис. 2 и 3). Опыты были проведены при максимальной скорости воздушного потока, соответствующей числу Рейнольдса, рассчитанного по размеру грани Re = 4,66-104. Все модели изготавливались из оргстекла толщиной 5 мм и устанавливались на подложку, выполненную из этого же материала. Все грани квадратных призм, в том числе и верхняя поверхность, а также подложка, на которой устанавливались модели, покрывались тонким слоем раствора черного цвета, полученного путем смешивания типографской краски и керосина. Исследуемые модели помещались в рабочую камеру аэродинамической трубы, где устанавливался необходимый гидродинамический режим воздушного потока. Эксперимент продолжался до тех пор, пока не появлялась четко выраженная картина обтекания модели потоком воздуха. Схема установки исследуемых моделей в аэродинамической трубе представлена на рис. 4.
Рис. 2. Общий вид ряда из двух моделей квадратных призм высотой Н = 150 мм, расположенных друг за другом для проведения сажемасляной визуализации: а - ф = 0°, L/H = 0; б - ф = 45°, L/H = 0
Рис. 3. Общий вид ряда из двух моделей квадратных призм высотой Н = 300 мм, расположенных друг за другом для проведения сажемасляной визуализации: а - ф = 0°, L/H = 0; б - ф = 45°, L/H = 0
Рис. 4. Схема расположения испытуемых моделей в аэродинамической трубе во время эксперимента
В результате удалось получить наглядную картину движения воздушных потоков вблизи призм при различных калибрах и углах атаки. Некоторые результаты приведены на рис. 5-8.
На рис. 5-8 видно, что пограничный слой, образующийся на пластине (подложке), приближаясь к впереди стоящей на плоскости квадратной призме выше по течению, подвергается трехмерному разделению. Область течения в непосредственной близости от передней призмы характеризуется наличием подковообразного вихря, возникающего на его передней грани и простирающегося вдоль боковых граней ряда призм в направлении ниже по течению. Течение вблизи подковообразного вихря имеет неустойчивый характер.
На боковых гранях наблюдаются зоны рециркуляции в виде вихрей. Эти вихри покрывают значительную часть боковых граней призмы. Такое явление, наблюдаемое во время сажемасляной визуализации воздушного потока, заключается в неустойчивом с периодическими пульсациями течении, возникающем ниже на боковых гранях призмы, где внутренний подковообразный вихрь испытывает разворот.
Рис. 5. Картина обтекания ряда квадратных призм потоком воздуха. Высота призм 300 мм, калибр 0,5, угол 0°
Рис. 6. Картина обтекания ряда квадратных призм потоком воздуха. Высота призм 300 мм, калибр 1, угол 45°
Рис. 7. Картина обтекания ряда квадратных призм потоком воздуха. Высота призм 150 мм, калибр 1, угол 0°
Рис. 8. Картина обтекания ряда квадратных призм потоком воздуха. Высота призм 150 мм, калибр 2, угол 45°
За второй призмой формируется сводообразный вихрь. Его отпечатки на подложке в результате визуализации обнаружены в виде двух противоположно вращающихся вихрей. Верхняя грань второй призмы находится в зоне отрыва течения, которая начинается около переднего края верхней грани на впереди стоящей квадратной призме.
Очевидно, что описанная выше качественная структура сложных отрывных течений будет непосредственно сказываться на характере изменения локальной теплоотдачи. Полученные данные моделируют теплообменные процессы зданий и сооружений высокой этажности в зависимости от температуры окружающей среды, скорости и угла атаки воздушного потока, а также местоположения здания. Подробный отчет представлен в [10, 12].
Измерение распределения давления
Для проведения экспериментов по определению полей статического давления было изготовлено две модели. Первая выполнена в виде квадратной призмы сечением 50^50 мм и высотой 300 мм. Две боковые и задняя грани модели были изготовлены из органического стекла толщиной 5 мм, передняя грань - из текстолита толщиной 8 мм. В передней грани модели для отбора статического давления были просверлены отверстия диаметром 0,8 мм по нормали к поверхности в вертикальном и горизонтальном направлениях с шагом 10 мм. Общее количество приёмников давления по высоте составило 38 шт. и по ширине - 10 шт. Для измерения полей статического давления в пограничной области до и после призмы, на которой установлена модель, были выполнены отверстия диаметром 0,8 мм с шагом 10 мм: до модели в количестве 9 шт., за моделью - 19 шт. Схема их размещения показана на рис. 9.
Вторая модель представляет собой прямоугольную призму сечением 50x50 мм, высотой 150 мм. Суммарное количество приемников давления также составляет 38 шт. (рис. 10).
Измерения перепада давления производились дифференциальным многоканальным микроманометром с ценой деления 1 Па. В качестве опорного берется статическое давление в канале для невозмущенного течения. Общий вид манометра представлен на рис. 11.
Скорость воздушного потока перед моделью измерялась с помощью трубки Пито - Прандтля, расположенной в центре канала выше по потоку от модели. Показания трубки Пито регистрируются наклонным микроманометром. Показания многоканального микроманометра снимались с помощью цифрового фотоаппарата Nikon, далее показания с фотографии оцифровывались по специальной программе Get Data Graph Digitizer.
Были проведены эксперименты на моделях высотой 150 и 300 мм при углах атаки набегающего потока воздуха 0 и 45°.
В результате проведенных экспериментов были рассчитаны коэффициенты давления Ср в каждой точке и построены графики его изменения по образующей модели (рис. 12). Подробный отчет представлен в [3].
Рис. 9. Схема модели квадратной призмы высотой 300 мм:
1 - подложка; 2 - квадратная призма; 3 - отверстия-приемники диаметром 0,8 мм; 4 - трубки ПВХ диаметром 3 мм
на многоканальный микроманометр
Рис. 10. Схема модели квадратной призмы высотой 150 мм:
1 - подложка; 2 - квадратная призма; 3 - отверстия-приемники диаметром 0,8 мм; 4 - трубки ПВХ диаметром 3 мм
Рис. 11. Общий вид дифференциального многоканального микроманометра
-ш-
X/H
X/H
(X/L)
Рис. 12. Распределение полей давления по граням квадратной призмы высотой 300 мм, а также в областях до и после призмы при угле атаки воздушного потока ф = 0° (Re = 4,66-104):
□ - отметки по продольной оси; ◊ - отметки по боковой грани. Данные I.P. Castro и A.G. Robins для куба 60x60x60 мм: x - куб, установленный на дне канала; Д - куб, установленный на подложке; Д' - для боковой грани
Отметим характерные особенности поведения коэффициента давления Ср. Как видно, на плоской поверхности канала перед призмой давление возрастает от значения в потоке до его величины на лобовой стенке призмы. При этом протяженность области торможения составляет примерно 2-3 калибра по размеру грани призмы. Размер же зоны восстановления потока за кормовой гранью призмы существенно больше, и он оценивается величиной в 5-6 калибров. Как и следовало ожидать, в кормовой зоне за счет мощного отрыва
потока имеет место разрежение, величина Ср превышает скоростной напор. В итоге, разность давлений на наветренной и подветренной сторонах призмы равна примерно двум скоростным напорам. Этот вывод является важным с точки зрения определения аэродинамических сил, действующих на строительную конструкцию.
На лобовой грани коэффициент давления не превышает Ср < 0,8, что говорит об отсутствии полного торможения набегающего потока. Это может быть вызвано растеканием заторможенных струек тока от центра к краям призмы. Действительно, если посмотреть на распределение коэффициента давления Ср поперек призмы, то отчетливо видно, что давление посредине призмы выше, чем на краях. Разность давлений вызывает растекание потока по малому размеру призмы. Это находится в соответствии с данными Castro (1979), что говорит о значительном влиянии растекания газа по поверхности призмы на характер распределения давления.
Для сопоставления и анализа данных результатов сравним картину, полученную при проведении сажемасляной визуализации, и данные, полученные при изучении полей статического давления на поверхностях перед призмой и после призмы. Перед призмой образуется область повышенного давления, что имеет явное визуальное подтверждение. За призмой формируется область пониженного давления с последующим его повышением с удалением от модели и формированием вихревого следа. Следовательно, можно сделать вывод, что данные, полученные в проведенных экспериментах методом визуализации и по измерениям статического давления, качественно совпадают между собой. Подробный отчет приведен в [11, 12].
Исследование теплообмена
Для исследования теплообмена моделей высотных зданий и сооружений было изготовлено 2 модели в виде квадратных призм сечением 50x50 мм и высотой 150 и 300 мм. Корпус каждой модели выполнен из текстолита толщиной 10 мм с проволочным омическим нагревателем из нихромовой проволоки диаметром 0,6 мм. Для обеспечения равномерности теплоподвода нихромо-вая проволока располагается с шагом 5 мм. Подогрев производится одновременно со всех четырех сторон модели, нагревательные секции подключаются последовательно, так что в опытах поддерживается тепловой режим qCT = const. Омические нагреватели через тонкий слой лакоткани (электроизоляция) накрываются стальными пластинами толщиной 3,0 мм.
Модель крепится к нижней стенке рабочей камеры так, чтобы ее продольная ось была перпендикулярна направлению воздушного потока. Локальная температура стенок опытной модели измеряется при помощи термопар, за-чеканенных на одной из граней каждой квадратной призмы на наружной поверхности и расположенных посредине грани в горизонтальном и вертикальном направлениях. Термопары были выполнены из хромель-алюмеля диаметром 0,2 мм и укладывались в выфрезерованные пазы заподлицо с наружной поверхностью грани. Места заделки термопар тщательно зачеканивались. Термопары соединяются с модулем Termo Lab-16v2.5, который в свою очередь соединен
с компьютером. Температура воздуха в рабочей камере трубы измеряется также хромель-аллюмелевой термопарой диаметром 0,2 мм. Скорость воздушного потока перед моделью измеряется с помощью трубки Пито - Прандтля, расположенной в центре канала выше по потоку от модели. Показания трубки Пито регистрируются наклонным микроманометром ММН-2400 (5)-1.0, а затем по известной зависимости рассчитывается скорость.
Общий вид моделей представлен на рис. 13. Места расположения термопар представлены на рис. 14. Блок-схема экспериментальной установки для измерения локальной и средней теплоотдачи представлена на рис. 15.
а б
Рис. 13. Общий вид испытуемых моделей квадратных призм сечением 50x50 мм: а - высотой 150 мм; б - высотой 300 мм
+ 30 + 29 + 28 + 27 + 26 + 25
З23> + 22 + 21 + 20
+ 17 + 16 + 15 + 14 + 13 + 12 + 11 + 10
31 *9 34
+ 6 + 5 + 4 + 3 + 2 + 1
і_В ,а
35
а
Ка"
Рис. 14. Схема расположения термопар на передней грани моделей квадратных призм: а - высотой 300 мм; б - высотой 150 мм
Рис. 15. Блок-схема экспериментальной установки для измерения теплоотдачи:
1 - входное сопло; 2 - испытуемая модель; 3 - трубка Пито - Прандтля; 4 -микроманометр; 5 - регулируемый однофазный автотрансформатор; 6 - амперметр; 7 - вольтметр; 8 - термоэлектрические преобразователи (термопары); 9 - адаптер термопар; 10 - теплограф; 11 - персональный компьютер; 12 -принтер; 13 - отчет; 14 - барометр; 15 - термометр технический
Каждая исследуемая модель в виде квадратной призмы помещается в рабочую камеру аэродинамической трубы, подключаются нагревательные элементы в электрическую цепь и поддерживается заданная температура поверхности. В аэродинамической трубе устанавливается необходимый гидродинамический режим воздушного потока и задается необходимый угол атаки взаимодействия модели с окружающей средой. Поскольку измерительные термопары располагаются только на одной из боковых поверхностей призмы, то измерение теплоотдачи от других поверхностей производится последовательным поворотом модели на 90°. В итоге, измерительная процедура при фиксированном угле атаки состоит из пяти последовательных измерений для каждой боковой поверхности. Опытами охвачен диапазон чисел Рейнольдса от 16 600 до 46 600. Угол атаки воздушного потока изменяется в диапазоне от 0 до 45° с интервалом в 5°, что позволяет получить полную информацию о влиянии на тепловые потери призм при любой их ориентации относительно ветра. После определения температуры стенки и введения необходимых поправок подсчитываются локальные значения коэффициента теплоотдачи по формуле
а = qw/(Tm - ^), (1)
где qW = qэл - qпoт, а qэn - суммарная плотность теплового потока, подводимого электрическим нагревателем; qПoT - тепловые потери за счет излучения, свободной конвекции и теплопроводности. Данные теплопотери определяются расчетным образом и учитываются при обработке опытных данных. Максимальная погрешность измерения коэффициента теплоотдачи по сделанным оценкам составляет 7,3 %.
Суммарный коэффициент теплоотдачи определяется интегрированием локальных распределений по всей поверхности грани
а = 1/L2■\\аlf6X■dY, (2)
где а*/ - локальное значение коэффициента теплоотдачи. При интегрировании (2) полагается, что относительные профили температур в горизонтальных сечениях подобны и совпадают с измеренным профилем в среднем по высоте квадратной призмы значением. Такое осреднение увеличивало погрешность определения <а> до 9,7 %.
Распределение локального коэффициента теплоотдачи по высоте квадратной призмы высотой 150 мм при нулевом угле атаки и вариации числа Яе демонстрируется на рис. 16 и 17. При нулевом угле атаки (рис. 17, а) на лобовой (А-В) и кормовой (С-П) сторонах квадратной призмы в области ее основания, где формируются подковообразный вихрь и вихревая дорожка, наблюдается сильное изменение теплообмена, где он заметно выше, чем на верхней ее части.
Рис. 16. Схема расположения рассматриваемых граней экспериментальной модели в виде квадратной призмы сечением 50x50 мм и высотой 150 мм
L,mm.
145
135
125
115
105
95
85
75
65
55
45
35
25-
15
5
L, мм.
145
• А-В B(D)-C(A)
■ C-D ір =0град. Rt = l,661(f
а,Вт/м20С
135
125
115
105
95
85
75
65
55
45
35
25
15
5
о ó
б
• А-В B(D)-C(A)
• C-D
f = Оград. R'=3,321(f
а,Вт/м20С
L,mm.
145 135 125 asios
95 85 7565 5545 35 25 15 5
. _а,Вт/м20С
Рис. 17. Распределение локального коэффициента теплоотдачи по высоте квадратной призмы высотой 150 мм при разных вариациях числа Рейнольдса: а - Яе = 16 600; б - Яе = 33 200; в - Яе = 46 600
а
в
Отмеченное подобие коэффициентов теплоотдачи наблюдается и при других числах Рейнольдса и может быть использовано в инженерных оценках локальных неоднородностей температур или тепловых потоков при турбулентном обтекании квадратной призмы.
Если призма повернута к потоку на ф = 45°, то максимум теплообмена достигается на передней кромке клина. Далее, по мере роста толщины пограничного слоя и стабилизирующего влияния дельтообразных вихрей теплоотдача сильно понижается и достигает значений, близких к кормовой зоне призмы, обтекаемой нормально направленным потоком воздуха. Слабая зависимость интенсивности теплообмена в кормовой зоне от формы плохообте-
каемых тел показана в многочисленных экспериментах, и в частности M.K. Chyu, V. Natarajan, Е. Спэрроу, К. Тьен и T. Ota [З]. Если обратиться к результатам визуализации, то можно видеть много сходственных черт вих-реобразования в кормовой зоне при изменении угла атаки, что в итоге приводит к выравниванию интенсивности теплообмена.
Полученные данные моделируют теплообменные процессы зданий и сооружений высокой этажности в зависимости от температуры окружающей среды, скорости и угла атаки воздушного потока, а также местоположения здания.
Данные результаты являются новыми. Их новизна определяется постановкой задач, а также широким диапазоном углов атаки и числом Рейнольдса. Они дополняют известные данные A.A. Жукаускаса, A.H. Лейзерона, T. Iga-rashi, I.P. Castro, P.W. Bearman, R.J. Goldstein, S.Y. Yoo, M.K. Chung, B.E. Lee, H. Reiher, H. Sakata и др. [З] по визуализации теплообмена квадратной призмы при вариации углов атаки и скорости набегающего воздушного потока.
Следует отметить, что все опытные данные получены в широком диапазоне углов атаки воздушного потока и чисел Рейнольдса для всех граней моделей, что позволяет обобщить и построить корреляционные зависимости для локальных и интегральных теплопотерь опытных моделей квадратных призм, а также сопоставить эти данные с имеющимися в литературе.
БИБЛИОГРAФИЧЕСКИЙ СПИСОК
1. Гныря, А.И. Результаты визуализации течения воздушного потока вдоль ряда из двух кубов, расположенных на плоскости друг за другом / A. И. Гныря, В. И. Терехов, С.В. Коробков // Вестник Т^СУ. - 1GG9. - № З. - С. 117-114.
1. Гныря, А.И. Результаты визуализации течения воздушного потока вдоль одиночного куба, установленного на плоскости / A. И. Гныря, В. И. Терехов, С. В. Коробков // Вестник Т^СУ. - 1GG9. - № 4. - С. 17G-178.
3. Жаркой, Р. А. Доступный ресурсо-энергосберегающий дом малоэтажной застройки / РА. Жаркой // Проектирование и строительство в Сибири. - Новосибирск, 1GG9. -№ 1 ^G). - С. 4б^.
4. Жаркой, Р. А. Малоэтажный ресурсо-энергосберегающий дом. Уменьшение тепловых потерь с вентиляционным воздухом / РА. Жаркой // Молодая мысль: Наука. Технологии. Инновации. Материалы I Межвузовской научной конференции студентов, магистрантов, аспирантов и молодых ученых 17-19 марта 1GG9 г. - Братск, 1GG9. - С. 14-18.
З. Терехое, В.И. Вихревая картина турбулентного обтекания и теплообмен одиночного куба на плоской поверхности при различных углах атаки / В.И. Терехов, A.K Гныря, С.В. Коробков // Теплофизика и аэромеханика. - 1G1G. - Т. 17. - № 4. - С. З11-ЗЗЗ.
6. Влияние угла атаки набегающего потока воздуха на процесс теплообмена тел, имеющих форму квадратной призмы / A.K Гныря, В.И. Терехов, С.В. Коробков [и др.] // Вестник Т^СУ. - 1G1G. - № З. - С. 1Зб-147.
7. Кошин, А.А. Aнализ существующих материалов по конвективному теплообмену плохообтекаемых тел [Электронный ресурс] / A.A. Кошин, С.В. Коробков, A.K Гныря // Перспективы развития фундаментальных наук : труды VII Международной конференции студентов и молодых учёных. Россия, Томск, lG-ІЗ апреля 1G1G г. / под ред. Г.В. Лями-ной, ЕА. Вайтулевич. - Томский политехнический университет, 1G1G. - С. б11-б14.
8. Сигимое, И.Г. Aнализ параметров ограждающих стен зданий в условиях Сибирского региона [Электронный ресурс] / И.Г. Сигимов, РА. Жаркой, A.K Гныря // Перспективы развития фундаментальных наук : труды VII Международной конференции студентов и молодых учёных. Россия, Томск, lG-ІЗ апреля 1G1G г. / под ред. Г.В. Ляминой, ЕА. Вайтулевич. - Томский политехнический университет, 1G1G. - С. ббІ-ббЗ.
9. Жаркой, Р.А. Натурные исследования теплозащитных свойств ограждающих конструкций малоэтажных зданий [Электронный ресурс] / Р.А. Жаркой, С.В. Коробков, А.И. Гныря // Перспективы развития фундаментальных наук : труды VII Международной конференции студентов и молодых учёных. Россия, Томск, 20-23 апреля 2010 г. / под ред. Г.В. Ляминой, Е.А. Вайтулевич. - Томский политехнический университет,
2010. - С. 592-594.
10. Экспериментальные исследования аэродинамики моделей зданий при помощи сажемасляной визуализации [Электронный ресурс] / А.И. Гныря, С.В. Коробков, А.А. Кошин [и др.] // Перспективы развития фундаментальных наук : труды VIII Международной конференции студентов и молодых учёных. Россия, Томск, 26-29 апреля 2011 г. / под ред. Г. В. Ляминой, Е. А. Вайтулевич. - Национальный исследовательский Томский политехнический университет, 2011. - С. 555-557.
11. Исследование структуры течения воздушного потока путем измерения полей статического давления [Электронный ресурс] / А.И. Гныря, С.В. Коробков, А.А. Кошин [и др.] // Перспективы развития фундаментальных наук : труды VIII Международной конференции студентов и молодых учёных. Россия, Томск, 26-29 апреля 2011 г. / под ред. Г. В. Ляминой, Е. А. Вайтулевич. - Национальный исследовательский Томский политехнический университет, 2011. - С. 570-572.
12. Аэродинамическая структура обтекания моделей зданий. Сажемасляная визуализация и поля давлений / А.А. Кошин, Д.И. Мокшин, С.В. Коробков [и др.] // Сборник трудов Международ. научн. конф. «Интеграция, партнерство и инновации в строительной науке и образовании», Москва, 19-21 октября 2011 г. - М. : Изд-во ФГБОУ ВПО «МГСУ». -
2011. - Т. 1. - С. 138-143.
