The direct strong-form collocation method with reproducing kernel approximation is introduced to efficiently and effectively solve the natural convection problem within a parallelogrammic enclosure. As the convection behavior in the fluid-saturated porous media involves phase coupling, the resulting system is highly nonlinear in nature. To this end, the local approximation is adopted in conjunction with Newton–Raphson method. Nevertheless, to unveil the performance of the method in the nonlinear analysis, only single thermal natural convection is of major concern herein. A unit square is designated as the model problem to investigate the parameters in the system related to the convergence; several benchmark problems are used to verify the accuracy of the approximation, in which the stability of the method is demonstrated by considering various initial conditions, disturbance of discretization, inclination, aspect ratio, and reproducing kernel support size. It is shown that a larger support size can be flexible in approximating highly irregular discretized problems. The derivation of explicit operators with two-phase variables in solving the nonlinear system using the direct collocation is carried out in detail.