Numerical modeling of Boussinesq equations by finite element method