Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
attachments_24-10-2012_21-31-48 / МСС-Lekt7(ldealnaya)PDF.pdf
Скачиваний:
310
Добавлен:
09.02.2015
Размер:
554.6 Кб
Скачать

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

1.5Задача о движении сферы в бесконечной идеальной несжимаемой жидкости.

Пусть абсолютно твердая сфера радиуса a движется поступательно со скоростью V (t) относительно неподвижной системы декартовых координат в неограниченном объёме идеальной несжимаемой жидкости. Массовые силы положим равными нулю, т.е.

~

F = 0.

Обозначим неподвижные координаты через x1; x2; x3. Свяжем с телом декартову систему координат x1; x2; x3, которая движется поступательно вместе со сферой в направлении оси x1. Оси подвижных координат направим по осям неподвижной системы координат, а центр расположим в центре сферы.

p

Обозначим через r = xixi = x21 + x22 + x23 расстояние от начала координат до некоторой точки подвижного пространства, а через # угол между лучем Or и осью Ox1, причем cos # = x1=r (рис. 5).

Рис. 5: Движение сферы в идеальной жидкости

При движении сферы частицы жидкости отодвигаются сферой в направлении нормали к поверхности тела со скоростью V cos #. Этим действием сфера приводит в движение всю жидкость. Пусть ~v(x) = [v1(x); v2(x); v3(x)] вектор скорости жидкости относительно подвижных координат. На границе сферы проекция вектора скорости ~v жидкости на нормаль к сфере будет равна скорости точки сферы в направлении нормали

 

 

 

vn = v~n r=a = vr r=a = V cos #

(1.66)

Граничное условие (1.66) называется условием непроницаемости (поток жидкости внутрь тела равен нулю). Условие непроницаемости диктует характер взаимодействия тела с жидкостью. При 0 6 # < =2 жидкость продвигается вперёд и в стороны. На экваторе, т.е. при # = =2 радиальная компонента вектора скорости vr равна нулю, тангенциальная компонента вектора скорости v# = V (t). При =2 < # 6 жидкость увлекается за сферой и, тем самым, обеспечивается безотрывность обтекания.

Другое условие на движение жидкости состоит в том, что на бесконечности жидкость покоится, т.е.

lim vr ! 0 (1.67)

r!1

14

В нашем случае массовые силы отсутствуют, несжимаемая однородная жидкость является баротропной жидкостью с функцией давления P = p=%0 + const:, движение жидкости начинается из состояния покоя, т.е. при t = 0 вихри в жидкости отсутствуют. Следовательно, все условия теоремы Лагранжа о сохраняемости потенциальности течения выполнены.

Течение жидкости вызванное движением в ней жесткой сферы будет потенциальным, т.е. ~v = grad ', при этом

vr = v~n = grad ' ~n =

@'

ni =

@'

;

v# = v~ =

1 @'

;

(1.68)

 

 

 

 

 

 

 

@r

 

r @#

 

@xi

 

 

 

 

 

Где ~n и ~ единичные векторы нормали и касательной к поверхности сферы. Потенциальное течение несжимаемой однородной жидкости описывается уравнением Лапласа (1.54)

' = 0

(1.69)

и интегралом Коши-Лагранжа в подвижных координатах (1.42), в котором следует положить U = 0, P = p=%0 + const:, ~vпер = [V (t); 0; 0], а функцию f(t) приписать к потенциалу. В результате интеграл Коши-Лагранжа примет вид:

@'

 

@'

 

1

 

p

 

(1.70)

 

 

 

 

V (t) +

 

 

jgrad 'j2

+

 

+ const: = 0

@t

@x1

2

%0

Константу выбираем из условия на бесконечности для давления p(x; t) ! p1. По-

r!1

сле этого из интеграла Коши-Лагранжа найдём распределение давления

@'

@'

1

(1.71)

p(x; t) = p1 %0

 

+ %0

 

V (t) %0

 

jgrad 'j2

@t

@x1

2

Искомый гармонический потенциал должен удовлетворять условию непроницаемо-

сти (1.66) и условию на бесконечности, т.е.

 

 

 

 

 

 

 

@'

r=a =

V

cos

# ;

lim '

!

0 ;

lim

@'

!

0

(1.72)

@r

 

 

r!1

 

r!1

@r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задача (1.69), (1.72) является классической задачей Неймана для внешности шара. Её решением, убывающем на бесконечности, является гармоническая функция

' = A@x1

r

= A r3 = A

r2

(1.73)

@

1

 

x1

cos #

 

Константу A найдём из условия непроницаемости

 

@r

r=a

= @r A

r2

r=a

@'

 

@

 

 

cos #

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

V a3

 

 

 

 

 

 

 

 

A =

 

 

 

 

 

 

 

 

 

и потенциал скорости принимает вид:

2

 

 

 

 

 

 

 

'(x; t) =

a3V (t) cos #

 

 

 

 

2r2

=

2A cos #

= V cos #

(1.74)

a3

 

 

 

=

a3V (t)x1

(1.75)

2r3

15

Далее найдем компоненты вектора скорости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a3V (t)

 

x1

 

 

a3V (t)

1 +

3x2

 

 

a3V (t)

 

 

v1

= ';1

=

 

 

 

 

;1

=

 

 

 

 

1

 

=

 

 

 

 

1 + 3 cos2

#

2

 

r3

 

2r3

 

r2

 

2r3

 

 

 

a3V (t)x

1

1

 

 

 

 

3a3V (t) x1x2

 

 

 

3a3V (t)

 

 

v2

= ';2

=

 

 

 

 

 

;2

=

 

 

 

 

 

=

 

 

 

 

 

sin # cos #

 

2

 

 

 

r3

2

 

 

r5

 

 

 

2r3

 

 

 

 

 

a3V (t)x

1

1

 

 

 

 

3a3V (t) x1x3

 

 

 

3a3V (t) x3

 

 

v3

= ';3

=

 

 

 

 

 

;3

=

 

 

 

 

 

=

 

 

 

 

 

 

cos #

 

2

 

 

 

r3

2

 

 

r5

 

 

 

2

 

 

r4

 

(1.76)

(1.77)

(1.78)

Найдём радиальную и тангенциальную составляющие вектора скорости в плоскости x1x2

@'

a3V (t) cos #

 

 

@'

 

a3V (t) sin #

 

) j~vj =

a3V (t)

vr(r; #) =

 

=

 

 

 

;

v#(r; #) =

 

=

 

 

 

 

 

 

@r

 

r3

 

r@#

 

r3

 

r3

На поверхности сферы, т.е. при r = a

 

 

 

 

 

 

 

 

(1.79)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j~v(a)j = q

 

= V (t)

(1.80)

vr(a; #) = V (t) cos # ;

 

v#(a; #) = V (t) sin #

)

 

vr2 + v#2

Отсюда видно, что линии тока начинаются и заканчиваются на поверхности сферы, поскольку

vr(a; #) = V (t) cos # ;

v#(a; #) = V (t) sin #

vr(a; #) = V (t) cos # ;

v#(a; #) = V (t) sin #

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 6: Скорость жидкости в симметричных точках и линии тока (помечены красным цветом)

Давление на поверхности сферы определяется по формуле (1.71), которая для случая постоянной скорости движения сферы принимает вид:

@'

 

% v2 a

)

 

% V 2

3 cos2 # 2

 

p(a; #) = p1 + %0

 

(a)V

0 (

= p1 +

0

(1.81)

@x1

2

 

2

16

1.5.1Задача об обтекании сферы идеальной несжимаемой жидкостью.

Пусть теперь жидкость обрекает неподвижную сферу. Систему координат оставим прежней Ox1x2x3 с центром O в центре сферы. Жидкость течет из бесконечности вдоль оси x1 против её направления, т.е. скорость жидкости на бесконечности равна V . Такое течение, также как и в предыдущем случае является потенциальным течением.

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

 

 

= 0 ;

 

 

 

(1.82)

условия непроницаемости

 

 

 

 

 

 

 

wn = w~n r=a = wr

r=a =

@

r=a

= 0

(1.83)

@r

 

 

 

 

 

 

 

 

и условия на бесконечности

 

 

 

~

 

 

(1.84)

 

rlim grad

 

 

 

 

 

 

! V k1

 

 

!1

 

 

 

 

 

 

Рассмотрим случай, когда сфера отсутствует, а течение есть. Тогда такое течение описывается потенциалом

'0 = V x1

(1.85)

Очевидно, что уравнение Лапласа и условие на бесконечности для этого потенциала удовлетворяются.

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

 

 

= ' + '0 =

a3V (t)x

1

V (t)x1 = V r +

 

a3

 

cos #

 

(1.86)

 

 

2r3

 

2r2

 

Вычислим компоненты вектора скорости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

= V 1

 

a3

cos # ;

 

 

 

 

 

@'

 

 

 

1 +

a3

sin #

(1.87)

wr(r; #) =

 

 

v#(r; #) =

 

= V

 

@r

r3

r@#

2r3

Проверяем условие непроницаемости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

 

 

 

a3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@r

 

r=a = V cos # 1 r3 r=a = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие на бесконечности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a3V x1

 

 

~

 

 

 

 

~

 

 

 

 

 

 

r!1

 

 

= r!1

 

 

 

2r3

 

 

V k

1 =

 

V k

1

 

 

 

 

 

lim grad

 

 

lim grad

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, все условия задачи (1.82)–(1.84) выполняются, а потенциал

, пред-

ставленный формулой (1.86) является решением задачи (1.82)–(1.84) об обтекании жесткой недеформируемой сферы идеальной несжимаемой жидкостью.

2по крайней мере очевидно, что скорость лобовой точки сферы станет равной нулю

17

Соседние файлы в папке attachments_24-10-2012_21-31-48