Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Gidromekhanika_1392-1_2 / Конспект лекций по гидромеханике.pdf
Скачиваний:
224
Добавлен:
27.02.2016
Размер:
2.25 Mб
Скачать

КАФЕДРА «СУДОСТРОИТЕЛЬНОЕ ПРОИЗВОДСТВО»

 

ТЕХНИЧЕСКАЯ ФИЗИКА

гидромеханика

ЛЕКЦИЯ 8.

ГИДРОДИНАМИКА ВЯЗКОЙ ЖИДКОСТИ

z

 

 

 

 

 

 

 

8.1 Модель вязкой жидкости

 

 

 

 

 

 

 

 

 

 

 

τzx

 

 

 

Приступая

к

рассмотрению

 

 

 

 

 

 

 

движения вязкой жидкости, необходимо

 

 

 

 

 

 

 

прежде всего уяснить терминологию , т. е.

 

 

 

 

 

 

 

 

 

 

 

 

ux

смысл, вкладываемый в понятие «вязкая

 

 

 

 

 

 

 

 

жидкость». С математических позиций

O

 

 

x

необходимо

установить

вид

 

 

функциональной

зависимости

для

y

 

 

 

 

напряжений, либо, другими словами,

 

 

 

 

сформировать модель вязкой жидкости. В

 

 

 

 

дальнейшем под вязкой мы будем понимать

жидкость, удовлетворяющую трем гипотезам: линейности, однородности и изотропности.

Рис. 8.1

 

8.1.1. Гипотеза линейности .

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

плоскости xOy (рис. 8.1), что дает

 

ux

 

τ

zx

= µ

 

z

 

 

 

 

 

Воспользуемся результатом, полученным при рассмотрении теоремы Гельмгольца о движении жидкой частицы. Согласно теореме, скорость угловой

деформации относительно оси y

 

 

 

1

 

ux

 

uz

 

γ

y

=

 

+

 

 

 

 

 

2

 

z

 

x

 

 

 

 

 

 

Так как движение происходит в плоскости xOy, то uz = 0 и

γ y =

1 ux

 

2 z

 

 

 

и, следовательно, касательное напряжение

(8.1)

τzx = 2µγ y

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

Page 68 из 87

1

Конспект лекций по гидромеханике

КАФЕДРА «СУДОСТРОИТЕЛЬНОЕ ПРОИЗВОДСТВО»

 

ТЕХНИЧЕСКАЯ ФИЗИКА

гидромеханика

деформаций, и связаны с ними линейной зависимостью. При этом коэффициент

пропорциональности остается неизменным и равным 2µ.

Кроме того, согласно закону Стокса касательные напряжения, как показано выше, пропорциональны скоростям угловой деформации, а

нормальные - скорости линейной деформации, т.е.

u

x

,

uy

,

u

z

.

 

 

 

 

 

 

 

x

y

 

z

Таким образом, можем записать

uy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

x

 

 

 

 

 

 

 

(8.2)

τxy =τyx = 2µγ z = µ

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

 

 

 

 

 

 

 

 

и т.д.

Рассмотрим теперь нормальные напряжения, возникающие от сил вязкости. Согласно закону Стокса, их можно записать в виде так называемых девиаторов напряжения, имеющих вид:

σxx = 2µ uxx

u

σ yy = 2µ yy (8.3)

σzz = 2µ uzz

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

pxx = −p + 2µ uxx pyy = −p + 2µ uyy pzz = −p + 2µ uzz

Выполним следующую операцию: из утроенной величины ( pxx + pyy + pzz ). Это дает:

3 pxx (pxx + pyy + pzz )= −3 p +6µ uxx

(8.4)

pxx вычтем сумму

− −3 p + 2µ uxx + uyy + uzz = 6µ uxx 2µ div ur

откуда найдем

 

u

x

 

2

r

 

pxx + pyy + pzz

pxx = 2µ

 

 

µ div u

+

 

 

 

3

3

 

x

 

 

Вкачестве давления в вязкой жидкости принимают среднее

арифметическое, т.е. p = − pxx + pyy + pzz . И, следовательно,

3

Page 69 из 87

2

Конспект лекций по гидромеханике

КАФЕДРА «СУДОСТРОИТЕЛЬНОЕ ПРОИЗВОДСТВО»

 

ТЕХНИЧЕСКАЯ ФИЗИКА

 

 

 

 

 

 

гидромеханика

pxx = −p +2µ

ux

 

2

µ div ur

 

x

3

 

 

 

 

 

 

 

uy

 

 

2

r

(8.5)

pyy = −p +2µ

 

 

 

 

µ div u

y

 

3

pzz = −p + 2µ

uz

 

 

2

µ div ur

 

 

 

3

 

 

z

 

 

 

 

Для несжимаемой жидкости divu = 0 , и выражения упрощаются.

8.1.2. Гипотеза однородности

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

8.1.3. Гипотеза изотропности

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

8.2 Уравнение движения вязкой жидкости. (уравнение Навье-Стокса)

18 марта 1822 года в докладе, представленном Французской академии наук, Клод Луи Навье писал о полученных им уравнениях: «Хотя уравнения основаны на гипотезе Ньютона о том, что касательные напряжения пропорциональны скорости деформации, никак нельзя сказать, что они не выражают ничего нового».

Уравнения движения вязкой жидкости можно получить из уравнений движения в напряжениях (2.16), выполнив некоторые преобразования. Рассмотрим лишь одну проекцию этих уравнений:

du

x

= X +

1 p

 

+

∂τyx

+

∂τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xx

 

 

 

 

zx

dt

 

ρ

x

 

y

 

z

 

Как было показано при рассмотрении модели вязкой жидкости, нормальные напряжения

pxx = −p + 2µ uxx 32 µ div ur

Для упрощения задачи будем считать жидкость несжимаемой ( divur = 0 ), тогда

 

 

pxx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

2

ux

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ux

 

 

 

 

 

 

 

(8.6)

 

 

 

 

 

=

 

 

 

 

p + 2µ

 

 

 

 

 

 

 

= −

 

 

 

 

 

+ 2µ

 

 

 

 

 

 

 

 

x

 

x

 

 

x

 

x

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Касательное напряжение

 

 

 

 

 

 

 

uy

 

 

 

u

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τyx

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= µ

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂τyx

 

 

 

 

 

 

uy

 

 

 

u

 

 

 

 

 

 

2uy

 

2u

 

 

(8.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

= µ

 

 

 

 

 

+

 

 

 

 

 

= µ

 

 

 

 

 

 

 

 

+ µ

 

 

 

 

2

 

 

 

y

 

y

 

 

x

 

y

 

yx

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂τzx

 

 

 

 

 

 

 

ux

 

 

 

 

 

 

 

 

 

 

 

 

 

2

ux

 

2

uz

 

 

 

 

 

 

 

 

 

 

 

 

uz

 

 

 

 

(8.8)

 

 

 

 

= µ

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

= µ

 

 

 

 

 

 

+ µ

 

 

 

 

 

 

 

 

 

z

 

y

 

z

 

 

x

 

z

2

zx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Суммируя (8.6), (8.7) и (8.8) и группируя члены, получаем:

Page 70 из 87

3

Конспект лекций по гидромеханике

 

 

КАФЕДРА «СУДОСТРОИТЕЛЬНОЕ ПРОИЗВОДСТВО»

ТЕХНИЧЕСКАЯ ФИЗИКА

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гидромеханика

 

p

 

2u

 

2u

 

2u

 

 

 

2u

 

2uy

 

2u

z

 

 

+ µ

 

 

x +

 

 

x +

 

 

x

+ µ

 

 

x +

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

2

y

2

z

2

 

 

x

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yx zx

Третий член можно записать в виде:

µ x uxx + uyy + uzz = µ x div ur

но жидкость несжимаема, и divu = 0 . Таким образом получаем:

dux

 

1 p

 

µ

 

2

 

 

 

2

 

 

 

2

 

 

 

 

= X

+

 

 

ux

+

 

ux

+

 

ux

(8.9)

 

 

 

 

 

 

 

2

 

 

 

2

 

 

dt

ρ ∂ x

ρ

 

x

y

z

2

 

 

 

 

 

 

 

 

 

 

 

Выражение в скобках есть ни что иное, как оператор Лапласа - 2ur , а µρ =ν . Окончательно получаем:

ux

ux

 

ux

 

ux

 

1 p

 

 

2

ux

 

2

ux

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

ux

(8.10)

 

+ux x

+uy

 

+uz

 

= X

 

+ ν

 

 

+

 

+

 

 

t

y

z

ρ

x

x2

y2

z2

Аналогично можно расписать и две другие проекции. Полученная система уравнений движения вязкой жидкости и носит название системы уравнений Навье-Стокса.

В векторной форме можно записать

r

r

 

1

2 r

(8.11)

a

= F

 

grad p +ν u

ρ

Как следует из (8.11), это уравнение отличается от уравнения движения идеальной жидкости дополнительным членом (ν 2ur), учитывающим действие сил вязкого трения.

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

величины: ux , uy , uz и p. Принципиально это оказывается возможным, так как

три уравнения Навье-Стокса (в проекциях) плюс уравнение неразрывности образуют замкнутую систему. Плотность и вязкость, входящие в них,

считаются известными, а проекции массовых сил (X, Y, Z) задаются условиями конкретной задачи.

С чисто математических позиций уравнения Навье-Стокса относится к классу нелинейных дифференциальных уравнений в частных производных второго порядка. Одно из наиболее неприятных из их свойств - нелинейность, обусловленная наличием конвективных членов ускорения. Следует отметить, что до настоящего времени вследствие практически непреодолимых математических трудностей не получено ни одного общего решения уравнений Навье-Стокса в их полном виде, т.е. при сохранении всех конвективных членов и всех членов, учитывающих вязкость. Известны лишь отдельные частные решения.

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

Page 71 из 87

4

Конспект лекций по гидромеханике