This paper implements topology optimization on two-dimensional manifolds. In this paper, the material interpolation is implemented on a material parameter in the partial differential equation used to describe a physical field, when this physical field is defined on a two-dimensional manifold; the material density is used to formulate a mixed boundary condition of the physical field and implement the penalization between two different types of boundary conditions, when this physical field is defined on a threedimensional domain with its boundary conditions defined on the two-dimensional manifold corresponding a surface or an interface of this three-dimensional domain. Based on the homeomorphic property of two-dimensional manifolds, typical two-dimensional manifolds, e.g., sphere, torus, Möbius strip and Klein bottle, are included in the numerical tests, which are provided for the problems on fluidic mechanics, heat transfer and electromagnetics.