Interstitial photodynamic therapy (iPDT) with 5-aminolevulinic acid (ALA) is a possible alternative treatment for malignant brain tumors. Further evaluation is, however, required before it can be clinically applied. Computational simulation of the photophysical process in ALA-iPDT can offer a quantitative tool for understanding treatment outcomes, which depend on various variables related to clinical treatment conditions. We propose a clinical simulation method of ALA-iPDT for malignant brain tumors using a singlet oxygen (O1₂) model and O1₂ threshold to induce cell death. In this method, the amount of O1₂ generated is calculated using a photosensitizer photobleaching coefficient and O1₂ quantum yield, which have been measured in several previous studies. Results of the simulation using clinical magnetic resonance imaging data show the need to specify the insertion positions of cylindrical light diffusers and the level of light fluence. Detailed analysis with a numerical brain tumor model demonstrates that ALA-iPDT treatment outcomes depend on combinations of photobleaching and threshold values. These results indicate that individual medical procedures, including pretreatment planning and treatment monitoring, will greatly benefit from simulation of ALA-iPDT outcomes.