We present a three-dimensional (3D) lattice model based on the differential adhesion hypothesis (DAH) to study self-assembly and fusion of multicellular aggregates, which forms the foundation for the scaffold-less biofabrication of tissues and organs, known as "bio-printing". In this new technology, live multicellular aggregates are used as bio-ink to make tissue or organ constructs via the layer-by-layer deposition technique in biocompatible hydrogels; the printed bio-constructs embedded in the hydrogels are then placed in bioreactors to undergo the fusion process to form the desired functional tissue or organ products. Our in-silico approach is an agent-based method, which uses the kinetic Monte Carlo (KMC) algorithm to evolve the cellular system on a 3D lattice. In this model, the cells and the hydrogel media, in which cells are embedded, are coarse-grained to material's points on the lattice, where the cellcell and cell-medium interaction are quantified by adhesion and cohesion energies. In a multicellular aggregate system with a fixed number of cells and fixed amount of hydrogel media, the effect of cell differentiation, proliferation, and death are tactically neglected in the current model, which can be readily included in the model, and the interaction s primarily dictated by the interfacial energy between cell and cell as well as between cells and medium particles on the lattice, respectively. The KMC method is applicable to transient simulations of fusion of multicellular aggregates at the time and length scale appropriate to biofabrication. Numerical experiments are presented to demonstrate fusion and cell sorting during tissue and organ maturation processes.