ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. “Машиностроение”. 2012
31
( )
;
2
d
f
d
ψ
ψ ψ
π
=
(12)
( , )
( )
sin ,
( )
ij
ij
ij
ij
g g
f
d
d
g g
σ
θ
θ θ
θ θ
σ
=
(13)
где
,
ψ θ
–
азимутный и зенитный углы сферической системы коор-
динат;
σ
–
полное сечение столкновений
Цикл из шагов 2–4 повторяется ровно
c
s
раз, т. е.
c
c
1
.
s
s
S t S
+
≤ Δ <
(14)
Этап II.
Алгоритм сдвига можно представить выражением сме-
щения каждой
i
-
й частицы:
(
) ( )
,
i
t
t
t
t
+ Δ = + Δ
r
r
c
(15)
где
r
–
радиус-вектор сферической системы координат.
На этом этапе также моделируется взаимодействие частиц с по-
верхностью канала.
Результаты моделирования.
Численный эксперимент дает боль-
шое количество информации об эволюции
( )
t
C
каждой частицы рас-
сматриваемой системы. На рис. 3 и 4 представлены изменения давления
Рис. 3. Изменение давления в канале по длине и радиусу