WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!

Pages:     | 1 |   ...   | 8 | 9 || 11 | 12 |   ...   | 27 |

jj Qij,0 = ln r0, uij,0 = ui0, j = 0,..., M, i =1,2, ( ) ( ) i В настоящей работе проведена реализация сформулированной выше конечно-разностной схемы.

Литература 1. Нигматулин Р. И. Динамика многофазных сред. – М.: Наука, 1987.

2. Rajagopal K. R., Tao L. Mechanics of mixtures. – London: World Scientific Publishing, 1995.

3. Кучер Н. А. Метод слабой аппроксимации и анализ схем расщепления в газовой динамике. – Кемерово, 1997.

4. Самарский А. А. Теория разностных схем. – М.: Наука, 1977.

БЕССЕТОЧНЫЙ МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ ДЛЯ МОДЕЛИРОВАНИЯ ВЗАИМОДЕЙСТВИЯ ГРАВИТАЦИОННЫХ ТЕЧЕНИЙ С ВЕРТИКАЛЬНОЙ ПРЕГРАДОЙ Т. С. Рейн, В. Н. Клепче Кемеровский государственный университет Распространение подводных гравитационных течений и их взаимодействие с подводными сооружениями и трубопроводами представляет собой актуальную задачу газо- и нефтедобычи. Исследование данной тематики носит фундаментальный характер. Современные газо- и нефтедобывающие комплексы зачастую располагаются в шельфовой зоне на значительном удалении от берега, вследствие чего подводные трубопроводы и линии коммуникаций подвержены опасному влиянию гравитационных и взвесенесущих течений. Подобные течения формируются, когда более «тяжелая» жидкость (преимущественно в горизонтальном направлении) распространяется внутрь более легкой жидкости. Последние экспериментальные исследования показали, что проблемы взаимодействия гравитационных течений с подводными структурами принципиально отличаются от хорошо изученных задач, в которых течения с постоянной плотностью взаимодействуют с твердыми телами.

В настоящей работе представлено исследование взаимодействия гравитационного (плотностного) течения с вертикальными препятствиями.

Подобные задачи имеют большую актуальность в приложениях. Примерами таких течений являются: снежные лавины, сели, пирокластические потоки, возникающие при извержениях вулканов, подводные переносящие взвеси течения, возникающие при перемещениях осадочных пород на наклонном дне. Многочисленные примеры техногенных катастроф, вызванных плотностными течениями, приведены в монографии [1].

Можно отметить, что задаче о накате гравитационного потока на прибрежные и донные конструкции посвящено сравнительно небольшое количество фундаментальных исследований. В работе [2] проводится теоретическое исследование перелива жидкости через препятствия; численное моделирование взаимодействия гравитационного течения с прямоугольным препятствием на дне канала выполнено авторами [3], там же проведено сопоставление с результатами [4]. Экспериментальные результаты сформировали некоторое представление о качественном характере эффектов, наблюдаемых при взаимодействии гравитационного течения с донными преградами (в частности, круговыми цилиндрами, моделирующими трубопроводы, расположенные на морском дне), но их сложно использовать для количественного сопоставления с расчетами. Кроме того, численное моделирование выявило наличие ряда новых динамических эффектов, для исследования которых необходимо проведение дополнительных расчетов.

В силу значительных деформаций свободных границ и границ раздела, а также возможного взаимного перемешивания слоев жидкости применение для обозначенных классов задач классических методов математического моделирования, таких как методов конечных и граничных элементов, не представляется возможным. Появление новых численных методик, базирующихся на способах построения функций формы, не нуждающихся в информации об узловой связности (сетке, покрывающей расчетную область), позволило как усовершенствовать классические методы математического моделирования, так и построить ряд новых методов, объединенных названием бессеточные. В данной работе для моделирования процессов взаимодействия жидкости с преградами, сопровождающихся нарушением связности расчетной области, используется бессеточный метод конечных элементов (Meshless Finite Element Method), который был адаптирован для численного моделирования динамики несжимаемой вязкой жидкости, включающего этапы сильно нелинейного поведения течения с образованием многосвязной расчетной области [5]. Характерной чертой метода является возможность с высокой степенью точности численно определять поле давления, которое необходимо для расчета гидродинамических нагрузок. Для устранения нефизических осцилляций функции давления в работе используется метод конечных приращений. Для аппроксимации функции давления используются линейные базисные функции (функции формы расширенной интерполяции Лапласа), а для аппроксимации функции скорости – квадратичные базисные функции (функции формы Сибсона) [6].

Выпишем постановку задачи. Пусть в некоторой области течения D происходит движение ньютоновской вязкой несжимаемой жидкости, описываемое системой уравнений Навье-Стокса (1)–(2).

Dui ¶¶ r =- p + tij +rfi, (1) Dt ¶xi ¶x j ¶ui = 0. (2) ¶xi Так как жидкость вязкая, то на твердых стенках выполняется условие прилипания:

% ui = ui, i =1,2. (3) На свободной границе заданы плотности поверхностных сил % pn =tijnj - pni =sni, (4) здесь ni, nj – компоненты вектора внешней нормали к свободной границе.

Для нестационарной задачи необходимо также задать начальное положение свободной границы и распределение неизвестных функций во всей области течения:

u(x,0) = u0(x), p(x,0) = p0(x). (5) В работе рассматривается задача об обрушении столба «тяжелой» жидкости с большей плотностью в горизонтальном направлении внутрь более легкой жидкости с меньшей плотностью (рис. 1). Представлено сравнение полученных значений гидродинамических нагрузок на вертикальные стенки расчетной области с экспериментальными данными.

Рис. 1. Расчетная область Литература 1. Simpson J. E. Gravity currents: In the environment and the Laboratory. – Cambridge University Press, 1997.

2. Lane-Serff G. F., Beal L. M., Hadfield T. D. Gravity current flow over obstacles // J. Fluid Mech. – 1995. – V. 292. – P. 39–53.

3. Gonzales-Juez E., Meiburg E., Constantinescu G. Gravity currents impinging on bottom-mounted square cylinders: flow fields and associated forces // J. Fluid mech. – 2009. – V. 631. – P. 65–102.

4. Ermanyuk E. V., Gavrilov N. V. On internal waves generated by large-amplitude circular and rectilinear oscillations of a circular cylinder in a uniformly stratified fluid // J.

Fluid Mech. – 2008. – V. 613. – P. 329–356.

4. Onate E. A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation // Comput. Meth. Appl. Mech. Eng. – 2000. – Vol.

182, №1–2. – Р. 355–370.

5. Шокин, Ю. И. Метод дифференциального приближения. – Новосибирск:

Наука, 1979. – 224 с.

5. Афанасьев К. Е., Рейн Т. С. Моделирование задач гидродинамики вязкой несжимаемой жидкости со свободными границами бессеточным методом естественных соседей // Вычислительные технологии. – 2008. – Т. 13, № 4. – С. 7–24.

6. Беликов В. В., Иванов В. Д., Конторович В. К., Корытник С. А., Семенов А. Ю. Несибсоновская интерполяция – новый метод интерполяции значений функции на произвольной системе точек // Вычислительная математика и математическая физика. – 1997. – Т. 37, №1. – С. 11–17.

ИССЛЕДОВАНИЕ СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ MMP(2) GI МЕТОДОМ ПРОСЕЯННОГО ПОТОКА Е. С. Ронжина, И. А. Ивановская Филиал Кемеровского государственного университета в г. Анжеро-Судженске Пусть на вход системы обслуживания поступает MMP -поток [2] сдвоенных заявок, то есть в момент наступления события в рассматриваемом потоке в систему одновременно поступают две заявки. Будем считать, что MMP -поток определяется эргодической цепью Маркова k(t), заданной матрицей Q – её инфинитезимальных характеристик qk k2, набором неотрицательных величин l 0.

k Дисциплина обслуживания определяется тем, что одна из этих заявок поступает в первый, а другая вторая – во второй блоки обслуживания и занимает любой из свободных приборов, на котором выполняется ее обслуживание. Продолжительности обслуживания различных заявок стохастически независимы и одинаково распределены для всех приборов соответствующего блока. Пусть B1 x, B2 x – функции распределения време( ) ( ) ни обслуживания для первого и второго блоков соответственно.

Обозначим i1(t),i2(t)– число приборов, занятых в момент времени t в первом или втором, соответственно, блоке обслуживания, а стационарное распределение вероятностей значений процесса {i1(t),i2(t)} обозначим П(i1,i2)= P{i1(t)= i1,i2(t)= i2}.

Для рассматриваемой системы ни двумерный случайный процесс {i1(t),i2(t)} изменения во времени состояний системы, ни трехмерный случайный процесс {k(t),i1(t),i2(t)} не являются марковскими.

Для исследования рассматриваемой системы MMP(2) GI применим метод просеянного потока [4]. Предлагаемый метод позволяет проблему исследования немарковской системы обслуживания с неограниченным числом приборов свести к задаче анализа нестационарного марковизируемого потока.

На оси времени t отметим моменты наступления событий этого потока. Выделим некоторый момент времени t1. Пусть t1 = 0. Будем полагать, что заявки входящего потока, поступившие в систему в момент времени t < t1 = 0, формируют события двумерного просеянного потока с вероятностями S1(t) = 1 - B1(-t), S2 (t) =1 - B2 (-t), а с вероятностью 1 - S1(t) и 1 - S2 (t) не рассматриваются.

Обозначим {n1(t),n2 (t)} – двумерный процесс, компоненты которого характеризуют число событий просеянных потоков, наступивших до момента времени t.

Если в некоторый начальный момент времени t0 < t1 система обслуживания свободна, то есть в ней нет обслуживаемых заявок, то для момента времени выполняется равенство t i1(t1) = n1(t1), i2(t1)= n2(t1), (1) то есть число ik (t1) приборов, занятых в k -м блоке обслуживания рассматриваемой системе обслуживания, равно числу nk (t1) событий просеянного потока, наступивших до момента времени t1.

Для рассматриваемой системы трехмерный случайный процесс {k(t),n1(t),n2(t)} является нестационарной двумерной цепью Маркова, для исследования которой естественно воспользоваться теорией марковских цепей. Здесь k(t) – значения цепи Маркова, управляющей входящим MMP -потоком, а n(t) – число событий просеянного потока, наступивших до момента времени t.

Для распределения вероятностей P(k,n1,n2,t) = P{k(t) = k,n1(t) = n1,n2(t) = n2} по формуле полной вероятности запишем следующее равенство:

P(k,n1, n2,t + Dt)= P(k, n1,n2,t)(1- lkDt)(1+ qkkDt)+ P(k,n1 -1,n2 -1,t)lkDtS1(t) S2(t)+ P(k, n1 -1, n2,t)lkDtS1(t)(1- S2(t))+ P(k, n1,n2 -1,t)lkDtS2(t)(1- S1(t))+ + P(k,n1, n2,t)lkDt(1- S1(t))(1- S2(t))+, n2,t)qnkDt + o(Dt).

P(n,nn Далее запишем систему дифференциальных уравнений Колмогорова [3] ¶P(k, n1,n2,t) = (- lk + qkk )P(k,n1,n2,t)+ lkS1(t)S2(t)P(k,n1 -1,n2 -1,t)+ ¶t + [P(k,n1 -1,n2,t)S1(t)(1 - S2 (t))+ P(k,n1,n2 -1,t)(1 - S1(t))S2 (t)]lk + + lk (1- S1(t))(1 - S2(t))P(k,n1,n2,t)+,n2,t)qnk. (2) P(n,nnk Начальное условие для решения P(k,n1,n2,t) в момент времени tопределим равенством R(k), если n1,n2 = 0, P(k,n1,n2,t0 ) = 0, если n1,n2 > 0.

Обозначив jun1 2 jun1(t ) (t ) H(k,u, w,t)= e eiwn P(k,n1, n2,t)= R(k)M{e eiwn k(t) = k} (3) n1 =0n2 =H (k,u, w,t0 ) = R(k), получаем следующую задачу Коши:

¶H(k,u,w,t) ju jw = lk[S1(t)(e -1) + S2 (t)(e -1) + ¶t ju jw + (e -1)(e -1)S1(t)S2 (t)] H(k,u, w,t)+ H(k,u, w,t)qnk. (4) n Систему (4) запишем в матричной форме, обозначив вектор-строку H (u,w,t) = {H (1,u, w,t,), H (2,u, w,t),L}, получим ¶H(u, w,t) = ¶t ju jw ju jw = H(u, w,t)[Q +(S1(t)(e -1) + S2(t)(e -1) + (e -1)(e -1)S1(t)S2(t))]L, H (u,w,t0 ) = R. (5) Решив эту задачу, найдём функцию H (u, w,t). Тогда, принимая во внимание равенство (3), характеристическую функцию величины n(t) определим в виде jun1(t ) jwn2 (t ) Me e = H (u,t)E. (6) Полагая t = t1 = 0, в силу равенства (1), характеристическая функция стационарного распределения вероятностей П(i1,i2 ) числа приборов, занятых в каждом блоке системы MMP(2) GI, будет иметь следующий вид:

jun1(t ) jwn2 (t ) h(u) = Me e = H (u,0)E, (7) что позволяет найти все вероятностно-временные характеристики системы обслуживания MMP(2) GI.

Уравнение (4), определяющее характеристики системы MMP(2) GI, предлагается решать в асимптотическом условии растущего времени, полагая, что среднее значение времени обслуживания b ®.

Литература 1. Гнеденко Б. В., Коваленко И. Н. Введение в теорию массового обслуживания. –3-е изд., испр. и доп. – М.: КомКнига, 2005. – 408 с.

2. Назаров А. А., Терпугов А. Ф. Теория массового обслуживания. – Томск:

Изд-во НТЛ, 2005. – 228 с.

3. Назаров А. А., А.Ф. Терпугов А. Ф. Теория вероятностей и случайных процессов. Томск: Изд-во НТЛ, 2006. – 204 с.

4. Назаров А. А., Моисеева С. П. Метод асимптотического анализа в теории массового обслуживания. – Томск: Изд-во НТЛ, 2006. – 112 с.

Работа выполнена при поддержке АВЦП «Развитие научного потенциала высшей школы (2009–2010 годы)», проект № 4761 «Разработка методов исследования немарковских систем массового обслуживания и их применение к сложным экономическим системам и компьютерным сетям связи».

ИССЛЕДОВАНИЕ n-ФАЗНОЙ СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ C НЕОГРАНИЧЕННЫМ ЧИСЛОМ ОБСЛУЖИВАЮЩИХ ПРИБОРОВ И ПРОИЗВОЛЬНЫМ РАСПРЕДЕЛЕНИЕМ ВРЕМЕНИ ОБСЛУЖИВАНИЯ М. В. Санникова, И. Р. Гарайшина Филиал Кемеровского государственного университета в г. Анжеро-Судженске Рассмотрим n-фазную систему массового обслуживания (СМО) с неограниченным числом обслуживающих приборов [1], на вход которой поступает простейший поток заявок с параметром l. Поступившая заявка занимает одну из свободных линий, последовательно выполняя обслуживание, начиная с первой фазы. Линия считается занятой, если занята одна из её фаз. Завершив обслуживание на k -й фазе k =1,n -1, заявка с веро( ) ятностью 1- rk покидает систему, а с вероятностью rk переходит на следующую фазу для продолжения обслуживания. Завершив обслуживание на п-ой фазе, заявка покидает систему. Продолжительности фаз обслуживания стохастически независимы друг от друга.

Считаем, что время обслуживания на k -м приборе k =1,n имеет ( ) произвольную функцию распределения Bk (x).

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

Разделим входящий поток на N независимых простейших потоков с l параметром. Заявки каждого потока направляются для обслуживания в N СМО на свободную линию. Таким образом, получаем совокупность N однолинейных n -фазных систем массового обслуживания, на вход которых l поступает простейший поток с параметром. Заявка, заставшая линию N занятой, теряется, но при N ® вероятность отказа в постановке заявки на обслуживание стремится к нулю.

Введем следующие обозначения:

Pages:     | 1 |   ...   | 8 | 9 || 11 | 12 |   ...   | 27 |



© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.