v1

5 2
4
1 v2 Introducción al Método
S
de los Elementos Finitos
Parte 2
Formulación abstracta del MEF para problemas elípticos

Alberto Cardona, Víctor Fachinotti

Cimec-Intec (UNL/Conicet), Santa Fe, Argentina
 Espacios de Hilbert
1. En formulaciones variacionales para solución de BVP, trabajaremos con
espacios V más grandes que los vistos hasta ahora, llamados espacios de
Hilbert. Se dotará a V de varios productos escalares determinados por el BVP.

2. Sucesión de Cauchy: sucesión v1 , v2 ,… , vn ∈ V que verifica

∀ε > 0, ∃N ∈ / vi − v j < ε si i, j > N

La sucesión de Cauchy converge a v si v − vi → 0 cuando i → ∞

3. Espacio de Hilbert: Es un espacio lineal V, equipado de producto escalar y

norma asociada, completo (o sea que las sucesiones de Cauchy convergen
respecto de la norma del espacio).

Introducción al Método de los Elementos Finitos 2

 Espacios de Hilbert: ejemplos 1D
• Espacio L2(I) de funciones de cuadrado integrable en I=(a,b)

{
L 2 (I) = v : v está definida en I y ∫ v 2 dx < ∞
I
}
dotado del producto interno: (v, w) = ∫ vw dx
I

( )
1

∫ dx
1
= = ( v, v ) 2
2 2
y la norma: v L 2 (I)
v
I

La desigualdad de Cauchy: (v, w) ≤ v L 2 (I)

w L 2 (I)

−β
Ejemplo: v( x) = x , x ∈ I = (0,1)
⎧ ln x 0 = ∞
1
si β = 1/ 2
1
⎪⎪
= ∫ x −2 β dx = ⎨ x1− 2 β 1 ⎧
2
v ∞ si β > 1/ 2
=⎨
L 2 (0,1)
⎪
⎪⎩1 − 2 β 0 ⎩(1 − 2 β ) si β < 1/ 2 ⇒ v ∈ L 2 (0,1) si β <1/2
0 −1

Introducción al Método de los Elementos Finitos 3

 Espacios de Hilbert: ejemplos 1D (cont.)
• Espacio H1(I) en I=(a,b): H1 (I) = {v : v, v′ ∈ L 2 (I)}

dotado del producto interno: (v, w) H1 (I) = ∫ ( vw + v′w′ ) dx

I

{∫ ⎡⎣v + ( v′) ⎤⎦ dx}

1 1

= ⎡(v, v) H1 (I) ⎤
2
= 2 2 2
y la norma: v H1 (I) ⎣ ⎦
I

• Espacio H10 (I) en I=(a,b): H10 (I) = {v : v ∈ H1 (I) y v(a) = v(b) = 0}

⎧ −u '' = f ( x), x∈I
• Dado el problema modelo: (D) ⎨
⎩u (a) = u (b) = 0
luego: (V ) Hallar u ∈ H10 (I) / ( u′, v′ ) = ( f , v ) ∀v ∈ H10 (I)

• Nota:
– H10 (I) es más grande que el espacio V de funciones lineales a trozos que
veníamos usando.
– En MEF, la norma del error es simplemente la norma en H1(I).

Introducción al Método de los Elementos Finitos 4

 Espacios de Hilbert en Ω∈Rd, d=2,3
{
• L 2 (Ω) = v : v está definida en Ω y ∫
Ω
v2d Ω < ∞ }
– Producto escalar: (v, w) = ∫ vwd Ω
Ω

( ∫ v dΩ)
1

v = 2 2
– Norma:
Ω
⎧ ∂v ⎫
• H (Ω) = ⎨v : v ∈ L 2 (Ω), ∈ L 2 (Ω), i = 1,… , d ⎬
1

⎩ ∂xi ⎭
– Producto escalar: (v, w) H1 ( Ω ) = ∫ ( vw + ∇v ⋅∇w ) d Ω
Ω

= ⎡ ∫ ( v + ∇v ⋅∇v ) d Ω ⎤ = ⎡(v, v) H1 ( Ω ) ⎤
1 1
2 2 2
– Norma: v H1 ( Ω ) ⎣ Ω ⎦ ⎣ ⎦
• H10 (Ω) = {v : v ∈ H1 (Ω), v Γ = 0}

Introducción al Método de los Elementos Finitos 5

 Espacios de Hilbert en Ω∈Rd, d=2,3 (cont.)
• Dado el problema modelo (ec. de Poisson + CB homogéneas)
⎧−Δu = f en Ω
(D) ⎨
⎩ u=0 sobre Γ
luego: (V ) Hallar u ∈ H10 (Ω) / ∫Ω
∇u ⋅∇v d Ω = ∫ fv d Ω
Ω
∀v ∈ H10 (Ω)

a(u , v) = ( f , v)
o, equivalentemente:
1 1
(M ) Hallar u ∈ H10 (Ω) / a(u , u ) − ( f , u ) ≤ a(v, v) − ( f , v) ∀v ∈ H10 (Ω)
2 2
F(u ) ≤ F (v )
• Nota: (V) se dice la formulación débil de (D), y la solución de (V) es la
solución débil de (D). Ésta no es necesariamente una solución clásica de (D).
Para que lo sea, u debe ser suficientemente regular de modo que Δu esté
definida en el sentido clásico.

Introducción al Método de los Elementos Finitos 6

 Interpretación geométrica del MEF
Consideremos el BVP de difusión-reacción con CB homogéneas:
⎧−Δu + u = f en Ω
(D) ⎨
⎩ u=0 sobre Γ
(V ) Hallar u ∈ H10 (Ω) / ∫ ∇u ⋅∇v d Ω + ∫ uv d Ω = ∫ fv d Ω ∀v ∈ H10 (Ω)
Ω Ω Ω

u, v = ( f , v)
Sea Vh un subespacio de dimensión finita de H10 (Ω) . El MEF
aplicado al problema de difusión-reacción da: u
(Vh ) Hallar uh ∈ Vh / uh , v = ( f , v) ∀v ∈ Vh recta ≡ Vh
Tomando v ∈ Vh ⊂ H10 (Ω) en (V), (V)−(Vh) resulta:
uh
u − uh , v = 0 ∀v ∈ Vh
La solución MEF uh es la proyección con resp a 〈.,.〉 de plano ≡ H10 (Ω)
la solución exacta u sobre Vh, o sea que uh es el
elemento de Vh más próximo u con resp a ||.||H10 ( Ω ), i.e.:
u − uh H1 ( Ω ) ≤ u − v H1 ( Ω ) ∀v ∈ Vh
0 0
Introducción al Método de los Elementos Finitos 7
 CB naturales y esenciales
Consideremos el BVP de difusión-reacción con CB Neumann:
⎧−Δu + u = f en Ω
⎪
(D) ⎨ ∂u
⎪⎩ =g sobre Γ
∂n
(V ) Hallar u ∈ H1 (Ω) /
∫Ω ( ∇u ⋅∇v + uv ) d Ω = ∫Ω fv d Ω + ∫Γ gv d Γ ∀v ∈ H1 (Ω)

u, v = ( f , v) + ( g , v)Γ
1
(M ) Hallar u ∈ H1 (Ω) / F(u ) = a(u , u ) − ( f , u ) − ( g , u )Γ ≤ F(v) ∀v ∈ H1 (Ω)
2

• La CB Neumann, que no tiene que ser impuesta explícitamente sobre u, se

denomina CB natural.
• La CB Dirichlet u=u0 sobre Γ, que debe ser satisfecha explícitamente por u, se
conoce como CB esencial.
Introducción al Método de los Elementos Finitos 8
 Problema continuo en forma abstracta
• Objetivos:
– Dar un tratamiento unificado a muchos problemas de la Mecánica y la Física,
a fin de no repetir el mismo argumento en distintos casos concretos.
• Entender la estructura básica del MEF.

Los ingredientes de una formulación abstracta son:

1. Un espacio de Hilbert V, con producto escalar (.,.)V y norma ||.||V, donde se busca
la solución.
2. Una forma bilineal a:V×V→R, determinada por el problema (elíptico) a resolver,
t.q. ∀u,v∈V: a(u, v) = a(v, u ) (simetría)
∃γ > 0 / a(u , v) ≤ γ u V
v V
(continuidad)
∃α > 0 / a(v, v) ≥ α v
2
V
(coercividad o V-elipticidad)
3. Una forma lineal L:V→R, determinada por los datos del problema, continua, i.e.,
∃λ > 0 / L(v) ≤ λ v V
∀v ∈ V
Introducción al Método de los Elementos Finitos 9
Formulación abstracta del MEF para problemas elípticos (cont)

1
(M ) Hallar u ∈ V / F(u ) = min F(v), con F(v) = a(v, v) − L(v)
v∈V 2
(V ) Hallar u ∈ V / a(u , v) = L(v), ∀v ∈ V

Teorema: (M) y (V) son equivalentes, i.e., u satisface (M) si y sólo

si u satisface (V). Además, ∃!u∈V, y se verifica la estimación de
estabilidad
λ
u ≤
V
α

Introducción al Método de los Elementos Finitos 10

 Formulación abstracta del MEF para problemas elípticos (cont)
Demo. (M )⇒(V ): sea v∈V y ε∈R arbitrarios. Luego, v+εu∈V, así que
F(u ) ≤ F(u + ε v) , ∀ε ∈
g(0) ≤ g(ε ) ⇒ g tiene un mínimo en ε = 0
1
Siendo g(ε ) = a(u + ε v, u + ε v) − L(u + ε v)
2
1 ε ε ε2
= a(u , u ) + a(u , v) + a(v, u ) + a(v, v) − L(u ) − ε L(v)
2 2 2 2
1 ε2
= a(u, u ) − L(u ) + ε a(u, v) − ε L(v) + a(v, v)
2 2
luego g′(ε ) = a(u , v) − L(v) + ε a(v, v)
Como g tiene un mínimo en 0, debe cumplirse
g′(0) = a(u, v) − L(v) = 0 (V ) (QED)

Introducción al Método de los Elementos Finitos 11

 Formulación abstracta del MEF para problemas elípticos (cont)
Demo. (V )⇒(M ): dados u (solución) y w (arbitrario) en V, u+w∈V, así que
1
F(u + w) = a(u + w, u + w) − L(u + w)
2
1 1 1 1
= a(u , u ) + a(u, w) + a( w, u ) + a( w, w) − L(u ) − L( w)
2 2 2 2
1 1
= a(u, u ) − L(u ) + a(u , w) − L( w) + a( w, w) ≥ F(u ), (M) (QED)
2 = 0 (V )
2 ≥0
= F( u )

Introducción al Método de los Elementos Finitos 12

 Formulación abstracta del MEF para problemas elípticos (cont)
Demo. estimación de estabilidad:

a(u , v) = L(v) ∀v ∈ V
λ
α ||u||2V ≤ a(u, u ) = L(u ) ≤ λ||u||V ⇒ || u ||V ≤ (QED)
α
V-elipticidad de a Continuidad de L

Demo. unicidad: supongamos que u1,u2∈V sean soluciones de (V)

a(u1 , v) = L(v) ∀v ∈ V
a(u2 , v) = L(v) ∀v ∈ V
⇒ a(u1 − u2 , v) = 0 ∀v ∈ V

Definiendo u=u1−u2, llego a un nuevo BVP con L(v)=0 ∀ v∈V, i.e. λ=0. Luego:
Estimación de estabilidad ||u||V = 0 ⇒ u = u1 − u2 = 0 ⇒ u1 = u2 (QED)

Introducción al Método de los Elementos Finitos 13

 Estimación del error de discretización
• Sea Vh un subespacio de V, con dimensión finita M, y sea {ϕ1, ϕ2,…, ϕM}
una base para Vh, de modo que toda v∈Vh puede representarse
M
v = ∑ ϕiηi , ηi ∈
i =1

• Usando Vh, obtenemos los problemas discretos análogos a (M) y (V):

(M h ) Hallar u h ∈ Vh / F(uh ) ≤ F(v), ∀v ∈ Vh
(Vh ) Hallar uh ∈ Vh / a(uh , v) = L(v), ∀v ∈ Vh

– Como ϕ j ∈ Vh ⇒ a(uh , ϕ j ) = L(ϕ j ), j = 1, 2,… , M

M M
– Como uh ∈ Vh ⇒ uh = ∑ ϕi ξi , ξi ∈ ⇒ ∑ a(ϕi , ϕ j )ξi = L(ϕ j ), j = 1, 2,… , M
i =1 i =1

– Forma matricial de (Vh): Aξ = b con Aij = a(ϕi , ϕ j ), b j = L(ϕ j ).

Introducción al Método de los Elementos Finitos 14

 Estimación del error de discretización (cont)
• Dado que a(.,.) es simétrica: a(ϕi , ϕ j ) = a(ϕ j , ϕi ) ⇒ Aij = Aji
⇒ la matriz A es simétrica
• Dado que a(.,.) es V-elíptica:
M M M M
a(v, v) = a(∑ ϕiηi , ∑ ϕ jη j ) =∑∑ηi a(ϕi , ϕ j )η j = η ⋅ Aη ≥ α v > 0 si v ≠ 0 (η ≠ 0)
2
V
i =1 j =1 i =1 j =1

⇒ la matriz A es positiva definida

⇒ la matriz A es no singular
⇒ ∃!ξ / Aξ = b ⇒ ∃!uh ∈ Vh solución de (Vh )
λ
Además, uh ∈ Vh verifica la estimación de estabilidad ||uh ||V ≤
α
Demo.: a(uh , v) = L(v) ∀v ∈ Vh
α ||uh ||2V ≤ a(uh , uh ) = L(uh ) ≤ λ||uh ||V
V-elipticidad de a Continuidad de L

Introducción al Método de los Elementos Finitos 15

 Estimación del error de discretización (cont)
• Teorema: Sea u∈V la solución de (V) y uh∈Vh⊂V. Entonces:
γ
u − uh ≤ u−v ∀v ∈ Vh
V
α V

• Demo.:
1. u solución de (V ) ⇒ a(u , w) = L( w), ∀w ∈ Vh ⊂ V
2. uh solución de (Vh ) ⇒ a(uh , w) = L( w), ∀w ∈ Vh
3. Restando: a(u − uh , w) = 0, ∀w ∈ Vh
4. Sea w = uh − v, con v ∈ Vh arbitrario.
5. Operando: α u − uh V ≤ a(u − uh , u − uh )
2
V-elipticidad de a
= a(u − uh , u − uh ) + a(u − uh , w)
= a(u − uh , u − uh + w)
Continuidad de a = a(u − uh , u − v) ≤ γ u − uh V
u−v V (QED)

Introducción al Método de los Elementos Finitos 16

 Norma de energía
Considerando
1. a(.,.) es continua, i.e., ∃γ > 0 / a(u, v) ≤ γ u V
v V
∀u, v ∈ V
2. a(.,.) es V-elíptica o coerciva, i.e., ∃α > 0 / a(v, v) ≥ α v ∀v ∈ V
2
V

podemos definir una nueva norma llamada norma de energía:

v a = [ a(v, v) ] , ∀v ∈ V
1
2

Esta norma es equivalente a ||.||V, i.e., ∃ constantes positivas c = α , C = γ , tal que

c v V
≤ v a ≤C v V
∀v ∈ V

El producto escalar asociado a ||.||a es (u, v)a = a(u, v)

Verifica la desigualdad de Cauchy (u , v)a ≤|| u ||a || v ||a

Introducción al Método de los Elementos Finitos 17

 Estimación del error de discretización en norma de energía
• Teorema: Sea u∈V la solución de (V) y uh∈Vh⊂V. Entonces:
u − uh a
≤ u −v a
∀v ∈ Vh
⇒ Medida en la norma de energía, uh es la mejor aproximación a u.

• Demo.:
1. u solución de (V ) ⇒ (u , w)a = L( w), ∀w ∈ Vh ⊂ V
2. uh solución de (Vh ) ⇒ (uh , w)a = L( w), ∀w ∈ Vh
3. Restando: (u − uh , w)a = 0, ∀w ∈ Vh
4. Sea w = uh − v, con v ∈ Vh arbitrario.
5. Operando: u − uh a = (u − uh , u − uh )a
2

= (u − uh , u − uh )a + (u − uh , w)a
= (u − uh , u − uh + w)a
Desigualdad
de Cauchy
= (u − uh , u − v)a ≤ u − uh a
u−v a (QED)

Introducción al Método de los Elementos Finitos 18

 Ejemplos concretos
Ejemplo 1: sea V=H (Ω), Ω ⊂
1 2
.
a(u , v) = ∫ ( uv + ∇u ⋅∇v ) dx
Ω

L(v) = ∫ fv dx, f ∈ L 2 (Ω )
Ω

En este caso, a(u , v) = L(v), ∀v ∈ H1 (Ω) es la forma débil del problema de Neumann
∂u
(D) − Δu + u = f en Ω, = 0 sobre Γ
∂n
Se verifica ∀v, w ∈ H (Ω)
1

• a(v, w) = a(w, v) ⇒ a(.,.) es una forma bilineal simétrica

• a(v, v) = v ⇒ a(.,.) es V-elíptica con α=1
2
H1 ( Ω )

• a(v, w) ≤ [ a(v, v) ] [ a(w, w) ] = v H1 ( Ω ) w H1 ( Ω ) ⇒ a(.,.) es continua con γ=1

1 1
2 2

• L(v) = ∫Ω fv dx ≤ f L2 ( Ω ) v L2 ( Ω ) ⇒ L(.) es continua con λ = f L ( Ω )

2

– Estimación de estabilidad: u H1 ( Ω )
≤ f L2 ( Ω )

– Estimación de error: u − uh H (Ω)

1 ≤ u−v H (Ω)
1 ∀v ∈ H1 (Ω)

Introducción al Método de los Elementos Finitos 19

 Ejemplos concretos (cont)

Ejemplo 2: sea V=H10 (I), I = (0,1), a(v, w)= ∫ v′w′dx, L(v)= ∫ fv dx.
I I

En este caso, a(u , w) = L(v), ∀v ∈ H10 (I) es la forma débil del problema:
(D) − u ′′ = f en I, u (a ) = u (b) = 0
• Se verifica ∀v, w ∈ H10 (I)
a(v, w) = a(w, v) ⇒ a(.,.) es simétrica
a(v, w) ≤ v′ L (I) w′ L (I) ≤ v H1 (I) w H1 (I) ⇒ a(.,.) es continua con γ = 1
2 2

a(v, v) ≥ α v ⇒ a(.,.) es V-elíptica con α = 1 2

2
H1 (I)

L(v) = ∫ fv dx ≤ f L 2 (I)
v L2 (I)
⇒ L(.) es continua con λ = f L2 (I)
I

• Luego:

– Estimación de error: u − uh 1
H (I)
≤ 2 u−v 1
H (I)
∀v ∈ H10 (I)
– Estimación de estabilidad: u H1 (I)
≤ f L 2 (I)

Introducción al Método de los Elementos Finitos 20

 Ejemplos concretos (cont. Ejemplo 2)
Demo. de V-elipticidad de a(.,.):
=0 x
v( x) = v(0) + ∫ v′( y ) dy Desigualdad de Cauchy
0 1 1
x
⎛ 1
⎞ ⎛ ⎞ 1 2 1 2

v( x) ≤ ∫ v′( y ) dy ≤ ∫ 1 ⋅ v′ dy ≤ ⎜ ∫ 12 dy ⎟ ⎜ ∫ ( v′ ) dy ⎟
2

0 0 ⎝0 ⎠ ⎝0 ⎠
1 1 1 1

∫ v dx ≤ ∫ ∫ ( v′) dydx = ∫ ( v′) dx

2 2 2

0 0 0 0
1 1 1

∫ dx + ∫ ( v′) dx ≤ 2 ∫ ( v′) dx
2 2 2
v (QED)
0 0 0
||v|| a ( v ,v )
H1 ( 0,1)

Teorema: ∀v∈H1(Ω), ∃ cte. C >0 dependiente de Ω, t.q.

v
2
L2 ( Ω ) ≤C v ( 2
L2 ( Γ ) + ∇v
2
L2 ( Ω ) ) Desigualdad de
Poincaré-Friedrichs
• Para Ω≡Ι=(0,1), v(0)=0, resulta C =1.

Introducción al Método de los Elementos Finitos 21

 Ejemplos concretos (cont)
Ejemplo 3: sea V=H10 (Ω), Ω ⊂ 2
, a(u, v)= ∫ ∇u ⋅∇v d Ω, L(v)= ∫ fv d Ω.
Ω Ω
En este caso, a (u , v)=L(v), ∀v ∈ H10 (Ω) es la forma débil del problema de Poisson
(D) − Δu = f en Ω, u = 0 sobre Γ
• Se verifica ∀v, w ∈ H10 (Ω)
a (v, w) = a (w, v) ⇒ a(.,.) es simétrica
a(v, w) ≤ ∇v L2 ( Ω )
∇w L2 ( Ω )
≤ v H1 ( Ω )
w H1 ( Ω )
⇒ a(.,.) es continua con γ = 1
a(v, v) = ∫ ∇v ⋅∇v d Ω ≥ C −1 ∫ v 2 d Ω, C ∈ +
← Desigualdad de Poincaré-Friedrichs
Ω Ω

( C + 1) a(v, v) ≥ ∫ v 2 d Ω + ∫ ∇v ⋅∇v d Ω ⇒ a(.,.) es V-elíptica, con α = (C + 1) −1

Ω Ω
=1/ α
||v||2 1
H (Ω)

L(v) = ∫ fv d Ω ≤ f L2 ( Ω )
v L2 ( Ω )
⇒ L(.) es continua con λ = f L2 ( Ω )
Ω

• Estimación de error: u − uh H (Ω)

1 ≤ ( C + 1) u − v H (Ω)
1 ∀v ∈ H10 (Ω)

• Estimación de estabilidad: u H1 ( Ω )
≤ ( C + 1) f L2 ( Ω )

Introducción al Método de los Elementos Finitos 22

 Ejemplos concretos (cont)
d 4u
Ejemplo 4: ( D) 4
= f en I=(0,1), u (0) = u ′(0) = u (1) = u ′(1)
dx
• Definimos los espacios H 2 (I)= {v : v, v′, v′′ ∈ L 2 (I)}
H 02 (I)= {v : v ∈ H 2 (I), v(0) = v′(0) = v(1) = v′(1) = 0}

{∫ ⎡⎣v + ( v′) + ( v′′) ⎤⎦ dx}

1
2 2 2
2
con norma v H 2 (I)
=
I

• La forma débil del problema (D) consiste en hallar u ∈ V=H 02 (I) tal que
∫ ′′ ′′ ∫ ∀ ∈ 2
u v dx = fv dx v H 0 (I)
I I

a(u , v) = L(v)
• a(.,.) es una forma bilineal y simétrica
• a(v, w) ≤ v′ L 2 (I)
w′ L 2 (I)
≤ v H 2 (I)
w H 2 (I) ⇒ a(.,.) es continua con γ=1
• a(v, v) ≥ α v = α ∫I ⎣ ( ) ( ) ⎤⎦ dx ⇒ a(.,.) es V-elíptica con α=1/3
⎡ + ′ + ′′
2 2 2 2
H 2 (I)
v v v

• L(.) es continua con λ = f L (I) 2

• Luego: u 2
H (I)
≤3 f L 2 (I)
u − uh 2
H (I)
≤3 u−v 2
H (I)
∀v ∈ H 02 (I)
Introducción al Método de los Elementos Finitos 23
 Ejemplos concretos (cont)
Ejemplo 5: Consideremos el problema bi-armónico:
∂u
(D) Δ 2u = f en Ω ⊂ d , d = 1,2,3, u=
∂n
= 0 sobre Γ ( Δ u = Δ ( Δu ) )
2

• Definimos el espacio H k (Ω)= {v : v ∈ L 2 (Ω), Dα v ∈ L 2 (Ω),| α |≤ k}

con norma v Hk ( Ω )
= ( D α )2
∑ ∫ v dΩ
|α |≤ k Ω

∂|α |v
donde D v = α1 α 2 , α = (α1 , α 2 );| α | = α1 + α 2 ;α1 , α 2 = 0,1, 2,…
α

∂x1 ∂x2

• La forma débil del problema (D) resulta:

⎧ ∂v ⎫
Hallar u ∈ H (Ω)= ⎨v : v ∈ H (Ω), v =
2 2
= 0 sobre Γ ⎬ , tal que
∂n
0
⎩ =0
⎭ =0

∂ ∂v
(V ) ∫ fv d Ω = ∫ Δ 2u v d Ω = − ∫ ∇ ( Δu ) ⋅∇v d Ω + ∫ ( Δu ) v d Γ = ∫ Δu ⋅ Δ v d Ω − ∫ Δ u d Γ
Ω Ω Ω Γ ∂n Ω Γ ∂n
=L( v ) =a ( u , v )

Introducción al Método de los Elementos Finitos 24

 Ejemplos concretos (cont. Ejemplo 5)
• a(.,.) es una forma bilineal simétrica, pues a(u , v) = ∫Ω Δu ⋅ Δv d Ω = a(v, u )

• a(.,.) es continua con γ=1, pues a(u , v) ≤ Δu L2 ( Ω ) Δv L2 ( Ω ) ≤ Δu H2 ( Ω ) Δv H2 ( Ω )

• a(.,.) es V-elíptica con α=1/(1+C+C2). Usando la identidad :

2
⎛ ∂ ⎞ ∂
2
⎛ ⎞
d 2
v 1 v
Δv L2 ( Ω ) = ∫ (Δv) 2 d Ω = ∑ ∫ ⎜ ⎟ d Ω + ∫ ⎜ ⎟ d Γ , ∀v suave t.q. v Γ = 0.
2

Ω i , j =1 Ω ⎝ ∂xi ∂x j ⎠ Γ
R ⎝ ∂n ⎠
=0
2
incluimos las derivadas cruzadas que faltan en Δv
2
L2 ( Ω )
para lograr v H2 (Ω)
:
2
⎛ ∂ 2v ⎞
d
Δv L2 ( Ω ) + v L2 ( Ω ) + ∇v L2 ( Ω ) = ∑ ∫ ⎜ ⎟ dΩ + v + ∇v ≡ v
2 2 2 2 2 2

i , j =1 Ω ⎝ ∂xi ∂x j ⎠
L2 ( Ω ) L2 ( Ω ) H2 ( Ω )

Luego usamos Poincaré-Friedrichs :

≤ Δv + (C + 1) ∇v ≤ Δv + (C + 1)C Δv = (1 + C + C 2 ) Δv
2 2 2 2 2 2
v H2 ( Ω ) L2 ( Ω ) L2 ( Ω ) L2 ( Ω ) L2 ( Ω ) L2 ( Ω )

(QED)
• L(.) es continua con
λ= f L2 ( Ω )

Introducción al Método de los Elementos Finitos 25

 Ejemplos concretos (cont)
Ejemplo 6: Consideremos el problema estacionario de convección-difusión:
(D) − μΔu + β ⋅∇u = f en Ω, μ ∈ +
,β ∈ 2
, u = 0 sobre Γ
Supongamos ||β||/μ moderado.
• La forma débil del problema (D) se obtiene haciendo =0

∂u
(V ) ∫ fv d Ω = ∫ ( − μΔu + β ⋅∇u ) v d Ω = ∫ ⎡⎣ μ∇u ⋅∇v + ( β ⋅∇u ) v ⎤⎦ d Ω + ∫ μ v d Γ
Ω Ω Ω Γ ∂n
=L( v ) =a ( u , v )

• L(.) es una forma lineal continua.

• a(.,.) es una forma bilineal, continua y V-elíptica, pero no simétrica.

Teorema: si L(.) es una forma lineal continua, y a(.,.) es una forma bilineal, continua
y V-elíptica, pero no simétrica, se puede demostrar que hay solución única a (V), y
está acotada. Sin embargo, en este caso no existe problema de minimización
asociado a (V).

Introducción al Método de los Elementos Finitos 26

 Ejemplos concretos (cont)
Ejemplo 7: Consideremos el problema estacionario de conducción de calor en Ω⊂R3
⎧ ⎡ k1 0 0⎤
⎪
⎪−∇ ⋅ ( k∇u ) = f en Ω, con k = ⎢⎢ 0 k2 0 ⎥⎥ , ki ∈ +
Ecuación del calor
⎪
(D) ⎨ ⎢⎣ 0 0 k3 ⎥⎦
⎪u = 0 en Γ1 CB Dirichlet
⎪
⎪⎩k∇u ⋅ n = g en Γ 2 CB Neumann

Definimos como espacio para la solución V = {v : v ∈ H1 (Ω), v = 0 sobre Γ1}

• La forma débil del problema (D) se obtiene haciendo

∫ fv d Ω = − ∫ ∇ ⋅ ( k∇u ) v d Ω = ∫ k∇u ⋅∇v d Ω − ∫ k∇u ⋅ n v d Γ = ∫ k∇u ⋅∇v d Ω − ∫ g v d Γ

Ω Ω Ω Γ Ω Γ2

⇒ ∫ k∇u ⋅∇v d Ω = ∫ fv d Ω + ∫ g v d Γ
Ω Ω Γ2
a ( u ,v ) L( v )

• L(.) es una forma lineal continua si f ∈L2(Ω), g∈L2(Γ2).

• a(.,.) es una forma bilineal, simétrica, continua y V-elíptica si
∃c, C ∈ + / c ≤ ki (x) ≤ C ∀x ∈ Ω
Introducción al Método de los Elementos Finitos 27