In previous work [21] , we studied the local well-posedness of weak solution to the 1-D full compressible Navier-Stokes equation with initial data of small total variation. Specifically, the local existence, the regularity, and the uniqueness in certain function space of the weak solution have been established. The basis for previous study is the precise construction of fundamental solution for heat equation with BV conductivity. In this paper, we continue to investigate the global stability and the time asymptotic behavior of the weak solution. The main step is to construct the "effective Green's function", which is the combination of the heat kernel with BV coefficient in short time and the Green's function around constant state in long time. The former one captures the quasi-linear nature of the system, while the latter one respects the dissipative structure. Then the weak solution is written into an integral system in terms of this "effective Green's function", and the time asymptotic behavior is established based on a priori estimate. Contents 1. Introduction 1 2. Preliminary 5 2.1. Weak solutions 5 2.2. Heat kernel estiamtes 6 3. Green's function 9 3.1. Analysis on eigenvalues 9 3.2. Singular parts of Green's function 13 3.3. Regular parts of Green's function 22 4. Global well-posedness 31 4.1. Representation by Green's function 32 4.2. Global existence, uniquenss and large time behavior 38 Appendix A. Coefficients of the approximated eigenvalues λ * j 61 Appendix B. The expression of M * ,k j in Lemma 3.4 62 Appendix C. The expression of M k j in (3.54) 63 Acknowledgments 64 References 64