Abstract
This paper describes numerical procedures of a new finite element method for thermal hydraulic analysis using equal-order velocity-pressure interpolation. The main feature of the proposed method lies in applying particular numerical schemes to the pressure gradient terms in the momentum conservation equations and the mass balance terms in the pressure equation for the sake of strong coupling of velocity and pressure fields. Some numerical results show the basic potential of the proposed method.