stringtranslate.com

Métodos sin malla

20 puntos y sus células Voronoi.

En el campo del análisis numérico , los métodos meshfree son aquellos que no requieren conexión entre nodos del dominio de simulación, es decir, una malla , sino que se basan en la interacción de cada nodo con todos sus vecinos. Como consecuencia, las propiedades extensivas originales, como la masa o la energía cinética, ya no se asignan a los elementos de la malla sino a los nodos individuales. Los métodos sin malla permiten la simulación de algunos tipos de problemas que de otro modo serían difíciles, a costa de tiempo de cálculo y esfuerzo de programación adicionales. La ausencia de una malla permite simulaciones lagrangianas , en las que los nodos pueden moverse según el campo de velocidades .

Motivación

Los métodos numéricos como el método de diferencias finitas , el método de volúmenes finitos y el método de elementos finitos se definieron originalmente en mallas de puntos de datos. En dicha malla, cada punto tiene un número fijo de vecinos predefinidos, y esta conectividad entre vecinos se puede utilizar para definir operadores matemáticos como la derivada . Luego, estos operadores se utilizan para construir las ecuaciones que se van a simular, como las ecuaciones de Euler o las ecuaciones de Navier-Stokes .

Pero en simulaciones donde el material que se simula puede moverse (como en la dinámica de fluidos computacional ) o donde pueden ocurrir grandes deformaciones del material (como en simulaciones de materiales plásticos ), la conectividad de la malla puede ser difícil de mantener sin introducir errores en la simulación. Si la malla se enreda o se degenera durante la simulación, es posible que los operadores definidos en ella ya no proporcionen valores correctos. La malla se puede recrear durante la simulación (un proceso llamado remallado), pero esto también puede introducir errores, ya que todos los puntos de datos existentes deben asignarse a un conjunto nuevo y diferente de puntos de datos. Los métodos sin malla están destinados a remediar estos problemas. Los métodos sin malla también son útiles para:

Ejemplo

En una simulación tradicional en diferencias finitas , el dominio de una simulación unidimensional sería alguna función , representada como una malla de valores de datos en puntos , donde

Podemos definir las derivadas que ocurren en la ecuación que se simula usando algunas fórmulas en diferencias finitas en este dominio, por ejemplo

y

Luego podemos usar estas definiciones y sus derivadas espaciales y temporales para escribir la ecuación que se está simulando en forma de diferencias finitas y luego simular la ecuación con uno de los muchos métodos de diferencias finitas .

En este ejemplo simple, los pasos (aquí el paso espacial y el paso de tiempo ) son constantes a lo largo de toda la malla, y los vecinos de malla izquierdo y derecho del valor de datos en son los valores en y , respectivamente. Generalmente en diferencias finitas se pueden permitir de manera muy simple pasos variables a lo largo de la malla, pero todos los nodos originales deben conservarse y pueden moverse de forma independiente sólo deformando los elementos originales. Si solo dos de todos los nodos cambian su orden, o incluso si solo se agrega o elimina un nodo de la simulación, eso crea un defecto en la malla original y la aproximación simple en diferencias finitas ya no puede mantenerse.

La hidrodinámica de partículas suavizadas (SPH), uno de los métodos sin malla más antiguos, resuelve este problema tratando los puntos de datos como partículas físicas con masa y densidad que pueden moverse con el tiempo y llevar consigo algún valor . SPH luego define el valor de entre las partículas por

donde es la masa de la partícula , es la densidad de la partícula y es una función central que opera en puntos de datos cercanos y se elige por su suavidad y otras cualidades útiles. Por linealidad, podemos escribir la derivada espacial como

Luego podemos usar estas definiciones de y sus derivadas espaciales para escribir la ecuación que se está simulando como una ecuación diferencial ordinaria y simular la ecuación con uno de muchos métodos numéricos . En términos físicos, esto significa calcular las fuerzas entre las partículas y luego integrar estas fuerzas a lo largo del tiempo para determinar su movimiento.

La ventaja de SPH en esta situación es que las fórmulas de y sus derivados no dependen de ninguna información de adyacencia sobre las partículas; pueden usar las partículas en cualquier orden, por lo que no importa si las partículas se mueven o incluso intercambian lugares.

Una desventaja de SPH es que requiere programación adicional para determinar los vecinos más cercanos de una partícula. Dado que la función kernel solo devuelve resultados distintos de cero para partículas cercanas dentro del doble de la "longitud de suavizado" (porque normalmente elegimos funciones kernel con soporte compacto ), sería una pérdida de esfuerzo calcular las sumas anteriores sobre cada partícula en una simulación grande. Por lo general, los simuladores SPH requieren algún código adicional para acelerar el cálculo del vecino más cercano.

Historia

Uno de los primeros métodos sin malla es la hidrodinámica de partículas suavizadas , presentada en 1977. [1] Libersky et al. [2] fueron los primeros en aplicar SPH en mecánica de sólidos. Los principales inconvenientes de SPH son los resultados inexactos cerca de los límites y la inestabilidad de la tensión que fue investigada por primera vez por Swegle. [3]

En la década de 1990 surgió una nueva clase de métodos sin malla basada en el método Galerkin . Este primer método llamado método del elemento difuso [4] (DEM), iniciado por Nayroles et al., utilizó la aproximación MLS en la solución de Galerkin de ecuaciones diferenciales parciales, con derivadas aproximadas de la función MLS. A partir de entonces, Belytschko fue pionero en el método Element Free Galerkin (EFG), [5] que empleaba MLS con multiplicadores de Lagrange para imponer condiciones de contorno, cuadratura numérica de orden superior en la forma débil y derivadas completas de la aproximación MLS que proporcionaba una mayor precisión. Casi al mismo tiempo, surgió el método de reproducción de partículas del núcleo [6] (RKPM), cuya aproximación motivó en parte a corregir la estimación del núcleo en SPH: para brindar precisión cerca de los límites, en discretizaciones no uniformes y precisión de orden superior en general. . En particular, en un desarrollo paralelo, los métodos de puntos materiales se desarrollaron casi al mismo tiempo [7] , que ofrecen capacidades similares. Los métodos de puntos materiales se utilizan ampliamente en la industria cinematográfica para simular grandes deformaciones mecánicas de sólidos, como la nieve en la película Frozen . [8] Chen, Liu y Li desarrollaron ampliamente el RKPM y otros métodos sin malla a finales de la década de 1990 para una variedad de aplicaciones y diversas clases de problemas. [9] Durante la década de 1990 y posteriormente se desarrollaron varias otras variedades, incluidas las que se enumeran a continuación.

Lista de métodos y acrónimos

Generalmente se considera que los siguientes métodos numéricos pertenecen a la clase general de métodos "sin malla". Las siglas se proporcionan entre paréntesis.

Métodos relacionados:

Desarrollo reciente

Las principales áreas de avance en los métodos sin malla son abordar problemas relacionados con el cumplimiento de límites esenciales, la cuadratura numérica y las deformaciones grandes y de contacto. [25] La forma débil común requiere una fuerte aplicación de las condiciones de contorno esenciales, sin embargo, los métodos sin malla en general carecen de la propiedad delta de Kronecker . Esto hace que la aplicación de las condiciones de contorno esenciales no sea trivial, al menos más difícil que el método de elementos finitos , donde se pueden imponer directamente. Se han desarrollado técnicas para superar esta dificultad e imponer condiciones con fuerza. Se han desarrollado varios métodos para imponer débilmente las condiciones de contorno esenciales , incluidos los multiplicadores de Lagrange , el método de Nitche y el método de penalización.

En cuanto a la cuadratura , generalmente se prefiere la integración nodal, que ofrece simplicidad, eficiencia y mantiene el método sin malla libre de cualquier malla (en lugar de usar la cuadratura de Gauss , que requiere una malla para generar puntos y pesos de cuadratura). Sin embargo, la integración nodal adolece de inestabilidad numérica debido a la subestimación de la energía de deformación asociada con los modos de longitud de onda corta [26] y también produce resultados inexactos y no convergentes debido a la subintegración de la forma débil. [27] Un avance importante en la integración numérica ha sido el desarrollo de una integración nodal conforme estabilizada (SCNI) que proporciona un método de integración nodal que no sufre ninguno de estos problemas. [27] El método se basa en el suavizado de deformaciones que satisface la prueba de parche de primer orden . Sin embargo, más tarde se descubrió que los modos de baja energía todavía estaban presentes en SCNI y se desarrollaron métodos de estabilización adicionales. Este método se ha aplicado a una variedad de problemas que incluyen placas delgadas y gruesas, poromecánica y problemas dominados por convección, entre otros. [25] Más recientemente, se ha desarrollado un marco para pasar pruebas de parche de orden arbitrario, basado en un método de Petrov-Galerkin . [28]

Un avance reciente en los métodos sin malla tiene como objetivo el desarrollo de herramientas computacionales para la automatización en modelado y simulaciones. Esto es posible gracias a la denominada formulación débil debilitada (W2), basada en la teoría del espacio G. [29] [30] La formulación W2 ofrece posibilidades para formular varios modelos (uniformemente) "suaves" que funcionan bien con mallas triangulares. Debido a que una malla triangular se puede generar automáticamente, resulta mucho más fácil volver a mallar y, por lo tanto, permite la automatización en el modelado y la simulación. Además, los modelos W2 se pueden hacer lo suficientemente suaves (de manera uniforme) para producir soluciones de límite superior (para problemas de conducción forzada). Junto con los modelos rígidos (como los modelos FEM totalmente compatibles), se puede unir cómodamente la solución por ambos lados. Esto permite una fácil estimación del error para problemas generalmente complicados, siempre que se pueda generar una malla triangular. Los modelos W2 típicos son los métodos de interpolación de puntos suavizados (o S-PIM). [17] El S-PIM puede estar basado en nodos (conocido como NS-PIM o LC-PIM), [31] basado en bordes (ES-PIM), [32] y basado en células (CS-PIM). [33] El NS-PIM se desarrolló utilizando la llamada técnica SCNI. [27] Luego se descubrió que NS-PIM es capaz de producir una solución de límite superior y bloqueo volumétrico libre. [34] El ES-PIM tiene una precisión superior y el CS-PIM se comporta entre el NS-PIM y el ES-PIM. Además, las formulaciones W2 permiten el uso de funciones de base polinómica y radial en la creación de funciones de forma (se adaptan a las funciones de desplazamiento discontinuo, siempre que estén en el espacio G1), lo que abre más espacios para futuros desarrollos. La formulación W2 también ha llevado al desarrollo de una combinación de técnicas sin malla con las técnicas FEM bien desarrolladas, y ahora se puede utilizar malla triangular con excelente precisión y suavidad deseada. Una formulación típica de este tipo es el llamado método suavizado de elementos finitos (o S-FEM). [35] El S-FEM es la versión lineal del S-PIM, pero con la mayoría de las propiedades del S-PIM y mucho más simple.

Es una percepción general que los métodos sin malla son mucho más caros que sus homólogos FEM. Sin embargo, un estudio reciente ha descubierto que algunos métodos sin malla, como S-PIM y S-FEM, pueden ser mucho más rápidos que sus homólogos FEM. [17] [35]

El S-PIM y el S-FEM funcionan bien para problemas de mecánica sólida. Para los problemas de CFD, la formulación puede ser más sencilla, mediante una formulación sólida. También se ha desarrollado recientemente un método de suavizado de gradiente (GSM) para problemas CFD, implementando la idea de suavizado de gradiente de forma sólida. [36] [37] El GSM es similar a [FVM], pero utiliza operaciones de suavizado de gradiente exclusivamente en modelos anidados, y es un método numérico general para PDE.

La integración nodal se ha propuesto como una técnica para utilizar elementos finitos para emular un comportamiento sin malla. [ cita necesaria ] Sin embargo, el obstáculo que debe superarse al utilizar elementos integrados nodales es que las cantidades en los puntos nodales no son continuas y los nodos se comparten entre múltiples elementos.

Ver también

Referencias

  1. ^ Gingold, RA; Monaghan, JJ (1 de diciembre de 1977). "Hidrodinámica de partículas suavizadas: teoría y aplicación a estrellas no esféricas". Avisos mensuales de la Real Sociedad Astronómica . 181 (3): 375–389. Código bibliográfico : 1977MNRAS.181..375G. doi : 10.1093/mnras/181.3.375 .
  2. ^ Libersky, Larry D.; Petschek, Albert G.; Carney, Theodore C.; Hipp, Jim R.; Allahdadi, Firooz A. (noviembre de 1993). "Hidrodinámica lagrangiana de alta tensión". Revista de Física Computacional . 109 (1): 67–75. doi :10.1006/jcph.1993.1199.
  3. ^ Swegle, JW; Hicks, DL; Attaway, SW (enero de 1995). "Análisis de estabilidad hidrodinámica de partículas suavizadas". Revista de Física Computacional . 116 (1): 123-134. Código Bib : 1995JCoPh.116..123S. doi :10.1006/jcph.1995.1010.
  4. ^ Nayroles, B.; Touzot, G.; Villon, P. (1992). "Generalizando el método de los elementos finitos: aproximación difusa y elementos difusos". Mecánica Computacional . 10 (5): 307–318. Código Bib : 1992CompM..10..307N. doi :10.1007/BF00364252. S2CID  121511161.
  5. ^ Belytschko, T.; Lu, YY; Gu, L. (30 de enero de 1994). "Métodos Galerkin libres de elementos". Revista internacional de métodos numéricos en ingeniería . 37 (2): 229–256. Código Bib : 1994IJNME..37..229B. doi :10.1002/nme.1620370205.
  6. ^ Liu, ala Kam; Jun, Sukky; Zhang, Yi Fei (30 de abril de 1995). "Métodos de reproducción de partículas del núcleo". Revista internacional de métodos numéricos en fluidos . 20 (8–9): 1081–1106. Código Bib : 1995IJNMF..20.1081L. doi :10.1002/fld.1650200824.
  7. ^ Sulsky, D.; Chen, Z.; Schreyer, HL (septiembre de 1994). "Un método de partículas para materiales dependientes de la historia". Métodos Informáticos en Mecánica e Ingeniería Aplicadas . 118 (1–2): 179–196. doi :10.1016/0045-7825(94)90112-0.
  8. ^ https://www.math.ucla.edu/~jteran/papers/SSCTS13.pdf [ URL básica PDF ]
  9. ^ Liu, WK; Chen, Y.; junio, S.; Chen, JS; Belytschko, T.; Pan, C.; Uras, RA; Chang, CT (marzo de 1996). "Descripción general y aplicaciones de los métodos de reproducción de partículas del núcleo". Archivos de métodos computacionales en ingeniería . 3 (1): 3–80. doi :10.1007/BF02736130. S2CID  122241092.
  10. ^ Atluri, SN; Zhu, T. (24 de agosto de 1998). "Un nuevo enfoque Meshless Local Petrov-Galerkin (MLPG) en mecánica computacional". Mecánica Computacional . 22 (2): 117–127. Código Bib : 1998CompM..22..117A. doi :10.1007/s004660050346. S2CID  3688083.
  11. ^ Oliveira, T.; Portela, A. (diciembre de 2016). "Colocación de forma débil: un método local sin malla en elasticidad lineal". Análisis de ingeniería con elementos de frontera . 73 : 144-160. doi :10.1016/j.enganabound.2016.09.010.
  12. ^ Chen, Shang-Ying; Hsu, Kuo-Chin; Fan, Chia-Ming (15 de marzo de 2021). "Mejora del método generalizado de diferencias finitas para el modelado estocástico de flujo subterráneo". Revista de Física Computacional . 429 : 110002. Código bibliográfico : 2021JCoPh.42910002C. doi :10.1016/J.JCP.2020.110002. S2CID  228828681.
  13. ^ Chen, Shang-Ying; Wei, Jian-Yu; Hsu, Kuo-Chin (1 de octubre de 2023). "Asimilación de datos para el modelado de flujo subterráneo en tiempo real con ajustes de nodos sin malla dinámicamente adaptables". Ingeniería con Computadoras . doi :10.1007/s00366-023-01897-6. ISSN  1435-5663.
  14. ^ WK Liu; S. Jun; YF Zhang (1995). "Métodos de reproducción de partículas del núcleo". En t. J. Número. Métodos Ing . 20 (8–9): 1081–1106. Código Bib : 1995IJNMF..20.1081L. doi :10.1002/fld.1650200824.
  15. ^ A. Behzadan; Su Majestad Shodja; M. Khezri (2011). "Un enfoque unificado para el análisis matemático de RKPM generalizado, RKPM de gradiente y GMLS". Computadora. Métodos. Aplica. Mec. Ing . 200 (5–8): 540–576. Código Bib : 2011CMAME.200..540B. doi :10.1016/j.cma.2010.07.017.
  16. ^ Medidor, Christoph; Leinen, Peter; Yserentant, Harry (enero de 2000). "El método de masa finita". Revista SIAM de Análisis Numérico . 37 (6): 1768-1799. doi :10.1137/S0036142999352564.
  17. ^ abcd Liu, GR 2.a ed.: 2009 Métodos sin malla , CRC Press. 978-1-4200-8209-9
  18. ^ Sarler B, Vertnik R. Meshfree
  19. ^ Li, B.; Habbal, F.; Ortiz, M. (17 de septiembre de 2010). "Esquemas de aproximación óptimos sin malla de transporte para flujos de fluidos y plásticos". Revista internacional de métodos numéricos en ingeniería . 83 (12): 1541-1579. Código Bib : 2010IJNME..83.1541L. doi :10.1002/nme.2869. S2CID  18225521.
  20. ^ Caminante, Wade A.; Langowski, Jörg (6 de julio de 2012). "El método de reemplazo repetido: un método lagrangiano puro sin malla para dinámica de fluidos computacional". MÁS UNO . 7 (7): e39999. Código Bib : 2012PLoSO...739999W. doi : 10.1371/journal.pone.0039999 . PMC 3391243 . PMID  22866175. 
  21. ^ Ooi, EH; Popov, V. (mayo de 2012). "Una implementación eficiente del método de ecuación integral de base radial". Análisis de ingeniería con elementos de frontera . 36 (5): 716–726. doi :10.1016/j.enganabound.2011.12.001. S2CID  122004658.
  22. ^ Zhang, Xiong; Liu, Xiao‐Hu; Canción, Kang-Zu; Lu, Ming‐Wan (30 de julio de 2001). "Método sin malla de colocación de mínimos cuadrados". Revista internacional de métodos numéricos en ingeniería . 51 (9): 1089-1100. Código Bib : 2001IJNME..51.1089Z. doi :10.1002/nme.200. S2CID  119952479.
  23. ^ Boroomand, B.; Soghrati, S.; Movahedian, B. (2009). "Funciones de base exponencial en la solución de problemas elásticos estáticos y armónicos de tiempo en un estilo sin malla". Revista internacional de métodos numéricos en ingeniería . 81 (8): 971–1018. doi :10.1002/nme.2718. S2CID  4943418.
  24. ^ Ghoneim, A. (marzo de 2015). "Un método de elementos finitos de interfaz sin malla para modelar la fusión y solidificación de solutos isotérmicos en sistemas binarios". Elementos finitos en análisis y diseño . 95 : 20–41. doi :10.1016/j.finel.2014.10.002.
  25. ^ ab Chen, Jiun-Shyan ; Hillman, Michael; Chi, Sheng-Wei (abril de 2017). "Métodos Meshfree: avances logrados después de 20 años". Revista de Ingeniería Mecánica . 143 (4): 04017001. doi :10.1061/(ASCE)EM.1943-7889.0001176.
  26. ^ Belytschko, Ted; Guo, Yong; Kam Liu, ala; Ping Xiao, Shao (30 de julio de 2000). "Un análisis de estabilidad unificado de métodos de partículas sin malla". Revista internacional de métodos numéricos en ingeniería . 48 (9): 1359-1400. Código bibliográfico : 2000IJNME..48.1359B. doi :10.1002/1097-0207(20000730)48:9<1359::AID-NME829>3.0.CO;2-U.
  27. ^ abc Chen, Jiun-Shyan ; Wu, Cheng-Tang; Yoon, Sangpil; Tú, Yang (20 de enero de 2001). "Una integración nodal conforme estabilizada para métodos sin malla Galerkin". Revista internacional de métodos numéricos en ingeniería . 50 (2): 435–466. Código Bib : 2001IJNME..50..435C. doi :10.1002/1097-0207(20010120)50:2<435::AID-NME32>3.0.CO;2-A.
  28. ^ Chen, Jiun-Shyan ; Hillman, Michael; Rüter, Marcus (3 de agosto de 2013). "Una integración variacionalmente consistente de orden arbitrario para los métodos sin malla de Galerkin". Revista internacional de métodos numéricos en ingeniería . 95 (5): 387–418. Código Bib : 2013IJNME..95..387C. doi :10.1002/nme.4512. S2CID  124640562.
  29. ^ ab Liu, GR (2009). "Teoría del espacio AG y una forma débil debilitada (W2) para una formulación unificada de métodos compatibles e incompatibles: teoría de la Parte I". Revista internacional de métodos numéricos en ingeniería . 81 (9): 1093-1126. doi :10.1002/nme.2719. S2CID  123009384.
  30. ^ ab Liu, GR (2009). "Teoría del espacio AG y una forma débil debilitada (W2) para una formulación unificada de métodos compatibles e incompatibles: aplicaciones de la Parte II a problemas de mecánica de sólidos". Revista internacional de métodos numéricos en ingeniería . 81 (9): 1127-1156. doi :10.1002/nme.2720. S2CID  119378545.
  31. ^ Liu GR, Zhang GY, Dai KY, Wang YY, Zhong ZH, Li GY y Han X, Un método de interpolación de puntos conforme linealmente (LC-PIM) para problemas de mecánica de sólidos 2D, Revista internacional de métodos computacionales , 2 (4): 645–665, 2005.
  32. ^ GR Liu, GR Zhang. Métodos de interpolación de puntos suavizados basados ​​en bordes. Revista internacional de métodos computacionales, 5 (4): 621–646, 2008
  33. ^ Liu, GR; Zhang, GY (20 de noviembre de 2011). "Un espacio G normalizado y una formulación débil debilitada (W2) de un método de interpolación de puntos suavizados basado en células". Revista internacional de métodos computacionales . 06 (1): 147-179. doi :10.1142/S0219876209001796.
  34. ^ Liu, GR; Zhang, GY (14 de mayo de 2008). "Solución de límite superior a problemas de elasticidad: una propiedad única del método de interpolación de puntos linealmente conformes (LC-PIM)". Revista internacional de métodos numéricos en ingeniería . 74 (7): 1128-1161. Código Bib : 2008IJNME..74.1128L. doi :10.1002/nme.2204. S2CID  54088894.
  35. ^ abc Liu, GR, 2010 Métodos de elementos finitos suavizados , CRC Press, ISBN 978-1-4398-2027-8 . [ página necesaria ] 
  36. ^ Liu, GR; Xu, George X. (10 de diciembre de 2008). "Un método de suavizado de gradiente (GSM) para problemas de dinámica de fluidos". Revista internacional de métodos numéricos en fluidos . 58 (10): 1101-1133. Código Bib : 2008IJNMF..58.1101L. doi :10.1002/fld.1788. S2CID  53983110.
  37. ^ Zhang, Jian; Liu, GR; Lam, KY; Li, Hua; Xu, G. (noviembre de 2008). "Un método de suavizado de gradiente (GSM) basado en una ecuación que rige la forma fuerte para el análisis adaptativo de problemas de mecánica de sólidos". Elementos finitos en análisis y diseño . 44 (15): 889–909. doi :10.1016/j.finel.2008.06.006.
  38. ^ Liu, GR (20 de noviembre de 2011). "Sobre la teoría del espacio G". Revista internacional de métodos computacionales . 06 (2): 257–289. doi :10.1142/S0219876209001863.

Otras lecturas

enlaces externos