|
|
|
|
|
| Remark:Currently, a total variation diminishing (TVD) scheme is applied for the famous double mach reflection test problem, a mach-10 shock wave impinging on a 30 degree wedge (refer to the 1984 Woodward and Colella's paper or Quirk's thesis). The simulation above uses the MinMod slope limiter and an exact Riemann Solver developed by E. F. Toro. As we see, the AMR yields high resolution data while only five levels of adaptively refined grids are used in the simulation. In order to generate iso-contours for density, we first reconstruct continuous piece-wise linear node-centered data from piece-wise constant cell-centered output for each level grid, then a method called marching cubes is employed to extract isocontours within each cell. The reconstructed node-centered data is discontinuous along coarse-fine level interfaces due to the grid inconsistency there. Oscillations are observed around the shock wave close to the top boundary. This is due to inaccurate treatment of boundary conditions there. Currently, transmissive/transparent BCs are imposed. To overcome the problem, we can apply Dirichlet boundary conditions by an exact solution along the top boundary. Click on each of the images above, you may see dynamic animations. If you like further to know how the adaptively refined grids are partitioned indistributed computing, please clickhere. |