Математическое моделирование комплексных процессов / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Вычисл. центр ; [редкол.: В. С. Мингалев (отв. ред.) и др.]. - Апатиты : Кольский филиал АН СССР, 1982.

Q .= r.l.T 1 1 x X I Для i -й грани ячейки, моделирующей горную породу, форму лы для вычисления потоков в случае ее контакта с такой же ячей кой, будут иметь вид- R . ( W . - N . ) 1 ' 1 г R . ( W . - N . ) U . - o f . P . 1 1 1 1 1 ri j 1 XI R . ( W . - N . ) v , - o / . * T . - f t . * P i' l i l ri 'i r . ( w . - n . ) e . - ( p . + t . ) u .*<*!.-( p . + t " . ) v.*fi. ( 6 ) iv 1 v 1 v n ri' 1 1 ' XI Xl' 1 J ) Б . ( W . - N . ) + R . * C 2 «U.‘ o ^ . +R . ( C 2 - 2 C 2 ) V.* B. n 4 1 l ' 1 О 1 1 l ' o s / i J x Б . ( W . - N . ) + R . ( C 2 - 2 C 2 ) U . V . + R . * C 2 *V.*ft. x i ' l V i v o s 7 1 1 і о 1 J 1 6,„. ( W . - N . ) + R , ( C 2 - 2 C 2 ) U /oZ .+R . ( C 2 - 2 C 2 )V * B. 4 >i 4 l i 7 i v о s 7 x i i 4 о s ' 1 J 1 T . ( W . - N . ) + 4 R . * C 2 ‘V.* o ^ . + 4 R / C 2 *V.* 3. 1 1 V * 1 0 - 1 1 * 1 s 1^1 где R . , W., N .. U . , V., P T P T Б Б ... BM., Т . — i ’ l i’ i’ i ’ ri! n ’ x i ’ x i ’ n ’ x i f i 1 параметры потока на і-й грани; R i - плотность; Wif N j - нормальные скорости движения ребра сетки и по тока; - радиальная и осевая компоненты скорости потока; - нормальные и касательные относи тельно ребра величины радиальных и осевых компонент напряжений; и ь Р Г 1 » Ѵі ^ х і » T r i , Т. гд. l i - -\f: Г :, X X I x !+ r ! I 1 ' ‘Pi* T i “ РаДиальная, осевая, азимутальная, сдвиговая компоненты напряжения; - средняя в интервале от по t, ребра; і» ~ проекции ребра 1 ,- на оси и д Чс+1 Длина х ; о^і= хі /Ц ? ~ синус и косинус углр наклона ребра К ОСИ г В плоскости ti k T : / 2 „ Формулы нахождения потоков через границу ячеек, когда обе моделируют газообразную среду или одна - газообразную, а друга* - горную породу, аналогичны приведенным в / 5 / . Шаг по времени выбирается из условия ( не г - I - + - I A и , I шіп( г . ) х г ______ ( U +C , с —и ) 4 о ’ о ' m i n ( x . ) ' V ( и + с , с - и ) Нахождение искомых величин осуществляется в последователь ности, представленной на блок-схеме (рис. 1 ). 124