11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

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.

Time-Dependent Problems 229<br />

We now wr<strong>it</strong>e Pn(x,t) = n−1 −x2<br />

m=0 cm(t)Hm(x)e , where the cm(t)’s are the<br />

Fourier coefficients of Pn at t ∈]0, ˆ T]. Subst<strong>it</strong>uting Pn in (10.2.22) and recalling<br />

(1.7.1), we obtain<br />

(10.2.24)<br />

d<br />

dt cm(t) = −2ζ(m + 1)cm(t), 0 ≤ m ≤ n − 1, ∀t ∈]0, ˆ T].<br />

Now, we can find the explic<strong>it</strong> expression of the coefficients. We note that, from (10.2.23),<br />

the in<strong>it</strong>ial datum is Pn(x,0) = (I∗ v,nU0)( √ 4ζx), v(x) := ex2, x ∈ R (see (6.7.6)).<br />

Another scheme is obtained w<strong>it</strong>h the in<strong>it</strong>ial cond<strong>it</strong>ion Pn(x,0) = (Π ∗ v,nU0)( √ 4ζx),<br />

x ∈ R (see (6.7.3)). This amounts to assigning the values cm(0), 0 ≤ m ≤ n −1, as the<br />

Fourier coefficients of U0. As the reader can easily check, after this modification, the<br />

solution to (10.2.22) turns out to be the approximation of V by the Galerkin method.<br />

Indeed, Pn also satisfies<br />

(10.2.25)<br />

<br />

R<br />

<br />

∂Pn<br />

∂(Pnv) ∂φ<br />

φv dx = −ζ<br />

∂t R ∂x ∂x dx, ∀φ ∈ Sn−1, ∀t ∈]0, ˆ T].<br />

10.3 Approximation of linear first-order problems<br />

We now focus the attention to another time-dependent partial differential equation. The<br />

unknown U : [−1,1] × [0,T] → R, satisfies<br />

(10.3.1)<br />

∂U ∂U<br />

(x,t) = −ζ (x,t), ∀x ∈] − 1,1], ∀t ∈]0,T],<br />

∂t ∂x<br />

(10.3.2) U(x,0) = U0(x), ∀x ∈ [−1,1],<br />

(10.3.3) U(−1,t) = σ(t), ∀t ∈]0,T],

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

Saved successfully!

Ooh no, something went wrong!