Compressible flow problems are frequently encountered in mechanical engineering. In fluid flow simulation at high Reynolds number, due to reduction of the viscous effects, compressibility can be neglected. By dropping the viscous term out of the Navier stokes equation, the Euler equation is obtained. To investigate supersonic and hypersonic flows the Euler equation can be used. Simulation of compressible flow is usually very time consuming hence its computational cost is high. Parallel processing with Graphical Processing Unit (GPU) has secondly been widely used to reduce such costs. In this work numerical solution of the governing equations for inviscid compressible flow using HLL, HLLC and WENO5 methods have been described and the accuracy of them for one dimensional problems (for which the analytical solution of Riemann is available) has been investigated to validate results. The explosion problem in two dimensions has been solved. The present solution has been compared with the Godunov method. The solver of the Riemann problems in the present work is almost the same as the exact solver, however an additional source term has been included, to implement the boundary condition for the 2D problem. The immersed boundary method has been used. In this method a ghost fluid is used for the extension of the boundary condition with high accuracy. Using this method one can use the rectangular Cartesian grid method instead of considering a boundary fitted grid. This will simplify the implementation costs. To validate the immersed boundary method, the triple points resulted from the trailing works will be located by an interpolation method and the result will be compared with experimental and numerical results reported by others. The computed codes for the simulation of the compressible fluid flow have been written using the CUDA programming language and may have been executed in parallel on the of graphical processing units. The GPU that has been use for this work was GerfoceGTX580, with 512 cores and 1536 MB accessory memory and 192 Gbyte/sec speed of data transfer. The speed up ratio is directly related to the number of cores, and it will be effectively increased as the number of computational cores is increased. The speed up ratio for the HLL code for the one dimensional problem was 30X and for the two dimensional case it was more than 80X. For the WENO5 code we achieved a 36X speed up ratio for the 1D case and we reached 174X speed up ratio for the 2D problem. The above mentioned speeds up ratio are very considerable and for compressible solvers have not been reported previously. In this research the effect of increasing the size of the block on the speed up ratio has been studied, and we showed that best size of the block for the used graphical card is 512x1. Using immersed boundary method, the speed up ratio was decreased however with increasing the domain size the speed up ratio will approach to the speed of ratio of the other method which does not use the immersed boundary technique. Keywords : Compressible flow, CUDA, HLL, HLLC, WENO, Immersed Boundary