Generally speaking, finite-element methods in computational fluid dynamics are universal, but uneconomical. In the present paper, in order to overcome this defect in FEM, we propose a new method using the discrete del operator which is a coordinate-free differential operator in discrete space. This operator in discrete space is defined as an element average of the gradient of the shape function, and it has three characteristics: orthogonality, identity and symmetry. Furthermore, the discrete del operator is useful in nonmemorizing and in easy coding. Since the analytical expression of the discrete del operator is a vector in two or three dimensions, the natural description of programing becomes objective and compact, which is more understandable for nonspecialists of CFD.