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.

1. Condiciones <strong>de</strong> bor<strong>de</strong> esenciales, en las que se especifican los valores que adopta la variable <strong>de</strong><br />

estado u.<br />

2. Condiciones <strong>de</strong> bor<strong>de</strong> naturales en las que:<br />

a) se especifica el valor <strong>de</strong>l flujo o<br />

b) se especifica una combinación lineal <strong>de</strong>l flujo y <strong>de</strong> la variable <strong>de</strong> estado.<br />

Por último notemos que en el caso <strong>de</strong> problemas en los cuales se especifiquen condiciones <strong>de</strong><br />

bor<strong>de</strong> naturales en ambos extremos, <strong>de</strong>berá satisfacerse una condición global <strong>de</strong> conservación que<br />

está dada por<br />

σ 0 + σ L =<br />

∫ L<br />

0<br />

f (x) _dx (3.13)<br />

3.3. Formulación variacional <strong>de</strong>l problema<br />

Resulta claro que el enunciado clásico <strong>de</strong>l problema <strong>de</strong> contorno (3.1) y (3.2) implica condiciones<br />

sobre la solución que pue<strong>de</strong>n resultar muy difíciles <strong>de</strong> satisfacer aun para situaciones en las que<br />

los datos resulten razonables. Lo que se preten<strong>de</strong> con una formulación variacional es obtener un<br />

enunciado propicio para po<strong>de</strong>r aplicar el método <strong>de</strong> elementos finitos. Es <strong>de</strong>seable a<strong>de</strong>más que<br />

la formulación posea simetría y preferiblemente que las funciones <strong>de</strong> <strong>prueba</strong> y <strong>de</strong> peso posean el<br />

mismo grado <strong>de</strong> suavidad, es <strong>de</strong>cir, que el or<strong>de</strong>n <strong>de</strong> <strong>de</strong>rivación, que presenten ambas funciones en<br />

el enunciado <strong>de</strong>l problema, sea el mismo. Reescribamos nuestro problema<br />

− d [<br />

k(x) du(x) ]<br />

+ c(x) du(x) + b(x)u(x) = f(x)<br />

dx dx dx<br />

α 0<br />

du(0)<br />

dx<br />

x ∈ Ω i , i = 1, 2, 3, 4<br />

[<br />

k(x) du(x) ]<br />

= 0 en x = x 1<br />

dx<br />

⎫⎪ ⎬<br />

[<br />

− k(x) du(x) ]<br />

(3.14)<br />

=<br />

dx<br />

ˆf en x = x 2<br />

[<br />

k(x) du(x) ]<br />

= 0 en x = x 3<br />

dx<br />

+ β du(L)<br />

0u(0) = γ 0 , α L<br />

dx + β ⎪<br />

Lu(L) = γ ⎭<br />

L<br />

en don<strong>de</strong> las condiciones <strong>de</strong> bor<strong>de</strong> son la expresión general <strong>de</strong> alguna <strong>de</strong> las siguientes situaciones<br />

• u (0) = ū 0 ; u (L) = ū L<br />

• k(0) du(0)<br />

dx<br />

= ¯σ 0 ; −k(L) du(L)<br />

dx<br />

= ¯σ L<br />

• k(0) du(0)<br />

dx<br />

= p 0 [u(0) − u ∞ ] ; −k(L) du(L)<br />

dx<br />

= p 0 [u(L) − u ∞ ]<br />

Recor<strong>de</strong>mos que el dominio Ω se ha dividido en cuatro subdominios como se indica en la fig.1.<br />

Como en cada subdominio la función incógnita <strong>de</strong>be ser suave, po<strong>de</strong>mos construir la función <strong>de</strong><br />

error para cada uno <strong>de</strong> ellos<br />

42<br />

R Ωi = − [k u ,x ] ,x<br />

+ c u ,x + b u − f

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

Saved successfully!

Ooh no, something went wrong!