
книги из ГПНТБ / Борисенко А.И. Аэродинамика и теплопередача в электрических машинах
.pdfполучаем уравнение энергии в виде, характеризующем процессы внутри рассматриваемого элемента,
|
dT |
dw |
д*Т |
(1-68) |
|
РсV dt |
Р Ж |
ді2 • |
|
|
d9 I Ddw — |
|||
Так как в |
силу уравнения неразрывности |
|||
|
dw |
V dS |
p dv , |
d t ' v d l ~ |
„ |
, I |
|||
= 0' T0 член |
|
~ТШ’==Ж 7Г (где ^=1/р-удель- |
кый объем) характеризует работу деформации элемента (здесь это работа только сил давления, так как силы вязкости не проявляются). Таким образом, тепло, подве денное к газу (Я д2Т/дР), идет на увеличение внутренней энергии рcydT/dt и на работу расширения р dw/dl объ ема элемента.
Если нет теплообмена (адиабатическое движение), нет сил трения и жидкость не совершает механической работы, то интеграл уравнения энергии для струйки (1-66) носит название уравнения Бернулли и имеет вид:
сѵТ + |
+ —- + 0 2 = const. |
(1-69) |
Для несжимаемой жидкости, когда изменения внут ренней энергии нет и термодинамические процессы, про исходящие внутри жидкости, не сказываются на ее дви жении, уравнение Бернулли часто пишут в виде
1 Г + ^ + у = С0П8+'- |
(1'70) |
В этой записи все члены имеют размерность длины. Смысл этого уравнения: для всех точек идеальной жид кой струйки сумма высот — скоростной (wzl2g), пьезо метрической (p/pg) и геодезической (у), отмериваемых от одной горизонтальной плоскости — есть величина по стоянная.
Давление р, входящее в уравнение (1-69), в гидрав лике называют гидростатическим, а величину р + рш2/2 —
гидродинамическим давлением.
Во многих задачах (в частности, при рассмотрении вентиляции электрических машин) влияние изменения величины у ничтожно; замечая еще, что p = pRT и R =
— сѵ—Сѵ, можно написать:
срТ + - ^ - = c o n st= < ? p7'0. |
(1-71) |
40
Температура Т0 = Т -\-~^г носит название темпера
туры заторможенного потока.
Уравнение Бернулли служит основой почти всех ги дродинамических расчетов. Сопротивление (потери) учи тываются введением дополнительного члена
2 |
2 |
|
|
0W\ |
0Wo |
|
|
Р . + - 2- = Р г + - 2^ + |
ЛРі,2- |
0 "7 2) |
|
Подчеркнем, что уравнение Бернулли справедливо для |
|||
одной и той же струйки |
(линии тока); константа в нем |
||
имеет различное значение для каждой линии тока. |
|||
В общем случае трехмерного |
движения |
уравнение |
|
энергии в прямоугольной |
системе |
координат, |
связанной |
с частицей, имеет вид: |
|
|
|
(1-73)
В векторной форме
pcH 4 r s p * v ( 4 r + wv:r ) + />(Ѵ«0 = ЯѵіТ + pD, (1-74)
где
(xD = p,
|
d w x |
|
dwn |
dwx |
d w t |
(div w)a -j- |
dy |
+ |
d x УЧ |
ö z + |
dx + |
|
dw |
|
dwv |
|
(1-75) |
+( dy |
+ |
dz } |
|
||
|
|
— так называемая функция рассеивания (диссипации); это работа сил трения, переходящая в тепло, которое рас сеивается в массе жидкости.
Физический смысл этого уравнения: тепло, подведен ное к жидкости ХѴ2Т и выделяющееся вследствие работы сил вязкости pD, идет на увеличение внутренней энергии рс dTjdt и на работу расширения p(Vw).
41
Динамический пограничный слой. Отсутствие сколь жения на твердых границах («прилипание») приводит к большим градиентам скорости вблизи этих границ и зна чительных касательных напряжений т = ц dw/dn внутри тонкого пограничного слоя. В этом слое влияние сил вязкости соизмеримо с влиянием сил инерции, тогда как вне его силы вязкости пренебрежимо малы по сравнению с другими силами (инерции, давления, массовыми) и те чение близко к потенциальному.
Расстояние от стенки до места, где скорость течения «почти равна» скорости внешнего течения, называется толщиной пограничного слоя б. Эта величина является несколько неопределенной.
Рассмотрим |
двухмерный установившийся поток |
(ш2= 0, d/dz=0; |
dfdt = 0 ) вязкой несжимаемой жидкости |
вдоль плоской границы твердого тела (хотя можно пока зать, что выводы справедливы и для случая криволиней ной стенки). Уравнения неразрывности и движения (при отсутствии массовых сил) при этом запишутся в виде
|
|
|
dwx |
|
I |
dwu |
|
|
|
|
|
|
|
|
|
|
дх |
|
* |
ду |
|
|
|
|
|
|
|
Wx |
dwx |
ШУ |
dwx |
__ |
|
1 |
др |
1 |
г ( |
d2Wx |
1 |
dswx |
(1-76) |
|
дх 1 |
ду |
|
|
р |
дх |
‘ |
|
у дх2 |
Т |
ду2 |
||
wx |
dwy |
dwy |
~~ |
|
1 |
др |
1 |
V ( д*щ |
1 |
d2Wy |
|
||
|
дх |
|
ду |
|
Р |
ду |
‘ |
|
1дх2 П |
|
ду2 |
|
|
|
Пусть |
ось |
X направлена |
вдоль |
границы, а |
ось у — |
по нормали к ней. Оценим величины отдельных членов, входящих в это уравнение. Для этого введем характер ный масштаб I длин вдоль оси х и в качестве масштаба вдоль у возьмем толщину пограничного слоя б, положив
x=lU У = цЬ,
где £ и г| — безразмерные величины.
Умножая уравнение неразрывности на dy и интегри руя его в пределах от 0 до б, получаем вертикальную составляющую скорости на границе пограничного слоя
S)=- j^ |
i f^f'-Л|. |
0 |
О |
42
Если скорость wx изменяется от значения нуль у стен ки до некоторого значения wox на внешней границе слоя, которое сравнимо с некоторой характерной для данного потока скоростью w (например, со скоростью невозму щенного потока), то производная dwx/dl мало зависит от толщины слоя и будет порядка w. Следовательно, попе речная скорость wy имеет внутри слоя величину порядка
W Ö / L
Тогда в левой части второго уравнения (1-76) полу чим: порядок величины первого члена wxdwxjdx, равный w2!t\ порядок величины второго члена wydwx/dy, также равный w2jl.
В правой части |
этого |
уравнения |
получим, |
принимая |
|||||||||
в качестве |
масштаба давления |
скоростной напор |
р&у2/2 , |
||||||||||
порядок |
величины первого члена |
~ ~ —j57-. |
равный |
w2[l; |
|||||||||
|
|
|
|
|
|
|
|
Р |
дх |
|
|
д |
dw- |
порядок |
величины |
второго |
члена ѵ |
d2w, |
= |
12 |
|||||||
|
|
дЬ |
|||||||||||
равный vw/F; порядок величины третьего члена |
d2w~ |
||||||||||||
ѵ - ~ = |
|||||||||||||
д2 |
д |
dwx , равный ѵда/82. |
|
|
|
|
|
|
|
||||
ârj |
|
dr; |
|
|
|
|
|
|
|
|
|
|
|
Рассматривая |
порядок |
величин всех членов |
уравне |
||||||||||
ния по отношению |
к порядку |
величины ш2//, получаем, |
|||||||||||
что член |
V• |
d2wx |
s |
|
|
|
ѵш> |
W2 |
|
V |
|
1 |
|
- дх2 |
будет порядка |
|
|
Т |
|
wl |
|
Re |
|||||
а порядок величины члена ѵ- |
d2w~ |
будет |
|
|
W2 |
||||||||
ду2 |
|
Т |
|||||||||||
___lj_ _ 1_ |
|
|
|
|
|
|
|
|
|
|
|
||
— ‘ ö 2 |
Re |
|
|
|
|
|
|
|
|
|
|
|
|
Сравнивая между собой члены, характеризующие вли-
d2wx
яние вязкости, можно сказать, что ѵ- ду2 будет в /2/82
d2w
раз больше ѵ—-^а—. Опыт показывает, что при больших
Re отношение |
Ь/1 мало. Поэтому при больших числах Ре |
||
,, |
|
d2wx |
• • |
можно пренебречь членом ѵ |
|
Предполагая далее, что все члены уравнения должны быть одного порядка, приходим к выводу, что толщина слоя б должна иметь порядок величины
Аналогичной оценкой порядка величины членов третье го уравнения получим:
dwy |
w2S |
wtl |
dwy |
w2S |
dp_ |
~T |
дх |
~~г~’ |
~~ду~ |
~ ~ ~ ; |
P dy |
||
|
d2wy |
|
vwS |
d2Wy |
\w |
|
|
дх2 |
|
|
öy2 |
|
|
Так как все члены второго уравнения очень малы по сравнению с членом др/ду, то можно принять при боль ших Re
дрІду = О, |
(1-78) |
т. е. давление внутри пограничного слоя одинаково вдоль нормали к контуру тела и равно давлению на внешней границе слоя на рассматриваемой нормали.
Таким образом, приближенные уравнения неразрыв ности и движения в пограничном слое могут быть запи-
с а н ы в ф о р м е |
dwx |
|
|
|
|
|
I |
ÖWy — п- |
|
)1 |
|
dwx I1 |
дх |
1 |
ду |
|
1 |
dwx |
|
1 dp 1 |
d2wx |
[ ( 1 - 7 9 ) |
|
а* Н |
äy - |
|
P дх |
д у • |
|
|
1 |
||||
|
dp |
- 0 . |
|
1 |
|
|
ду |
|
) |
Результаты численного решения этих уравнений для плоской пластинки в предположении, что профили ско ростей в пограничном слое подобны во всех сечениях x = const, т. е. wx/w0=f(y/Ö) , и сопоставление их с опыт ными данными приведены на
рис. 1-15.
Хотя найденное реше ние асимптотическое, но, как видно из графика, уже при
т V
т- е- ПРИу =
= 4 х /У Я ех, скорость'практически равна скорости невоз-
Рис. 1-15. Распределение скоро стей в пограничном слое на пла стинке при ламинарном режиме течения (сплошная линия — теория).
44
мущенного потока іѵ0. Эту величину можно принять за толщину пограничного слоя для сечения х, в котором сосредоточено все влияние вязкости при ламинарном об текании пластинки.
Соответственно двум режимам течения вязкой жидко сти различают ламинарный и турбулентный пограничные слои. В ламинарном пограничном слое движение жидко сти происходит слоями без пульсаций скорости; в тур булентном имеются хаотические, подчиненные только статистической закономерности пульсации скорости и связанное с ними молярное перемешивание жидкости.
Напряжение трения на стенке при этом течении
dw — ? *=о = Г* ----PW’xW'y,
где второй член характеризует молекулярную вязкость, связанную с молярным объемом импульсов.
В непосредственной близости к стенке основную роль играет первое слагаемое, а вдали от нее — второе. Для определения расстояния уо, на котором существенно влияние молекулярной вязкости р, заметим, что в этой
области порядок величины скорости |
равен ш * = ]/ то/р |
|
и соответствующее число Re= y0w jv |
должно быть |
по |
рядка единицы; отсюда y0=v/w *. |
|
|
Таким образом (рис. 1-16), в непосредственной |
бли |
зости от стенки |
к ней прилегает тонкая прослойка жидко |
|
сти, в которой |
усредненная скорость мала и |
меняется |
по линейному закону от нуля до значения ю* = |
|А 0/Р- |
Эту прослойку, как было сказано, называют вязким подслоем. При турбулентном течении продольные состав ляющие пульсационных скоростей значительно больше поперечных и течение в вязком подслое также турбулент-
Рис. 1-16. Развитие пограничного слоя на пластине.
45
ное, но с линейным законом усредненных скоростей. Вне области вязкого подслоя (резкой границы между вязким подслоем и остальным потоком жидкости провести, ко нечно, нельзя) основную роль играют турбулентные на пряжения pw'xw'y.
Из величин, определяющих движение жидкости с плот ностью р, создающее напряжение то, можно составить единственную комбинацию, имеющую размерность dwjdy,
а именно У Xolply, и поэтому должно иметь место равен ство
d w ___L |
i / " _I°_, |
|
d y ~ x y |
V |
p |
где x — так называемая универсальная постоянная, кото рая определяется только опытом (по результатам испы таний в трубах).
Интегрирование дает логарифмический закон распре
деления скоростей: |
|
да = ^ -1 п ^ * -+ С , |
(1-80) |
где г/ш*/ѵ — так называемое безразмерное расстояние от стенки.
Из опыта для труб: х = 1/2,4 и С=5,84.
Для плоской пластинки толщина турбулентного по граничного слоя
8 = 0,37//R£T.
Переход от ламинарной формы движения к турбу лентной зависит от числа Rex=T2j0x/v. С увеличением ско рости при всех прочих одинаковых условиях область пе рехода приближается к переднему краю хорошо обте каемого тела. Структура пограничного слоя конечной толщины дана на рис. 1-16.
Интегральные соотношения. Помимо нелинейности приближенных уравнений пограничного слоя трудности решения их заключаются в необходимости удовлетворить граничным условиям на стенке, где скорость жидкости равна нулю, и на бесконечности, где исчезает влияние вязкости. Это решение предполагает, строго говоря, без граничную толщину пограничного слоя, в которой ско рость асимптотически приближается к скорости невозму щенного потока. Однако теоретические оценки, приведен-
46
Рис. 1-17. Интегральные соотношения получаются при менением законов сохранения к элементу пограничного
ABCD.
ные выше и подтвержденные опытом, показывают, что при ламинарном движении, например, уже на безразмер
ном расстоянии порядка 1 /|/R e разница между скоро стью потенциального обтекания и скоростью в погранич ном слое практически исчезает. Эта идея пограничного слоя оказалась весьма плодотворной. Она привела к за мене дифференциальных уравнений для каждой точки потока в пограничном слое интегральными соотношения ми, которые получаются в результате применения зако нов сохранения к элементу длины слоя. Метод интеграль ных соотношений по своему замыслу, оперирующий с не которым образом усредненными параметрами, оказывает ся сравнительно малочувствительным к точности приня того в расчет закона распределения продольной скорости по нормали. Расчеты показывают, что во многих случаях достаточно, чтобы принятое распределение было только более или менее правдоподобным.
Изменение количества движения в объеме ABCD (рис. 1-17) должно быть равно секундному импульсу сил,
складывающихся |
из силы трения r0dx |
на |
грани |
АС и |
результирующей силы давления, равной |
р |
dS |
, |
|
ах, как |
||||
это видно из рис. |
1-17. |
|
|
|
Это дает: |
|
|
|
|
гг
d Г |
, |
d Г 2 , |
dd |
|
Р З Г |
J 9 * ° J y = - % - р -3F ’ |
|
о |
|
9 |
' |
47
или, замечая, что при р = const получаем p-f-рДОд /2 =const
и, следовательно, dpfdx — — pwa dwjdx,
5 |
|
5 |
|
|
- ^ - J p wx{w0 — wx) d y + - ^ ~ [ р(ш0 — ts»,)dj/ = |
v |
(1-81) |
||
о |
|
b |
|
|
В дальнейших вычислениях для несжимаемой |
жидко |
|||
сти удобно первый |
интеграл записать в виде |
хю~0 §**, а |
||
второй — в виде ш05*, |
где |
|
|
|
6* |
= |
|
|
|
|
5, |
оо |
|
(1-82) |
|
|
|
dy.
О
Знак 6, оо означает, что с некоторого значения б прак тически безразлично, до какого расстояния от стенки про изводят интегрирование.
Величины б* и б** имеют размерность длины и носят название толщины вытеснения и толщины потери импуль са соответственно.
При этом вместо (1-82) получаем: |
|
|
||
dè* |
w\ |
(23**-1-8*) = |
to |
(1-83) |
dx |
w0 |
Po^o |
||
|
|
|
|
Это основное интегральное соотношение теории по граничного слоя было выведено Карманом и называется уравнением импульсов.
Уравнение энергии Л. С. Лейбензон получил анало гичным путем, применяя закон сохранения энергии к эле менту ламинарного пограничного слоя. В ней фигуриру ет толщина потери энергии
5, оо |
|
|
wx \ |
(1-84) |
|
l - - f |
К |
|
wo |
) |
|
О |
|
|
Условные толщины пограничного слоя. Толщины вы теснения б*, потерь импульса б**, энергии б*** значи-
48
тельно точнее характеризуют область влияния вязкости, чем расплывчатое понятие толщины пограничного слоя как расстояния от стенки до места, где скорость стано вится близкой к скорости в ядре потока.
Толщина вытеснения — расстояние, на которое нужно отодвинуть поток газа от стенок тела, чтобы не изменил ся массовый расход его при обтекании тела потенциаль ным потоком несжимаемой жидкости. Толщина потери импульса — дополнительное к толщине вытеснения рас стояние, на которое должен быть отодвинут от стенки те ла поток, чтобы не изменилось количество движения при обтекании тела потенциальным потоком несжимаемого газа.
Отметим, что распределение давления по телу хорошо совпадает с опытом, если расчет вести по теории потен циального обтекания тела невязким газом для тела, на контур которого наращена толщина вытеснения.
При наиболее грубом приближении — линейном рас пределении скорости (wx= w0yl6) в пограничном слое не сжимаемой жидкости имеем:
6* = 6/2; б** = 6/6; 6*** = 6/4.
Опыт показывает, что при движении жидкости с поло жительным градиентом давления в некоторой точке, точ нее говоря области, у поверхности тела возникает воз вратное движение жидкости. Физически это можно объ яснить тем, что при движении жидкости в направлении возрастающего давления частицы жидкости, лежащие ближе к поверхности тела и имеющие под действием тре ния наименьшую скорость, затормозятся и затем под влиянием давления, направленного против основного те чения, начнут двигаться обратно, создавая таким обра зом отрыв пограничного слоя. Из графика распределения
скоростей (рис. |
1-18) в пограничном слое |
следует, что |
в месте отрыва |
производная от скорости |
по нормали |
к стенке у самой стенки обращается в нуль. За местом отрыва вследствие возвратного движения образуется об ласть возникновения вихрей конечных размеров, которые проникают в область внешнего течения и изменяют фор му движения. __
Так как для ламинарного слоя S//~ 1 / у Re; W yj W x ^ 1 / У Re
и при изменении Re изменяются только поперечные мас штабы длины и скорости, а движение сохраняется подоб-
4—233 |
49 |