Shock-focusing concave reflectors can have parabolic, circular or elliptic cross-sections. They produce effectively a very high pressure at the focusing point. In the past, many optical images have been obtained on shock focusing via experiments. Measurement of field variables is, however, difficult in the experiment. Using the wave propagation algorithm and the Cartesian embedded boundary method, we have successfully obtained numerical Schlieren images that appear very much like the experimental results. In addition, we obtained the detailed field variables such as pressure, velocity, density and vorticity in the unsteady domain. The present numerical results have made it possible to understand the shock focusing phenomenon in more detail than before.