Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
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
- Page 37 and 38: Discontinuous Galerkin Methods 25 [
- Page 39 and 40: Mixed Finite Element Methods on Pol
- Page 41 and 42: Mixed FE Methods on Polyhedral Mesh
- Page 43 and 44: Mixed FE Methods on Polyhedral Mesh
- Page 45 and 46: Mixed FE Methods on Polyhedral Mesh
- Page 47 and 48: 4 Hybridization and Condensation Mi
- Page 49 and 50: Mixed FE Methods on Polyhedral Mesh
- Page 51 and 52: is symmetric and positive definite,
- Page 53 and 54: with some coefficient α ∈ R wher
- Page 55 and 56: 44 E.J. Dean and R. Glowinski so fa
- Page 57 and 58: 46 E.J. Dean and R. Glowinski 2 A L
- Page 59 and 60: 48 E.J. Dean and R. Glowinski S:T=
- Page 61 and 62: 50 E.J. Dean and R. Glowinski minim
- Page 63 and 64: 52 E.J. Dean and R. Glowinski 6 On
- Page 65 and 66: 54 E.J. Dean and R. Glowinski Fig.
- Page 67 and 68: 56 E.J. Dean and R. Glowinski and
- Page 69 and 70: 58 E.J. Dean and R. Glowinski 7 Num
- Page 71 and 72: 60 E.J. Dean and R. Glowinski Fig.
- Page 73 and 74: 62 E.J. Dean and R. Glowinski Assum
- Page 75 and 76: Higher Order Time Stepping for Seco
- Page 77 and 78: u n+1 h Optimal Higher Order Time D
- Page 79 and 80: Optimal Higher Order Time Discretiz
- Page 81 and 82: Optimal Higher Order Time Discretiz
- Page 83 and 84: Optimal Higher Order Time Discretiz
- Page 85 and 86: Optimal Higher Order Time Discretiz
- Page 87: Optimal Higher Order Time Discretiz
- Page 91 and 92: Optimal Higher Order Time Discretiz
- Page 93 and 94: Optimal Higher Order Time Discretiz
- Page 95 and 96: Optimal Higher Order Time Discretiz
- Page 97 and 98: Optimal Higher Order Time Discretiz
- Page 99 and 100: Optimal Higher Order Time Discretiz
- Page 101 and 102: Optimal Higher Order Time Discretiz
- Page 103 and 104: 96 I. Sazonov et al. To provide a p
- Page 105 and 106: 98 I. Sazonov et al. In the first s
- Page 107 and 108: 100 I. Sazonov et al. Fig. 1. An ex
- Page 109 and 110: 102 I. Sazonov et al. H z 1 exact F
- Page 111 and 112: 104 I. Sazonov et al. (a) (b) Fig.
- Page 113 and 114: 106 I. Sazonov et al. Scattering Wi
- Page 115 and 116: 108 I. Sazonov et al. 6.4 Scatterin
- Page 117 and 118: 110 I. Sazonov et al. (a) (b) Fig.
- Page 119 and 120: 112 I. Sazonov et al. [MHP96] K. Mo
- Page 121 and 122: 114 R. Sanders and A.M. Tesdall I R
- Page 123 and 124: 116 R. Sanders and A.M. Tesdall imp
- Page 125 and 126: 118 R. Sanders and A.M. Tesdall alo
- Page 127 and 128: 120 R. Sanders and A.M. Tesdall (a)
- Page 129 and 130: 122 R. Sanders and A.M. Tesdall D C
- Page 131 and 132: 124 R. Sanders and A.M. Tesdall 8.6
- Page 133 and 134: 126 R. Sanders and A.M. Tesdall 0.3
- Page 135 and 136: 128 R. Sanders and A.M. Tesdall [TR
- Page 137 and 138: 132 S. Lapin et al. Ω R γ Ω 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