Bone fracture healing is one of the important topics in biomechanics, demanding computation simulations due to the difficulty of obtaining experimental or clinical results. In this study, we adopt the design space optimization method which was established by the authors as a tool for the simulation of bone growth using its evolutionary characteristics. As the mechanical stimulus, strain energy density is used. We assume that bone tissues over a threshold strain energy density will be differentiated and bone tissues below another threshold will be resorbed. Under compression and torsion as loadings, the filling process of the defect is well illustrated following the given mechanical criterion. It is shown that the design space optimization is an excellent tool for simulating the evolutionary process of bone growth, which has not been possible otherwise.