11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Domain-Decompos<strong>it</strong>ion Methods 253<br />

Theorem 11.2.1 - Let Uk, 1 ≤ k ≤ m, be the solutions of problem (11.2.2) w<strong>it</strong>h<br />

σ1 = σ2 = 0. Let η (n)<br />

j , 0 ≤ j ≤ n, in (11.2.3), be the nodes of the Legendre (w ≡ 1)<br />

Gauss-Lobatto formula (3.5.1). Let pn,k, 1 ≤ k ≤ m, be the polynomials satisfying<br />

(11.2.4), (11.2.5), (11.2.7) w<strong>it</strong>h σ1 = σ2 = 0, and (11.2.8) w<strong>it</strong>h γk := n(n+1)<br />

sk+1−sk−1 ,<br />

1 ≤ k ≤ m − 1. Then, we have<br />

(11.2.9) lim<br />

n→+∞<br />

sk<br />

sk−1<br />

(Uk − pn,k) 2 dx +<br />

sk<br />

sk−1<br />

Proof - For any 1 ≤ k ≤ m, define the weights ˜w (n,k)<br />

j<br />

(U ′ k − p ′ n,k) 2 dx<br />

1<br />

2<br />

= 0, 1 ≤ k ≤ m.<br />

:= 1<br />

2 (sk −sk−1) ˜w (n)<br />

j , 0 ≤ j ≤ n<br />

(see (3.5.6)). Therefore, by (3.5.1) w<strong>it</strong>h w ≡ 1, we get the integration formula<br />

(11.2.10)<br />

sk<br />

sk−1<br />

p dx =<br />

n<br />

j=0<br />

p(θ (n,k)<br />

j ) ˜w (n,k)<br />

j<br />

∀p ∈ P2n−1, 1 ≤ k ≤ m.<br />

Let Xn denote the space of continuous functions φ in I, such that φ(s0) = φ(sm) = 0<br />

and such that the restriction of φ to any interval ¯ Sk is a polynomial of degree n. Then,<br />

by noting that γk[ ˜w (n,k)<br />

n<br />

satisfies<br />

(11.2.11)<br />

= −<br />

= −<br />

m<br />

<br />

n<br />

k=1<br />

i=1<br />

m<br />

<br />

k=1<br />

= −<br />

Sk<br />

m<br />

k=1<br />

+ ˜w (n,k+1)<br />

0 ] = 1, 1 ≤ k ≤ m − 1, the function πn ∈ Xn<br />

<br />

I<br />

π ′′<br />

nφ dx +<br />

sk<br />

sk−1<br />

[p ′′ n,kφ](θ (n,k)<br />

i ) ˜w (n,k)<br />

i<br />

=<br />

m<br />

<br />

n<br />

k=1<br />

i=1<br />

π ′ nφ ′ dx =<br />

m−1 <br />

k=1<br />

<br />

p ′′ n,kφ dx +<br />

<br />

+<br />

m−1 <br />

k=1<br />

[fφ](θ (n,k)<br />

i ) ˜w (n,k)<br />

i<br />

m<br />

<br />

k=1<br />

lim<br />

x→s −<br />

k<br />

m−1 <br />

k=1<br />

Sk<br />

π ′ nφ ′ dx<br />

(π ′ nφ)(x) − lim<br />

x→s +<br />

k<br />

(π ′ <br />

nφ)(x)<br />

<br />

p ′ n,kφ − p ′ <br />

n,k+1φ (sk)<br />

γk( ˜w (n,k)<br />

n<br />

<br />

+ ˜w (n,k+1)<br />

<br />

0 ) p ′ n,kφ − p ′ <br />

n,k+1φ (sk)<br />

=: Fn(φ), ∀φ ∈ Xn.

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

Saved successfully!

Ooh no, something went wrong!