Modeling groundwater flow by lattice Boltzmann method in curvilinear coordinates