Generally speaking, the finite-element method in computational fluid dynamics is a universal numerical method, but in the computation of a flow phenomenon with unstructured grids, it takes much CPU time to construct computational grids. Therefore it is important to develop an analytical formulation for the purpose of reducing the time required to construct computational grids. In the present paper, in order to overcome the defect of the finite-element method, we propose a new mixed- element method using the discrete del operator defined as an element average of the gradient of the shape function in discrete space. The analytical expression of the discrete del operator for the mixed-element method is a vector in two or three dimensions, and becomes a nonmemorizing method. Furthermore when we use the discrete del operator, the natural description of programing becomes objective and compact.