  A lattice Boltzmann model with interacting particles was developed in order
to simulate the magneto-rheological characteristics of magnetic fluids. In the
frame of this model, $6\, +\,1$ species of particles are allowed to move across
a $2D$ triangular lattice. Among these species, $6$ of them carry an individual
magnetic dipole moment and interact themselves not only as a result oflocal
collisions, as in current Lattice Boltzmann models, but also as a result of
nearest neighbours magnetic dipole-dipole interaction. The relative
distribution of the individual magnetic moments is determined by the intensity
of an external static magnetic field acting on the whole system. This model
exhibits some relevant characteristics of real magnetic fluids, i.e.,
anisotropic structure formation as a result of magnetic field induced
gas-liquid phase transition and magnetic field dependence of the sound velocity
and the attenuation coefficient.
