Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

6854

.pdf
Скачиваний:
0
Добавлен:
23.11.2023
Размер:
901.92 Кб
Скачать

Министерство образования и науки Российской Федерации

Федеральное государственное бюджетное образовательное учреждение высшего образования

«Нижегородский государственный архитектурно-строительный университет»

Кочев А.Г.

МОДЕЛИРОВАНИЕ ПРОЦЕССОВ В СИСТЕМАХ ТЕПЛОГАЗОСНАБЖЕНИЯ И ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТАЛЬНЫХ ИССЛЕДОВАНИЙ

Учебно-методическое пособие по подготовке к практическим занятиям

(включая рекомендации по выполнению самостоятельной работы)

для обучающихся по дисциплине «Моделирование процессов в системах теплогазоснабжения и планирование экспериментальных исследований» по направлению подготовки 08.04.01 Строительство

профиль Возобновляемые источники энергии и энергоэффективность в зданиях

Нижний Новгород ННГАСУ

2016 г

2

УДК 726.5+536.24 (075.8)

Кочев А.Г. / Моделирование процессов в системах теплогазоснабжения и планирование экспериментальных исследований . [Электронный ресурс]: учеб.-метод. пос. / А.Г. Кочев; Нижегор. гос. архитектур. - строит. ун-т – Н. Новгород: ННГАСУ, 2016. – 53 с.–

1электрон. опт. диск (CD-RW).

Внастоящем учебно-методическом пособии по дисциплине «Моделирование процессов в системах теплогазоснабжения и планирование экспериментальных исследований» даются конкретные рекомендации учащимся для освоения основного и дополнительного материала дисциплины для достижения целей, обозначенных в учебной программе дисциплины. Цель учебно-методического пособия — помощь в подготовке к практическим занятиям, в освоении требуемого объёма знаний самостоятельной работы студентов.

Учебно-методическое пособие предназначено для обучающихся в ННГАСУ по дисциплине «Моделирование процессов в системах теплогазоснабжения и планирование экспериментальных исследований» по направлению подготовки 08.04.01 Строительство, профиль Возобновляемые источники энергии и энергоэффективность в зданиях.

Учебно-методическое пособие ориентировано на обучение в соответствии с календарным учебным графиком и учебным планом по основной профессиональной образовательной программе направления 08.04.01 Строительство, профиль Возобновляемые источники энергии и энергоэффективность в зданиях, утверждённым решением научно-технического совета (НТС) ННГАСУ от 9.06.2015г. (протокол № 2).

©

А.Г. Кочев, 2016

©

ННГАСУ, 2016

 

3

 

 

Содержание

 

 

Введение

5

1.

1. Методы термодинамического и физико-математического описа-

 

 

ния процессов в системах ТГС

5

 

1.1. Уравнение Лапласа

6

 

1.2. Решения задач на вдув. Уравнение Фурье

12

 

1.3. Численное моделирование теплообменных процессов

15

2.

Основы теории подобия. Моделирование тепловых и гидро-

 

 

аэродинамических процессов

30

 

2.1. Основы теории подобия

30

 

2.2. Критерии (числа) подобия

32

 

2.3. Степень турбулентности потока. Миделевое сечение

35

 

2.4. Аэролинамические трубы и гидролотки

36

3.

Статистический анализ эксперимента. Планирование эксперимен-

 

 

тальных исследований

41

 

Библиографический список

47

4.

Вопросы для контрольных работ

48

4

ВВЕДЕНИЕ

В процессе передачи тепла через светопропускающее ограждение от внутреннего воздуха к наружному участвуют одновременно три вида теплообмена (теплопроводность, конвекция, излучение), которые влияют друг на друга. Для математического описания процесса переноса тепла при стационарном режиме можно использовать дифференциальные уравнения теплопроводности и конвекции (уравнение Лапласа).

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

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

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

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

1. Методы термодинамического и физико-математического описания процессов в системах ТГС.

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

Ro = R в

+ ΣRст + R вп + R н.

(1.1)

Теплозащитные качества

светопропускающих

ограждений в

5

основном зависят от термического сопротивления воздуха у внутренней поверхности остекления Rв и термического сопротивления воздушных прослоек Rвп. В окнах с двойным остеклением термическое сопротивление Rв составляет 30÷35 % от общего сопротивления теплопередаче, термическое сопротивление Rвп составляет 45÷55 %. При тройном остеклении Rв составляет 20÷25 %, Rвп составляет 60÷70 %. Наибольшую трудность при определении Rв и Rвп представляет конвективная составляющая теплового потока.

 

t

x

 

 

 

 

 

 

 

 

 

 

 

 

 

t В

 

 

 

 

 

 

tВ

 

1

 

 

 

 

 

 

 

 

 

 

FН

 

 

В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τВ

τ1

F2

 

 

 

 

tСР1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FН

 

 

 

 

 

 

 

 

 

 

 

 

 

 

FВ

 

 

 

 

 

 

FВ

2

 

3

 

 

 

F1

tСР ОК

 

Н

 

 

 

 

 

tСР2

 

 

Н

 

q

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ 2

 

τН

 

 

B

F2

 

 

 

 

 

 

 

В

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н

 

vy

 

 

 

 

t Н

 

 

F1

 

F3

F4

 

 

 

 

vx

 

 

 

 

 

V

 

 

 

 

 

 

 

y

R1

R*ВП

R**ВП

R2

RН

 

 

 

 

 

 

 

 

 

 

R В

y

Rв

R 1

R''R

 

R''' R'''' R

 

R

 

 

 

 

 

 

 

2

3

н

 

 

 

 

 

 

 

 

 

вп

вп

вп

 

 

 

 

 

 

 

 

 

R'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вп

 

 

 

 

 

 

Рис. 1.1 Схема передачи теплоты через остекление оконного проема

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

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

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

1.1. Уравнение Лапласа

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

6

t н= tж2

2 = н t пов н

tж( ), ( )

 

ОК

о

t( ;)

tв= tж1

1 = в t пов в

Рис. 1.2. Схема расчетной области исследований

Задача решается для двумерного поля температур в стационарных условиях, т. е. требуется решить дифференциальное уравнение Лапласа:

2

t

+

2

t

= 0 .

(1.2)

 

 

2

 

 

2

 

 

x

y

 

Задача решается в следующих пределах: x > 0; 0 < y < δо .

Чтобы найти численные значения температур (температурное поле) сразу внутри и на поверхности исследуемой области, необходимо знать краевые условия – размеры и геометрию исследуемой области, теплофизические характеристики применяемых материалов и закон взаимодействия между окружающей средой и поверхностью области (граничные условия).

Записываем граничные условия:

 

t(

 

 

 

 

 

 

= α × [t

ж1 - t(

 

,0)]

 

 

 

 

- λ ×

 

 

x,0)

 

 

 

 

 

 

x

;

 

(1.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

y

2 × [t(

 

 

o )- t

ж2 ] ;

 

 

t(

 

 

o ) = α

 

 

- λ ×

x

 

x

(1.4)

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t(0,

 

) = α(

y) × [t

ж (y) - t(0,

 

)] ;

 

- λ ×

y

 

y

(1.5)

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

t ж1 - t ж2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

t(¥, y) = t ж1

-

 

 

 

 

× (

 

+

)

.

(1.6)

 

 

 

δо

 

 

 

λ

 

 

 

 

 

1

+

+

1

 

 

 

 

α1

 

 

 

 

 

α1

λ

α2

 

 

 

 

 

 

Граничное условие (1.6) отвечает стационарному одномерному полю температуры вдали от оконного проема.

Принимается, что температура внутреннего воздуха tв, коэффициент теплоотдачи αв отк и градиент температуры в поперечном сечении откоса являются постоянными величинами. Потерями тепла через верхние и нижние части откоса пренебрегаем.

Температура tж (y) является известной функцией изменения температуры воздуха, ограничивающего поверхность оконного откоса и

 

 

 

 

 

стены, и определяется

 

 

 

 

 

образом: при

 

 

= 0 :

наружной

следующим

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tж (y) = tж1

= tв ; при y = δо : tж (y) = t ж2

= tн ; при x = 0

и 0 < y < δ: t ж (y) = tв ;

 

 

 

= 0

и δ <

 

 

 

 

 

 

 

 

при

x

y

< (δ + δок ) : tж (y) = tср ок ; при x = 0 и (δ + δок ) < y < δо :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tж (y) = tн .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент теплоотдачи α(y) является известной функцией на поверхностях оконного откоса и наружной стены; определяется следующим образом: при y = 0 : α(y) = α1 = αв ; при y = δо : α(y) = α2 = αн ;

при x = 0 и 0 < y < δ: α(y) = αв отк ; при x = 0 и δ < y < (δ + δок ) : α(y) = αок отк ; при x = 0 и (δ + δок ) < y < δо : α(y) = αн отк .

Численные значения соответствующих коэффициентов принимаются с учетом исследований К.Ф.Фокина и В.Н.Богословского для вынужденной и естественной конвекции.

Величина αн для конкретного оконного проема может определяться через значение аэродинамического коэффициента Сv, для чего находится значение скорости в расчетный период vo, а затем определяется скорость обтекания стены вблизи оконного проема v. Зная эти скорости, можно определить величину αн для условий вынужденной конвекции.

Выбрав в качестве масштаба длины толщину стены δо, а масштаба температуры – температурный напор (tж1-tж2), приведем уравнение Лапласа (1.2) и граничные условия (1.3)-(1.6) к безразмерному виду:

2θ

+

2θ

= 0 .

(1.7)

x 2

y2

 

 

 

При этом x>0, 0<y<1.

Граничные условия запишутся в виде:

-

∂θ(x,0) = Bi1 × [1 - θ(x,0)] ;

(1.8)

 

y

 

 

8

 

 

-

∂θ(x,1) = Bi

2 × θ(x,1)

;

 

y

 

 

∂θ(0, y) = Bi(y) × [θ(0, y) - θ

ж (y)] ;

x

Bi1 × y + 1

 

 

θ(¥, y) =1 -

= F(y) .

1 + Bi1

+

Bi1

 

 

 

 

Bi2

 

 

 

 

 

Здесь

(1.9)

(1.10)

(1.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ(x, y) =

t(x,

 

y) - t ж2

;

 

 

 

 

(1.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t ж1 - t ж2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θж (y) =

t ж (y) - t ж2

; 0 £ θж (y) £1 ;

(1.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t ж1 - t ж2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α1 × δo

 

 

α2 × δo

 

α(

 

 

 

 

x =

 

x

; y =

 

y

; Bi1

=

; Bi2

=

; Bi(y) =

y) × δo

.

(1.14)

 

 

 

 

 

 

 

 

 

δo

δo

 

λ

 

 

 

 

λ

 

λ

 

 

Таким образом, безразмерное температурное поле зависит от

следующих величин

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ = θ(x, y,Bi1 , Bi2 , Bi(y),θж (y)) .

 

 

 

 

(1.15)

 

Будем искать решение уравнения (1.15) в виде суммы двух

функций:

 

 

 

 

 

θ(x, y) = u(x, y) + v(x, y) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.16)

Подставим (1.16) в уравнение Лапласа (1.7) и граничные условия

(1.8)-(1.11):

 

 

 

 

 

2 v +

2 v = -

2 u -

2 u

= 0; x > 0; 0 < y <1;

x 2

y2

x 2

y2

 

 

∂v(x,0) - Bi1 × v(x,0) = -

∂u(x,0) + Bi1

× u(x,0) ;

y

 

 

 

y

 

∂v(x,1) + Bi2 × v(x,1) = -

∂u(x,1) - Bi2

× u(x,1) ;

y

 

 

 

y

 

∂v(0, y) - Bi(y) × v(0, y) = -

∂u(0, y) + Bi(y) ×[u(0, y) - θ

ж (y)];

x

x

 

v(∞, y) = −u(∞, y) + F(y) .

(1.17)

(1.18)

(1.19)

(1.20)

(1.21)

Подберем функцию u(x,y) такой, чтобы она удовлетворяла уравнению Лапласа (1.7) и доставляла функции v(x,y) на всех границах, кроме границы х=0, однородные граничные условия. Этому условию удовлетворяет функция u(x,y)=F(y). Действительно, подставив ее в (1.17)- (1.21), получим

9

2 v

+ 2 v = 0 ;

 

x 2

y2

 

 

∂v(x,0) = Bi1

× v(x,0)

;

y

 

 

 

∂v(x,1) = -Bi2

× v(x,1)

;

y

 

 

 

∂v(0, y) = Bi(y) ×[v(0, y) + F(y) - θ

ж (y)] ;

x

 

v(¥, y) = 0 .

 

(1.22)

(1.23)

(1.24)

(1.25)

(1.26)

Таким образом, неизвестная функция v(x,y) представляет собой отклонение от известного стационарного неодномерного поля температуры вдали от оконного проема. Будем искать решение задачи (1.22)-(1.26) способом Фурье. Представим искомое решение в виде произведения двух функций, одна из которых зависит только от координаты х, а другая – от координаты у

v(x, y) = X(x) × Y(y) .

(1.27)

Подставив (1.27) в (1.22), получим:

··

 

 

 

 

Y(y)

= -

X // (x)

= -ε2 .

(1.28)

Y(y)

 

 

X(x)

 

Вданном случае точки означают дифференцирование по х, штрихи

по у; e2=const, так как левая часть равенства (1.28) зависит только от у, а правая – только от х.

Из равенства (1.28) следует

X // (x) - ε2

× X(x) = 0 ;

(1.29)

··

× Y(y)

= 0 .

 

Y(y) + ε2

(1.30)

Решением (1.29) будет функция

 

 

 

X(x) = с1 × еε×x + c2

× e-ε×x .

(1.31)

Решение (1.30) будет иметь вид

 

 

 

Y(y)= с3 ×sin(ε× y)+ c4 ×cos (ε× y).

(1.32)

Подставим (1.31) и (1.32) в (1.27); получим

 

v(x, y) = (c1 × еε×x + c2 × e-ε×x ) × (c3 × sin (ε × y) + c4 × cos (ε × y)) .

(1.33)

Из условия ограниченности (1.26) функции v(x,y) на бесконечности

постоянная с1=0. Тогда (1.33) перепишется в виде

 

v(x, y) = (A × sin (ε × y) + B × cos (ε × y)) × e-ε×x .

(1.34)

В данном случае А и В – постоянные величины: А=с2×с3; В=с2×с4. Тем самым задача полностью решена, и безразмерное

температурное поле (1.15) имеет вид:

 

 

10

 

 

 

 

 

θ(x, y) =1 -

Bi1 × y + 1

 

+ ¥

An

×ϕn (y) × e-εn ×x

. (1.35)

 

 

Bi1

 

 

 

1 + Bi1

+

 

n=1 Bi1

 

 

Bi2

 

 

 

 

 

 

 

 

 

 

 

 

 

Анализ температурного поля (1.35) показывает, что при малых значениях безразмерной продольной координаты х от х=0 до х=х1 распределение температуры и скорости ее изменения по х в отдельных точках стены зависит от особенностей граничных условий при х=0. Математически в этой области поле температуры в стене будет определяться не только первым, но и последующими членами ряда (1.35). Эту область стены, непосредственно примыкающую к оконному проему, в которой скорость изменения температуры в направлении координаты х зависит от вида граничных условий при х=0, назовем по аналогии нестационарными температурными полями неупорядоченной областью теплопередачи через стену.

Начиная с некоторого расстояния х>х1 от оконного проема, граничные условия при х=0 начинают играть второстепенную роль, и процесс теплопередачи через стену полностью определяется только граничными условиями при у=0 и у=1, физическими свойствами материала стены и ее размерами и формой, хотя температурное поле вдоль стены и изменяется. Температурное поле в этой области х>х1 математически описывается первым членом ряда (1.35)

θ(x, y) =

Bi1 × y + 1

 

+

A1

×ϕ1 (y) × e-ε1 ×x .

(1.36)

 

 

Bi1

 

 

 

1 + Bi1

+

 

 

Bi1

 

 

Bi2

 

 

 

 

 

 

 

 

 

 

 

 

 

Безразмерный поток тепла, подводимый к стене на этом участке в предположении одномерного температурного поля, равен

 

 

 

 

Qo = -

 

Bi1 × x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bi

.

 

 

 

 

(1.37)

 

 

 

 

 

1 + Bi1

+

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Bi2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Их отношение соответственно равно

 

 

 

 

 

 

 

 

 

 

 

 

(1 + Bi1

+ Bi1 ) × ∑ An

/

(0) × (1 - e

-εn

×x

) - εn

× ∫ϕn (y) × dy]

 

 

×[ϕ n

 

 

 

¥

 

 

 

 

 

 

 

2

 

2

y1

 

 

Qд =1 +

 

 

 

εn

 

 

 

 

 

 

 

0

.

 

Bi2

 

 

 

 

 

 

 

 

 

 

 

 

n =1

 

 

 

 

 

 

 

 

 

 

 

 

Qo

 

 

 

 

Bi12 × x 2

 

 

 

 

 

 

(1.38)

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]