Implementar el método de elemento finito a la ecuación ADR. aplicación al estado de morelos.

MONICA ANDREA NAVARRO RAMIREZ

En las últimas décadas la contaminación del aire ha sido objeto de investigación a nivel mundial, debido a los daños que provoca, principalmente, a la humanidad y al ambiente. Es por esto que los modelos computacionales son cada vez más indispensables debido de su potencia y flexibilidad para monitorear y predecir el movimiento y transformación de la contaminación. Hasta ahora, hay pocos estudios sobre estos aspectos de contaminación en el estado [1]. Actualmente para obtener datos en Morelos se utiliza un programa llamadoHYSPLIT [2] que accesa a datos de una sola estación meteorológica, la del aeropuerto, aproximadamente a 20 km del centro de la ciudad de Cuernavaca. En este proyecto, se trabaja con la ecuación de Advección-Difusión-Reacción (ADR) Ut = �����vUx + Ux x +F para modelar los fenómenos de transporte, especialmente para la transferencia de calor y energía [3]. El proyecto consistió, primero, en la derivación del algoritmo matemático del método de elemento finito (FEM) en una dimensión, empezando con la forma estacionaria, Ut = 0, y posteriormente la forma no estacionaria. Para cualquier problema, el uso de FEM consiste en cuatro conceptos básicos [4]: la forma fuerte, forma débil, aproximación de elemento finito y el sistema de ecuaciones lineales. Se utilizó la aproximación de Galerkin/Ritz y el método de Runge-Kutta 4 (RK4) para la resolución del sistema de ecuaciones lineales.Una vez derivadas las ecuaciones con las cuales se iba a hacer la implementación del código en el lenguaje de programación FORTRAN- 90. Después, se llevó a cabo un estudio numérico de estabilidad, probando diferentes valores de velocidad, coeficiente de difusión, x y t , esto para diferentes números de nodos. Para la parte de implementación al estado de Morelos, se construyó un perfil de viento usando datos de 27 estaciones meteorológicas del estado de los meses de octubre y noviembre de 2010. Y se aproximaron algunas trayectorias para probar el código con valores reales. Se tomaron unos pasos preliminares para la implementación de la ecuación en dos dimensiones. Una de las ventajas de FEM es la apliación de una malla irregular, la cual da mayor felixibildad y mejor adaptación a la topología. Se llevó a cabo la triangulación de la malla de tipo Delaunay, se construyó vía el algoritmo de Bowyer-Watson [5, 6] y se implementó el código.

Tipo de documento: Tesis de maestría

Formato: Adobe PDF

Audiencia: Investigadores

Idioma: Español

Área de conocimiento: INGENIERÍA Y TECNOLOGÍA

Campo disciplinar: CIENCIAS TECNOLÓGICAS

Nivel de acceso: Acceso Abierto