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

efn.uncor.edu
from efn.uncor.edu More from this publisher
15.07.2014 Views

2. Las condiciones de borde no pueden construirse en forma arbitraria, sino que deben ser compatibles con la ecuación diferencial que gobierna el problema. En este sentido la formulación variacional posee la ventaja de que cualquier condición de borde que pueda incorporarse naturalmente al problema, como consecuencia de la integración por partes, será automaticamente compatible con la ecuación diferencial. 3. Recordemos que las condiciones de borde entran en la definición del problema variacional (3.18) en dos formas: las condiciones de borde esenciales, que son la especificación del valor de la solución, entran en el problema a partir de la definición del espacio admisible de las funciones del problema, mientras que las condiciones naturales, que conllevan la especificación de las derivadas de la solución, determinan la forma de la ecuación variacional. Ejemplos: i) −u(x), xx +u(x) = f(x) 0 < x < L u(0) = 0, u(L) = 0 ii) −u(x), xx +u(x) = f(x) 0 < x < L u(0), x = γ 0 , Los enunciados variacionales de estos problemas son u(L), x = γ L V-i) Encontrar u en H0 1 tal que (es decir que u satisface las condiciones (3.19) y además u(0) = u(L) = 0) L∫ L∫ (W ,x u ,x + W u) dx = W fdx, ∀ W ∈ H0 1 0 0 V-ii) Encontrar u ∈ H 1 tal que L∫ (W ,x u ,x + W u) dx = L∫ W fdx − γ 0 W (0) + γ L W (L), ∀ W ∈ H 1 0 0 Estos ejemplos muestran que las condiciones esenciales se introducen en V-i) a partir de la definición de la clase de funciones admisibles para proponer las funciones de prueba y las naturales como datos en el segundo miembro de la ecuación en V-ii). 4. Por último, en analogía con la energía U definida anteriormente, podemos escribir la norma energética para el problema (3.18) (con b(x) > 0) [∫ L ‖e‖ E = También puede utilizarse la norma H 1 equivalente 0 [∫ L ‖e‖ 1 = 0 ( k e 2 ,x + b e 2) ] 1/2 dx (3.20) ( e 2 ,x + e 2) ] 1/2 dx (3.21) Ejercicio N ◦ 26: Construir la forma débil indicando en cada caso el espacio de funciones de prueba admisibles. 44

a) b) c) − (k(x) u (x) , x ) , x +u(x), x +u (x) = 0 en 0 < x < 1; 1 < x < 2; 2 < x < 3; 3 < x < 4; con [k (1) u ,x (1)] = [k (2) u ,x (2)] = 0; [k (3) u ,x (3)] = 10 y ⎧ ⎨ k(x) ⎩ 1 0 ≤ x < 1 2 1 ≤ x < 2 1 2 < x < 4 ; ⎧ ⎨ ⎩ u (0) = 0 u (4) , x = 3 −u (x) ,xx + u (x) = δ (x − 1) 0 < x < 2 u (0) ,x = 2, u (2) ,x + u (2) = 3 u (x) ,xx + u (x) = 0 0 < x < 1; 1 < x < 2 [u ,x (1)] = 1; u (0) = u (2) = 0 Ejercicio N ◦ 27: los problemas de valores de contorno siguientes están mal condicionados o no satisfacen las hipótesis descriptas en las notas. ?‘Qué está mal en ellos? a) 1 (x − 1) u (x) ,xx + u (x) = sen (x) 0 < x < 2 b) c) d) y u (0) = u (2) = 0 (1 − x 2 ) u (x) ,xx + u (x) ,x = 3 0 < x < 2 y u (0) = 0 ; u (2) ,x = 1 u (x) ,xx + x u (x) ,x = 3 0 < x < 1 y u (0) ,xx = 0 ; u (1) = 0 −u (x) ,xx + u (x) ,x = x 0 < x < 1 y 1 2 u (0) ,x + u (1) ,x = u (1) ; u (0) ,x + 2 u (1) ,x = u (0) 3.4. Funciones de aproximación basadas en la interpolación de Lagrange En esta sección extenderemos las ideas respecto a las funciones de interpolación elementales introduciendo una técnica que conduce a lo que se conoce como elementos finitos Lagrangeanos, nombre que deriva del conceptos de interpolación utilizando polinomios de Lagrange: 1. Consideremos un elemento típico Ω e , aislado de la malla, y establezcamos un sistema local de coordenadas ξ, cuyo origen se ubica en el centro del elemento, escalado de tal forma de satisfacer que ξ = −1 en el extremo izquierdo y ξ = 1 en el derecho, Fig.2.a. Esto se realiza mediante la transformación general que para el caso de un elemento de k + 1 nodos se escribe ξ(x) = 2 x − xe 1 − x e k+1 h e (3.22) 45

a)<br />

b)<br />

c)<br />

− (k(x) u (x) , x ) , x +u(x), x +u (x) = 0<br />

en<br />

0 < x < 1; 1 < x < 2; 2 < x < 3; 3 < x < 4;<br />

con<br />

[k (1) u ,x (1)] = [k (2) u ,x (2)] = 0; [k (3) u ,x (3)] = 10<br />

y<br />

⎧<br />

⎨<br />

k(x)<br />

⎩<br />

1 0 ≤ x < 1<br />

2 1 ≤ x < 2<br />

1 2 < x < 4<br />

;<br />

⎧<br />

⎨<br />

⎩<br />

u (0) = 0<br />

u (4) , x = 3<br />

−u (x) ,xx<br />

+ u (x) = δ (x − 1) 0 < x < 2<br />

u (0) ,x<br />

= 2, u (2) ,x<br />

+ u (2) = 3<br />

u (x) ,xx<br />

+ u (x) = 0 0 < x < 1; 1 < x < 2<br />

[u ,x (1)] = 1; u (0) = u (2) = 0<br />

Ejercicio N ◦ 27: los problemas <strong>de</strong> valores <strong>de</strong> contorno siguientes están mal condicionados o no<br />

satisfacen las hipótesis <strong>de</strong>scriptas en las notas. ?‘Qué está mal en ellos?<br />

a)<br />

1<br />

(x − 1) u (x) ,xx<br />

+ u (x) = sen (x) 0 < x < 2<br />

b)<br />

c)<br />

d)<br />

y u (0) = u (2) = 0<br />

(1 − x 2 ) u (x) ,xx<br />

+ u (x) ,x<br />

= 3 0 < x < 2<br />

y u (0) = 0 ; u (2) ,x<br />

= 1<br />

u (x) ,xx<br />

+ x u (x) ,x<br />

= 3 0 < x < 1<br />

y u (0) ,xx<br />

= 0 ; u (1) = 0<br />

−u (x) ,xx<br />

+ u (x) ,x<br />

= x 0 < x < 1<br />

y<br />

1<br />

2 u (0) ,x + u (1) ,x = u (1) ; u (0) ,x + 2 u (1) ,x = u (0)<br />

3.4. <strong>Funciones</strong> <strong>de</strong> aproximación basadas en la interpolación <strong>de</strong><br />

Lagrange<br />

En esta sección exten<strong>de</strong>remos las i<strong>de</strong>as respecto a las funciones <strong>de</strong> interpolación elementales<br />

introduciendo una técnica que conduce a lo que se conoce como elementos finitos Lagrangeanos,<br />

nombre que <strong>de</strong>riva <strong>de</strong>l conceptos <strong>de</strong> interpolación utilizando polinomios <strong>de</strong> Lagrange:<br />

1. Consi<strong>de</strong>remos un elemento típico Ω e , aislado <strong>de</strong> la malla, y establezcamos un sistema local <strong>de</strong><br />

coor<strong>de</strong>nadas ξ, cuyo origen se ubica en el centro <strong>de</strong>l elemento, escalado <strong>de</strong> tal forma <strong>de</strong> satisfacer<br />

que ξ = −1 en el extremo izquierdo y ξ = 1 en el <strong>de</strong>recho, Fig.2.a. Esto se realiza mediante la<br />

transformación general que para el caso <strong>de</strong> un elemento <strong>de</strong> k + 1 nodos se escribe<br />

ξ(x) = 2 x − xe 1 − x e k+1<br />

h e (3.22)<br />

45

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

Saved successfully!

Ooh no, something went wrong!