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 ...
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