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

Гусев / Методы научных исследований

.pdf
Скачиваний:
172
Добавлен:
22.05.2015
Размер:
1.43 Mб
Скачать

Vo - объемное количество воздуха, поступающего удаляемого из помещения в единицу времени (воздухообмен);

tух и tпр - температуры воздуха соответственно удаляемого (уходящего) из помещения и поступающего (приточного) в него.

2. Уравнения расходов воздуха и потерь давления p

ваэрационных проемах:

-объемный расход:

Vo = ωo Fo ;

(7.16)

- массовый расход:

 

Go = ρωo Fo = μFo 2 pρ ;

(7.17)

 

ω 2

 

p =ξ

0

ρ ,

(7.18)

2

 

 

 

где ωo - скорость воздуха в проеме;

Fo - площадь аэрационного проема; μ - коэффициент расхода проема;

ξ - коэффициент местного сопротивления проема.

3. Закономерности развития осесимметричных конвективных (при соотношении сторон источников тепла в плане не более 1:3) и плоских приточных (при соотношении сторон приточных отверстий более 1:10) струй на основных участках.

Для конвективных потоков в направлении координаты развития струи Z :

 

 

 

βg 1/ 3

1/ 3

 

5 / 3

 

Vz

 

 

 

 

 

Q

z

 

(7.19)

 

 

 

 

c

 

ρ

к

 

 

 

 

 

 

p

 

в

 

 

 

 

 

 

 

βg 1/ 3

1/ 3

 

1/ 3

 

ω

 

 

 

 

 

 

 

Q

z

 

(7.20)

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

z

cp ρ в

 

 

 

 

161

tz = (tz tВ) (βgc2 p ρ2 )в1/ 3 Qк2 / 3 z5 / 3

(7.21)

Для плоских (приточных) струй в направлении развития по оси x по отношению к соответствующему параметру в живом сечении отверстия (индекс o ):

 

 

Vх

=

1,41

 

 

 

x

 

 

 

(7.22)

 

Vo

m

bo

 

 

 

 

 

 

 

 

 

 

ωх

= m

bo

 

 

 

 

 

(7.23)

 

 

 

 

 

 

 

 

 

ωo

 

 

 

 

x

 

 

 

 

 

 

 

tх

 

tх tВ

= n

bo

 

 

to =

 

 

 

,

(7.24)

 

tпр tВ

x

где Qк

- тепловая мощность конвективного источника тепла;

bo

- размер (ширина)

 

плоского отверстия для выпуска

приточной струи;

 

 

 

 

 

 

 

 

 

 

 

 

m -

коэффициент изменения скорости по оси струи;

n - коэффициент изменения избыточной температуры в струе.

4. Уравнения физических констант воздуха в зависимости от определяющей температуры:

 

β =

 

1

 

 

;

 

 

 

 

 

 

 

 

 

(7.25)

 

273 + t

 

 

 

 

 

 

 

 

 

 

 

λ = (2,45 + 0,0076t)102 ;

 

 

(7.26)

 

v = 0,415 108

(237 + t)2,5

;

 

(7.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

387 + t

 

 

 

 

 

353

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ =

 

,

 

 

 

 

 

 

(7.28)

 

 

273 + t

 

 

 

 

 

 

5. Уравнения для определения аэрационного воздухо-

обмена

Vo

и

избыточной

температуры

 

уходящего

воздуха

tух :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2gβ 1/ 3

1/ 3

 

 

 

1/ 3

 

 

(μF)пр2 (μF)выт2

1/ 3

 

V

=

 

 

 

Q

 

Н

 

 

 

 

 

 

;

(7.29)

 

 

 

 

(μF)2

+(μF)2

o

c ρ

 

изб

 

 

 

o

 

 

 

 

p

в

 

 

 

 

 

 

 

 

пр

выт

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

162

 

 

 

 

 

 

 

 

 

 

 

 

2

2

1/ 3

 

t

= t

t

= (2gβc2

ρ2 )1/ 3Q2 / 3

Н1/ 3

 

(μF)пр(μF)выт

 

,(

(μF)2 +(μF)2

ух

yx

пр

p

в

изб

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр

выт

 

7.30)

где Qизб - теплоизбытки в помещении; Ho - размер аэрационных проемов;

«пр» и «выт» - параметры, относящиеся соответственно

кприточным и вытяжным аэрационным проемам.

Вусловиях установившихся тепловых и аэродинамиче-

ских процессов и несжимаемой среды имеем:

τt = ωτ = ρτ = 0

В этом случае анализ уравнений (7.19-7.30) методом теории подобия позволяет получить критерии, определяющие процесс аэрации:

- критерий Релея Ra

Ram (Gr Pr)m = β glλ3cvρ ρ m (tпн tв)

-критерий Архимеда Ar

Arв = (β)в ωgl2 (tпн tв )

-критерий Рейнольдса Re Re в = (ωv)lв

-критерий Эйлера Eu

Euв = (ρ)вpω2

- критерий Нуссельта Nu

 

akl

n

Num =

 

= c Ram

(λ)m

Анализ полученных выражений показывает, что выполнить подобие всех условий однозначности и равенство всех

163

определяющих критериев в натуре и модели, в данном случае физически невозможно.

Например, согласно требованиям теории полного подобия необходимо соблюсти равенство критериев Релея в натуральном объекте и в модели - [Ram ]Н = [Ram ]М :

 

 

gc

ρ

ρ

l3 (t

t

 

 

 

 

 

gc

ρ

ρ

l3 (t

 

t

 

 

 

β

 

 

 

)

 

=

 

β

 

 

 

 

)

 

λv

λv

 

 

 

 

 

пн в

 

 

 

 

 

 

пн

в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

М

 

 

 

 

 

 

 

 

 

 

Н

Отсюда определим масштабы подобия Сi для всех однородных физических величин:

 

С

С С С

 

 

3

 

 

β

g c

ρ

 

 

 

 

 

 

Cl Ctt =1 ,

(7.31)

 

С С

 

 

 

λ v

 

m

 

 

где Ctt - масштаб разности температур (tпн tв ) .

Если допустить, что при одинаковой рабочей среде в модели и натуре, имеет место равенство коэффициентов тем-

пературного расширения βН = βМ , плотностей

среды

ρН = ρМ , коэффициентов теплопроводности λМ = λН

и кине-

матической вязкости νН =νМ , а так же ускорений и теплоемкостей воздуха (Сg = Cc =1), то из ( 7.31) следует:

Ctt = 1Сl3 или (tпн tв)М = (tпн tв)Н Сl3

А это значит, что если принять линейный масштаб модели хотя бы равным Сl = 5 , то для того чтобы обеспечить условие подобия в модели следует увеличить, разность температур поверхности источника тепла и окружающего воздуха (t- tв) в 125 раз по сравнению с натурой, что технически невыполнимо.

То же самое можно сказать и относительно скоростей воздуха в модели. При том же линейном масштабе они должны быть больше, чем в объекте исследований, в 5 раз.

Для выхода из создавшейся ситуации в теории моделирования вводится понятие автомодельности.

164

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

Так, например, если в натуре режимы естественных и вынужденных течений турбулентные ([Ram ]Н [Ram ]Кр ;

[ReВ ]Н [ReВ ]Кр ), то в модели их подобие обеспечивается условием:

[Ram ]М [Ram ]Кр

[ReВ]М [ReВ]Кр .

Это дает возможность ввести понятие о масштабах этих критериев:

СRa

=

[Ra

m

]

М

[Ram ]Кр

 

 

 

[Ram ]Н

 

[Ram ]Н

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

СRem

=

 

[Re

в

]

М

 

[Reв ]Кр

 

.

 

[Reв ]Н

[Reв ]Н

 

 

 

 

 

 

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

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

165

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

7.3. Особенности моделирования стохастических систем

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

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

Принципы подобия в статистическом смысле требуют выполнение следующих условий [14]:

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

166

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

Масштаб подобия в стохастических натуральном и модельном объектах, характеризующиxся случайными величи-

нами, соответственно X Н.1, X Н.2 ,....., X Н.n

и X М.1, X М.2 ,....., X М.n

с плотностями распределения

f (X Н ) и

f (X М ) рассматрива-

ется как отношение математических ожиданий Xv соответст-

венных случайных величин:

 

 

CvX .i.r =

XwН.i

,

i =1,2,....., n .

 

X М.i

 

 

Подобие в поведении систем будет соблюдаться при подобии математических ожиданий (начальных моментов). Этот вид подобия называется статистическим подобием первого порядка [15].

Условия подобия более высоких порядков выражаются через подобие центральных моментов. Для обеспечения подобия в статистическом смысле необходимо, что бы каждая константа удовлетворяла условию:

 

μh1 ,......,hn

( X Н )

n

1/ hk

CvX .i.r =

CvXh.ii.r

 

 

 

( X М )

 

1

 

n

 

 

 

μh

,......,h

i =1

 

 

 

 

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

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

167

Из теории размерности следует следующее соотноше-

ние:

v

K

l

v

 

 

I j =

v

= CXaij.i.r =1

j =1,2,3......m

(7.32)

 

K

i=1

 

где K - критерий подобия; j - номер критерия;

a - показатели степени, устанавливаемые на основе анализа размерностей;

l- число величин, определяющих явление.

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

вие:

 

 

Ivj (1

)I j I j (1

),

 

 

 

 

 

(7.33)

где

- некоторая вероятность.

 

 

 

 

 

Вероятность условия (7.33) определяется выражением:

P(

 

I j Ivj

 

< ,.........j =1,2,.....,m)=

 

 

 

 

 

..... f1,2,...,m (I1, I2,....Im )dI1,dI2,...,dIm

(где

f1,2,...m (I1, I2 ,...Im )

- плотность

совместного распределения

индикаторов подобия I j ) и называется вероятностью подобия при погрешности подобия . Вероятность подобия по (7.33) определяет подобие поведения одной наугад взятой модели и одного взятого наугад натурального образца.

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

168

Для некоррелированных индикаторов подобия можно записать:

 

 

 

 

f1,2,...,m (I1, I2 ,..., Im )= f1 (I1 ),

f2 (I2 ),...., fv (Im )

 

 

 

 

 

и вычисление вероятности подобия проводится по выраже-

нию:

 

 

 

 

 

 

 

 

 

 

P(

 

I j Ivj

 

< ,......... j =1,2,....., m)= m

P(

 

I j Ivj

 

< )

= m f j (I j )dI j

 

 

 

 

 

 

 

 

j=1

 

 

 

 

j

j=1

 

I j I j

 

<

 

 

 

 

 

 

 

 

 

 

Числовые характеристики плотности распределения в этом случае представляют собой следующее:

- для математического ожидания:

I j = 1

- для дисперсии

 

 

2l

D(vX2 i

)

2l

 

 

 

D(I j )= aij2 *

= aij2ω2

(Xi

).

 

 

 

 

i =1

Xi

 

i =1

 

 

где ω(Xi )=

Dv(Xi )

- коэффициент вариации величины X i .

 

Xi

 

 

 

 

 

 

Вычисление вероятности подобия при некоррелированных индикаторах в предположении о логарифмическом нормальном распределении плотности вероятности f (I ) можно провести по следующему выражению:

P(I j Ivj < ,.. j =1,2,..,m)= hj = Ф* lg(1 +σuj)lg I0 Ф* lg(1 σuj)lg I0 lg I0 = −1,1513σuj2

Среднеквадратичное отклонение находится из решения уравнения:

 

 

 

σ2

 

 

ω2 (I

j

)= exp

uj

 

1

 

 

 

0,1886

 

 

 

 

 

 

 

Значение функции Ф* = Ф + 0,5 определяется по стандартным таблицам, например [4].

169

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

I j 1

< ,

j = 1,2 ,..., m

предполагает многозначность значений индикатора подобия в интервале (1 ;1 + ).

Определение возможного интервала (X Нmax X Нmin ) значений какой-либо величины X Н.k для натурального образца при условии, что значение этой же величины X М.k измерено на модели и известен интервал возможных значений множителя преобразования (масштаба) CX .r при заданной вероятно-

сти подобия h и погрешности

 

определяется следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

1

 

 

aij1

v

X Н.k

> X M .k

 

 

 

 

 

 

CX .r

1 +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.34)

 

min

 

 

 

 

1

 

aij1

v

 

X

 

> X

 

 

 

 

 

 

C

X .r .

 

 

 

 

 

 

Н.k

 

M .k 1

 

 

 

При этом максимальное и минимальное значение множителя перехода представляют отношения:

X max Ckrmax = Xkmin

k

Если принять экспериментально найденное (осредненное по некоторому количеству параллельных опытов) значение величины в модели X&&&М.k в качестве математического

170