The evolution of the Universeˇlled with dark energy (DE) with or without perfect uid is discussed. In doing so, we consider a number of cosmological models, namely Bianchi type I, III, V, VI 0 , VI, and FRW one. For the anisotropic cosmological models we have used proportionality condition as an additional constraint. The exact solutions to theˇeld equations in quadrature are found in the case of a BVI model. It was found that the proportionality condition used here imposed severe restriction on the energyÄmomentum tensor, namely, it leads to isotropic distribution of matter.Anisotropic BVI 0 , BV, BIII, and BI DE models with variable EoS parameter ω have been investigated by using a law of variation for the Hubble parameter. In this case the matter distribution remains anisotropic, though, depending on the concrete model, there appear different restrictions on the components of energyÄmomentum tensor. That is why we need an extra assumption such as a variational law for the Hubble parameter. It is observed that, at the early stage, the EoS parameter ω is positive, i.e., the Universe was matter-dominated at the early stage but at later time the Universe is evolving with negative values, i.e., the present epoch. DE model presents the dynamics of EoS parameter ω whose range is in good agreement with the acceptable range by the recent observations.A spatially homogeneous and anisotropic locally rotationally symmetric Bianchi-I spaceÄtimě lled with perfect uid and anisotropic DE possessing dynamical energy density is studied. In the derived model, the EoS parameter of DE (ω (de) ) is obtained as time varying and it is evolving with negative sign which may be attributed to the current accelerated expansion of Universe. The distance modulus curve of derived model is in good agreement with SNLS type Ia supernovae for high red-shift value which in turn implies that the derived model is physically realistic.A system of two uids within the scope of a spatially at and isotropic FRW model is studied. The role of the two uids, either minimally or directly coupled in the evolution of the dark energy parameter, has been investigated. In doing so, we have used three different ansatzs regarding the scale factor that gives rise to a variable decelerating parameter. It is observed that, in the noninteracting case, both the open and the at Universes can cross the phantom region whereas in the interacting case only the open Universe can cross the phantom region. The stability and acceptability of the obtained solution are also investigated.