We numerically investigate the thermally unstable accretion disks around black holes. We adopt an evolutionary viscous stress equation to replace the standard alpha-prescription based on the results of two MHD simulations. We find a kind of interesting oscillations on some running models in limit-cycle outburst state. The oscillations arise near the inner boundary and propagate radially outwards. We deem that they are the trapped $p$-mode oscillations excited by sonic-point instability. We direc...