  We obtain two-sided heat kernel estimates for Riemannian manifolds with ends
with mixed boundary condition, provided that the heat kernels for the ends are
well understood. These results extend previous results of Grigor'yan and
Saloff-Coste by allowing for Dirichlet boundary condition. The proof requires
the construction of a global harmonic function which is then used in the
$h$-transform technique.
