Abstract
The present paper describes a numerical solution to topology optimization problems of domains in which boundary value problems of partial differential equations are defined. Density raised to a power is used instead of the characteristic function of the domain. A design variable is set by a function on a fixed domain which is converted to the density by a sigmoidal function. Evaluation of derivatives of cost functions with respect to the design variable appear as stationary conditions of the Lagrangians. A numerical solution is constructed by a gradient method in a design space for the design variable.