PDE-toolbox de Matlab
(Primeros Auxilios)
Curso 2004-2005
Irene Peral Walias
Explicacin sobre un ejemplo.
Sea el problema de contorno
u = 4 en B(0,1). u = 1 en B(0,1).
2 2 u ( x, y ) = x + y
Vamos a utilizar la herramienta pdetool de MATLAB para resolverlo, mediante el Mtodo de Elementos Finitos. Procedemos de la siguiente manera: 1. PDE Toolbox. En la pantalla de trabajo de MATLAB escribimos pdetool y damos enter.
2. Dominio. Introducimos el dominio del problema (B(0,1)).
- Pulsamos en la tecla subrayada en la pantalla y pintamos una circunferencia arbitraria.
- Con el ratn, hacemos doble clic en la figura e introducimos el radio y el centro de la circunferencia de forma precisa.
- Finalmente ajustamos los ejes. Para ello, pinchamos en Options Axes limits e indicamos los limites que deseemos.
Una vez hecho esto, pinchamos en Options Axes equal, obteniendo el dibujo del dominio deseado.
3. Condiciones de Contorno. Introducimos los valores en el contorno. Para ello pulsamos, con el ratn, sobre la tecla , sombreada en amarillo en la figura. En este momento, aparece dibujado (como se muestra en la figura inferior) la frontera de nuestro dominio.
Pulsamos en cada una de las regiones de nuestra frontera y nos aparece la pantalla que mostramos, donde introducimos los valores de contorno. (Obsrvese que en el problema que nos ocupa son condiciones de Dirichlet).
4. Ecuacin. Seguidamente pasamos a escribir la ecuacin de nuestro problema. Para ello, pulsamos, con el ratn, sobre la tecla , y aparece la siguiente pantalla donde escribimos nuestra ecuacin. (Obsrvese que el problema que nos ocupa es elptico).
5. Mallado. Pulsamos en la tecla
, para realizar una triangulacin. Obtenemos:
Pulsando la tecla: hacemos una triangulacin ms fina, es decir, cuantas ms veces la pulsamos ms fino queda el mallado. 6. Solucin. Para obtener la solucin del problema hemos de pulsar en la tecla pantalla: y obtenemos la
7. Grfica de la Solucin. Para dibujar la solucin pulsamos la tecla y obtenemos la siguiente pantalla:
Como se puede observar, hemos marcado que nos dibuje la solucin en tres dimensiones, en color y que muestre el mallado. Una vez seleccionado esto, pulsamos y obtenemos el dibujo deseado.
8. Exportacin de la Triangulacin y la Solucin. Exportamos a MATLAB la malla. Para ello pulsamos obteniendo la siguiente pantalla:
y ah Export Mesh,
Donde: p es una matriz de tamao 2xN que contiene las coordenadas de los nodos de la malla, es decir,
x = p(1,:) y = p(2, :)
e guarda informacin sobre las aristas. t almacena informacin de los tringulos (vrtices de los tringulos, subdominos). Para exportar la solucin obtenida mediante la herramienta pdetool pulsamos en y ah Export Solution, obteniendo la pantalla:
Donde: u guarda la solucin en cada nodo. Matriz de dimensin Nx1. Una vez hecho esto, en la pantalla de trabajo de MATLAB, como conocemos la solucin exacta, hacemos un estudio del error cometido en norma infinito y dibujamos la solucin exacta con el comando pdesurf. >> whos Name ans e p t u Size 1x8 7x32 2x146 4x258 146x1 Bytes Class 16 char array 1792 double array 2336 double array 8256 double array 1168 double array
Grand total is 1702 elements using 13568 bytes >> x=p(1,:); >> y=p(2,:); >> uexact=(x.^2+y.^2)'; >> norm(uexact-u,inf) ans = 0.0047 >> pdesurf(p,t,uexact)
Observacin.
La herramienta pdetool de MATLAB, solo sirve para problemas con dominios bidimensionales.
