A numerical method is developed to calculate the shape of a liquid aspherical lens modulated by electrostatic force by solving a multiphysics problem that consists of electric field module and fluid flow module. The electric field module is used to calculate the electrical force on the lens surface and the fluid flow module is used to describe the velocity field inside the liquid lens by solving incompressible Navier-Stokes equation with the electrical force as a boundary condition. The arbitrary Lagrangian Eulerian (ALE) method is adopted to get the lens shape evolution since it is volume conservative. Lens shapes at different voltages and different kinds of electrodes are discussed by simulation results. Optical properties of the lens are then studied.