книги из ГПНТБ / Борисенко А.И. Аэродинамика и теплопередача в электрических машинах
.pdfРейнольдс в 1883 г. показал, что значение коэффици ента £ зависит от безразмерной комбинации средней ско рости течения wm, диаметра трубы сі, плотности р и вязкости ,р текущей жидкости. Безразмерная комбина ция из этих .параметров может быть .построена единст венным способом; она называется числом Рейнольдса
Ре = і ^ . |
(1-42) |
Рейнольдс показал также, что от величины этого чис ла зависит режим течения жидкости в трубе. Пока число Re меньше некоторого критического значения, движение является прямолинейным слоистым, что хорошо видно с помощью окрашенной струйки. В потоке она отчетливо видна в виде резко выделяющейся нити (рис. 1-9). Не большие колебания могут иметь место, однако целост ность нити не нарушается. Если же число Re превышает критическое значение, то нить расплывается и равномер но окрашивает жидкость в трубе. Сопротивление дви жению при этом резко возрастает. Критическое значение числа ReKp приблизительно равно 2 300. Это значение ReKp называют нижним критическим числом Рейнольдса, так как рядом позднейших опытов было показано, что принимая специальные меры против возникновения на чальных возмущений, критическое число Re можно по высить до 25 000 и более. Это значение называют верх ним критическим числом.
Движение жидкости при докритических значениях числа Рейнольдса называют ламинарным, при этом тече ние .происходит слоями. При превышении критического значения числа Рейнольдса появляющиеся возмущения начинают развиваться, движение частиц (достаточно больших по сравнению с молекулой) принимает беспо рядочный характер, скорости в каждой точке становятся различными как по величине, так и по направлению, колеблясь около некоторого среднего значения. Такой
Рис. 1-9. Прибор Рейнольд са для наблюдения за ре жимом течения жидкости в трубе.
режим |
течения |
называется |
|
||
турбулентным. |
|
Беспорядоч |
|
||
ность и хаотичность пульсаций |
|
||||
скорости ведет к усилению пе |
|
||||
ремешивания жидкости, к вы |
|
||||
равниванию скоростей |
(рис. |
|
|||
1-10), но вместе с тем это пе |
|
||||
ремешивание требует и боль |
|
||||
шей затраты энергии, так что |
|
||||
сопротивление |
движению при |
Рис. 1-10. Распределение |
|||
турбулентном |
режиме |
возра |
скоростей в трубе при ла |
||
стает (рис. 1-11). |
|
течении |
минарном (/) и турбу |
||
При |
ламинарном |
лентном (2) режимах тече |
|||
теплопередача |
в |
направле |
ния. |
||
|
нии, .нормальном к движению, осуществляется молеку лярным обменом энергией — теплопроводностью. При турбулентном режиме течения в очень тонком слое, не посредственно прилегающем к стенке, поперечные пуль сации затухают и распределение скоростей линейное (рис. 1-12), поэтому его называют ламинарным или вяз ким подслоем. В этом подслое перенос тепла соверша ется только теплопроводностью. Вне этого слоя наиболь шее значение имеет турбулентное перемешивание жидко сти, беспорядочное поперечное движение отдельных, до статочно больших (по сравнению с молекулами) частиц жидкости.
Пульсационные скорости, наблюдающиеся при тур булентном течении, являются причиной возникновения дополнительных турбулентных напряжений. Например, за счет пульсационных скоростей w 'x, w 'y и w 'z в направ лении осей X, у, z соответственно масса жидкости, теку-
Рис. 1-11. Зависимость коэффициента сопротив ления в круглой трубе от числа Re.
I — ламинарное течение; 2 — турбулентное течение.
31
Вязкий подслой |
|
у |
|
Рис. 1-12. Ламинарный под- |
Рис. 1-13. Турбулентные напряжения |
||
слой в турбулентном по- |
на трех взаимно перпендикулярных |
||
токе. |
птлощадках |
возникают |
вследствие |
|
пульсационных скоростей |
w'x, w'y, |
W ' z .
щая в направлении оси х, приобретет дополнительный импульс рw'xw'y, что эквивалентно касательному напря жению хху (рис. 1-13). Аналогичные соотношения можно определить и для других составляющих
Так как симметричные составляющие равны между собой Хху==Хух\ Xxz~^zxl Xyz==Xzyt то различных касательных напряжений шесть.
Пограничные слои. Название «вязкий подслой» свя зано с тем, что при внешнем обтекании тел вблизи его поверхности образуется так называемый пограничный слой, в формировании которого основную роль играют касательные напряжения. Опыт показывает, что измене ние скорости вдоль нормали от значения нуль на поверх ности тела до наибольшей происходит в относительно тонком пограничном слое. Толщиной его называют рас стояние от стенки до той точки, где скорость практиче ски не отличается от скорости вдали от стенки.
Тщательно поставленные опыты с пластинкой показа ли, что если X — расстояние от ребра входа и шм— скорость невозмущенного потока, то при представлении числа Рейнольдса в виде Ре*=ргіу0х/|л,, то толщина ла
минарного пограничного слоя равна 8„ = 5,5A:/'|/Re3Cи при
32
Ре8=рдо0 8л/}х = 5,5|/^Г?еж |
3000 |
ламинарное течение пе. |
|||
реходит в турбулентное. |
|
|
|
|
|
Одинаковый порядок величины критического числа |
|||||
для внутреннего |
(около |
2 300) |
и |
внешнего |
обтекания |
(3 000) послужил |
основанием для |
построения |
аналогии |
этих течений, в которой радиус трубы соответствует тол щине пограничного слоя, а скорость на оси — скорости вне пограничного слоя. Установлено, что и структура пограничного слоя в обоих случаях одинакова — имеется вязкий подслой с линейным распределением скоростей.
Толщина турбулентного слоя 8Т— 0,37/VRe*, а вяз'
кого подслоя
При теплообмене наряду с описанным выше вязким пограничным слоем (динамическим) есть и тепловой по граничный слой, также небольшой толщины, в котором температура жидкости (газа) Тт резко изменяется до значения, равного температуре стенки Тс:
(1-44)
В целях удобства и сохранения единой формы записи принимают и при турбулентном течении
тт= Ах dwfdy; <7= — cpAq дТж!ду, |
(1-45) |
где Ах и Ад — коэффициенты турбулентного обмена, ко торые рассматривают как коэффициент (турбулентной) вязкости Ах и коэффициент теплопроводности срАч соот ветственно. При этом в общем случае
Т:= Ф + |
AJ dw/dy, <7= |
— (Я + |
£рЛ7) дТж/ду. |
(1-46) |
Следует иметь в виду, |
что в области развитой тур |
|||
булентности |
(вне вязкого |
подслоя) |
Ах > ц; срАч |
и в |
ней первыми слагаемыми обычно пренебрегают.
Законы сохранения для струйки. Пусть параметры по сечению струйки не изменяются, сечение ее изменяется плавно, средняя линия изгибается по большому радиусу. Рассмотрим сохранение массы, импульса и энергии для
жидкости, ограниченной сечениями 11—22 |
и Г1'—2'2' |
в моменты времени t и t + Ы соответственно |
(рис. 1-14). |
3 -233 |
33 |
Рис. 1-14. Переход объема струйки 11—22 и в по ложение Г Г —2'2' за время бt.
Если ki, kn, km — величины, '.переносимые жидкостью (вещества, импульса или энергии) в областях I, II, III, то в силу законов сохранения можно написать:
(£j + kn)t = (kn -f- km )i+lt
или
^ I I I ( f + 8/) |
^lt ~ k nt |
^ U ( t + U ) ’ |
T. e. изменение величины k при перемещении за время бt должно быть компенсировано изменением ее внутри объема. Если движение установившееся, то в дважды заштрихованном объеме ничто не изменится, т. е.
Äii(/+80=Â m »и вм ест0 (1_47) будет:
|
k \ u ( t + b t ) |
, |
(1-48) |
Закон сохранения массы |
для установившегося |
дви |
|
жения (d / d t |
— 0 ) в силу (1-48) |
имеет вид |
|
т. е. |
= 'р852ку28^. |
|
|
34
Сокращая на 8t, получаем: |
|
piW1Si = ß2.WzS2 или pwS = const. |
(1-49) |
Если сечения 11 и 22 взяты достаточно |
близко друг |
к другу на расстоянии 81, то в общем случае д/дІфО
можно написать (приписывая знак |
«+ » |
вытекшей мас |
|||
се) : |
|
|
|
|
|
mu\{t+bt) |
m\t ==тш |
|
т щі+ыу |
||
Но |
|
|
|
|
|
m1ІІ(<+М) mit - |
^ род St SS -j- |
|
(рм> |
8S-f-...) |
|
pw Ы SS = |
8/ SS |
д (pw) |
|
||
1 |
T "SI; |
|
|||
mi\t m U(t+lt) |
Лdt- Ы bs ы. |
Отсюда после сокращения на öS 81 Ы получим диф ференциальное уравнение неразрывности для одномер ного течения (для струйки) — здесь приведены три его формы, вытекающие одна из другой:
др |
_|_б) (рш) |
|
|
■О (б); |
dt |
öl |
|
|
|
|
|
|
||
|
|
|
-Р-5Г |
|
В общем случае трехмерного движения |
||||
|
dp |
I |
д (Рwx) I д (ршц) |
д (pmz) _ _ Q |
|
dt |
' |
dx * dy * |
dz |
dp dt
(1-50)
''
или в векторной форме
| f + VPW = 0 . |
С1'52) |
Если жидкость несжимаема (р = const), то
yw = 0.
Уравнение движения. Если через 1 обозначить им пульс (количество движения), то .по теореме импульсов изменение импульса за единицу времени равно равнодей ствующей внешних сил
а « + іін )і+м— а і + і ц ) * = 2л»- |
а-53) |
Для установившегося движения d/dt — Q и -поэтому |
|
1цг — Ь і(/+»г |
|
3* |
35 |
Поток импульса (векторная величина) за время бt через нормальное сечение струйки 5і определится как произведение массы объема (1 -1 ) (р^Ш і бt) на ско рость Wi, т. е. piSitOj бt Wi. То же для сечения 22: р2w2S2 hi w2. Приписывая знак плюс потоку импульса, выходя щему из контрольного объема И —22 (совпадающему с направлением внешней нормали п), напишем измене ние потока импульса:
р2S2W2 бt Wz—piSiWi бt wi, |
(1-54) |
которое должно быть равно импульсу результирующей всех сил (массовых и поверхностных) за тот же проме жуток времени Ы, приложенных к жидкости в контроль ном объеме 11—22. Это дает после сокращения на Ы
p2S2ay2w2 — |
= R. |
(1-55) |
Принимая во внимание уравнение непрерывности, бу |
||
дем иметь: |
|
(1-56) |
R=m (w 2—Wi) |
|
|
(или в проекциях на прямоугольные |
оси |
координат х, |
У, z) |
|
|
Rx = m{w2X — wlXy, Ry ~ m { w 2y — wiy); Rz =m {w 2Z — wlZ),
где |
|
m = p(wn)5. |
(1-57) |
Чтобы получить дифференциальное уравнение дви жения для струйки в общем случае d / d t ^ 0, перепишем левую часть, опуская знаки векторных величин, так как теперь сечения 11 и 22 должны быть взяты сколь угодно близко друг к другу и направления векторов скорости в сечениях 11 и 22 можно считать совпадающими. (При составлении уравнения следует иметь в виду, что при Ы— Ч) объем стремится совпасть с исходным объемом.)
I ІШ+М |
ІІ(< + 8<) |
■ІИ ' рw U 85 -f- |
|
+ -fif- (PW btbSw )bl-\-... |
— (pay 5185 w) -j- |
||
~ (p 85 bl w) bt = |
pay 85 ~ |
blbf~[-p bS bl |
« (i_58) |
(имея в виду, что вдоль трубки тока рw 65 = const, а во времени рб5 6/= р dV масса постоянна).
36
Скорости по сечению струйки постоянны (5ш/<3г=0) и силы трения -будут отсутствовать. Поэтому, обозначая проекцию вектора напряжения массовых сил f на ось трубки через //, правую часть последнего выражения можно записать в виде
р/, SS 8/ Ы+ jp S S - |
|
|
^ - U b S + . ы = |
|
||||||||
|
|
|
= |
[p/j — 'dläp_ |
SV Ы. |
|
|
(1-59) |
||||
Это дает дифференциальное уравнение одномерного |
||||||||||||
движения |
|
|
|
|
|
|
|
|
|
|
|
|
|
dw |
__ |
dw |
, |
|
dw |
|
t |
1 |
|
dp |
(1-60) |
|
dt |
|
dt |
' |
W |
dl |
|
' 1 |
p |
dl |
||
|
|
|
|
|||||||||
В общем случае трехмерного движения вязкой не |
||||||||||||
сжимаемой жидкости |
|
|
|
|
|
|
|
|
||||
dwx |
dwx |
, |
|
dwx |
I |
|
dw. |
|
|
dwx |
|
|
dt |
dt |
|
|
|
|
|
|
dy |
|
■® *-ar = |
|
|
|
1 |
dp |
|
|
d2wx |
|
d2wx |
fI . |
d2wx |
|
||
■ ■ и - |
p |
dx |
- И (■ dx2 |
Г |
dy2 |
dz2 |
|
|||||
dwy |
dwy |
, |
|
dwy |
I |
. |
ISUÜTldw» |
I |
dwy |
|
||
dt |
dt |
" г |
wx |
|
3~T""d |
+ |
-■ |
|
|
|
||
|
dx |
|
Wy |
d yV- |
■Wz ~ ~ - = |
(1-61) |
||||||
|
1 |
dp |
|
|
2i% |
|
|
|
|
d2wv |
||
= / « |
- f |
¥ |
+ |
v |
dx2 |
+ |
1 |
^ |
+ |
dz2 |
|
|
dw* |
dw* |
|
• |
dw- |
г |
|
dw, |
I |
|
док |
|
|
dt |
dt |
|
~Tx |
+ |
«о* -Шdy - + |
®* dz |
|
|||||
|
|
dp |
|
|
ä2Wf |
I |
d2wz |
|
â2wz |
|
||
■ ■ fz — |
|
d y |
|
|
dx2 |
”1 |
dy2 |
|
d z 2 |
|
||
|
|
|
|
|
|
и в векторной |
форме ( [A = VV |
d2 |
â2 |
||
dx2 + |
dy2 |
||||
dw |
dw |
(wy)w = |
f -------vP-i-'vAw. |
||
dt |
dt |
||||
|
P |
|
d2 dz2
(1-62)
Умножая уравнение (1-59) векторно на радиус-вектор относительно произвольно выбранного начала, получа ем уравнение моментов им-пульсов
М= от(r2 Xw2—ГіX w4) |
(1-63) |
37
или |
(1-64) |
M = m(r2W2n—riWin), |
|
где Win и w2n — проекции относительной |
скорости w на |
радиусы гі и г2 соответственно. В этой форме уравнение движения называется уравнением Эйлера. Оно широко попользуется при расчетах лопастных машин (компрес соров п турбин).
Закон сохранения энергии. Для составления уравне
ния энергии одномерного течения (все параметры |
меня |
|
ются только .вдоль одной координаты — длины /, |
отсчи |
|
тываемой по оси струйки |
от какой-либо точки на ней) |
|
запишем закон сохранения энергии в виде |
|
|
(Е1+ En)t + SL + |
= (Еи + Ет ),+м . |
(1-65) |
где Е — суммы внутренней и кинетической энергии; 6Е — работа сил (массовых и поверхностных) за время Ы', SQ — тепло, подведенное за время St.
Изменение полной энергии равно сумме изменен кинетической m (w \—w\)j2 (где т = р,®, Ы S x= p2w2 Ы S 2)
и внутренней энергии тсѵ (Т2— Ті). Работа массовых сил (тяжести) равна mg(y2—yi); работа поверхностных сил равна работе сил дазлеішя p2S2w2 St—piSiWi St. Если, кроме того, обозначить через &L„nem работу, совершен ную жидкостью, например, в турбине и вынесенную на ружу, то, обозначая через SQ подведенное тепло, можно написать:
wö |
■m f |
ih |
Pl |
4" |
пг |
||||
|
V |
P2 |
Pl |
|
+ |
m g (y » -y i)- ■3L„ |
|
|
|
Сокращая на m и обозначая SQ/m = Sq и 8LBmum/rn= |
||||
= 6/внеш работу |
трения и работу турбины |
на |
единицу |
|
массы, получаем: |
|
|
|
|
|
= f - -*- + * & - * , ) + |
|||
|
^внеш 4 " ty- |
|
|
(1 -66) |
Для вывода дифференциального уравнения энергии обозначим сумму внутренней и кинетической энергии
38
CvT + zv2/2 |
единицы массы через |
е и вычислим |
сначала |
||
(предположив d/dt=/=0 ) |
|
|
|
|
|
(^ut+u |
в ш ) + (£„,«, - |
Еи) = Л Щdtр п . |
ы + |
||
V »p8S t» S < + iM n ? > 8/ + |
.... |
s p bSwbt |
|
||
|
де |
de |
|
|
|
|
-pbSbl Ы ^ dt |
■ w - W |
|
|
|
так как p 65 6/ = p 6У не меняется во времени, |
а р öS w |
||||
не меняется вдоль трубки тока. |
|
|
время 61, |
||
Далее вычислим работу |
массовых |
сил за |
т. е. 5LM.c = p//6S 6/ wht и работу поверхностных сил (так как dw/dr = 0 , то силы вязкости отсутствуют!)
SL = р SS wit — pbSw bt + d- ^ - U b S U + ...J =
- (l± p -blbS5t. dl
Так «ак движение одномерное, то тепло может быть подведено только через торцевые сечения. Поэтому мож но написать (имея в виду направление потока тепла):
SQ: |
|
дТді bS Ы |
дТді bS ы - |
|
г 5/+ - ) 858/ |
дТ |
|
|
ді - |
||
Подставляя |
все |
это в (1-65), |
после сокращения на |
б/ 65 6/ (считая А = |
const), получаем: |
d |
д (pw) . |
^ dt |
= Рfi® ~ dl I |
- |
д*Т |
(1-67) |
|
' |
dl2 |
||
|
Умножение уравнения движения (1-61) на w дает:
dw |
с |
dp |
?w ~ d f ~ ? wh ~ ~ w |
• |
Следовательно, изменение кинетической энергии про исходит за счет работы массовых (рw f i ) и поверхност ных (wdp/dl) силВычитая это произведение из (1-67),
39