28.10.2014 Views

Cálculo vectorial discreto para esquemas en diferencias sobre ...

Cálculo vectorial discreto para esquemas en diferencias sobre ...

Cálculo vectorial discreto para esquemas en diferencias sobre ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>en</strong> la id<strong>en</strong>tidad (2) y consideramos los escalares<br />

a 1 jj = h2<br />

2 (α j + α 3+j + 2(β j + β 3+j ) +<br />

3∑<br />

(γ ji + γ ij ), j = 1, 2, 3;<br />

i=1<br />

i≠j<br />

b 1 = h 2 (α 4 − α 1 + 2(β 4 − β 1 ) + γ 13 − γ 31 + γ 21 − γ 12 ,<br />

b 2 = h 2 (α 5 − α 2 + 2(β 5 − β 2 ) + γ 21 − γ 12 + γ 32 − γ 23 ,<br />

b 3 = h 2 (α 6 − α 3 + 2(β 6 − β 3 ) + γ 31 − γ 13 + γ 32 − γ 23 ,<br />

y las matrices de ord<strong>en</strong> 3, A 1 = h 2 ⎡<br />

⎡<br />

y A 3 = −h 2 ⎢<br />

⎣<br />

⎢<br />

⎣<br />

a 1 11 0 −γ 13<br />

0 a 1 22 0<br />

−γ 31 0 a 1 33<br />

⎤<br />

⎡<br />

⎥<br />

⎦, A 2 = −h 2 ⎢<br />

⎣<br />

⎤<br />

2β 4 γ 21 0<br />

⎥<br />

γ 21 2β 5 γ 32 ⎦<br />

0 γ 32 2β 6<br />

⎤<br />

2β 1 γ 12 0<br />

⎥<br />

γ 12 2β 2 γ 23 ⎦, <strong>en</strong>tonces L h (u) = −div (A∇u) + 〈b, ∇u〉 + qu es un<br />

0 γ 23 2β 3<br />

operador <strong>en</strong> difer<strong>en</strong>cias que satisface la propiedad requerida.<br />

El sigui<strong>en</strong>te resultado determina propiedades estructurales del esquema <strong>en</strong> términos<br />

de los coefici<strong>en</strong>tes del operador L h . Para ello mant<strong>en</strong>dremos las notaciones anteriores<br />

y <strong>para</strong> cada j = 1, . . . , 6 d<strong>en</strong>otaremos por r j a la suma de la j-ésima fila de la matriz<br />

A.<br />

Proposición. Se satisfac<strong>en</strong> las sigui<strong>en</strong>tes propiedades:<br />

i) L h es un esquema casi-simétrico sii A 1 = A t 1 y A 2 = A 3 . Además, L h es un<br />

esquema simétrico sii se satisfac<strong>en</strong> las igualdades anteriores y b = 0.<br />

ii) L h es un esquema no negativo sii q ≥ 0, A es una Z-matriz y r j ≥ −hb j ,<br />

r 3+j ≥ hb j , j = 1, 2, 3, <strong>para</strong> h sufici<strong>en</strong>tem<strong>en</strong>te pequeño. En particular, cuando L h<br />

es un esquema no negativo se debe verificar que a 1 jj ≥ h|b j |, j = 1, 2, 3, <strong>para</strong> h<br />

sufici<strong>en</strong>tem<strong>en</strong>te pequeño.<br />

iii) Si lím<br />

h→0<br />

hb = 0, <strong>en</strong>tonces L h es de tipo positivo sii q ≥ 0 y A es una M-matriz<br />

diagonalm<strong>en</strong>te dominante de forma estricta.<br />

iv) Si b = 0, <strong>en</strong>tonces L h es un esquema de tipo no negativo sii q ≥ 0 y A es una<br />

M-matriz diagonalm<strong>en</strong>te dominante.<br />

Agradecimi<strong>en</strong>tos<br />

Trabajo parcialm<strong>en</strong>te financiado por la Comisión Interministerial de Ci<strong>en</strong>cia y Tecnología,<br />

mediante el proyecto BFM2000-1063 y por la ETSECCPB.<br />

Refer<strong>en</strong>cias<br />

[1] E. B<strong>en</strong>dito, A. Carmona and A.M. Encinas,“Solving Boundary Value Problems on<br />

Networks Using Equilibrium Measures”, J. Func. Anal., 171 (2000), 155-176.

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

Saved successfully!

Ooh no, something went wrong!