Sonoporation is a non-invasive method that uses ultrasound for drug and gene delivery for therapeutic purposes. Here, both Finite Element Method (FEM) and Lattice Boltzmann Method (LBM) are applied to study the interaction physics of microbubble oscillation and collapse near flexible tissue. After validating the Finite Element Method with the nonlinear excited lipid-coated microbubble as well as the Lattice Boltzmann Method with experimental results, we have studied the behavior of a three-dimensional compressible microbubble in the vicinity of tissue. In the FEM phase, the oscillation microbubble with a lipid shell interacts with the boundary. The range of pressure and ultrasound frequency have been considered in the field of therapeutic applications of sonoporation. The viscoelastic and interfacial tension as the coating properties of the microbubble shell have been investigated. The presence of an elastic boundary increases the resonance frequency of the microbubble compared to that of a free microbubble. The increase in pressure leads to an expansion in the range of the microbubble’s motion, the velocity induced in the fluid, and the shear stress on the boundary walls of tissue. An enhancement in the surface tension of the microbubble can influence fluid flow and reduce the shear stress on the boundary. The multi-pseudo-potential interaction LBM is used to reduce thermodynamic inconsistency and high-density ratio in a two-phase system for modeling the cavitation process. The three-dimensional shape of the microbubble during the collapse stages and the counter of pressure are displayed. There is a time difference between the occurrence of maximum velocity and pressure. All results in detail are presented in the article bodies.