For three‐dimensional fluid flows in complex geometries, it is convenient to make predictions using a non‐orthogonal boundary‐fitted mesh. The present paper describes an economical method of solving the equations of motion for two and three dimensional problems using such meshes. The locations on the mesh at which the depenent variables are calculated, and the methods used to solve the equations, are key issues in the development of a successful algorithm; these are discussed in the present paper. Results obtained when the proposed method is applied to several problems are also described. The method is intended for flows in which compressibility effects do not dominate.