Bibliografía de la Clase:
Una ecuación diferencial parcial (EDP) es una ecuación que contiene una función multivariable desconocida (ej. ) y sus derivadas parciales. Este tipo de ecuaciones describen varios fenómenos físicos como el calor, el sonido, dinámica de fluidos, etc.
El orden de la EDP está determinado por el grado de la derivada parcial más alto.
Una EDP de segundo orden genérica tiene la forma:
en una cierta región . Se dice que la EDP es de tipo:
- Hiperbólica: si en .
- Parabólica: si en .
- Elíptica: si en .
Ejemplos de EDP elípticas son la ecuación de Poisson y Laplace:
Ecuación de Poisson
Usada en el estudio de problemas físicos como distribución de calor, energía potencial sobre un punto y dinámica de fluidos.
(1)
donde f(x,y) es una función definida sobre una región con borde .
La ecuación de Poison es del tipo elíptica ya que y , por lo tanto .
Esta ecuación es equivalente a:
donde es el operador de Laplace o Laplaciano:
Ecuación de Laplace
es un caso particular de la ecuación de Poisson donde que significa que la región es plana.
Si el valor de en la región es determinada por el valor de la función en el borde se les denomina condiciones de borde de Dirichlet que viene dado por:
En la figura, la región está definida como:
Para resolver numéricamente este tipo de ecuaciones se utiliza una adaptación del método de diferencias finitas. El método consiste en generar una grilla equiespaciada en la región para obtener aproximaciones de la función en los puntos de la grilla iterativamente.
El número de puntos de la grilla para e son y respectivamente, por lo tanto el espaciado entre los puntos es:
Para expresar numéricamente e se utiliza la serie de Taylor alrededor de . De manera aproximada se tiene:
y luego se suman para obtener una expresión para :
Podemos usar las expansiones de Taylor mostradas anteriormente para resolver la ecuación de Poisson expresando:
ahora podemos expresar la ecuación (1) en el punto (x_i,y_j) aproximadamente como:
Agrupando términos, y expresando se tiene:
(2)
Para definir las condiciones de borde, se debe especificar que el valor de en la frontera está definido por la función .
(3)
Las ecuaciones (2) y (3) definen un sistema de ecuaciones que puede ser representado matricialmente.
Determine la distribución del calor en estado estable en una superficie metálica de dimensiones utilizando una grilla de tamaño . Las condiciones de borde en los ejes e son 0 y el calor en los otros dos bordes crece linealmente desde 0 a 100 grados.
Estudie cómo se utiliza el método de Gauss Seidel para encontrar los y encuentre la solución de la ecuación (5). Encuentre el error asociado al método numérico considerando que la solucíón analítica es .
En un recinto cuadrado de lado 1 considere la ecuación de Poisson:
con condición de borde igual a cero en los cuatro costados. Encuentre la solución de este problema usando el método de diferencias finitas, considerando .