stringtranslate.com

Método de elementos discretos

Un método de elementos discretos ( MDE ), también llamado método de elementos distintos , pertenece a una familia de métodos numéricos para calcular el movimiento y el efecto de una gran cantidad de partículas pequeñas. Aunque DEM está muy estrechamente relacionado con la dinámica molecular , el método generalmente se distingue por su inclusión de grados de libertad de rotación , así como contacto con estado y, a menudo, geometrías complicadas (incluidos los poliedros). Con los avances en la potencia informática y los algoritmos numéricos para la clasificación del vecino más cercano, se ha hecho posible simular numéricamente millones de partículas en un solo procesador. Hoy en día, DEM está ganando amplia aceptación como un método eficaz para abordar problemas de ingeniería en materiales granulares y discontinuos, especialmente en flujos granulares, mecánica de polvos y mecánica de rocas. DEM se ha ampliado al método de elementos discretos extendido teniendo en cuenta la transferencia de calor , [1] la reacción química [2] y el acoplamiento a CFD [3] y FEM [4] .

Los métodos de elementos discretos son relativamente intensivos en términos computacionales, lo que limita la duración de una simulación o el número de partículas. Varios códigos DEM, al igual que los códigos de dinámica molecular, aprovechan las capacidades de procesamiento paralelo (sistemas compartidos o distribuidos) para aumentar el número de partículas o la longitud de la simulación. Una alternativa a tratar todas las partículas por separado es promediar la física de muchas partículas y tratar así el material como un continuo . En el caso del comportamiento granular similar a un sólido , como en la mecánica de suelos , el enfoque continuo generalmente trata el material como elástico o elastoplástico y lo modela con el método de elementos finitos o un método sin malla . En el caso de un flujo granular similar a un líquido o un gas, el enfoque continuo puede tratar el material como un fluido y utilizar dinámica de fluidos computacional . Sin embargo, los inconvenientes de la homogeneización de la física a escala granular están bien documentados y deben considerarse cuidadosamente antes de intentar utilizar un enfoque continuo.

La familia DEM

Las diversas ramas de la familia DEM son el método de elementos distintos propuesto por Peter A. Cundall y Otto DL Strack en 1979, [5] el método de elementos discretos generalizado (Williams, Hocking & Mustoe 1985), el análisis de deformación discontinua (DDA) ( Shi 1992) y el método de elementos finitos discretos desarrollado simultáneamente por varios grupos (p. ej., Munjiza y Owen). El método general fue desarrollado originalmente por Cundall en 1971 para problemas de mecánica de rocas. Williams, Hocking y Mustoe (1985) demostraron que DEM podría verse como un método generalizado de elementos finitos. Su aplicación a problemas de geomecánica se describe en el libro Numerical Methods in Rock Mechanics (Williams, Pande & Beer 1990). Las 1ª, 2ª y 3ª Conferencias Internacionales sobre Métodos de Elementos Discretos han sido un punto común para que los investigadores publiquen avances en el método y sus aplicaciones. Williams, Bicanic y Bobet et al. han publicado artículos de revistas que revisan el estado de la técnica. (vea abajo). El libro The Combined Finite-Discrete Element Method contiene un tratamiento completo del método combinado de elementos finitos y elementos discretos . [6]

Simulación de elementos discretos con partículas dispuestas según una fotografía de Peter A. Cundall . Como proponen Cundall y Strack (1979), los granos interactúan con fuerzas elásticas lineales y fricción de Coulomb. La cinemática de los granos evoluciona a través del tiempo mediante la integración temporal de su equilibrio de fuerza y ​​torque. El comportamiento colectivo es autoorganizado con zonas de corte discretas y ángulos de reposo, característicos de los materiales granulares no cohesivos.

Aplicaciones

El supuesto fundamental del método es que el material consta de partículas separadas y discretas. Estas partículas pueden tener diferentes formas y propiedades que influyen en el contacto entre partículas. Algunos ejemplos son:

Las industrias típicas que utilizan DEM son:

Esquema del método

Una simulación DEM se inicia generando primero un modelo, lo que da como resultado la orientación espacial de todas las partículas y la asignación de una velocidad inicial . Las fuerzas que actúan sobre cada partícula se calculan a partir de los datos iniciales y las leyes físicas y modelos de contacto relevantes. Generalmente, una simulación consta de tres partes: la inicialización, el paso de tiempo explícito y el posprocesamiento. El paso de tiempo generalmente requiere un paso de clasificación del vecino más cercano para reducir el número de posibles pares de contactos y disminuir los requisitos computacionales; A menudo, esto sólo se realiza periódicamente.

Es posible que sea necesario considerar las siguientes fuerzas en simulaciones macroscópicas:

A nivel molecular, podemos considerar:

Todas estas fuerzas se suman para encontrar la fuerza total que actúa sobre cada partícula. Se emplea un método de integración para calcular el cambio en la posición y la velocidad de cada partícula durante un determinado paso de tiempo a partir de las leyes del movimiento de Newton . Luego, las nuevas posiciones se utilizan para calcular las fuerzas durante el siguiente paso y este ciclo se repite hasta que finaliza la simulación.

Los métodos de integración típicos utilizados en un método de elementos discretos son:

DEM térmico

El método de elementos discretos se aplica ampliamente para la consideración de interacciones mecánicas en problemas de muchos cuerpos, particularmente materiales granulares. Entre las diversas extensiones de DEM, la consideración del flujo de calor es particularmente útil. En términos generales, en los métodos DEM térmico, se considera el acoplamiento termomecánico, mediante el cual se consideran las propiedades térmicas de un elemento individual para modelar el flujo de calor a través de un medio macroscópico granular o de múltiples elementos sujeto a una carga mecánica. [9] Las fuerzas entre partículas, calculadas como parte del DEM clásico, se utilizan para determinar áreas de verdadero contacto entre partículas y así modelar la transferencia conductiva de calor de un elemento sólido a otro. Otro aspecto que se considera en DEM es la conducción, radiación y convección del calor en fase gaseosa en los espacios entre partículas. Para facilitar esto, es necesario considerar las propiedades de la fase gaseosa entre elementos en términos de presión, conductividad del gas y camino libre medio de las moléculas de gas. [10]

Fuerzas de largo alcance

Cuando se tienen en cuenta las fuerzas de largo alcance (normalmente la gravedad o la fuerza de Coulomb), es necesario calcular la interacción entre cada par de partículas. Tanto el número de interacciones como el costo de cálculo aumentan cuadráticamente con el número de partículas. Esto no es aceptable para simulaciones con una gran cantidad de partículas. Una posible forma de evitar este problema es combinar algunas partículas, que están lejos de la partícula en cuestión, en una pseudopartícula. Consideremos como ejemplo la interacción entre una estrella y una galaxia distante : el error que surge al combinar todas las estrellas de la galaxia distante en una masa puntual es insignificante. Para decidir qué partículas se pueden combinar en una pseudopartícula se utilizan los llamados algoritmos de árbol . Estos algoritmos organizan todas las partículas en un árbol, un árbol cuádruple en el caso bidimensional y un ocárbol en el caso tridimensional .

Sin embargo, las simulaciones en dinámica molecular dividen el espacio en el que se realiza la simulación en células. Las partículas que salen por un lado de una celda simplemente se insertan en el otro lado ( condiciones de contorno periódicas ); Lo mismo ocurre con las fuerzas. La fuerza ya no se tiene en cuenta después de la llamada distancia de corte (normalmente la mitad de la longitud de una celda), de modo que una partícula no se ve influenciada por la imagen especular de la misma partícula en el otro lado de la celda. Ahora se puede aumentar el número de partículas simplemente copiando las células.

Los algoritmos para hacer frente a la fuerza de largo alcance incluyen:

Método combinado de elementos finitos y discretos

Siguiendo el trabajo de Munjiza y Owen, el método combinado de elementos finitos discretos se ha desarrollado aún más para diversas partículas irregulares y deformables en muchas aplicaciones, incluidas la formación de tabletas farmacéuticas, [11] simulaciones de flujo y embalaje, [12] y el análisis de impacto. [13]

Ventajas y limitaciones

Ventajas

Desventajas

Ver también

Referencias

  1. ^ Peng, Z.; Doroodchi, E.; Moghtaderi, B. (2020). "Modelado de transferencia de calor en simulaciones de procesos térmicos basadas en el método de elementos discretos (DEM): desarrollo de teoría y modelo". Progresos en Ciencias de la Energía y la Combustión . 79, 100847: 100847. doi : 10.1016/j.pecs.2020.100847. S2CID  218967044.
  2. ^ Papadikis, K.; Gu, S.; Bridgewater, AV (2009). "Modelado CFD de la pirólisis rápida de biomasa en reactores de lecho fluidizado: modelado del impacto de la contracción de la biomasa" (PDF) . Revista de Ingeniería Química . 149 (1–3): 417–427. doi :10.1016/j.cej.2009.01.036.
  3. ^ Kafui, KD; Thornton, C.; Adams, MJ (2002). "Modelado de fluido continuo de partículas discretas de lechos fluidizados gas-sólidos". Ciencias de la Ingeniería Química . 57 (13): 2395–2410. Código Bib : 2002ChEnS..57.2395K. doi :10.1016/S0009-2509(02)00140-9.
  4. ^ Triviño, LF; Mohanty, B. (2015). "Evaluación de la iniciación y propagación de grietas en roca a partir de ondas de tensión inducidas por explosiones y expansión de gas mediante sismometría de agujeros cruzados y método FEM-DEM". Revista Internacional de Mecánica de Rocas y Ciencias Mineras . 77 : 287–299. Código Bib : 2015IJRMM..77..287T. doi :10.1016/j.ijrmms.2015.03.036.
  5. ^ Cundall, Peter. A.; Strack, Otto DL (1979). "Modelo numérico discreto para ensamblajes granulares" (PDF) . Geotécnica . 29 (1): 47–65. doi :10.1680/geot.1979.29.1.47.
  6. ^ Munjiza, Ante (2004). El método combinado de elementos finitos y discretos . Chichester: Wiley. ISBN 978-0-470-84199-0.
  7. ^ Alizadeh, Mohammadreza; Hassanpour, Ali; Pasha, Mehrdad; Ghadiri, Mojtaba; Bayly, Andrés (1 de septiembre de 2017). "El efecto de la forma de las partículas sobre la segregación prevista en mezclas de polvo binario" (PDF) . Tecnología en polvo . 319 : 313–322. doi :10.1016/j.powtec.2017.06.059. ISSN  0032-5910.
  8. ^ Behjani, Mohammadreza Alizadeh; Motlagh, Yousef Ghaffari; Bayly, Andrés; Hassanpour, Ali (7 de noviembre de 2019). "Evaluación del rendimiento de la mezcla de mezclas de polvos farmacéuticos en un mezclador continuo utilizando el método de elementos discretos (DEM)". Tecnología en polvo . 366 : 73–81. doi :10.1016/j.powtec.2019.10.102. ISSN  0032-5910. S2CID  209718900. Archivado desde el original el 21 de febrero de 2020.
  9. ^ Gan, Yixiang; Hernández, Francisco; Hanaor, Dorian; Annabattula, Ratna; Kamlah, Marc; Pereslavtsev, Pavel (2014). "Análisis de elementos térmicos discretos de una manta reproductora sólida de la UE sometida a irradiación de neutrones". Ciencia y Tecnología de Fusión . 66 (1): 83–90. arXiv : 1406.4199 . Código Bib : 2014FuST...66...83G. doi :10.13182/FST13-727. S2CID  51903434.
  10. ^ Tsory, Tal; Ben-Jacob, Nir; Brosh, Tamir; Levy, Avi (2013). "Modelado térmico DEM-CFD y simulación de transferencia de calor a través de lecho empacado". Tecnología en polvo . 244 : 52–60. doi :10.1016/j.powtec.2013.04.013.
  11. ^ Lewis, RW; Gethin, DT; Yang, XS; Rowe, RC (2005). "Un método combinado de elementos finitos y discretos para simular la formación de tabletas en polvo farmacéutico". Revista internacional de métodos numéricos en ingeniería . 62 (7): 853. arXiv : 0706.4406 . Código Bib : 2005IJNME..62..853L. doi :10.1002/nme.1287. S2CID  122962022.
  12. ^ Gethin, DT; Yang, XS; Lewis, RW (2006). "Un esquema bidimensional combinado de elementos finitos y discretos para simular el flujo y la compactación de sistemas que comprenden partículas irregulares". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 195 (41–43): 5552. Código bibliográfico : 2006CMAME.195.5552G. doi :10.1016/j.cma.2005.10.025.
  13. ^ Chen, Y.; Mayo, IM (2009). "Miembros de hormigón armado sometidos a impactos por caída de peso". Actas del ICE - Estructuras y Edificaciones . 162 : 45–56. doi :10.1680/stbu.2009.162.1.45.
  14. ^ SJ Antonio (2007). "Vínculo entre las propiedades de una sola partícula y las propiedades macroscópicas en conjuntos de partículas: papel de las estructuras dentro de las estructuras". Transacciones filosóficas de la Royal Society de Londres, Serie: A. 356 : 2879-2891.
  15. ^ SJ Antony, D. Chapman, J. Sujatha y T. Barakat (2015). "La interacción entre las inclusiones de diferentes tamaños y su proximidad a los límites de las paredes influye en la naturaleza de su distribución de tensiones dentro de las inclusiones dentro del embalaje de partículas". Tecnología en polvo . 286 : 286, 98-106.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  16. ^ Él, Yi; Bayly, Andrew E.; Hassanpour, Ali; Müller, Frans; Wu, Ke; Yang, Dongmin (1 de octubre de 2018). "Un método SPH-DEM acoplado basado en GPU para flujo de partículas y fluido con superficies libres". Tecnología en polvo . 338 : 548–562. doi : 10.1016/j.powtec.2018.07.043 . ISSN  0032-5910.

Bibliografía

Libro

Periódico

Actas