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.

Time-Dependent Problems 225<br />

The next step is to discretize (10.2.11) w<strong>it</strong>h respect to the variable t. This will be<br />

investigated in section 10.6. Let us first analyze the convergence of the solution of the<br />

semi-discrete problem (10.2.8), (10.2.9), (10.2.10) to the solution of problem (10.2.5),<br />

(10.2.6), (10.2.7). Arguing as in section 9.3, we can propose a variational formulation<br />

for equation (10.2.5) as follows:<br />

(10.2.12)<br />

1<br />

−1<br />

1<br />

∂V<br />

∂V ∂(φw)<br />

φw dx = −ζ<br />

∂t −1 ∂x ∂x<br />

dx +<br />

1<br />

−1<br />

fφw dx,<br />

∀φ ∈ X, ∀t ∈]0,T],<br />

where the solution and the test functions both belong to the space X ≡ H 1 0,w(I) (see<br />

(5.7.6)). Existence and uniqueness of weak solutions of (10.2.12) are discussed in lions<br />

and magenes(1972) for the case w ≡ 1.<br />

W<strong>it</strong>h the help of the quadrature formula (3.5.1), the approximating polynomials satisfy<br />

(10.2.13)<br />

n<br />

j=0<br />

+<br />

<br />

∂pn<br />

∂t φ<br />

<br />

(η (n)<br />

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

j<br />

n<br />

(fφ)(η (n)<br />

j=0<br />

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

1<br />

∂pn ∂(φw)<br />

= −ζ<br />

−1 ∂x ∂x<br />

j , ∀φ ∈ P 0 n, ∀t ∈]0,T].<br />

Interesting properties are obtained from this equation. Consider ν := α = β w<strong>it</strong>h<br />

−1 < ν ≤ 1. For any t ∈]0,T], setting φ(x) := pn(x,t), x ∈ I, one gets<br />

(10.2.14)<br />

=<br />

n<br />

(fpn)(η (n)<br />

j=0<br />

⎛<br />

1 d<br />

⎝<br />

2 dt<br />

n<br />

[pn(η (n)<br />

j=0<br />

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

j<br />

1<br />

≤<br />

2<br />

j ,t)] 2 ˜w (n)<br />

j<br />

⎞<br />

n<br />

[pn(η (n)<br />

j=0<br />

⎠ + ζ<br />

1<br />

−1<br />

j ,t)] 2 ˜w (n)<br />

j<br />

∂pn<br />

∂x<br />

1<br />

+<br />

2<br />

∂(pnw)<br />

∂x<br />

dx<br />

dx<br />

n<br />

[f(η (n)<br />

j=0<br />

j ,t)] 2 ˜w (n)<br />

The last inequal<strong>it</strong>y in (10.2.14) is derived from the relation ab ≤ 1<br />

2 (a2 +b 2 ), a,b ∈ R.<br />

j .

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

Saved successfully!

Ooh no, something went wrong!