To simplify the problem, we assume the shock wave will not reach the
boundary of our simulation region(this can be controlled by looking at the
profile of gas parameters at certain times). Thus we can simply select fixed
boundary conditions.
*f*(*x*_{b},*t*)=*f*(*x*_{b},0).