A numerical scheme based on Boundary Domain Integral Method (BDIM) is being
developed to model the bubbling behaviour of fluidized beds by using the E
ulerian approach. Each phase is described by a modified Navier-Stokes equat
ion. While modified Helmholtz fundamental solution has been used for the ki
nematics, the elliptic diffusion-convective one has been used for the kinet
ics of the both components motion. A gas-phase is described by its physical
properties, but the solid-phase needs additional physical models. The most
promising set of the equations in the sense of fast numerical convergence
is presented in present work.