15.07.2014 Views

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

I(a 1 , a 2 , a 3 , ..., a M ) = 1 2<br />

es <strong>de</strong>cir, que se <strong>de</strong>ben satisfacer las siguientes igualda<strong>de</strong>s<br />

Cumpliéndose según la ec. (1.12) que<br />

la ec.(1.25) conduce a<br />

∫<br />

Ω<br />

(u − û) 2 dΩ (1.24)<br />

∂I<br />

∂a l<br />

= 0, l = 1, 2, ...M (1.25)<br />

∫<br />

Ω<br />

∂û<br />

∂a l<br />

= φ l (1.26)<br />

φ l (u − û) dΩ = 0 (1.27)<br />

que en este caso resulta idéntica a la que se obtiene mediante el método <strong>de</strong> Galerkin.<br />

1.3.2. Ejercicios:<br />

Ejercicio N ◦ 1: Dada la función<br />

u(x) = 0,05 − 1,5x + 4x 2 − 1,5x 3 − 0,7x 4 con 0 ≤ x ≤ 1<br />

se <strong>de</strong>sea aproximarla utilizando los siguientes conjuntos <strong>de</strong> funciones <strong>de</strong> aproximación:<br />

a) φ m = x m (1 − x); m = 1, 2, ...<br />

b) φ m =sin(mπx); m = 1, 2, ...<br />

En ambos casos, la solución particular ψ(x) queda <strong>de</strong>finida por la recta que une los extremos<br />

<strong>de</strong> u(x).<br />

Para minimizar el error utilice:<br />

1. El método <strong>de</strong> colocación con un punto (x = 0,5) y con dos puntos (x = 0,3 y x = 0,7).<br />

2. El método <strong>de</strong> Galerkin con uno y dos términos.<br />

Saque conclusiones comparado los resultados obtenidos al usar las distintas aproximaciones.<br />

Ejercicio N ◦ 2: Un problema unidimensional <strong>de</strong> calor produjo los siguientes resultados:<br />

Distancia 0,0 0,2 0,4 0,6 0,8 1,0<br />

Temperatura 20 30 50 65 40 30<br />

Ajuste una curva suave a estos datos utilizando el método <strong>de</strong> Galerkin y un conjunto admisible<br />

<strong>de</strong> funciones <strong>de</strong> <strong>prueba</strong>.<br />

Ejercicio N ◦ 3: Resolver el ejercicio 1 utilizando el método <strong>de</strong> colocación en subdominios. Utilizar<br />

como funciones <strong>de</strong> peso W l dos subregiones ([0,0.5] y [0.5,1]) y como funciones <strong>de</strong> <strong>prueba</strong> las<br />

propuestas en 1.a).<br />

1.4. Aproximación a la solución <strong>de</strong> ecuaciones diferenciales<br />

1.4.1. Condiciones <strong>de</strong> bor<strong>de</strong> satisfechas por elección <strong>de</strong> las funciones <strong>de</strong> <strong>prueba</strong><br />

Consi<strong>de</strong>remos la ecuación general<br />

A (u) = L (u) + p = 0 en Ω (1.28)<br />

en la que L es un operador diferencial lineal y p es in<strong>de</strong>pendiente <strong>de</strong> u. Por ejemplo, en el caso<br />

<strong>de</strong> la ecuación diferencial <strong>de</strong>l calor (1.9) que representa el flujo bidimensional se tiene<br />

L (u) = ∂ (<br />

k ∂u )<br />

+ ∂ (<br />

k ∂u )<br />

∂x ∂x ∂y ∂y<br />

6

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!