PLUTO Test Problems
4.4-patch2
|
Double Mach reflection test problem. More...
Functions | |
void | InitDomain (Data *d, Grid *grid) |
void | UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid) |
Sets the initial condition for a planar shock front making an angle of with a reflecting wall:
where is the shock position. The ideal equation of state with
is used. The wedge is represented by a reflecting boundary starting at x=1/6 along the lower y-boundary. As the shock reflects off the lower wall, a complex flow structure develops with two curved reflected shocks propagating at directions almost orthogonal to each other and a tangential discontinuity separating them. At the wall, a pressure gradient sets up a denser fluid jet propagating along the wall. Kelvin-Helmholtz instability patterns may be identified with the rolls developing at the slip line. This feature is very sensitive to numerical diffusion and it becomes visible at high resolution and/or low dissipative schemes.
In the frames below we show configuration # 02 using a high-order finite difference scheme with 5-th order WENO-Z reconstruction and RK3
time stepping. The resolution is 960 x 240.
References
void InitDomain | ( | Data * | d, |
Grid * | grid | ||
) |
void UserDefBoundary | ( | const Data * | d, |
RBox * | box, | ||
int | side, | ||
Grid * | grid | ||
) |
Assign user-defined boundary conditions:
x < 1/6
and reflective boundary otherwise.