CC BY
15
УДК 621.321
МОДЕЛЬ ВЯЗКОГО ТРЕНИЯ В МАГНИТОЖИДКОСТНОМ ГЕРМЕТИЗАТОРЕ
НА ОСНОВЕ МЕТОДА МОНТЕ-КАРЛО
А.И. ТИХОНОВ1, Ю.Б. КАЗАКОВ1, К.М. ПИРОГОВ2, А.С. МАРТЫНОВ1 1ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
Иваново, Российская Федерация E-mail: ait@dsn.ru
2ФГБОУВПО «Ивановский государственный университет», Иваново, Российская Федерация
E-mail: piroqov K.M.@.ru
Авторское резюме
Состояние вопроса: Основным недостатком аналитических и численных методов расчета процессов вязкого трения является сложность математического аппарата. В настоящее время большую популярность при решении подобных задач приобретает метод Монте-Карло, отличающийся простотой в алгоритмизации, в том числе, с учетом реализации этих алгоритмов в технологии параллельных вычислений.
Материалы и методы: Для расчета положения магнитной жидкости использовался метод конечных элементов. Для расчета поля скоростей жидкости использовался метод Монте-Карло.
Результаты: Разработан метод расчета момента вязкого трения магнитожидкостного герметизатора. Сделаны выводы о перспективах использования метода Монте-Карло при решении задач гидродинамики. Выводы: Использование комбинированного метода позволяет существенно упростить математический аппарат расчета вязкого трения.
Ключевые слова: магнитная жидкость, метод конечных элементов, метод Монте-Карло, вязкое трение, магни-тожидкостный герметизатор.
MODEL OF VISCOUS FRICTION IN MAGNETIC LIQUID SEALER ON THE BASIS
OF MONTE CARLO METHOD
A.I. TIKHONOV1, YU.B. KAZAKOV1, K.M. PIROQOV2, AS. MARTYNOV1 Ivanovo State Power Engineering University, Ivanovo, Russian Federation
E-mail: ait@dsn.ru 2Ivanovo State University, Ivanovo, Russian Federation E-mail: piroqov K.M.@.ru
Abstract
Background: The main drawback of analytical and numerical methods for calculating the processes of viscous friction is the complexity of the mathematical apparatus. At present, the Monte Carlo method becomes more and more popular in solving such problems. It is simple in algorithmization, including the implementation of these algorithms in parallel computing technology.
Materials and methods: For the calculation of the magnetic fluid finite element method is used. The Monte Carlo method is also applied when the velocity field of fluid is calculated.
Results: The method for calculating the moment of viscous friction of magnetic liquid stealer is developed. The conclusions about the prospects of using the Monte Carlo method for solving hydrodynamics problems are made. Conclusions: The usage of the combined method allows substantially to simplify the mathematical apparatus of viscous friction calculation.
Key words: magnetic fluid, finite element method, Monte Carlo method, viscous friction, magnetic liquid sealer.
Традиционно при решении задач моделирования физических процессов предпочтение отдается аналитическим методам, позволяющим достаточно быстро получить точное решение. Частным случаем данного подхода является использование численных методов моделирования, сводящих задачу к численному решению систем уравнений, построенных на основе аналитических методов. Основным недостатком этих методов моделирования является сложность математического аппарата, которая экспоненциально возрастает, например, по мере уменьше-
ния количества допущений, при возрастании мерности задачи и т.п. В частности, одними из математически наиболее сложных считаются задачи гидродинамики с учетом вязкости моделируемой жидкости, к числу которых относится и расчет момента вязкого трения магнитожидкостного герметизатора (МЖГ).
Вместе с тем в математике существует направление, связанное с методом Монте-Карло, которое основано на использовании теории вероятностей для решения различных задач. Из физики известно, что любому закону природы
можно дать статистическую интерпретацию, особенно если исследуемый процесс можно интерпретировать с позиций молекулярно-кинетической модели вещества. Главное преимущество метода Монте-Карло состоит в простоте математического аппарата, который базируется на понятных физических представлениях о природе моделируемых процессов.
Основной недостаток недетерминированных методов, в частности метода Монте-Карло, состоит в их медленной сходимости, что предъявляет повышенные требования к компьютерным ресурсам. Однако их простота и наглядность обусловливают возрастающий к ним интерес. Особую актуальность эти методы приобретают в связи с развитием многопроцессорной вычислительной техники, так как недетерминированные методы легко распараллеливаются, в отличие о детерминированных, требующих, как правило, последовательных вычислений. Так как развитие компьютерных технологий в настоящее время напрямую связывается с многопроцессорными вычислительными системами, то можно утверждать, что метод Монте-Карло является одним из наиболее перспективных.
Для примера рассмотрим версию метода Монте-Карло для моделирования процессов вязкого трения в МЖГ. Исходными данными для расчета является параметрически формируемая конечно-элементная модель МЖГ (рис. 1), позволяющая рассчитать магнитное поле с учетом положения магнитной жидкости, а также само это положение с учетом заданного перепада давления (рис. 2) [1].
вторяет логику поведения реальных структурных элементов (частиц) моделируемой жидкости, подчиняющуюся статистическим законам.
Рис. 1. Конечно-элементная модель магнитожидкостного герметизатора
Расчет методом Монте-Карло реализуется в форме имитационного эксперимента, который строится на совокупности правил, регламентирующих поведение отдельных элементов модели во времени. Это поведение в общих чертах по-
Рис. 2. Положение магнитной жидкости в МЖГ при заданном перепаде давления, рассчитанное с помощью системы конечно-элементного моделирования магнитного поля
Примем ряд допущений:
1. Расчетная область обладает осевой симметрией, что позволяет свести задачу к плоскомеридианной постановке.
2. Магнитная жидкость является несжимаемой ньютоновской жидкостью с заданным динамическим коэффициентом вязкости, величина которого является константой.
3. Положение магнитной жидкости не зависит от скорости вращения вала и однозначно определяется картиной магнитного поля (рис. 2).
4. Скорость вращения вала постоянна.
5. Скорость слоев магнитной жидкости, непосредственно примыкающих к твердым поверхностям, равна линейной скорости движения данной поверхности, которая не достигает критических значений.
6. В течении жидкости отсутствуют вихревые составляющие.
Имитационная модель магнитной жидкости строится из множества частиц, обладающих одинаковой массой m0, совершающих броуновское движение с одинаковыми скоростями v0 в пределах подобласти плоской конечно-элементной модели, занятой магнитной жидкостью (рис. 2). Изначально каждая i-я частица обладает нулевым моментом импульса относительно оси вращения вала, но при соударении с j-й границей раздела сред она увлекается валом, приобретая момент импульса относительно оси вращения
Loi = mo02<» j, (1)
где r - расстояние от i-й частицы до оси вращения; - угловая скорость j-й границы относительно оси вращения.
При свободном движении частицы ее момент импульса относительно оси вращения вала сохраняется. Изменение этого момента импульса может происходить только при взаимодействии частицы с твердой поверхностью, что никак не
отражается на скорости частицы в осевой плоскости модели.
Так как массы частиц одинаковы, то при столкновении двух частиц, они обмениваются импульсами и моментами импульса. Это значит, что траектории двух столкнувшихся частиц представляют собой две пересекающиеся прямые, движение вдоль которых осуществляется с постоянной скоростью. Так как частицы не отличимы друг от друга, то с учетом сказанного сам факт взаимодействия частиц можно в учет не принимать, считая их движение броуновским.
Численный эксперимент состоит в имитации броуновского движения множества частиц, численностью от одного до нескольких миллионов. Частицы переносят момент импульса по плоскости расчетной области в пределах подобласти, заполненной магнитной жидкостью. После определенного времени релаксации система частиц придет в состояние термодинамического равновесия, характеризующегося неизменным во времени значением среднего момента импульса относительно оси вращения вала в каждой элементарной области, занятой магнитной жидкостью. При этом распределение величины среднего момента импульса будет соответствовать уравнению Лапласа с граничными условиями Дирихле на поверхностях твердых тел (-0 = -0;) и граничными условиями Неймана на свободной границе жидкости.
Плоскость модели покрыта конечно-элементной сеткой. Каждый треугольный элемент является сечением кольца, ось которого совпадает с осью вращения вала. Это кольцо с определенным приближением можно считать обособленной струей, угловая скорость вращения которой определяется средним моментом импульса всех частиц в пределах элемента
ык
X X Гк ®к/ X - 'к
¡=1
¡=1
¡=1
1Чкт0гк
(2)
где -й, гк, и юк,- - момент импульса, расстояние до оси и угловая скорость ¡-й частицы, находящейся в пределах к-го треугольного элемента, относительно оси вращения вала; N1,: — количество частиц в к-м треугольном элементе; гк — расстояние от центра тяжести к-го элемента до оси вращения.
В выражении (2) масса частицы т0 сократилась, то есть эту величину из модели можно исключить, оперируя величиной приведенного момента импульса частицы, связанного с ее моментом импульса относительно оси вращения вала отношением
-01
-
0¡ = Г0/ю0/ - "
т0
(3)
Согласно закону Ньютона, сила вязкого трения между слоями жидкости, приходящаяся
на единицу площади поверхности раздела слоев, находится как
1 = -Ч-У, (4)
дх
где п — динамическая вязкость жидкости; координатные оси х и у перпендикулярны друг другу, струя движется в направлении оси х.
Соответственно, момент вязкого трения между двумя слоями жидкости с площадью контакта дв относительно оси вращения вала МЖГ для осесимметричной задачи рассчитывается по формуле
дм = -цг2 — дв. (5)
дг
При этом на каждую кольцевую струю, соответствующую к-му треугольному элементу конечно-элементной сетки, действует момент трения
¡=3/
мк =-2плГк3 X
¡=1 V Лгм
(6)
где гк — расстояние от центра тяжести к-го элемента до оси вращения вала; Лши — разность угловых скоростей к-й струи и струи, примыкающей к ¡-й стороне к-го треугольного элемента; Лгк— расстояние между центрами тяжести к-го элемента и элемента, примыкающего к -й стороне к-го элемента; Iк— длина ¡-й стороны к-го треугольного элемента.
Общий момент трения находится путем суммирования моментов трения всех струй.
Таким образом, задача расчета момента вязкого трения МЖГ решена в три этапа: определение положения магнитной жидкости на основе расчета магнитного поля методом конечных элементов; расчет распределения угловых скоростей элементарных струй методом Монте-Карло; и наконец, расчет моментов трения каждой элементарной струи на конечно-элементной сетке на основе конечно-разностной аппроксимации закона Ньютона для вязкого трения. Задача решена с учетом ряда допущений. Однако главная особенность метода Монте-Карло состоит в простоте модификации алгоритмов расчета. Перспектива развития метода Моте-Карло в плане решаемой задачи состоит в обосновании и программной реализации таких правил поведения частиц, которые смогут снять принятые допущения. Разрабатывается алгоритм, позволяющий распараллелить расчетный процесс в технологии СЫРА.
Список литературы
1. Система расчета магнитожидкостного герметизатора с параметрически формируемой полевой моделью и определением критического положения магнитной жидкости / Ю.Я. Щелыкалов, Ю.Б. Казаков, А.И. Тихонов, С.М. Перминов: сб. науч. тр. Х Междунар. Плесской конф. по магнитным жидкостям. Плес. 9—12 сентября 2002. — Иваново, 2002. — С. 374—381.
References
1. Shchelykalov, Yu.Ya, Kazakov, Yu.B., Tikhonov, A.I., Perminov, S.M. Sistema rascheta magnitozhidkostnogo germeti-zatora s parametricheski formiruemoy polevoy model'yu i opre-
deleniem kriticheskogo polozheniya magnitnoy zhidkosti [Calculation System of Magnetic Liquid Sealer with the Parameter - depended Field Model and Determination of Critical Position of Magnetic Liquid], in \ Mezhdunarodnaya konferentsiya po mag-nitnym zhidkostyam. Ples, Rossiya [The Xth International Ples Conference on Magnetic Liquids], Ivanovo, 2002, pp. 374-381.
Тихонов Андрей Ильич,
ФГБОУВПО «Ивановский государственный энергетический университет имени В. И. Ленина», доктор технических наук, доцент, заведующий кафедрой физики, e-mail: ait@dsn.ru
Казаков Юрий Борисович,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор, заведующий кафедрой электромеханики, e-mail: elmash@em.ispu.ru
Пирогов Константин Михайлович,
ФГБОУВПО «Ивановский государственный университет»,
доктор технических наук, профессор кафедры информационных технологий в экономике и организации производства,
e-mail: piroqov K.M.@.ru
Мартынов Александр Станиславович,
ФГБОУВПО «Ивановский государственный университет»,
магистрант,
e-mail: alexmart@mail.ru
