Специальный курс "Математические модели в гидродинамике". Глушко А.В - 12 стр.

UptoLike

Составители: 

12
Распределение скоростей в произвольно малой частице сплошной
среды
Рассмотрим совокупность точек среды с координатами
iiii
xdxx
ρ
+=+
,
удаленных от точки
O
(
)
123
,,
на расстояние
0,
,,
ii
dd
ρξξρ
=≤ где
ρ
0
0
>
достаточно малое число.
Поле скоростей
(
)
(
)
(
)
1
,0,
vvxtC
=×∞
будем считать непрерывно
дифференцируемым. Обозначим скорость точки
O
:
0
(,)
vxtv
=
, а скорость
произвольно выбранной точки
1
O
(
)
xdx
+
- через
1
v
. За время
t
вектор
1
OO
ρ
=
перейдет в
///
1
О O
ρ
=
и
(
)
//
;
ρρρρ
==
очевидно, что
/
ρ
+
0
v
t
∆=
ρ
+
1
vt
или
/
ρ
=
ρ
(
)
10
vvt
−∆
.
Разложим
1
v
в окрестности точки
O
:
()()
(
)
0
3
100
0
0
1
1lim11
j
j
j
v
vv
x
ρ
ρροο
=
=++=
(1)
Отсюда
3
/
0
0
1
(1)
j
j
j
dv
ot
dx
ρρρρ
=

=++∆



(2)
Пусть
123
,,
eee
- орты осей
123.
,,
OxOxOx
. Тогда
()
3
112233
1
,
kk
k
vxtveveveve
=
=++=
.
Из (1) имеем
()
33
100
11
1
k
kj
jk
j
v
vve
x
ρρο
==
=++
∑∑
Методом «добавить и отнять» получим
()
3333
100
1111
11
1
22
jj
kk
kjkj
jkjk
jkjk
vv
vv
vvee
xxxx
ρρρο
====

∂∂
∂∂
=+++++


∂∂

∑∑
(3)
Обозначим
1
;
2
jj
kk
jkkj
jkjk
vv
vv
e
xxxx
ω

∂∂
∂∂
=+=−


∂∂

;
;;,1,3
kjkj
eekjωω

===

;
                                                              12


     Распределение скоростей в произвольно малой частице сплошной
среды
     Рассмотрим совокупность точек среды с координатами xi +dxi =xi +ρi ,
удаленных от точки O ( x1 , x2 , x3 ) на расстояние ρi =d ξ i , dξ ≤ρ0, , где ρ0 >0
– достаточно малое число.
Поле     скоростей      v =v ( x, t )∈C1 (Ω ×(0, ∞))                      будем         считать     непрерывно
дифференцируемым. Обозначим скорость точки O : v ( x, t ) =v0 , а скорость
произвольно выбранной точки O1                              ( x +dx )   - через v1 . За время ∆t вектор

O O1 = ρ перейдет              в О / O1/ = ρ /                     и    ( ρ =ρ; ρ   /
                                                                                        =ρ /   ) очевидно,   что

                                                        (          )
  /                                     /
ρ +v0 ∆t = ρ + v1 ∆t или ρ = ρ v1 −v0 ∆t .

        Разложим v1 в окрестности точки O :
                                     ∂v
                                                                          ( lim ο (1)=1)
                                3
                   v1 =v0 +∑                       ⋅ ρ j +ρ0 ο (1)                                           (1)
                                j =1 ∂x j     0                            ρ0 → 0


Отсюда
                               �            3
                                                       dv                 �
                        ρ =ρ +�
                          /
                                 �          ∑ dx               ρ j +ρ0 o(1)� ∆t
                                                                            �
                                                                                                             (2)
                                   �        j =1        j    0
                                                                              �
Пусть e1 , e2 , e3 - орты осей Ox1 , Ox2 , Ox3. . Тогда
                                                                            3
                              v ( x, t ) =v1 e1 +v2 e2 +v3 e3 =∑ vk ek .
                                                                           k =1

Из (1) имеем
                                                   3    3
                                                          ∂vk
                                v1 =v0 +∑ ∑                    ek ρ j +ρ0 ο (1)
                                                j =1 k =1 ∂x j

Методом «добавить и отнять» получим


              1 3 3 � ∂v          ∂v� j 1 3 3 � ∂v ∂v� j
      v1 =v0 + ∑ ∑ � k + � ek ρ j + ∑ ∑ � k + � ek ρ j +ρ0 ο (1)                                             (3)
              2 j =1 k =1 �� ∂x j ∂xk�� 2 j =1 k =1 �� ∂x j ∂xk��

                1 � ∂v ∂v �                               � ∂vk ∂v� j
Обозначим e jk = � k + j � ;ωkj                         =�       −      ; e =�� e��kj ;ω�� =�� ωkj ; k , j =1,3 ;
                2 �� ∂x j ∂x�k �                           � ∂x j ∂x��k
                                                            �        �