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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Se propone una aproximación <strong>de</strong> la forma<br />

u ≃ û = ψ +<br />

M∑<br />

a m φ m (1.54)<br />

en la que se han elegido a las funciones ψ y φ m para que se satisfagan las condiciones <strong>de</strong> bor<strong>de</strong><br />

esenciales, es <strong>de</strong>cir, ψ = u _ y φ m = 0, (m = 1, 2...M) en Γ u . El problema <strong>de</strong> <strong>residuos</strong> pon<strong>de</strong>rados<br />

correspondiente resulta<br />

∫<br />

Ω W l<br />

[ ( ∂<br />

k ∂û )<br />

∂x ∂x<br />

+ ∂ (<br />

∂y<br />

k ∂û<br />

∂y<br />

m=1<br />

)]<br />

_dx dy + ∫ Ω W l F dx dy+<br />

+ ∫ (<br />

__<br />

Γ σ<br />

W l k ∂û )<br />

∂n + σ<br />

_ dΓ + ∫ (<br />

Γ c<br />

˜W l k ∂û<br />

)<br />

∂n + p (u − u ∞) dΓ = 0<br />

l = 1, 2...M<br />

(1.55)<br />

que para M → ∞ asegura la satisfacción <strong>de</strong> la ecuación diferencial en Ω y <strong>de</strong> las condiciones <strong>de</strong><br />

bor<strong>de</strong> naturales. La primer integral en la ec.(1.55) pue<strong>de</strong> ser ”<strong>de</strong>bilizada” utilizando el lema <strong>de</strong><br />

Green que establece las siguientes i<strong>de</strong>ntida<strong>de</strong>s para las funciones α y β<br />

∫<br />

α ∂β<br />

Ω ∂x<br />

∫Ω<br />

dx dy = − ∂α<br />

∂x<br />

∫Γ<br />

β dx dy + α β n x dΓ (1.56)<br />

∫<br />

Ω<br />

α ∂β<br />

∂y<br />

∫Ω<br />

dx dy = − ∂α<br />

∂y<br />

∫Γ<br />

β dx dy + α β n y dΓ (1.57)<br />

en la que n x y n y son los cosenos directores <strong>de</strong> la normal n (saliente) al contorno cerrado Γ que<br />

ro<strong>de</strong>a al dominio Ω en el plano x y. La integración sobre Γ se realiza en sentido antihorario.<br />

Utilizando estas i<strong>de</strong>ntida<strong>de</strong>s y notando que<br />

∂α<br />

∂n = ∂α<br />

∂x n x + ∂α<br />

∂y n y (1.58)<br />

la ec. (1.55) pue<strong>de</strong> reescribirse<br />

∫ ( ∂Wl<br />

−<br />

Ω ∂x k ∂û<br />

∂x + ∂W l<br />

∂y k ∂û ) ∫<br />

_dx dy + W l F dx dy+<br />

∂y<br />

Ω<br />

(<br />

+ W l k<br />

∫Γ ∂û )<br />

(<br />

__<br />

dΓ + W l k<br />

σ +Γ u +Γ c<br />

∂n<br />

∫Γ ∂û )<br />

σ<br />

∂n + σ<br />

_ dΓ+<br />

(<br />

+ ˜W l k<br />

∫Γ ∂û<br />

)<br />

c<br />

∂n + p (u − u ∞) dΓ = 0 (1.59)<br />

Limitando ahora la elección <strong>de</strong> las funciones <strong>de</strong> peso <strong>de</strong> tal modo que<br />

W l = 0 en Γ u ;<br />

__<br />

W l = −W l en Γ σ ; y ˜W l = −W l en Γ c (1.60)<br />

se observa que el término que contiene al gradiente <strong>de</strong> û <strong>de</strong>saparece y la ec.(1.59) queda<br />

∫ ( ∂Wl<br />

∂x k ∂û<br />

∂x + ∂W l<br />

∂y k ∂û ) ∫<br />

dx dy − W l F dx dy+<br />

∂y<br />

16<br />

Ω<br />

Ω<br />

∫<br />

_<br />

+ W l σdΓ + W l p (u − u ∞ ) dΓ = 0 (1.61)<br />

∫Γ σ Γ c

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

Saved successfully!

Ooh no, something went wrong!