6 sep. 2020 20:35:49 - TP4-Ej.2.

sm

Trabajo práctico 4 - Ejercicio 2

l 1 Q 1

Ecuacion de poisson unidimensional: d d

u x Q 0
dx dx

u 0 0
Condiciones de contorno:
u l 0

1) Resolución con aproximación de primer orden

Se propone la siguiente función de aproximación de primer orden:

2
ü1 x a0 a1 x a2 x

Aplicando las condiciones de contorno:

ü1 0 0 --> a0 0

ü1 l 0 --> a1 a2 l

Entonces:

ü1 x a2 x 1 x

El residuo es:

d d d
R1 x ü x Q ü x a2 1 2 x
dx dx 1 dx 1

d
R1 x a 1 2 x Q d
dx 2 a 1 2 x 2 a2
dx 2
R1 x 2 a2 Q

R1 x 2 a2 1

El funcional de residuos ponderados:

l l
R1 x Wi d x 0 2 a2 1 Wi d x 0
0 0

1/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

Ahora por los diferentes métodos, se calcula el valor de a2 y la función asociada:

*Método de colocación puntual:

se considera x1 0,5

1
R1 x1 2 a2 1 0 --> a2
2

1 1 2 x 1 x
ü1CP x l x x
2 2 2

*Método de colocación por subdominios:

1
1
R1 x d x 0 --> 2 a2 1 0 --> a2
2
0

1 1 2 x 1 x
ü1CS x l x x
2 2 2

*Método de cuadrados mínimos:

d
Wi x R x 2
d a2 1

1
2 a2 1 2d x 0 1
--> a2
2
0

1 1 2 x 1 x
ü1CM x l x x
2 2 2

*Método de Garlekin:

d 2
Wi x ü x l x x
d ai

1
2 1
2 a2 1 l x x dx 0 --> a2
2
0

1 1 2 x 1 x
ü1G x l x x
2 2 2

2/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

Para todos los métodos analizados, se obtiene el mismo valor de a2 y de función:

1
a2
2

x 1 x
ü1 x
2

2
x x
uan1 x Solución exacta
2 2
2
x x
ucp1 x
2 2
2
x x
ucs1 x
2 2

2
x x
umc1 x
2 2

2
x x
ug1 x
2 2

0,3
analítico
0,2 colocación putual
colocación por subdominios
0,1
mínimos cuadrados

0 galerkin

-0,1
y

-0,2

-0,3

-0,4

-0,5

-0,6
-1,2
-1,1-1-0,9
-0,8
-0,7
-0,6
-0,5
-0,4
-0,3
-0,2
-0,10 0,10,20,30,40,50,60,70,80,9 1 1,11,21,31,41,51,61,71,81,9 2 2,12,22,32,42,52,62,72,82,9 3 3,13,23,3

x
uan1 x
ucp1 x
ucs1 x
umc1 x
ug1 x

3/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

2) Resolución con aproximación de segundo orden

Se propone la siguiente función de aproximación de segundo orden:

2 3
ü2 x b0 b1 x b2 x b3 x

Aplicando las condiciones de contorno:

ü2 0 0 --> b0 0

ü2 l 0 --> 2
b1 b2 l b3 l

Entonces:

2 2 3
ü2 x b2 l x b3 l x b2 x b3 x

El residuo es:

d
ü x b3 1 x2 b2 1 x x b2 2 x b3
d x 2
d d
R2 x ü x Q
d x d x 2
d
b3 1 x2 b2 1 x x b2 2 x b3 2 b2 2 x b3 b3 x
d x

d
R2 x b3 1 x2 b2 1 x x b2 2 x b3 Q
d x

R2 x 2 b2 2 x b3 b3 x Q

R2 x 2 b2 6 b3 x 1

El funcional de residuos ponderados:

l l
R2 x Wi d x 0 2 b2 6 b3 x 1 Wi d x 0
0 0

Ahora por los diferentes métodos, se calcula el valor de a2 y la función asociada:

*Método de colocación puntual:

se considera: 1 2
x1 x2
3 3

1
R2 x1 2 b2 6 b3 x 1 0 --> b2 b3
2

R2 x2 ; b2 2 b2 6 b3 x 1 0 --> b3 0

4/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

1
Entonces: b2 b3 0
2

1 1 2 x 1 x
ü2CP x l x x
2 2 2

*Método de colocación por subdominios:

1
3
1 1
R2 x d x 0 --> b2 b
2 2 3
0

Ahora resolviendo la siguiente integral, considerando el valor de b2 obtenido antes:

2
3
R2 x d x 0 --> b3 0

1
Entonces: b2 b3 0
2

1 1 2 x 1 x
ü2CS x l x x
2 2 2

*Método de cuadrados mínimos:

d
Wi_b2 x R x 2
d b2 2

1
3
1 1
2 b2 6 b3 x 1 2d x 0 --> b2 b
2 2 3
0

Ahora resolviendo la siguiente integral, considerando el valor de b2 obtenido antes:

d
Wi_b3 x R x 6 x
d b3 2

2
3
2 b2 6 b3 x 1 6 xd x 0 --> b3 0

1
Entonces: b2 b3 0
2

5/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

1 1 2 x 1 x
ü2CM x l x x
2 2 2

*Método de Garlekin:

d 2
Wi_b2 x ü x l x x
d b2

1
3
2 1 9
2 b2 6 b3 x 1 l x x dx 0 --> b2 b
2 14 3
0

Ahora resolviendo la siguiente integral, considerando el valor de b2 obtenido antes:

d 2 3
Wi_b3 x ü x l x x
d b3

1
3
2 3 b3 0
2 b2 6 b3 x 1 l x x dx 0 -->
0

1
Entonces: b2 b3 0
2

1 1 2 x 1 x
ü2G x l x x
2 2 2

Para todos los métodos analizados, se obtiene el mismo valor de a2 y de función:

1 b3 0
b2
2

x 1 x
ü2 x
2

2
x x
uan2 x Solución exacta
2 2

2
x x
ucp2 x
2 2

2
x x
ucs2 x
2 2

2
x x
umc2 x
2 2

2
x x
ug2 x
2 2

6/7
 6 sep. 2020 20:35:49 - TP4-Ej.2.sm

0,2
analítico
0,1
colocación putual
0 colocación por subdominios
mínimos cuadrados
-0,1
galerkin
-0,2

-0,3
y

-0,4

-0,5

-0,6

-0,7

-0,8

-1,7
-1,6
-1,5
-1,4
-1,3
-1,2
-1,1-1-0,9
-0,8
-0,7
-0,6
-0,5
-0,4
-0,3
-0,2
-0,100,10,20,30,40,50,60,70,80,911,11,21,31,41,51,61,71,81,922,12,22,32,42,52,62,72,82,933,13,23,33,43,53,63,7

x
uan2 x
ucp2 x
ucs2 x
umc2 x
ug2 x

7/7