In this study we describe an efficient numerical method, based on the semi-Lagrangian (SL) semi-implicit (SI) method and Newton's method, for obtaining steady-state (SS) solutions of equations arising in distributed-loop models of the urine concentrating mechanism. Dynamic formulations of these models contain large systems of coupled hyperbolic partial differential equations (PDEs). The SL method advances the solutions of these PDEs in time by integrating backward along flow trajectories, thus allowing large time steps while maintaining stability. The SI approach controls stiffness arising from transtubular transport terms by averaging these terms in time along flow trajectories. An approximate SS solution of a dynamic formulation obtained via the SLSI method can be used as an initial guess for a Newton-type solver, which rapidly converges to a highly accurate numerical approximation to the solution of the ordinary differential equations that arise in the corresponding SS model formulation. In general, it is difficult to specify a priori for a Newton-type solver an initial guess that falls within the radius of convergence; however, the initial guess generated by solving the dynamic formulation via the SLSI method can be made sufficiently close to the SS solution to avoid numerical instability. The combination of the SLSI method and the Newton-type solver generates stable and accurate solutions with substantially reduced computation times, when compared to previously applied dynamic methods.