Next: 5.1.6 Resumé: Conservative-hyperbolic DE
Up: 5.1 Initial Value Problems
Previous: 5.1.4 Lax-Wendroff Scheme (LW)
5.1.5 Lax and Lax-Wendroff in Two Dimensions
(advective case: and
)
Lax scheme:
Figure 5.2:
Lax method in two dimensions
|
Lax-Wendroff:
For the second stage (half-step leapfrog) we need
etc., which requires
, which must be determined from
,
etc.
But: quantities with half-step spatial indices
(
,
etc.)
are given at half-step times
() only.
Modifying the LW scheme to allow for this, we have
Lax-Wendroff in 2 dimensions:
- Lax method to determine the -values at half-step time
:
etc.
- Evaluation at half-step time:
- leapfrog with half-step:
Lax-Wendroff in two dimensions
Figure:
First stage (= Lax) in the 2-dimensional LW method:
...
, ...
|
For only the points
(at ) are
used;
for
we use the points .
Problem: Drift between subgrids and .
Solution: If the given PDE contains a diffusive
term, this guarantees coupling. Otherwise, artificially add a small diffusive
term.
Stability analysis:
Fourier modes are now 2-dimensional:
Assuming
we find the CFL condition
Next: 5.1.6 Resumé: Conservative-hyperbolic DE
Up: 5.1 Initial Value Problems
Previous: 5.1.4 Lax-Wendroff Scheme (LW)
Franz J. Vesely Oct 2005
See also: "Computational Physics - An Introduction," Kluwer-Plenum 2001