This paper proposed meshless Radial Point Interpolation Meshless Methods (RPIM) method for numerical
solution of natural convection in Darcy porous square cavity. It is assumed that Boussinesq approximation is valid to characteristic the buoyancy effect as the driving force of the fluid flow. The Galerkin global weak form is used to
discretize the system equations. The multiquadratic radial basis function (RBF) is chosen as the shape and test function.
Comparing the numerical results obtained using the proposed method with those obtained using the conventional
methods shows very good agreement.