2D shallow water equations with depth-averaged k−εmodel is considered. A meshless method based on multiquadric radial basis functions isdescribed. This methods is based on the collocation formulation and does not require thegeneration of a grid and any integral evaluation. The application of this method to a flowin horizontal channel, taken as an experimental device, is presented. The results ofcomputations are compared with experimental data and are found to be satisfactory