Вычислительный эксперимент в аэродинамике вентиляции. Аверкова О.А. - 21 стр.

UptoLike

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

21
1. По заданным начальным условиям (при t=0
0
, 1,2,3,4
j j
y y j= = ) вычисляются коэффициенты:
1 0 1
2 0 2
3 0 3
4
; 0,5 ;
; 0,5 ;
; ;
.
j j j j j
j j j j j
j j j j j
j j
k h f y y k
k h f y y k
k h f y y k
k h f
= = +
= = +
= = +
=
При этом производится вычисление по одной формуле при
1,2,3,4
j
=
, а лишь затем необходимо перейти к следующей.
2. Определяется положение пылевой частицы и составляющие ее
скорости через время h:
0 1 2 3 4
j j j j j j
y y k k k k= + + + + ,
1,2,3,4.
j
=
3. Переопределяются начальные условия для следующего шага по
времени:
0
, 1,2,3,4
j j
y y j= =
и осуществляется переход к пункту 1. Вычисления производятся до
тех пор, пока пылевая частица не достигнет заданного сечения.
При расчете траекторий пылевой частицы в области с изменяющи-
мися во времени граничными условиями необходимо в каждый момент
времени для определения скорости воздуха пересчитывать интенсив-
ность источников (стоков), распределенных по границе течения.
1.2. Расчет максимального диаметра пылевых частиц в местном
отсосе от укрытия с одинарными стенками
Исследуемое укрытие (рис.4) используется для локализации пыле-
выделений узла загрузки конвейера на рудоподготовительных фабри-
ках.
                                                 21

   1.     ��       ��������            ���������                 ��������     (���   t=0
y j = y j0 ,   j = 1, 2,3, 4 ) ����������� ������������:

                           k1 j = h ⋅ f j ;      y j = y j 0 + 0,5 ⋅ k1 j ;
                          k2 j = h ⋅ f j ;       y j = y j 0 + 0,5 ⋅ k2 j ;
                              k3 j = h ⋅ f j ;        y j = y j 0 + k3 j ;
                                          k4 j = h ⋅ f j .

   ��� ���� ������������ ���������� �� ����� ������� ���
j = 1, 2,3, 4 , � ���� ����� ���������� ������� � ���������.
   2. ������������ ��������� ������� ������� � ������������ ��
�������� ����� ����� h:

               y j = y j 0 + (k1 j + 2k2 j + 2k3 j + k4 j ) / 6 , j = 1, 2,3, 4.

   3. ���������������� ��������� ������� ��� ���������� ���� ��
�������:

                                  y j0 = y j ,         j = 1, 2,3, 4

   � �������������� ������� � ������ 1. ���������� ������������ ��
��� ���, ���� ������� ������� �� ��������� ��������� �������.
   ��� ������� ���������� ������� ������� � ������� � ���������-
���� �� ������� ���������� ��������� ���������� � ������ ������
������� ��� ����������� �������� ������� ������������� ��������-
����� ���������� (������), �������������� �� ������� �������.


1.2. ������ ������������� �������� ������� ������ � �������
������ �� ������� � ���������� ��������

     ����������� ������� (���.4) ������������ ��� ����������� ����-
��������� ���� �������� ��������� �� �������������������� �����-
���.