We show that an explicit method for solving hyperbolic partial differential equations can be applied to a model of a renal tubule to obtain both dynamic and steady-state solutions. Appropriate implementation of this method eliminates numerical instability arising from reversal of intratubular flow direction. To obtain second-order convergence in space and time, we employ the recently developed ENO (Essentially Non-Oscillatory) methodology. We present examples of computed flows and concentration profiles in representative model contexts. Finally, we indicate briefly how model tubules may be coupled to construct large-scale simulations of the renal counterflow system.