Estabilización del Medio de Capa Perfectamente Acoplada Mediante Funciones de Estiramiento Multiaxiales

Abstract

En los últimos años se ha demostrado la efectividad de los bordes de capas perfectamente acopladas (PML) para absorber ondas elásticas salientes al dominio del modelo numérico. Bécache et al. (J. Comput. Phys., 188, 399-433, 2003) advierten que el medio PML es intrínsecamente inestable, es decir, inestable independientemente de la discretización que se realice. Una forma de estabilizar el medio PML es presentada por Meza-Fajardo y Papageorgiou (Bull. Seismol. Soc. Am., 98, 1811-1836, 2008), quienes proponen utilizar funciones de estiramiento multiaxiales. Sin embargo, su desarrollo es cuestionado por Fathi et al. (Int. J. Numer. Meth. Engng., DOI: 10.1002/nme.4780, 2014), quienes observan que se desestiman algunos términos del jacobiano de la transformación de coordenadas. Por este motivo la capa absorbente obtenida es denominada M-PML “no corregida”. En este trabajo se plantea la formulación fuerte de capa absorbente M-PML “corregida” mediante funciones de estiramiento multiaxiales utilizando todos los términos de la transformación de coordenadas. Se comparan los resultados del modelo híbrido estable con soluciones de referencia obtenidas a partir de modelos de dominio extendido.