stringtranslate.com

Método de elementos analíticos

El método de elementos analíticos ( AEM ) es un método numérico utilizado para la solución de ecuaciones diferenciales parciales . [1] [2] [3] Fue desarrollado inicialmente por ODL Strack en la Universidad de Minnesota . Es similar en naturaleza al método de elementos de contorno (BEM), ya que no se basa en la discretización de volúmenes o áreas en el sistema modelado; solo se discretizan los límites internos y externos. Una de las principales distinciones entre AEM y BEM es que las integrales de contorno se calculan analíticamente. Aunque originalmente se desarrolló para modelar el flujo de agua subterránea, [4] AEM se ha aplicado posteriormente a otros campos de estudio, incluidos estudios de flujo y conducción de calor, ondas periódicas y deformación por fuerza. [5]

Flujo alrededor de cilindros impermeables. Resuelto con el AEM utilizando 20 coeficientes en las expansiones en serie.

Base matemática

La premisa básica del método de elementos analíticos es que, para ecuaciones diferenciales lineales , las soluciones elementales pueden superponerse para obtener soluciones más complejas. Existe un conjunto de soluciones analíticas 2D y 3D ("elementos") disponibles para diferentes ecuaciones de gobierno. Estos elementos corresponden típicamente a una discontinuidad en la variable dependiente o su gradiente a lo largo de un límite geométrico (por ejemplo, punto, línea, elipse, círculo, esfera, etc.). Esta discontinuidad tiene una forma funcional específica (generalmente un polinomio en 2D) y puede manipularse para satisfacer las condiciones de contorno de Dirichlet, Neumann o Robin (mixtas). Cada solución analítica es infinita en el espacio y/o el tiempo.

Comúnmente, cada solución analítica contiene grados de libertad (coeficientes) que pueden calcularse para cumplir con las condiciones de contorno prescritas a lo largo del borde del elemento. Para obtener una solución global (es decir, los coeficientes de elemento correctos), se resuelve un sistema de ecuaciones de manera que se satisfagan las condiciones de contorno a lo largo de todos los elementos (utilizando la colocación , la minimización de mínimos cuadrados o un enfoque similar). En particular, la solución global proporciona una descripción espacialmente continua de la variable dependiente en todas partes del dominio infinito, y la ecuación gobernante se satisface exactamente en todas partes excepto a lo largo del borde del elemento, donde la ecuación gobernante no es estrictamente aplicable debido a la discontinuidad.

La capacidad de superponer numerosos elementos en una única solución significa que se pueden obtener soluciones analíticas para condiciones de contorno de cualquier complejidad. Es decir, se pueden resolver modelos que tienen geometrías complejas, contornos rectos o curvos, contornos múltiples, condiciones de contorno transitorias, capas de acuíferos múltiples, propiedades que varían por partes y propiedades que varían continuamente. Los elementos se pueden implementar utilizando expansiones de campo lejano de modo que los modelos que contienen muchos miles de elementos se puedan resolver de manera eficiente y con alta precisión.

El método de elementos analíticos se ha aplicado a problemas de flujo de agua subterránea regidos por una variedad de ecuaciones diferenciales parciales lineales, incluidas la ecuación de Laplace , la ecuación de Poisson , la ecuación de Helmholtz modificada, [6] la ecuación del calor y las ecuaciones biarmónicas . A menudo, estas ecuaciones se resuelven utilizando variables complejas que permiten utilizar técnicas matemáticas disponibles en la teoría de variables complejas. Una técnica útil para resolver problemas complejos es utilizar el mapeo conforme que mapea el límite de una geometría, por ejemplo, una elipse, sobre el límite del círculo unitario donde se conoce la solución.

En el método analítico de elementos se utilizan el potencial de descarga y la función de caudal , o el potencial complejo combinado. Este potencial vincula las propiedades físicas del sistema de agua subterránea, la carga hidráulica o los límites de flujo, a una representación matemática de un potencial. Esta representación matemática se puede utilizar para calcular el potencial en términos de posición y, por lo tanto, también para resolver problemas de flujo de agua subterránea. Los elementos se desarrollan resolviendo las condiciones de contorno para cualquiera de estas dos propiedades, la carga hidráulica o el límite de flujo, lo que da como resultado soluciones analíticas capaces de abordar numerosas condiciones de contorno.

Comparación con otros métodos

Como se mencionó, el método de elementos analíticos no se basa en la discretización del volumen o el área en el modelo, como en los métodos de elementos finitos o de diferencias finitas . Por lo tanto, puede modelar problemas complejos con un error del orden de la precisión de la máquina. Esto se ilustra en un estudio que modeló un acuífero isotrópico altamente heterogéneo al incluir 100.000 heterogeneidades esféricas con una conductividad aleatoria y rastrear 40.000 partículas. [7] El método de elementos analíticos se puede utilizar de manera eficiente como verificación o como una herramienta de detección en proyectos más grandes, ya que puede calcular de manera rápida y precisa el flujo de agua subterránea para muchos problemas complejos. [8] [9]

A diferencia de otros métodos de modelado de aguas subterráneas de uso común, por ejemplo, el método de elementos finitos o el método de elementos finitos diferentes , el AEM no divide el dominio del modelo en celdas. Esto ofrece la ventaja de que el modelo es válido para cualquier punto dado en el dominio del modelo. Sin embargo, también impone que el dominio no se divida tan fácilmente en regiones de, por ejemplo, conductividad hidráulica diferente, como cuando se modela con una cuadrícula de celdas; sin embargo, una solución a este problema es incluir subdominios en el modelo AEM. [10] También existen soluciones para implementar propiedades o estructuras que varían verticalmente en un acuífero en un modelo AEM. [11] [12] [13]

Véase también

Referencias

  1. ^ Strack, Otto DL, 1943- (1989). Mecánica de aguas subterráneas . Englewood Cliffs, NJ: Prentice Hall. ISBN 0-13-365412-5.OCLC 16276592  .{{cite book}}: CS1 maint: nombres múltiples: lista de autores ( enlace ) CS1 maint: nombres numéricos: lista de autores ( enlace )
  2. ^ Strack, Otto DL (agosto de 2017). Mecánica analítica de aguas subterráneas. doi :10.1017/9781316563144. ISBN 9781316563144. Recuperado el 20 de abril de 2020 . {{cite book}}: |website=ignorado ( ayuda )
  3. ^ Haitjema, HM (Henk M.) (1995). Modelado analítico de elementos del flujo de agua subterránea . San Diego: Academic Press. ISBN 978-0-08-049910-9.OCLC 162129095  .
  4. ^ Haitjema, HM (1995), "Modelado analítico de elementos", Modelado analítico de elementos del flujo de agua subterránea , Elsevier, págs. 203-305, doi :10.1016/b978-012316550-3/50005-2, ISBN 978-0-12-316550-3, consultado el 11 de septiembre de 2023
  5. ^ Steward, David R. (17 de septiembre de 2020). Método de elementos analíticos. Oxford University PressOxford. doi :10.1093/oso/9780198856788.001.0001. ISBN 978-0-19-885678-8.
  6. ^ Strack, ODL; Namazi, T (octubre de 2014). "Una nueva formulación para el flujo constante en múltiples acuíferos: un elemento analítico para la infiltración constante por partes". Water Resources Research . 50 (10): 7939–7956. doi :10.1002/2014WR015479.
  7. ^ Janković, I.; Fiori, A.; Dagan, G. (2006). "Modelado del flujo y transporte en acuíferos tridimensionales altamente heterogéneos: ergodicidad, gaussianidad y comportamiento anómalo—1. Cuestiones conceptuales y simulaciones numéricas". Investigación en recursos hídricos . 42 (6): W06D12. Código Bibliográfico :2006WRR....42.6D12J. doi : 10.1029/2005WR004734 . ISSN  1944-7973.
  8. ^ Hunt, Randall J. (2006). "Aplicaciones de modelado de aguas subterráneas utilizando el método de elementos analíticos". Aguas subterráneas . 44 (1): 5–15. doi : 10.1111/j.1745-6584.2005.00143.x . ISSN  1745-6584. PMID  16405461. S2CID  24530553.
  9. ^ Kraemer, Stephen R. (2007). "Modelado analítico de elementos de agua subterránea como programa de investigación (1980 a 2006)". Aguas subterráneas . 45 (4): 402–408. doi :10.1111/j.1745-6584.2007.00314.x. ISSN  1745-6584. PMID  17600570. S2CID  26319150.
  10. ^ Fitts, CR (2010). "Modelado de sistemas acuíferos con elementos analíticos y subdominios: MODELADO DE ACUÍFEROS CON SUBDOMINIOS ANALÍTICOS". Investigación en Recursos Hídricos . 46 (7). doi :10.1029/2009WR008331.
  11. ^ Bakker, Mark; Strack, Otto DL (10 de febrero de 2003). "Elementos analíticos para el flujo multiacuífero". Revista de hidrología . 271 (1): 119–129. Código Bibliográfico :2003JHyd..271..119B. doi :10.1016/S0022-1694(02)00319-0. ISSN  0022-1694.
  12. ^ Strack, ODL; Ausk, BK (agosto de 2015). "Una formulación para el flujo de agua subterránea integrado verticalmente en un acuífero costero estratificado: FLUJO DEL ACUÍFERO COSTERO ESTRATIFICADO". Water Resources Research . 51 (8): 6756–6775. doi : 10.1002/2015WR016887 .
  13. ^ Toller, Erik AL; Strack, Otto DL (2019). "Flujo de interfaz con conductividad hidráulica que varía verticalmente". Investigación de recursos hídricos . 55 (11): 8514–8525. Código Bibliográfico :2019WRR....55.8514T. doi :10.1029/2019WR024927. ISSN  1944-7973. S2CID  202924261.

Más información

Enlaces externos