30 de mayo de 2023 a 2 de junio de 2023 Ciencias Naturales, Exactas y Ténicas
America/Havana zona horaria

Solución numérica de la ecuación de Helmholtz basada en el uso de precondicionadores. Experimentos en 1D.

No programado
23h 59m

Ponente

Sofía Albizu-Campos Rodríguez (acultad de Matemática y Computación de la Universidad de La Habana)

Descripción

La obtención de soluciones numéricas de las distintas variantes de la ecuación de Helmholtz,Δϕ+k^2 ϕ=fdondek∈R_+ y f,ϕ:R^n→C, es un problema muy estudiado. En particular, la matriz de coeficientes del sistema de ecuaciones lineales derivado de la discretización de la ecuación, tantopor el método de Diferencias Finitas (FDM) como por el de Elementos Finitos (FEM), espor lo general mal condicionada. Incluso, los métodos iterativos basados en subespacios de Krylovfallan en ocasiones en que la matriz de coeficientes posee valores propios muy pequeños, por lo que se ha estudiado el uso de precondicionadores para acelerar la convergencia de estos métodos. En 2006, Y.A. Erlangga, C.W. Oosterlee y C. Vuik introdujeron el precondicionadorLaplaciano con Desplazamiento Complejo (CSLP, en inglés ComplexShiftedLaplacianPreconditioner), cuyo uso se ha estandarizado aunque no es infalible. La ecuación de Helmholtz en una dimensión (1D)homogénea (f≡0) con condiciones de frontera complejas de Dirichlet y Robin, tiene solución analítica conocida. Este trabajo se propone obtener aproximaciones numéricas de calidad de esta solución para grandes valores de k. Se utilizará en la discretización de la ecuación tanto FDM como FEM. Para la solución del sistema de ecuaciones lineales se empleará el método basado en subespacios de Krylov de Mínimo Residuo Generalizado (GMRES por sus siglas en inglés), adicionando el uso del CSLP y de su factorización incompleta LU. Finalmente, se proponen experimentos numéricos con comparaciones entre los métodos utilizados en cuanto a su eficiencia y la calidad de las aproximaciones obtenidas.

Autor primario

Sofía Albizu-Campos Rodríguez (acultad de Matemática y Computación de la Universidad de La Habana)

Coautor

Victoria Hernández Mederos (Instituto de Cibernética, Matemática y Física (ICIMAF))

Materiales de la presentación

Todavía no hay materiales.