Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate

24.02.2014 Views

80 J.Ch. Gilbert and P. Joly Therefore, since ψ m is an affine function, for any t ∈ [0, 1], there holds ∀x ≤ α m,k , 0 ≤ [ ψ m ( tR1 +(1−t)R 2 )] (x) ≤ 4. Hence ( ) α m tR1 +(1−t)R 2 = αm,k . In other words, any point of the segment [R 1 ,R 2 ] is a solution of (20), i.e., the set of solutions of (20) is convex. ⊓⊔ As a consequence of Lemmas 2 and 3, we know that any solution R of (20) is such that T R ≡{τ ∈ ]0,α m,k [ | [ψ m (R)] (τ) =0or4} is nonempty. Let us call tangent point an element of T R . Theorem 4 below is more precise, since it claims that there is at least M ≥ k tangent points τ j at which ψ m (R) takesalternatively the values 0 and 4. For any R, itisconvenient to construct and enumerate these tangent points in decreasing order: τ M+1 =0

Optimal Higher Order Time Discretizations 81 Proof. We proceed by contradiction, assuming that M ≤ k − 1. For j = 0,...,M − 1, one can find a point ( ) τ j+ 1 ∈ ]τ j+1,τ 2 j [ such that [ψ m (R)] ]τ j+1 ,τ j+ 1 ] ⊂ ]0, 4[. (26) 2 Consider the polynomial ϕ defined at x ∈ R by M−1 ∏ ( ) ϕ(x) =s 0 x − τj+ 1 . j=0 Hence ϕ ≡ s 0 if M = 0. This polynomial is of degree M ≤ k − 1, so that it is a possible increment to R in P k−1 .Fort>0, consider the polynomial p t = ψ m (R + tϕ), which verifies for all x ∈ R: p t (x) =[ψ m (R)](x)+tx m+1 ϕ(x). We shall get a contradiction and conclude the proof if we show that, for any small t>0, p t (x) ∈ ]0, 4[ for x ∈ ]0,α m,k ] (since then α m (R + tϕ) >α m,k and R would not be a local maximum). We shall only consider the case when [ψ m (R)] (α m,k ) = 4, since the reasoning is similar when [ψ m (R)] (α m,k )=0.Thens 0 = −1 by (25). • On the interval ]τ 1/2 ,α m,k ], ψ m (R) is greater than a positive constant (since it is positive on ]τ 1 ,α m,k ] by the definition of τ 1 and τ 1 0 small enough, p t ∈ ]0, 4[ on this interval. • Since ϕ(τ 1/2 ) = 0, p t (τ 1/2 ) = [ψ m (R)] (τ 1/2 ), which is in ]0, 4[ by the definition of τ 1/2 in (26). • On the interval ]τ 3/2 ,τ 1/2 [, ψ m (R) is less than a constant < 4(sinceit is < 4on]τ 2 ,τ 1/2 ] by the definition of τ 2 and τ 1/2 , see 2 and (26), and τ 2 0 small enough, p t (·) ∈ ]0, 4[ on the interval. • If s M < 0 then, on the considered interval, the tangent points are all at y = 4, ψ m (R) is positive, and ϕ is negative. Since the map x ↦→ [ψ m (R)] (x)/x =1+c 1 x + ... is greater than a positive constant on the considered interval, the map x ↦→ [ψ m (R)] (x)/x + tx m ϕ(x) =p t (x)/x is also positive on the interval for t>0 sufficiently small. It results that, for t>0 small enough, p t (·) ∈ ]0, 4[ on the considered interval. ⊓⊔ 2

Optimal Higher Order Time Discretizations 81<br />

Proof. We proceed by contradiction, assuming that M ≤ k − 1. For j =<br />

0,...,M − 1, one can find a point<br />

( )<br />

τ j+ 1 ∈ ]τ j+1,τ<br />

2 j [ such that [ψ m (R)] ]τ j+1 ,τ j+ 1 ] ⊂ ]0, 4[. (26)<br />

2<br />

Consider the polynomial ϕ defined at x ∈ R by<br />

M−1<br />

∏ ( )<br />

ϕ(x) =s 0 x − τj+ 1 .<br />

j=0<br />

Hence ϕ ≡ s 0 if M = 0. This polynomial is of degree M ≤ k − 1, so that<br />

it is a possible increment to R in P k−1 .Fort>0, consider the polynomial<br />

p t = ψ m (R + tϕ), which verifies for all x ∈ R:<br />

p t (x) =[ψ m (R)](x)+tx m+1 ϕ(x).<br />

We shall get a contradiction <strong>and</strong> conclude the proof if we show that, for any<br />

small t>0, p t (x) ∈ ]0, 4[ for x ∈ ]0,α m,k ] (since then α m (R + tϕ) >α m,k <strong>and</strong><br />

R would not be a local maximum).<br />

We shall only consider the case when [ψ m (R)] (α m,k ) = 4, since the reasoning<br />

is similar when [ψ m (R)] (α m,k )=0.Thens 0 = −1 by (25).<br />

• On the interval ]τ 1/2 ,α m,k ], ψ m (R) is greater than a positive constant<br />

(since it is positive on ]τ 1 ,α m,k ] by the definition of τ 1 <strong>and</strong> τ 1 0 small enough, p t ∈ ]0, 4[ on this interval.<br />

• Since ϕ(τ 1/2 ) = 0, p t (τ 1/2 ) = [ψ m (R)] (τ 1/2 ), which is in ]0, 4[ by the<br />

definition of τ 1/2 in (26).<br />

• On the interval ]τ 3/2 ,τ 1/2 [, ψ m (R) is less than a constant < 4(sinceit<br />

is < 4on]τ 2 ,τ 1/2 ] by the definition of τ 2 <strong>and</strong> τ 1/2 , see 2 <strong>and</strong> (26), <strong>and</strong><br />

τ 2 0 small enough, p t (·) ∈ ]0, 4[ on the interval.<br />

• If s M < 0 then, on the considered interval, the tangent points are all<br />

at y = 4, ψ m (R) is positive, <strong>and</strong> ϕ is negative. Since the map x ↦→<br />

[ψ m (R)] (x)/x =1+c 1 x + ... is greater than a positive constant on the<br />

considered interval, the map x ↦→ [ψ m (R)] (x)/x + tx m ϕ(x) =p t (x)/x is<br />

also positive on the interval for t>0 sufficiently small. It results that, for<br />

t>0 small enough, p t (·) ∈ ]0, 4[ on the considered interval. ⊓⊔<br />

2

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

Saved successfully!

Ooh no, something went wrong!