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

1

1. Математические модели жидкости

Сплошная среда, заполняющая некоторую пространственную об-

ласть и рассматриваемая как изолированная термодинамическая система, в соответствии с нулевым законом термодинамики имеет хотя бы одно естественное состояние. Характерным свойством вещества в жидком и газообразном агрегатных состояниях является то, что оно имеет несчетное множество естественных состояний. В качестве постулата считают, что все состояния, для которых плотность ρ неподвижной среды при постоянной температуре совпадает с начальным значением ρ0, являются естественными. Поэтому в качестве одного из аргументов активных переменных принимают якобиан J = ρ0/ρ, опре-

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

1.1. Жидкость как сплошная среда скоростного типа

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

σ и вектора q плотности теплового потока — реактивные перемен-

b

 

 

 

 

 

 

 

ρ,

 

среды

 

-

ные:

якобиан J = ρ0/ρ, равный отношению плотности ρ0

в

естественном состоянии к ее текущему значению

 

тензор скоро

 

стей

Vb ,

абсолютную температуру T и ее градиент ϑ, т. е.

 

 

 

 

 

 

 

 

A = A(J ,Vkl,T,ϑk), h = h(J ,Vkl,T,ϑk),

 

 

 

 

 

 

 

 

 

 

σij = σij(J ,Vkl,T,ϑk),

qi = qi(J ,Vkl,T,ϑk),

 

(1)

 

 

 

i, j, k, l = 1, 2, 3,

 

 

 

 

 

 

; ϑk

 

 

где Vkl и σij — компоненты тензоров V и σ

 

 

 

 

и

системыi

координат Ox1x2x3.

b

b

соответственно

 

 

 

q —

проекции векторов соответственно ϑ и q на оси прямоугольной

Подставив первое и третье соотношения (1) в уравнение ρ

du

 

=

dt

 

 

 

∂qi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= σjiVij

 

+ qV закона сохранения энергии (уравнение

переноса

∂xi

энергии), где u — массовая плотность внутренней энергии, с учетом равенства u = A + T h запишем

 

∂A dJ

 

∂A dVij

 

∂A dϑi

 

 

 

 

 

ρ

 

 

 

− σijVij + ρ

 

 

 

 

+ ρ

 

 

 

 

 

+

 

 

 

 

∂J dt

∂Vij

dt

∂ϑi dt

 

 

 

 

 

 

 

 

 

 

∂A

 

 

dT

 

dh

 

∂q

 

 

 

 

 

 

+ ρ

 

+ h

 

+ ρT

 

+

i

− qV = 0, (2)

 

 

 

 

 

∂T

dt

dt

∂xi

2

где t — время; qV

— объемная плотность мощности внутренних

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

 

 

 

 

 

 

 

 

 

 

 

 

dJ

=

∂J

 

d

 

∂xi

 

= J

∂ak

 

d

 

∂xi

 

= J

∂vk

= J Vkk,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

∂(∂xi/∂ak)

dt

∂ak

∂xi

 

dt

∂ak

∂xk

где ak — материальные координаты, определяющие положение частицы жидкости в начальной конфигурации, а vk — проекции вектора скорости на оси Oxk системы пространственных координат, то

ρ

∂A dJ

= ρ

∂A

J Vijδij = ρ0

∂A

Vijδij.

 

 

 

 

 

∂J dt

 

 

 

 

∂J

∂J

Используя это равенство и вычитая (2) из неравенства Клаузиуса —

 

 

 

dh

 

∂qi

 

qi

∂T

 

 

 

 

 

 

 

 

 

 

 

 

Дюгема в виде ρT

 

 

+

 

 

 

 

 

 

− qV >

0, получаем

 

dt

∂xi

T

∂xi

 

 

∂A

 

 

 

 

 

 

∂A

dT

 

 

 

 

 

 

 

 

 

− ρ0δij

 

− σij Vij

ρ

 

 

+ h

 

 

 

 

 

 

 

 

 

 

∂J

∂T

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂A dVij

 

∂A dϑi

qi ∂T

> 0, (3)

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

− ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Vij

 

dt

∂ϑi dt

T

∂xi

откуда в силу того, что скорости dT/dt, dVij/dt и dϑi/dt изменения реактивных переменных произвольны, следуют достаточные условия реализуемости рассматриваемого процесса:

 

∂A

∂A

 

∂A

 

∂A

 

 

 

qi ∂T

>

 

 

h = −

 

,

 

= 0,

 

= 0,

σij − ρ0

 

δij

 

Vij

 

 

 

0.

(4)

∂T

∂Vij

∂ϑi

∂J

 

T

∂xi

Представим каждую функцию из (1) в виде

 

 

 

 

 

 

 

 

A = A+ A(D),

h = h+ h(D),

σij = σ+ σ(D)

,

qi = q+ q(D)

,

(5)

 

 

 

 

 

 

 

ij

ij

 

 

i

i

 

 

где первые слагаемые в правых частях равенств, отмеченные индексом (·), не зависят от Vkl, т. е. представляют собой значения активных переменных при Vkl = 0, а вторые слагаемые, отмеченные индексом

(·)(D), зависят от Vkl и обращаются в нуль при Vij = 0. Тогда с учетом (4) получим, что A(D) = 0, т. е. A = A, а также h(D) = 0, т. е.

h = h. Введем диссипативную функцию δD = σij(D)Vij и представим неравенство (4) в виде

σij

∂A

Vij + δD

qi ∂T

> 0.

− ρ0δij

 

 

 

 

∂J

T

∂xi

Так как первое слагаемое в левой части этого неравенства линейно зависит от Vij, то

σ

= ρ0

∂A

δij.

(6)

 

ij

 

∂J

 

 

 

 

3

Поэтому второй закон термодинамики в данном случае принимает

вид δD

qi ∂T

> 0. Если не учитывать взаимного влияния процесса

T

 

∂xi

 

переноса теплоты и диссипации энергии, то получим два неравенства:

(D)

 

 

 

 

 

 

∂T

δD = σij

Vij > 0

и −qi

 

> 0, причем в случае δD > 0 говорят о вязкой

∂xi

жидкости.

Жидкость, как правило, изотропна, в связи с этим в линейном приближении тензор вязких напряжений имеет компоненты σij(D) =

= λDVkkδij + 2µDVij,

где λD

и µD — коэффициенты, аналогичные

константам Ламе, причем µD принято

называть

динамической

(или сдвиговой) вязкостью.

Основной

единицей измерения λD и

 

 

·

с. Если принять,

что p =

ρ0∂A/∂J — давление

µD является Па

 

жидкости, то с учетом (6) получим реологическое уравнение

 

σ

ij

= σ

+ σ(D) =

ij

+ λ

D

V

kk

δ

ij

+ 2µ

V

ij

,

(7)

 

 

ij

ij

 

 

 

 

D

 

 

 

характеризующее линейную вязкую (или ньютоновскую) жид-

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

Если положить A = A(J ,T ) и определить зависимость вектора плотности теплового потока от реактивных переменных, например, в виде закона Био — Фурье, то, учитывая (2), (4) и (7), после несложных преобразований можно получить уравнение теплопереноса в

жидкости

 

+ 2µDVijVij + ∂xi

λij

∂xj

+ qV − T ∂T Vkk, (8)

ρcε dt = λDVkk2

 

dT

 

 

 

 

 

(T ) ∂T

 

∂p

 

 

 

2A(J ,T )

 

 

 

 

 

 

 

 

где cε = −T

— теплоемкость

при

постоянной деформа-

∂T 2

 

 

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

удельная теплоемкость при постоянном объеме или изохорная теплоемкость, соответствующая изохорному (при постоян-

ном объеме среды) термодинамическому процессу, и обозначение cv.

Подставляя (7) в уравнения движения среды ρ

dvi

=

∂σji

+ bi, полу-

 

 

 

dt

∂xj

чаем уравнения Навье — Стокса — Дюгема для вязкой сжимаемой жидкости

ρ

dvi

+

 

∂p

=

∂(λDVkk)

+ 2

∂(µDVij)

+ bi,

(9)

dt

 

 

 

 

 

∂xi

∂xi

∂xj

 

4

где bi — проекции вектора b плотности объемных сил на оси Oxi. Полная система уравнений, описывающих течение жидкости, в рассматриваемом варианте математической модели содержит шесть неизвестных: ρ, T , p и три проекции vi вектора v скорости. Для их определения необходимо совместно решить уравнения (8), (9) и уравне-

ние неразрывности

+ ρ

∂vj

= 0, используя уравнение состояния

 

dt

 

 

 

 

 

∂xj

 

 

 

 

 

 

ρ = ρ(p,T ).

(10)

Согласно (7) среднее нормальное напряжение в жидкости σkk/3 = = −p + 3κDVkk, где κD = λD + 2µD/3 — объемная вязкость. Если принять, что это напряжение в движущейся жидкости определяется (как и в покоящейся) лишь давлением, то приходим к условию Сток-

са κD = 0. При этом с учетом равенства Vij =

1

 

∂vi

+

 

∂vj

из (9)

 

 

 

 

следуют уравнения Навье — Стокса

2

∂xj

 

∂xi

 

dvi

 

∂p

 

µD

∂vi

 

∂vj

 

 

 

 

∂vj

+ bi.

 

ρ

 

+

 

=

 

 

+

 

 

D

 

 

(11)

dt

∂xi

∂xj

∂xj

∂xi

∂xi

3

 

∂xj

Если жидкость несжимаемая и однородная (ρ = const, ∂vj/∂xj = = 0), а также µD = const, то из (9) получим

ρ

dvi

= −

∂p

+ µD

2vi

+ bi,

или ρ

dv

= −rxp + µDrx2 v + b, (12)

dt

∂xi

∂xj∂xj

dt

где rx и r2x — дифференциальные операторы Гамильтона и Лапласа, вычисляемые в системе пространственных координат с радиус-век-

тором x. Вместе с уравнением неразрывности эти уравнения образуют замкнутую систему относительно неизвестных функций v(x,t) и p(x,t), однако не удается корректно задать граничные условия для давления на непроницаемых границах области V , в которой рассматривается движение жидкости.

При сравнительно малой скорости течения жидкости и выполнении неравенства |vj∂vi/∂xj| |∂vi/∂t| говорят о ползучем движении жидкости, описываемом уравнениями

ρ

∂vi

= −

∂p

+ µD

2vi

+ bi,

∂t

∂xi

∂xj∂xj

которые отличаются от (12) отсутствием инерционных сил, вызванных переносным ускорением частиц жидкости. При bi = 0 и установившемся

движении (∂vi/∂t = 0) имеем

∂p

= µD

2vi

. Дифференцируя обе части

 

 

 

∂xi

∂xj∂xj

5

этого равенства по xi и суммируя по индексу i, с учетом равенства

∂vi = 0 получаем

∂xi

2p

 

2

 

∂v

= 0,

 

= µD

 

i

∂xi∂xi

∂xj∂xj

∂xi

т. е. для такого движения давление является гармонической функцией: r2xp = 0.

Если существует потенциал B0 векторного поля объемных

сил,

2.

е

. b = −rxB0,

то векторную форму

(12)

с учетом равенства

 

т

 

 

rx|v| = 2(v · rx)v + 2v × (rx × v) и выражения для полной производ-

2

dv/dt

2

∂v/∂t −v ×(rx ×v) = −rx(p/ρ +

ной

 

можно представить в виде

 

|v| /2 + B0/ρ) + (µD/ρ)rxv. Затем вычислим ротор левой и правой частей этого равенства и с учетом равенства rx × (b × c) = (c · rx)b − c(rx · b) − (b · rx)c + b(rx · c) запишем

∂(rx × v)

r

 

(

r

 

 

v) v + (

v)

v + (v

 

)

 

2

v

∂t

x

·

 

∂W

 

rx × rx ·

 

 

µD

 

 

 

x ×

 

 

 

· rx

rx × −

− rx · (rx × v) v =

 

 

+ (v · rx)W − (W · rx)v =

 

rxW , (13)

 

∂t

ρ

поскольку rx×v = W — вектор завихренности, rx · (rx×v) = 0 как смешанное произведение векторов с двумя одинаковыми сомножителя-

ми

для несжимаемой жидкости

, rx× (rxf(x) = 0

для любой

, rx ·v = 0

 

 

 

 

 

2

дважды дифференцируемой скалярной

функции

f(x)

и

x

(

 

xv) =

 

 

r

2

2

 

 

 

 

 

r

×

 

= rx(rx×v) = rxW . Так как ∂W /∂t + (v · rx)W = dW /dt, в итоге

получим уравнение переноса завихренности

 

 

 

 

 

 

 

 

dW

 

 

 

 

 

 

 

 

 

 

= (W · rx)v + νDrx2 W ,

 

 

 

 

 

(14)

 

 

dt

 

 

 

 

 

где νD = µD/ρ — кинематическая вязкость, м2/с.

Исключение давления из системы уравнений снимает и проблему формулировки граничных условий для функции p(x,t) на непроницаемых границах области V , но порождает аналогичную проблему по отношению к функции W (x,t). Равенством v = rx×ψ можно ввести векторную функцию тока ψ, тождественно удовлетворяющую ра-

венству ∂vj = 0. Тогда W = rx×(rx×ψ) = rx(rx · ψ) − r2 ψ. После

∂xj x

подстановки этого равенства в (13) получим единственное уравнение относительно векторной функции ψ(x,t), содержащее ее производные по пространственным координатам до четвертого порядка включительно. Преимущество математической модели, включающей это уравнение, состоит в том, что для функции ψ(x,t) удается корректно сформулировать по два граничных условия в каждой точке поверхности S, ограничивающей область V .

6

1.2. Идеальная жидкость

Гипотетическую жидкость, полностью лишенную свойства вязкости, принято называть идеальной. Для получения уравнений, описывающих ее движение, достаточно в (11) положить µD = 0. В итоге получим уравнения Эйлера

ρ

dvi

= −

∂p

+ bi, i = 1, 2, 3,

или ρ

dv

= −rxp + b, (15)

dt

∂xi

dt

где ρ — плотность жидкости; vi и bi — проекции векторов v скорости и b плотности объемных сил на оси Oxi системы пространственных координат; t — время; p — давление; rx — дифференциальный оператор Гамильтона, вычисляемый в системе пространственных координат с радиус-вектором x. Для неподвижной жидкости из (15) следуют уравнения ∂p/∂xi = bi, или rxp = b, описывающие состояние

гидростатического равновесия.

Если ρ зависит лишь от p, т. е. ρ = ρ(p), то жидкость называют

баротропной и вводят функцию давления

p

ρ(p),

(16)

P (p) = Z

 

dp

 

p0

 

 

 

равную работе, совершаемой при движении единицы массы баротропной жидкости при изменении давления от p0 до p. Во многих приложениях векторное поле плотности массовых сил не зависит от времени и обладает потенциалом B(x), поэтому b/ρ = −rxB. При перечисленных допущениях, учитывая rx|v|2 = 2(v · rx)v + 2v × (rx × v) и выражения для полной производной dv/dt, запишем векторную форму (15) в виде

∂v

+ W × v = −rx B + P +

 

v 2

,

 

 

|

|

(17)

∂t

 

2

где W = rx×v — вектор завихренности.

Если из массовых сил действует только сила тяжести, то B = = gh, где g — постоянное ускорение свободного падения, а h —

высота, отсчитываемая от некоторого уровня. Величины h, hp = P/g и hv = |v|2/(2g) называют напором соответственно геометрическим,

пьезометрическим и скоростным. Сумма значений этих напоров

составляет полный напор H = h + hp + hv. С

учетом введенных

обозначений вместо (17) получим

 

 

∂v

 

 

 

+ W × v = −grxH.

(18)

 

∂t

7

Поскольку rx ·(rx×v) = rx ·W = 0, согласно теореме Остроградского — Гаусса для области V пространства, ограниченной поверхностью S и лежащей в области определения поля скоростей, имеем

Z

W · ndS = 0,

(19)

S

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

хревой нитью.

Пусть dS1 и dS2 — площадки соседних нормальных сечений вихревой нити, а вектор завихренности направлен от dS1 к dS2. Так как поток вихря через боковую поверхность вихревой трубки равен нулю, то, применяя (19) к участку вихревой нити, ограниченному рассматриваемыми сечениями, получаем W1 dS1 = W2 dS2, где W1 и W2 — модули вектора W в сечениях dS1 и dS2 соответственно. Таким образом, |W | изменяется вдоль вихревой нити обратно пропорционально площади ее поперечного сечения. Вихревая нить (или трубка) не может оканчиваться внутри жидкости, она или замкнута или оканчивается на поверхности, ограничивающей жидкость. Справедлива теорема Кельвина (Томсона) о том, что при движении идеальной баротропной жидкости в потенциальном поле массовых сил циркуляция вектора скорости по замкнутому контуру остается постоянной во времени. Следовательно, такая жидкость обладает свойством сохранять безвихревое движение, определяемое выполнением условия W = rx×v = 0 в каждой точке области, занятой этой жидкостью. В этом случае векторное поле скоростей обладает потенциалом Φ(x,t):

v(x,t) = rxΦ(x,t).

(20)

Отсюда для несжимаемой жидкости получим уравнение Лапласа

rx2 Φ(x,t) = 0,

(21)

т. е. Φ(x,t) является гармонической функцией.

Для безвихревого движения сжимаемой жидкости из (17) и (20)

следует интеграл Лагранжа — Коши

∂Φ

+ B + P +

|rxΦ|2

= f(t),

(22)

∂t

2

 

 

 

8

где f(t) — одинаковая для всей области движения жидкости функция времени, определяемая обычно или из граничных условий, или по зависимости от t левой части (22) в какой-либо одной точке этой области. Если движение жидкости установившееся, то ∂Φ/∂t = 0, f(t) =

= const и (22) переходит в интеграл Бернулли B + P + |v|2/2 = const,

который в случае несжимаемой жидкости плотностью ρ0 имеет вид

 

 

 

 

 

 

ρ0B + p +

ρ0|v|2

= const.

(23)

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Дифференциал функции давления можно представить в виде dP =

=

 

1

dp =

a2

dρ, где a2

=

dp

> 0, или

dP

=

a2

. Тогда, используя (20) и

 

 

 

dt

 

 

ρ

 

ρ

 

 

 

 

 

 

 

ρ

dt

 

уравнение неразрывности (3.31), получаем

 

 

 

 

 

 

 

 

 

1

 

dP

+ rx2 Φ = 0.

(24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Система уравнений (22), (24) замкнута относительно неизвестных функций P и Φ, поскольку a2 можно также представить как функцию P .

В ряде прикладных задач аэродинамики и акустики течение среды можно рассматривать как возмущенное относительно известного движения или состояния покоя жидкости. Так, при малых возмущениях относительно состояния покоя, в котором ρ = ρ0 = const, систему (22), (24) можно линеаризовать, если в (22) пренебречь слагаемым |rxΦ|2/2

и положить f(t) ≡ 0, а в (24) принять a2 = a02 =

dp

ρ=ρ0

. Тогда при отсут-

 

 

∂Φ

 

 

1 ∂P

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ствии массовых сил получаем ∂t + P = 0 и

0

 

 

 

 

 

xΦ = 0. Исключая

a2 ∂t +

 

 

отсюда P или Φ, приходим к волновым уравнениям

 

 

2Φ

2P

 

 

 

 

 

 

 

 

 

 

 

 

 

= a02rx2 Φ или

 

 

= a02rx2 P,

(25)

 

∂t2

∂t2

 

относящимся к дифференциальным уравнениям гиперболического типа. Значение a0 в (25) является скоростью распространения малых возмущений в среде и носит название скорости звука.

При установившихся колебаниях с некоторой частотой ω функцию Φ можно представить произведением зависящей только от времени t периодической функции вида sinωt или cos ωt (или линейной комбинации этих функций) и функции Φ(x), зависящей лишь от пространственных координат и описывающей форму волны. Тогда первое уравнение (25) переходит в уравнение эллиптического типа r2xΦ+ (ω2/a20= 0,

называемое уравнением Гельмгольца.

dr
dR(r)

9

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

т. е. в отличие от вязкой жидкости отсутствует эффект прилипания

частиц к границе. Поэтому на непроницаемых границах совпадают проекции векторов скорости идеальной жидкости и заданной скорости vграницы на направление нормали, т. е. v · n = v· n, или (в случае безвихревого движения жидкости) (rxΦ) · n = v· n. При построении математической модели (ММ) течения идеальной жидкости на ее свободной поверхности должно быть задано давление, а на проницаемых границах области — вектор скорости или давление.

Пусть твердый шар радиусом r0 и массой m движется поступательно вдоль оси Ox3 со скоростью v0 в идеальной несжимаемой жидкости плотностью ρ0. Приняв обтекание шара безвихревым и осесимметричным относительно этой оси, представим (21) в подвижной системе сферических координат с началом в центре шара:

 

∂Φ

 

 

1 ∂

 

∂Φ

 

 

r2

 

 

 

+

 

 

 

 

 

sinϑ = 0,

(26)

∂r

∂r

sinϑ

∂ϑ

∂ϑ

где r — радиальная координата; ϑ [0, π] — угол, отсчитываемый от положительного направления оси Ox3 (рис. 0.1).

Рис. 0.1

На поверхности S шара (при r = r0) ∂Φ/∂r = v0 cos ϑ, что дает основание искать решение (26) в виде Φ(r,ϑ) = R(r)cos ϑ, удовлетворяющем

(26) при условии drd r2 − 2R(r) = 0. Подставляя сюда R(r) =

= rk, находим k(k + 1) − 2 = 0, т. е. k1 = 1 и k2 = −2. Таким образом, R(r) = C1r + C2/r2. В случае C1 6= 0 при r → ∞ имеем R(r) → ∞, что

10

противоречит физическому смыслу, поскольку вдали от движущегося шара жидкость неподвижна. Поэтому C1 = 0, а вторую константу

можно найти из

 

 

 

3

 

 

 

r = r0

3

/2.

 

-

 

 

 

 

 

 

: C2 = −v0r0

В ито

 

граничного условия при

 

 

 

 

 

 

 

 

 

 

ге получим Φ(r,ϑ) = −v0

 

r0

cos ϑ. Составляющие вектора v скорости

2r2

жидкости при обтекании движущегося шара будут равны

 

 

 

 

∂Φ

 

r3

1 ∂Φ v0 r3

 

 

 

vr =

 

 

= v0

 

0

cos ϑ, vϑ =

 

 

 

=

 

 

 

0

sinϑ.

 

(27)

 

 

 

 

 

 

 

 

 

 

∂r

 

r3

r ∂ϑ

 

2 r3

 

 

 

Кинетическая энергия жидкости при движении шара с учетом (27) составит

∞ π

+ vϑ2)r2 sinϑdr dϑ =

 

 

 

 

 

K = 2π Z Z

 

2 (vr2

 

 

 

 

 

 

 

ρ0

 

 

 

 

 

 

 

 

 

r0 0

 

 

 

 

dr

π

 

 

 

 

 

 

πρ0v2r6

1

 

 

v2

 

=

0

0

Z

 

Z (3cos2 ϑ + 1) sinϑdϑ =

 

πρ0v02r03

= m

0

,

4

 

r4

3

2

 

 

 

 

 

r0

 

0

 

 

e

 

 

где m = 2πρ0r3/3 — половина массы жидкости, вытесненной шаром. e 0

Суммарная кинетическая энергия системы «шар — жидкость» при поступательном движении шара будет KΣ = (m + m)v02/2. Изменение KΣ за время dt, согласно закону сохранения энергии, равно работе, со-

вершаемой силой

P ,

приложенной к шару

,

на перемещении

т е

.

 

 

 

 

 

(m + m)v02

 

 

 

 

e

 

v0 dt, .

 

 

или

 

d

 

 

 

 

 

 

dv0

 

Отсюда полу

 

 

 

 

 

 

2e

= (m + m)v0 dt = P v0.

 

dKΣ = P v0 dt, dv0

 

 

dt

 

 

 

-

чаем

(m + m)

 

= P .

Таким образом

,

ММ поступательного движения

 

dt

 

 

e

 

 

 

 

 

 

 

 

шара

массой

m

в идеальной несжимаемой жидкости плотностью

ρ0

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

эквивалентна ММ поступательного движения в вакууме материальной точки массой m + me . В связи с этим величину me называют присоеди-

ненной массой шара.

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

тензора второго ранга коэффициентов присоединенных масс.

Если на поле скоростей, определяемое соотношениями (27), наложить поле с составляющими vr = −v0 cos ϑ и vϑ = v0 sinϑ вектора ско-

рости, то получим поле скоростей с составляющими

 

 

vr = −v0

r3

cos ϑ,

vϑ = v0

1 +

r3

(28)

 

1 − r3

2r3 sinϑ,

 

 

0

 

 

 

0

 

 

которое

возникает при установившемся обтекании неподвижного шара

b

 

 

b

 

 

 

 

идеальной несжимаемой жидкостью, имеющей вдали от шара скорость

Соседние файлы в папке МММСС (От Кувыркина Г.Н.)