This preview shows page 1. Sign up to view the full content.
Unformatted text preview: MATH 572
Assignment #7 Numerical Methods for Scientiﬁc Computing II Winter 2005 due : Tuesday, April 19 LF = LaxFriedrichs , LW = LaxWendroﬀ
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/centereddiﬀ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 2norm?
Note
The ﬁnal exam is on Wednesday, April 27, 46pm 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
This note was uploaded on 12/08/2011 for the course MATH 623 taught by Professor Conlon during the Fall '08 term at University of Michigan.
 Fall '08
 CONLON
 Math, matlab, Scalar

Click to edit the document details