Optimization of variational boussinesq models door