BEARCLAW Application description

../bearclaw/examples/classicbear/2d/Advection/example1

The uniform velocity advection problem is solved with a number of algorithms and input parameters. The refinement ratio between succesive grid levels is 2 in all the computations presented below.

1. Lax-Friedrichs

The Lax-Friedrichs scheme introduces a large amount of numerical diffusion. When using a single level of grids the initial square pulse is rapidly diffused.

Single grid, CFL=0.5

When using smaller step sizes the numerical diffusion is even larger

Single grid, CFL=0.25

The effect is reduced when using multiple grid levels, but still noticeable

Three grid levels, CFL=0.5

Three grid levels, CFL=0.25

 

2. Upwind

Upwind also introduces significant numerical diffusion, but less than Lax-Friedrichs. One can notice that the error estimation procedure generates grids that conform more closely to the initial square discontinuity than in the Lax-Friedrichs computation, a reflection of the smaller values of numerical diffusion.

Single grid, CFL=0.5

Three grid levels, CFL=0.5

Four grid levels, CFL=0.5

 

3. Lax-Wendroff

The second order Lax-Wendroff scheme has a dispersive leading order error term. Numerical diffusion is less than in the previous schemes but the dispersive error leads to over/undershoots in the solution. The dispersive error subsides quickly with increasing resolution, a reflection of the scheme's second order of accuracy.

Single grid, CFL=0.5

Three grid levels, CFL=0.5

Four grid levels, CFL=0.5

 

Directory listing

Download directory in compressed tar file format (Note: Please unpack the tar file in your CLAW home directory.)


CLAWPACK Home page