Search for the development of a reliable mathematical model for understanding bubble dynamics behavior is an ongoing endeavor. A long list of complex phenomena underlies physics of this problem. In the past decades, the lattice Boltzmann (LB) method has emerged as a promising tool to address such complexities. In this regard, we have applied a 121-velocity multiphase lattice Boltzmann model (LBM) to an asymmetric cluster of bubbles in an acoustic field. A problem as a benchmark is studied to check the consistency and applicability of the model. The problem of interest is to study the deformation and coalescence phenomena in bubble cluster dynamics, and the screening effect on an acoustic multi-bubble medium. It has been observed that the LB model is able to simulate the combination of the three aforementioned phenomena for a bubble cluster as a whole and for every individual bubble in the cluster.