In the present work, a Boundary Domain Integral Method, which has been already established for the solution of viscous incompressible fluid flow through porous media, is extended to capture compressible fluid flow in porous media. The presented numerical scheme was used for solving the problem of double diffusive natural convection in a square porous cavity heated from a side, while the horizontal walls are maintained at different concentrations. The Brinkman extension of Darcy equation is used to model the flow through porous medium. The velocity‐vorticity formulation is employed enabeling the computation scheme to be partitioned into kinematic and kinetic parts. The results of double diffusive natural convection in porous cavity are presented in terms of velocity, temperature and concentration redistributions.