В то же время алгоритм должен хранить  историю изменений ячеек. Поэтому  рассматриваются дополнительные события, называемые событиями ячеек. Их обработка  происходит так же, как и обработка  столкновений частиц, но партнером  при столкновении считается граница  ячейки. В отличие от столкновения двух частиц, здесь только одна частица  участвует в изменениях в ячейке и в момент наступления события  скорость частицы не изменяется. Изменения  ячеек на границе области моделирования  рассматриваются следующим образом: вектор положения перепрыгивает  к противоположной стороне области  благодаря периодическим граничным  условиям.

5. Оптимальное количество ячеек

     В среднем для каждой ячейки в соседних для нее находится  частиц. Если , это количество сильно возрастает и после каждого столкновения приходится пересчитывать большое количество возможных событий. С другой стороны, если , то каждой частице приходится пересекать большое количество граничных ячеек в промежутке между столкновениями двух частиц, что приводит к большему количеству обрабатываемых событий для полного моделирования. Эти два замечания противопоставлены друг другу. Первое пропорционально , а второе - . На рис. 3 (слева) показано, что существует нечеткий минимум для времени симуляции при низкой плотности при оптимальном количестве ячеек в двумерном случае и в трехмерном. Программа может выбрать оптимальное количество ячеек до начала вычислений. Заметим, что в пределе для высоких плотностей, особенно в трехмерном случае, оптимальное количество ячеек не может быть достигнуто, так как размеры ячейки должны быть больше, чем диаметры частиц (рис. 3 справа). В этом случае должно быть выбрано как можно большим.

     Рис. 1.3 График зависимости времени симуляции от количества ячеек 

     В примере (рис. 3) использованы частиц (двумерный случай) и частиц (трехмерный случай) с относительными объемами (двумерный), (трехмерный случай на левом рисунке) и (трехмерный случай на правом рисунке).

2. МОДЕЛЬ  ТЕЧЕНИЯ ГАЗА ВСРЕДЕ С ФИЛЬТРАМИ

2.1 Описание  рабочей области

     Рассмотрим  двумерный случай течения жидкости в некоторой прямоугольной области. Определим, что вектор движения жидкости направлен слева направо. Будем  считать, что рассматриваемая область  является замкнутой, то есть молекулы жидкости, вылетая из области на правом срезе, передаются на вход (левый  срез). Появление на входе частиц осуществляется по принципу очереди  «Первый зашел, первый вышел». Границы  области будем считать отражающими. В рабочей области расположены фильтры. Каждый фильтр – это набор одинаковых прямоугольников с отражающими границами, расположенных параллельно горизонтальным границам области. Расстояния между прямоугольниками постоянны. Пример области с фильтрами показан на рис. 2.1

     Рис. 2.1 – Рабочая область с фильтрами

     Кроме основного вещества в области  также движутся включения – непроницаемые  тела формы, близкой к прямоугольной, обладающими свойствами упругости  и гибкости в той или иной степени. Эти включения перемещаются слева  направо вдоль рабочей области, так как на них действуют силы со стороны основного потока. Они могут быть остановлены при столкновении с фильтром. При большом количестве включений движение основного потока может прекратиться в силу закупорки фильтров.

     При решении задачи  используем событийно-дискретный подход. Включения будет представлять в виде совокупности окружностей, собранных  в слои. Каждый слой представляет собой  нить, составленную из окружностей, причем  расстояние между центрами соседних окружностей меньше  суммы радиусов этих сфер. Каждый слой сдвинут относительно соседних слоев так, что в целом окружности плотно упакованы, то есть их центры в начальный момент времени расположены в вершинах правильных треугольников.  При такой организации  представления включений нетрудно обеспечить разнообразие их форм с помощью регулирования количества слоев или длин нитей. Количество слоев влияет на гибкость включений. При количестве слоев, превышающем 3, можно говорить об упругости включения, мало изменяющими форму во время движения.

     Частицы жидкости при столкновении упруго отражаются. В зависимости от типа столкновения будем различать внутренние и внешние столкновения.

     Внешние отражения наблюдаются также при столкновении модельных частиц вещества с включениями, а также между отдельными молекулами включений, не являющиеся изначально соседним. Отражения для частиц включений, которые являются соседними, считаются внутренними. Но к ним применимы те же формулы, что и для внешних отражений. Полная энергия  всех модельных частиц при таком подходе сохраняется.

2.2 Инициализация  входных параметров

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

     На  этапе инициализации множество  частиц, границ ячеек и области  генерируется автоматически.

2.3 Реализация  событийного подхода

     Согласно  рассмотренному в первом разделе  алгоритму (рис. 1.1) на каждом шаге мы получаем из хипа («кучи») время ближайшего события  и номера частиц, для которых оно  должно произойти. Хип – это дерево, которое формируется из локальных времен всех частиц по следующему правилу: значение отцовского узла не превышает минимального времени до ближайшего события рассматриваемых частиц. Его инициализация происходит до первого прохождения алгоритма и в дальнейшем алгоритм лишь добавляет новые значения, убирая старые. В случае если количество модельных частиц отлично от некоторой степени двойки, то удобно дополнить дерево фиктивными значениями, заведомо большими любого из времен наступления события. В таком случае они не повлияют на результат.

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

     Работа  алгоритма продолжается до тех пор, пока не будет оставнолена. 
 

Выводы

     Метод молекулярной динамики является одним  из доминирующих методов компьютерного  моделирования физических процессов  в мировой компьютерной науке. Дополненный  методами теоретической физики, теории устойчивости, случайных процессов  и др., метод МД является основным инструментом исследований систем с  большим числом степеней свободы. Метод  молекулярной основан на численном решении классических уравнений движения частиц в некотором выделенном объеме (ячейке). Все частицы, находящиеся в одной ячейке, взаимодействуют друг с другом посредством заданного потенциала взаимодействия. Выбор потенциала взаимодействия частиц является наиболее важным этапом построения численной модели и, как правило, вызывает немало дискуссий. Для определения микроскопических характеристик среды, как правило, используется расчет корреляционных функций.

     В данной работе рассмотрен процесс моделирования  потока жидкости в замкнутой области с фильтрами. Для большей наглядности область была представлена прямоугольной формы. Для моделирования был использован событийный подход, основанный на пересчете столкновений (событий) для модельных частиц и обработке лишь тех, для которых момент наступления событий – минимален.

Список  использованной литературы

1. Allen P., Introduction to Molecular Dynamics Simulation – Computational Soft Matter: From Synthetic Polymers To Proteins, Lecture Notes, 2004 – p. 1-24.
2. Miller S., Luding S, Event-driven molecular dynamics in parallel / S. Miller, S. Luding – Journal of Computational Physics, 193, 2004 – p. 306–316.
3. Kuksin A., Morozov I., Norman G., Stegailov V., Valuev I., Standard of molecular dynamics modeling and simulation of relaxation -  Molecular simulation, 2005. V. 31. - p. 1005-1017.
4. De Michele C., Scala A., Schilling R., Sciortino F., Molecular correlation functions for uniaxial ellipsoids in the isotropic state, - J. Phys. Chem., 124:104509, 2004.
5. Scala A, De Michele C., Voightmann Th., Event-driven Brownian dynamics for hard spheres – ArXiv Condensed Matter e-prints, 2006.
6. S. Luding, Computer Physics Communications 147, 134, 2002.
7. B. D. Lubachevsky, Int. J. Comput. Simul. 2, 372, 1992.
8. D. E. Knuth, The Art of Computer Programming - Addison-Wesley, Reading, MA, 1968.
9. Чернышев Ю.К. Событийное программирование. Применение к решению некоторых задач физики / Ю.К. Чернышев // НАУ им. Н.Е. Жуковского «ХАИ». –– Х., 2008.