Project 7
Wave Equation in one dimension on accelerating grid
(beta = ax)
Implicit Method (9 point)
with Radiation Boundary Conditions
General Results
- The shift vector is beta = ax
- The code works when a < 0
- The code does not work well when a > 0.5 (approximately). The
code consistantly begins to blow up at the boundaries.
- The code blows up faster when as a gets bigger. When a = 0.6,
the code blows up after approximately 170 iterations. When a = 1.1,
the code lasts for about 20 iterations. When a = 1.6
the code blows up after 15 iterations.
Mpeg movies to be added in the future!