Based on the stress-strain relations (Matsumoto et al., 2004; Yuu and Umekage, 2004) obtained by our Distinct Element Method calculation, particulate matter can be considered as a continuum. In this study, we simulated the flow fields of the two-dimensional bubbling fluidized bed using the Smoothed Particle Hydrodynamics (SPH) method, which gives the imaginary particle trajectory of the continuum, to calculate the equation of motion for the particulate continuum and the finite difference method for Navier-Stokes equations of the gas phase. The calculated results agree well with the experimental data. This indicates that the continuum model proposed by us and our simulation method presented in this paper including the application of the SPH method make it possible to predict correctly the flow fields in a high particle concentration gas system represented by the bubbling fluidized bed in which particulate matter exists in both quiescent and moving states.