math572_hw7

# math572_hw7 - MATH 572 Assignment#7 Numerical Methods for...

This preview shows page 1. Sign up to view the full content.

This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: MATH 572 Assignment #7 Numerical Methods for Scientiﬁc Computing II Winter 2005 due : Tuesday, April 19 LF = Lax-Friedrichs , LW = Lax-Wendroﬀ 1. Consider the scalar 1D wave equation vt + cvx = 0 with c > 0. a) Plot (using Matlab) the ampliﬁcation factor ρ(ξh) in the complex plane for the upwind, LF, and LW schemes, for λ = 1/4, 1/2, 3/4, 5/4. Make a plot for each value of λ, showing all three methods on the same plot. Include a dotted line showing the unit circle. How are the methods similar? diﬀerent? b) Find the model equation for the LF and LW schemes. Compare these to the model equation for the upwind scheme (derived in class). Which scheme has the most artiﬁcial viscosity - upwind, LF, or LW? Which scheme has the least? c) The numerical wave speed has an expansion of the form c = c(1 + γ1 ξh + γ2 (ξh)2 ) + · · ·) in the long wave limit ξh → 0. Find γ1 , γ2 for the upwind and LF schemes (LW was done in class). Compare the size of the phase error in the three schemes. 2. Consider the LF scheme for vt + cvx = 0, where c may be positive or negative. a) Show that the scheme is stable in the ∞-norm if |c|λ ≤ 1. b) Show that the scheme converges in the ∞-norm if |c|λ ≤ 1. 3. Problem 3 on hw6 asked you to compute the solution of vt + vx = 0 using the upwind scheme for two initial data. Repeat now using LF and LW. Are the present results better or worse than those obtained using the upwind scheme? 4. Consider the equation φt + φxxx = 0 and the forward Euler/centered-diﬀerence scheme, un+1 = un + kD0 D+ D− un . What condition on the parameters h, k is necessary to ensure j j j stability of the scheme in the 2-norm? Note The ﬁnal exam is on Wednesday, April 27, 4-6pm in room 130 Dennison. The exam will cover the entire course, but PDEs will be emphasized. You may use two pages of notes (i.e. 2 × 8.5 in × 11 in). Calculators are not allowed. ...
View Full Document

{[ snackBarMessage ]}

Ask a homework question - tutors are online