|
|
|
|
|
| two sin waves | triangle pulse | one sine wave | rectangular pulse |
These are HW4 animations.
The following are animated GIFs showing the finite difference numerical solution to problem 1 as described in the above HW. The scheme used is Lax-Wendroff.
Clicking on an image will start the animation in a new window.
These simulations only show the pressure wave, \(p(x,t)\) and not the acoustic perturbation velocity \(u(x,t)\).
This solution was run with boundary conditions which caused the sound wave to reflect from both boundaries. This is what would happen inside a room with reflective walls such as concerete or wood.
|
|
|
|
|
| two sin waves | triangle pulse | one sine wave | rectangular pulse |
This solution was run with boundary conditions which caused the sound wave to reflect from only the left boundary but absorbed into the right boundary. This is what would happen inside a room with one wall treated with material to absorbe the sound waves reaching it.
|
| |||
|
|
|
|
|
| two sin waves | triangle pulse | one sine wave | rectangular pulse
|
The following are animations of the numerical solution to \(u_t+a u_x=0\). The solution used the finite volume method using 7 different numerical flux limiter functions to compare performance.
These 7 methods are defined in the problem statement in the report above.
The methods are
The following tables show the results of the simulations. 4 tables are given. Each table is for a different initial conditions. In all of these results, the maximum run time was \(5\) seconds. In order to reduce the size of the animation file, not every frame was captured from the simulation run.
Courant number used was \(0.9\), the advection speed was set at \(a=1\) and grid spacing was \(h=0.005\) meters. The domain is \([0,1]\) using cell centered grid.
These animations will run only once and stop at 5 seconds. To run it again, simply reload the web page using the browser reload button, this will cause the animation to start from the beginning again.
The wave packet is defined as \(u(x,0)=cos(16 \pi x)exp(-50(x-0.5)^2)\)
|
| |||
|
|
|
|
|
| Upwinding: Large diffusion at wave crest and trough. but no shift. Grade \(F\) |
Lax-Wendroff: Some diffusion at wave crest and trough, in addition of significant shift to the left relative to exact solution, grade: \(B-\) |
Beam-Warming: Similar to Lax-Wendroff, but shift is to the right. Grade \(B\) |
Minmod: diffusion at wave crest and trough, as with upwinding, but not as bad. No shifting. Grade \(C\)
|
|
|
|
|
|
| Superbee: NO shifting, and very small amount of diffusion at the crest and trough. Grade \(B+\) |
MC Limited: Similar to superbee, but a little more diffusion at the crest and trough. Grade \(B\) |
Van Leer: Similar to MC limited, but even more diffusion at the crest and trough. Grade \(B-\) |
|
The wave packet is defined as \(u(x,0)=sin(2 \pi x)sin(4 \pi x)\)
|
| |||
|
|
|
|
|
| Upwinding: No shifting, but large amount of diffusion at the crest and trough of the wave. Grade \(C\) |
Lax-Wendroff: Did well, no shifting nor diffusion seen. Grade \(A\) |
Beam-Warming: Similar to Lax-Wendroff. Grade \(A\) |
Minmod: No shifting, but small amount of diffusion near the crest and trough. Grade \(B\)
|
|
|
|
|
|
| Superbee: No shift and no diffusion. But at the crest and trough, the numerical solution appeared to less smooth than with Lax-Wendroff making a small straight edge. Grade \(A-\) | MC Limited: Similar to Lax-Wendroff, better than Superbee around the crest and trough, no straight edge appeared. Grade \(A\) | Van Leer: No diffusion and no shifting. Grade \(A\) |
|
A step function from \(x=0.25\) to \(x=0.75\).
|
| |||
|
|
|
|
|
| Upwinding: No ripples seen, follows the general form of the step function, but large amount of diffusion near the corners. Grade \(C\) |
Lax-Wendroff: Large ripples around the corners on the left of the step function. Less diffusion than upwinding. Grade \(C\) |
Beam-Warming: The ripples are larger and have a larger extent than Lax-Wendroff. Grade \(C-\) |
Minmod: No ripples and little diffusion. An improved version of upwinding. Grade \(C+\)
|
|
|
|
|
|
| Superbee: The best scheme for the step function. No ripples, very closely followed the exact solution, but very small diffusion is still there. Grade \(A-\) |
MC Limited: Similar to supperbee, but more diffusion. Grade \(B\) |
Van Leer: Similar to MC limited. Grade \(B+\) |
|
The initial condition used for this test is \(u(x,0)= (X > 0.1)(X < 0.3) +exp(-200(X-0.75)^2)\)
This test just combines the step function with the low frequency smooth test done above. Hence, the same comments will apply as above.
|
| |||
|
|
|
|
|
| Upwinding |
Lax-Wendroff |
Beam-Warming |
Minmod
|
|
|
|
|
|
| Superbee |
MC Limited |
Van Leer |
|
