  A new method is presented for mesoscopic simulations of particle dispersions
in nematic liquid crystal solvents. It allows efficient first-principle
simulations of the dispersions involving many particles with many-body
interactions mediated by the solvents. A simple demonstration is shown for the
aggregation process of a two dimentional dispersion.
