The purpose of this paper is to establish bilinear estimates in Besov spaces generated by the Dirichlet Laplacian on a domain of Euclidian spaces. These estimates are proved by using the gradient estimates for heat semigroup together with the Bony paraproduct formula and the boundedness of spectral multipliers.