In the framework of the Polyakov quark-meson model with two flavors, the bubble dynamics of a first-order phase transition in the region of high density and low temperature are investigated by using the homogeneous thermal nucleation theory. In mean-field approximation, after obtaining the effective potential with the inclusion of the fermionic vacuum term, we build a geometric method to search two existing minima, which can be actually connected by a bounce interpolated between a local minimum to an adjacent global one. For both weak and strong first-order hadron quark phase transitions, as fixing the chemical potentials at μ=306MeV and μ=310MeV, the bubble profiles, the surface tension, the typical radius of the bounce, and the saddle-point action as a function of temperature are numerically calculated in the presence of a nucleation bubble. It is found that the surface tension remains at a very small value even when the density is high. It is also noticed that the deconfinement phase transition does not change the chiral phase transition dramatically for light quarks and phase boundaries for hadron and quark matter should be resized properly according to the saddle-point action evaluated on the bounce solution.