Метод дискретного элемента

Метод дискретного элемента (DEM, от англ. Discrete element method) — это семейство численных методов предназначенных для расчёта движения большого количества частиц, таких как молекулы, песчинки, гравий, галька и прочих гранулированных сред. Метод был первоначально применён Cundall в 1971 для решения задач механики горных пород. Williams, Hocking и Mustoe детализировали теоретические основа метода. В 1985 они показали, что DEM может быть рассмотрен как обобщение метода конечных элементов (МКЭ, FEM). В книге Numerical Modeling in Rock Mechanics, by Pande, G., Beer, G. and Williams, J.R. описано применение этого метода для решения геомеханических задач. Теоретические основы метода и возможности его применения неоднократно рассматривалось на 1-й, 2-й и 3-й Международной Конференции по Методам Дискретного Элемента. Williams, и Bicanic (см. ниже) опубликовали ряд журнальных статей описывающих современные тенденции в области DEM. В книге The Combined Finite-Discrete Element Method, Munjiza детально описано комбинирование Метода Конечного Элемента и Метода Дискретного Элемента.

Этот метод иногда называют молекулярной динамикой (MD), даже когда частицы не являются молекулами. Однако, в противоположность молекулярной динамике, этот метод может быть использован для моделирования частиц с несферичной поверхностью. Различными ответвлениями семейства DEM являются метод отдельных элементов (distinct element method), предложенный Cundall в 1971, обобщенный метод дискретного элемента (generalized discrete element method), предложенный Williams, Hocking и Mustoe в 1985, дискретный деформационный анализ (discontinuous deformation analysis) (DDA) предложенный Shi в 1988, и метод конечных дискретных элементов (finite-discrete element method), предложенный Munjiza и Owen в 2004.

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

Применение

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

• жидкости и растворы, например сахар или белок;
• сыпучие вещества в элеваторе, такие как крупа;
• гранулированный материал, такой как песок;
• порошки, такие как тонер.

Типичные отрасли промышленности использующие DEM:

• Горнодобывающая
• Фармацевтическая
• Нефтегазовая
• Сельскохозяйственная
• Химическая

Основные принципы метода

Моделирование МДЭ начинается c помещения всех частиц в конкретное положение и придания им начальной скорости. Затем силы, воздействующие на каждую частицу, рассчитываются, исходя из начальных данных и соответствующих физических законов.

Следующие силы могут иметь влияние в макроскопических моделях:

• трение, когда две частицы касаются друг друга;
• отскакивание, когда две частицы сталкиваются;
• гравитация (сила притяжения между частицами из-за их массы), которая имеет отношение только при астрономическом моделировании.

На молекулярном уровне, мы можем рассматривать

• Силу Кулона, электростатическое притяжение или отталкивание частиц, несущих электрический заряд;
• Отталкивание Паули, когда два атома находятся вблизи друг от друга;
• Силу Ван дер Ваальса.

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

Типичные методы интеграции используемые в методе дискретного элемента:

• алгоритм Верлета,
• скорость Верлета,
• метод прыжка.

Дальнодействующие силы

Когда во внимание принимаются дальнодействующие силы (гравитация, сила Кулона), взаимодействия каждой пары частиц необходимо рассчитывать. Число взаимодействий, а следовательно, ресурсоёмкость расчёта, возрастает с увеличением количества частиц квадратично, что не приемлемо для моделей с большим числом частиц. Возможный путь решить эту проблему — объединить некоторые частицы, которые находятся на расстоянии от рассматриваемой частицы, в одну псевдочастицу. Рассмотрим, например, взаимодействие между звездой и отдаленнной галактикой: ошибка, возникающая из-за объединения массы всех звезд в удалённой галактике в одну точку, незначительна. Для того, чтобы определить, какие частицы могут быть объединены в одну псевдочастицу, используются так называемые древесные алгоритмы. Эти алгоритмы распределяют все частицы в виде дерева, квадрадерева в случае двухмерной модели и октадерева в случае трехмерной модели.

Модели в молекулярной динамике делят пространство, в котором происходит моделируемый процесс, на ячейки. Частицы, уходящие через одну сторону ячейки просто вставляются с другой стороны (периодические граничные условия); так же происходит и с силами. Силы перестают приниматься в расчёт после так называемой дистанции отсечения (обычно половина длины ячейки), так что на частицу не воздействует зеркальное расположение той же частицы на другой стороне ячейки. Таким образом, можно увеличивать количество частиц простым копированием ячеек.

Алгоритмы для обработки долговременных сил:

• Barnes-Hut,
• быстрый мультипольный метод.

Библиография

• P.A. Cundall, O.D.L. Strack, A distinct element model for granular assemblies. Geotechnique, 29:47—65, 1979.
• Williams, J.R., Hocking, G., and Mustoe, G.G.W., "The Theoretical Basis of the Discrete Element Method, " NUMETA 1985, Numerical Methods of Engineering, Theory and Applications, A.A. Balkema, Rotterdam, January 1985
• Shi, G, Discontinuous deformation analysis — A new numerical model for the statics and dynamics of deformable block structures, 16pp. In 1st U.S. Conf. on Discrete Element Methods, Golden. CSM Press: Golden, CO, 1989.
• Williams, J.R. and Pentland, A.P., "Superquadric and Modal Dynamics for Discrete Elements in Concurrent Design, " National Science Foundation Sponsored 1st U.S. Conference of Discrete Element Methods, Golden, CO, October 19-20, 1989.
• Pande, G., Beer, G. and Williams, J.R., Numerical Modeling in Rock Mechanics, John Wiley and Sons, 1990.
• Kawaguchi, T., Tanaka, T. and Tsuji, Y., Numerical simulation of two-dimensional fluidized beds using the discrete element method (comparison between the two- and three-dimensional models) Powder Technology, 96(2):129—138, 1998.
• Griebel, Knapek, Zumbusch, Caglar: Numerische Simulation in der Molekulardynamik. Springer, 2004. ISBN 3-540-41856-3.
• Bicanic, Ninad, Discrete Element Methods in Stein, de Borst, Hughes Encyclopedia of Computational Mechanics, Vol. 1. Wiley, 2004. ISBN 0-470-84699-2.
• 2nd International Conference on Discrete Element Methods, Editors Williams, J.R. and Mustoe, G.G.W., IESL Press, 1992 ISBN 0-918062-88-8
• Williams, J.R. and O’Connor, R., Discrete Element Simulation and the Contact Problem, Archives of Computational Methods in Engineering, Vol. 6, 4, 279—304, 1999
• Ante Munjiza, The Combined Finite-Discrete Element Method Wiley, 2004, ISBN 0-470-84199-0

Программное обеспечение

Открытые источники и некоммерческое программное обеспечение:

• BALL & TRUBAL (1979—1980) distinct element method (FORTRAN code), originally written by P.Cundall and currently maintained by C.Thornton.
• SDEC Spherical Discrete Element Code.
• YADE Yet Another Dynamic Engine, second incarnation of SDEC written from ground-up, GPL license.
• LIGGGHTS Open Source Discrete Element Method Particle Simulation Code, распространяется и поддерживается Кжиштофом Клоссом

Доступные за дополнительную плату пакеты DEM-программ, включая PFC3D,EDEM и Passage/DEM:

• Chute Maven (Hustrulid Technologies Inc.)
• PFC2D и PFC3D (Particle Flow Code in 2 Dimensions; Particle Flow Code in 3 Dimensions), PFC2D использует исходный код BALL, PFC3D использует исходный код TRUBAL.
• EDEM (DEM Solutions Ltd.)
• GROMOS 96
• ELFEN
• MIMES
• PASSAGE®/DEM (PASSAGE®/DEM Программа для расчёта частиц под воздействие разнообразных сил)
• Bulk Flow Analyst (Bulk Flow Analyst Программа для расчёта сыпучих материалов общего назначения)

См. также

• Метод конечных разностей
• Метод конечных элементов
• Метод подвижных клеточных автоматов