Mathematical modeling of cold cap: Effect of bubbling on melting rate