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.

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!