We show that an explicit method for solving hyperbolic partial differe
ntial 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 sp
ace and time, we employ the recently developed ENO (Essentially Non-Os
cillatory) methodology. We present examples of computed flows and conc
entration profiles in representative model contexts. Finally, we indic
ate briefly how model tubules may be coupled to construct large-scale
simulations of the renal counterflow system.