With the advent of 5G and the need for energy-efficient massively connected wireless networks, in this work, we consider an energy harvesting (EH) based multi-relay downlink cooperative non-orthogonal multiple access (NOMA) system with practical constraints. The base station serves NOMA users with the help of decode-and-forward based multiple EH relays, where relays harvest the energy from the base station's radio frequency. A relay is selected from the multiple K-relay by using a partial relay selection protocol. The system is considered to operate in half-duplex mode over a generalized independent and identical Nakagami$-m$ fading channel. The closed-form expression of outage probability and ergodic rate are derived for users, under the assumption of imperfect channel state information (CSI) and imperfect successive interference cancellation (SIC) at the receiver node. Expression of outage probability and ergodic rate for both users under the assumption of perfect CSI and perfect SIC are also presented. Further, the asymptotic expression for the outage probability is also shown. The derived analytical expressions are verified through Monte-Carlo simulations.