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.

212 Polynomial Approximation of Differential Equations<br />

When ǫ tends to zero, Uǫ converges to the discontinuous function: U0 ≡ 0 in [−1,1[,<br />

U0(1) = 1. For very small values of ǫ, various physical phenomena are related to the<br />

function Uǫ, which displays sharp derivatives in a neighborhood of the point x = 1. Due<br />

to this behavior, approximations of Uǫ by low degree global polynomials, are affected<br />

by oscillations. Nevertheless, we are not in the s<strong>it</strong>uation that characterizes the Gibbs<br />

phenomenon (see sections 6.2 and 6.8, as well gottlieb and orszag(1977), p.40). Ac-<br />

tually, since the solution is analytic for ǫ > 0, uniform convergence w<strong>it</strong>h an exponential<br />

rate is obtained w<strong>it</strong>h the usual techniques. Tau, Galerkin and collocation methods have<br />

been studied in canuto(1988). In particular, the tau method is equivalent to finding<br />

pǫ,n ∈ Pn, n ≥ 2, ǫ > 0, satisfying problem (7.3.5), where: A ≡ 0, B ≡ 1<br />

ǫ , q ≡ 0,<br />

σ1 = 0, σ2 = 1. Thus, w<strong>it</strong>h the notations of section 7.1, the corresponding linear system<br />

is<br />

(9.7.3)<br />

⎧<br />

⎨<br />

⎩<br />

−ǫ c (2)<br />

k + ck = 0 0 ≤ k ≤ n − 2,<br />

n<br />

k=0 ckuk(−1) = 0, n<br />

k=0 ckuk(1) = 1.<br />

The coefficients c (2)<br />

k , 0 ≤ k ≤ n, are given for ν := α = β in karageorgis and<br />

phillips(1989) (see also (7.1.9) and (7.1.10)). These are all pos<strong>it</strong>ive for ν > −1.<br />

Therefore, <strong>it</strong> is easy to deduce that ck > 0, 0 ≤ k ≤ n. This property is used in<br />

canuto(1988) to estimate the maximum norm of pǫ,n. For instance, from (1.3.9), we<br />

have<br />

(9.7.4) |pǫ,n(x)| =<br />

≤<br />

n<br />

k=0<br />

ck<br />

<br />

n + ν<br />

n<br />

=<br />

<br />

<br />

<br />

<br />

<br />

n<br />

<br />

<br />

<br />

ckuk(x) <br />

≤<br />

k=0<br />

n<br />

ck|uk(x)|<br />

k=0<br />

n<br />

ckuk(1) = pǫ,n(1) = 1, ∀x ∈ Ī, ∀ǫ > 0, ∀n ≥ 2.<br />

k=0<br />

In add<strong>it</strong>ion, using relation (1.4.8), in the Legendre case one gets<br />

(9.7.5) |pǫ,n(x)| ≤<br />

n<br />

k=0<br />

ck|Pk(x)| ≤ 1<br />

2 (1 + x2 )<br />

n<br />

k=0<br />

ckPk(1) = 1<br />

2 (1 + x2 ),<br />

∀x ∈ Ī, ∀ǫ > 0, ∀n ≥ 2.

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

Saved successfully!

Ooh no, something went wrong!