Abstract
A Eulerian-Lagurangian numerical method is very useful for groundwater flow with mass transport analysis. In this study, we improved the method for two-dimensional saturated-unsaturated flow in consideration of density dependent characteristic. Validity of this method was proved in comparison with the theorical solutions and the experimental results. Also, this paper show that the method enable to simulate precisely groundwater flow with mass transport without dependence on mesh size.