In RIPPLE, surface tension at free surfaces is replaced by a volume force derived from the continuum surface force model. A two-step projection method is used to solve for incompressible flow, with the pressure Poisson equation solved via a robust incomplete Cholesky conjugate gradient technique. Momentum advection is estimated with a technique very similar to the weakly monotonic upwind method of van Leer. Flow obstacles and curved boundaries interior to the mesh are representable with a partial cell treatment that borrows from two-phase flow physics. The numerical scheme of RIPPLE is based on a finite difference solution of a coupled set of partial differential equations governing incompressible fluid flow. Finite difference solutions to the incompressible Navier-Stokes equations are obtained on an Eularian, rectilinear mesh in Cartesian or cylindrical geometries. Free surfaces are represented with volume-of-fluid data on the mesh. Surface tension is modeled as a volume force derived from the continuum surface force model. The continuum surface force model for surface tension gives RIPPLE both robustness and accuracy in modeling surface tension effects at the free surface. RIPPLE can also model wall adhesion effects.