A surface capturing method is developed for steady water-air flow with gravity. Second-order accuracy is obtained with flux limiting and turbulence is modeled with Menter's model. The model is solved efficiently with a combination of multigrid and defect correction. Results for two test cases confirm the efficiency and accuracy of the method.