<p>This paper considers the meshfree collocation method for the flow model governed by Stokes and Darcy equations with the Beavers–Joseph–Saffman interface condition. This model is discretized by using hierarchical radial basis functions (H-RBFs) collocation method. In terms of trial, the hierarchical radial basis functions trial spaces are constructed by employing successive refinement scattered data sets and scaled compactly supported radial basis functions with varying support radii. It allows the easy construction of approximation spaces in arbitrary dimensions with arbitrary smoothness and avoids the huge workload caused by mesh-based methods at the same time. In terms of testing, it just needs to take values directly on the collocation points and greatly simplifies the difficulties caused by variation and integration. H-RBFs collocation method can not only solve the present problem on scattered points in complex domains with higher accuracy and lower computational cost, but also produce a highly sparse discrete algebraic system. Finally, these observations are obtained by taking the direct approach of numerical simulation.</p>