Orientational ordering in two-dimensional polymer solutions: Monte Carlo simulations of a bond fluctuation model