Теплопередача учебное пособие
..pdfслой заторможенной жидкости, в пределах которой скорость изменяется от нуля на поверхности тела до скорости невозмущенного потока вдали от тела. Этот слой заторможенной жидкости получил название гидродинамического пограничного слоя. Под толщиной слоя δ подразумевается такое расстояние от поверхности тела, на котором скорость будет отличной от v0 на опреде-
ленную заранее заданную малую величину, например на 1 %.
Таким образом, поток жидкости, проходящий вблизи поверхности пластины, условно можно разделить на пограничный слой и на внешний поток. В пограничном слое силы инерции соизмеримы с силами вязкого трения, а во внешнем потоке преобладают силы инерции.
Тепловой пограничный слой.
По аналогии введено понятие теплового пограничного слоя. Тепловой пограничный слой – это слой жидкости у стенки, в пределах которого температура изменяется от значения, равного температуре стенки, до значения, равного температуре жидкости вдали от тела (рис. 2.4).
Таким образом, все изменения температуры жидкости сосредоточиваются в сравнительнотонком слое, непосредственно прилегающем кповерхности тела.
101
2.4. ПОДОБИЕ И МОДЕЛИРОВАНИЕ ПРОЦЕССОВ КОНВЕКТИВНОГО ТЕПЛООБМЕНА
2.4.1. Приведение математической формулировки краевой задачи к записи в безразмерном виде
Рассмотрим продольное обтекание плоской поверхности твердого тела размером l0 безграничным потоком жидкости, температура и скорость которого вдали от тела постоянны и равны соответственно t0 , v0 (рис. 2.5) [1]. Температура поверхности тела равна tc . Будем считать, что tc > t0 , физиче-
ские параметры жидкости постоянны (учтем только подъемную силу, возникающую в результате зависимости плотности от температуры).
Диссипативный источник тепла не учитывается. Процесс стационарный.
Рис. 2.5. К постановке краевой задачи конвективного теплообмена
102
Расположим оси координат так, как показано на рис. 2.5. Будем считать, что вектор ускорения свободного падения совпадает с направлением оси 0x , в связи с этим компоненты
gx = g , g y = gz = 0 .
Размер тела по направлению оси 0z намного больше l0 . Введем переменную ϑ = t − t0 , где t – температура жидко-
сти; t0 = const – температура невозмущенного потока жидкости. При этом dt = dϑ .
Впотоке жидкости учтем подъемную силу ρgβϑ , считая
еесоизмеримой с силой вязкого трения.
Вданной постановке поля скоростей и температур можно описать дифференциальными уравнениями в приближении пограничного слоя:
|
|
|
∂ ϑ |
|
|
∂ϑ |
|
|
∂ ϑ 2 |
|
|
|
|
||||||
|
vx |
|
|
+vy ∂ y |
|
|
= a∂ |
|
|
; |
|
|
|||||||
∂ x |
y2 |
|
|
||||||||||||||||
|
∂ vx |
∂ vx |
|
|
∂ 2vx |
+ gβϑ ; |
|
|
|||||||||||
vx |
|
+vy |
|
|
=v∂ |
|
|
|
|
||||||||||
∂ x |
∂ y |
|
y2 |
|
|
||||||||||||||
|
|
|
|
∂ vx |
+ |
∂ vy |
|
= 0. |
|
|
|
|
|
||||||
|
|
|
|
∂ |
x |
∂ y |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Граничные условия: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
1) вдали от тела (y = ∞ |
, 0≤ |
|
≤x |
|
|
l0 ,− ∞ ≤ |
≤ +∞z |
) |
|
||||||||||
ϑ =ϑ 0 = t0 − t0 = 0; vх |
= v0 ; vу |
= 0; |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|||||||||||||||
2) на поверхности тела (y = 0, 0 ≤ x≤ |
|
l0 , − ∞ ≤ |
≤ +∞z |
|
|||||||||||||||
|
) |
||||||||||||||||||
ϑ =ϑ c = tc − t0 = const; |
|
vх = vy |
= 0. |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
В уравнениях и условиях однозначности можно различить три вида величин:
103
♦независимые переменные – координаты x, y;
♦зависимые переменные – ϑ , vx , vy ; которые однознач-
но определяются значениями независимых переменных при заданных краевых условиях;
♦ постоянные величины– это v0 , ϑ 0 , l0 , ϑ c , ν, a , g , β
и др.; они определяются краевыми условиями и для данной задачи являютсяпостоянными, независящимиотдругихпеременных.
Таким образом, искомые зависимые переменные ϑ , vx , vy
являются функцией независимых переменных и постоянных величин, входящих в краевые условия.
Приведем дифференциальные уравнения тепло- и массопереноса потока жидкости к записи в безразмерном виде. Для этого выберем масштабы приведения, в качестве которых выберем величины, входящие в краевые условия. Для пространственных координатвозьмемдлину l0 , для скорости – v0 , длятемпературы – ϑ c .
Введем обозначения безразмерных величин:
|
x |
|
|
y |
|
vx |
|
|||||
X = |
|
; Y = |
|
|
; Vx |
= |
|
|
; |
|||
|
|
|
v0 |
|||||||||
|
l0 |
|
|
l0 |
|
|
||||||
|
vy |
|
|
|
|
ϑ |
|
|
|
|
||
Vy = |
|
Θ = |
|
|
|
|
|
|
||||
|
|
; |
|
|
, |
|
|
|
|
|||
v0 |
|
|
|
|
||||||||
|
|
|
ϑ c |
|
|
|
|
тогда
x = l X ; |
y = l Y ; v |
x |
=v V ; |
0 |
0 |
0 x |
|
v y =v0Vy ; ϑ =ϑ cΘ . |
|
|
|
|
|
||
|
|
|
|
(а)
(б)
Подставим выражения (б) в дифференциальные уравнения тепло- и массопереноса:
(v0Vx ) |
|
( |
c |
) |
+ (v0Vy ) |
( |
c ) |
= a |
∂ |
|
|
|
( |
c ) |
|
|
0 |
; |
|
∂ |
ϑ |
Θ |
∂ |
ϑΘ |
|
|
|
∂ |
|
Θϑ |
|
|
|
|
l 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ (l0 X ) |
|
∂ (l0Y ) |
∂ |
|
|
|
(l0Y ) |
|
aϑ c |
|
|||||||
|
|
(l0Y ) ∂ |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
104
(v |
0Vx ) |
∂ (v0Vx ) |
|
+ (v0Vy ) |
∂ (v0Vy ) |
|
= |
∂ 2 (v0Vx ) |
|
+ gβϑ cΘ |
|
|
l02 |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ν |
|
|
|
|
|
|
|
|
; |
||||||||||||||||||||||||||||
∂ (l0 X ) |
∂ (l0Y ) |
|
(∂ (l0Y ))2 |
|
νv0 |
||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
∂ (v0Vx ) |
+ |
|
∂ |
(v0Vy ) |
= 0 . |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ (l0 X ) |
|
|
|
∂ (l0Y ) |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
После преобразования получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
v0l0 |
|
|
|
∂Θ |
|
|
|
|
∂Θ |
|
|
|
|
∂ |
|
Θ |
2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Vx |
|
|
|
|
|
+Vy |
|
|
|
|
|
|
= |
|
|
|
|
; |
|
|
|
|
|
(2.31) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
∂ X |
∂ |
|
|
|
|
Y |
2 |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y |
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
v0l0 |
|
|
|
∂ Vx |
|
|
|
|
|
|
|
∂ Vy |
|
∂ 2Vx |
|
|
gβϑ cl03 |
|
ν |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
Vx |
|
|
|
+Vy |
|
|
|
|
|
= |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
Θ ; |
|
|
(2.32) |
|||||||||||||
|
|
|
|
ν |
|
∂ X |
|
∂ Y |
|
|
|
|
Y |
2 |
|
|
ν |
2 |
|
|
|
v0l0 |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ Vx |
|
+ |
|
∂ Vy |
= 0 |
|
, |
|
|
|
|
|
|
|
|
|
|
|
(2.33) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ X |
∂ |
Y |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
здесь |
|
gβϑ cl02 |
|
Θ = |
|
gβϑ cl03 |
|
|
νΘ |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
νv0 |
|
|
|
ν2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
v0l0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Постоянные величины, содержащиеся в дифференциальных уравнениях (2.31)–(2.33), сгруппированы в безразмерные комплексы, число которых меньше числа размерных величин.
Граничные условия в безразмерном виде:
1) вдали от тела (Y = ∞ |
, 0≤ X≤ 1,− ∞ ≤ |
≤ +∞z |
) |
|
||
Θ = Θ =0 |
0; V=x |
1; V=y |
0; |
|
|
|
|
|
|
||||
2) на поверхности тела |
(Y = 0, 0 ≤ X≤ |
1,∞ ≤ |
≤ +∞z |
|
||
) |
||||||
Θ = Θ =с |
1; V=x |
V=у 0. |
|
|
|
|
|
|
|
|
Несмотря на то, что величины v0 , ϑ 0 , l0 и др., входящие в размерные граничные условия, могут иметь различные числовые значения, каждая из безразмерных величин Θ 0 , Θ c и др. имеет в рассматриваемом случае вполне конкретное значение.
105
Коэффициент теплоотдачи по известному распределению температуры определяется из уравнения теплоотдачи (2.30), которое для рассматриваемой задачи примет вид
|
|
λ |
|
∂ ϑ |
|
α = − |
|
|
|
|
. |
tc |
− t0 |
|
|||
|
|
∂ y y =0 |
Уравнение теплоотдачи в безразмерном виде:
α l0 |
∂Θ |
|
|
|
|
= − |
|
|
. |
λ |
∂ |
|
||
|
Y Y =0 |
2.4.2. Числа подобия и уравнения подобия
Помимо безразмерных зависимых переменных Θ , Vx , Vy
и безразмерных координат в дифференциальные уравнения входят также безразмерные комплексы, называемые числами подобия [1]:
|
αl0 |
, |
v0l0 |
, |
v0l0 |
, |
gβϑ cl03 |
. |
||
|
λ |
|
ν |
a |
ν2 |
|||||
Числам подобия присвоены имена ученых, внесших значи- |
||||||||||
тельный вклад в развитие гидродинамики или теплопередачи. |
||||||||||
Безразмерный комплекс |
|
|
|
|
|
|||||
|
|
|
Nu ≡ |
|
α l0 |
(2.34) |
||||
|
|
|
|
λ |
||||||
|
|
|
|
|
|
|
|
|
называется числом Нуссельта или безразмерным коэффициентом теплоотдачи. Число Нуссельта характеризует теплообмен между поверхностью твердого тела и подвижной окружающей средой. Одной из основных задач конвективного теплообмена является нахождение числа Nu, и затем коэффициента теплоотдачи α . Число Bi в задачах теплопроводности рассчитывается по формуле, схожей с формулой (2.34) для числа Nu. При этом
106
между этими безразмерными комплексами имеются принципиальные отличия. В формуле для числа Bi присутствует коэффициент теплопроводности λ твердого тела, и коэффициент теплоотдачи α является наперед заданной величиной, а в формуле для числа Nu – коэффициент теплопроводности λ жидкости, и коэффициент теплоотдачи α – искомая величина.
Второй безразмерный комплекс
Rе ≡ |
v0l0 |
(2.35) |
|
ν |
|||
|
|
называется числом Рейнольдса, которое характеризует отношение сил инерции к силам вязкого трения. Число Рейнольдса определяет характер как изотермического, так и неизотермического процессов течения подвижных сред.
Безразмерный комплекс
Ре ≡ |
v0l0 |
(2.36) |
|
a |
|||
|
|
называется числом Пекле. После его преобразования можно получить
|
v0l0 |
= |
ρсv0ϑ |
, |
||
|
|
|||||
|
a |
|
λ |
ϑ |
||
|
|
|
|
|
||
|
|
|
|
l0 |
здесь выражение в числителе отвечает за конвективный перенос тепла, а в знаменателе – за перенос тепла за счет процессов теплопроводности.
Четвертый безразмерный комплекс
Gr ≡ |
gβϑ cl03 |
(2.37) |
|
ν2 |
|||
|
|
называется числом Грасгофа. Оно характеризует подъемную силу, возникающую в жидкости вследствие разности плотностей нагретых и холодных областей подвижной среды.
107
После подстановки чисел подобия в безразмерные дифференциальные уравнения, получим:
|
|
|
∂Θ |
|
|
∂Θ |
|
|
|
|
|
|
∂ Θ |
2 |
|
|
|
|
|
|
||||
|
Pe Vx |
|
+Vy |
|
|
|
|
|
|
|
= |
|
|
|
; |
|
|
|
(2.38) |
|||||
|
∂ X |
∂ |
|
|
|
|
Y |
2 |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
Y |
∂ |
|
|
|
|
|
|
|||||||||
|
|
∂ Vx |
+Vy |
∂ Vx |
|
= |
|
|
Gr |
Θ + |
∂ |
2Vx |
|
|||||||||||
Rе Vx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
(2.39) |
|||||
∂ X |
∂ |
Y |
|
|
|
|
Rе |
∂ |
Y |
2 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
∂Θ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
Nu = − |
|
|
|
|
|
. |
|
|
|
|
|
(2.40) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
∂ |
Y Y =0 |
|
|
|
|
|
|
|
Уравнение (2.33) остается неизменным.
Безразмерные величины Θ , Vx , Vy , X , Y , Nu, Re, Pe, Gr можно разделить на три группы:
1)независимые переменные – это безразмерные координаты X и Y ;
2)зависимые переменные – это Θ , Vx , Vy , Nu;
3)постоянные величины – это Re , Pe , Gr ; они заданы краевыми условиями.
В результате имеем:
Nu = f1 ( X ,Y , Pe, Rе,Gr ) ; |
(2.41) |
|
Θ = |
f2 ( X ,Y , Pe, Rе,Gr ) ; |
(2.42) |
Vx |
= f3 ( X ,Y , Pe, Rе,Gr ) ; |
(2.43) |
Vy |
= f4 ( X ,Y , Pe, Rе, Gr ). |
(2.44) |
Уравнения(2.41)–(2.44) носят названияуравнений подобия.
Если учесть градиент давления ρ1 ∂∂ Рх в уравнении дви-
жения, при переходе к записи в безразмерном виде можно по-
|
l |
2 |
|
∂ P |
|
|
∂ |
|
|
P |
v |
l |
|
= ∂ |
|
(Eu Rе) . |
|
лучить |
0 |
|
|
= |
|
|
|
0 0 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
∂ |
|
2 |
|
|
ν |
|||||||||
|
νρv0 ∂ x |
X |
ρv0 |
|
|
|
∂ |
X |
|||||||||
108 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Величина
Eu = |
P |
(2.45) |
|
ρv02 |
|||
|
|
называется числом Эйлера. Оно характеризует отношение сил давления к силам инерции.
Если в числе Pe числитель и знаменатель умножить на коэффициент кинематической вязкости, то можно получить следующее соотношение:
Ре = |
v0l0 |
|
ν |
|
= Rе Рr . |
(2.46) |
||||||
ν |
|
|
||||||||||
|
|
a |
|
|
|
|||||||
Безразмернаявеличина Рr ≡ |
ν |
называетсячисломПрандтля: |
||||||||||
a |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|||
Рr = |
ν |
= |
µ с |
. |
(2.47) |
|||||||
|
|
|||||||||||
|
|
a |
λ |
|
Число Рr – это теплофизическая характеристика теплоносителя.
Для неметаллических капельных жидкостей Рr > 1 ; для неметаллических капельных жидкостей при больших температурах и газов Рr ≈ 1; дляжидких металлов Рr << 1 .
Числа подобия, составленные из наперед заданных параметров (постоянных) математического описания процесса, называют критериями подобия.
2.4.3. Моделирование процессов конвективного теплообмена и построение эмпирических критериальных уравнений
Изучение моделирования процесса в объекте заменяется исследованием этого процесса на модели [1].
Исследование процесса на модели осуществляется так, чтобы результаты его изучения можно было перенести на объект.
109
Моделирование включает в себя две задачи:
1)в модели необходимо осуществить процесс, подобный процессу, происходящему в объекте;
2)выполнить на модели все требуемые измерения и наблюдения.
Условия подобия физических процессов:
1. Подобные процессы должны быть качественно одинаковыми, т.е. они должны иметь одинаковую физическую природу и описываться одинаковыми по форме записи дифференциальными уравнениями.
2. Условия однозначности подобных процессов должны быть одинаковыми во всем, кроме числовых значений размерных постоянных, содержащихся в этих условиях.
3. Одноименные определяющие безразмерные переменные подобныхпроцессовдолжныиметьодинаковыечисловыезначения.
Сформулированные условия являются определением подобия физических процессов.
Подобные процессы можно рассматривать как один и тот же процесс, но взятый в различном масштабе, причем масштабы размерных величин могут быть неодинаковыми.
Знание математического описания является необходимой предпосылкой теории подобия.
Получение эмпирических формул. Прежде чем обрабаты-
вать опытные данные в числах подобия, нужно установить, от каких чисел зависит определенное значение. Это можно получить из системы дифференциальных уравнений, описывающих процесс [1].
Пусть было получено, что Nu = f (Re) .
По данным измерения подсчитывается значение Re и соответствующие значения Nu . Зависимость между числами подобия обычно представляется в виде степенных функций, например Nu = cRеn , где c, n являются постоянными безразмерными числами (рис. 2.6).
110