The planar 2D k-ε double equations turbulence model was adopted and transformed into non-orthogonalcurvilinear coordinates. The concentration convection-diffusion was introduced to planar 2D SIMPLEC algorithm offlow in non-orthogonal curvilinear coordinates. The numerical model of pollutant transportation in non-orthogonalcurvilinear coordinates was constructed. The model was applied to simulate the flow and pollutant concentrationfields. In the testing concentration field, two optimal operations of contamination discharging both along bank and inthe centerline at the first bend of the meandering channel were adopted. Comparison with available data showed themodel developed was successful, was valuable to engineering application.