Although convection-enhanced delivery can successfully bypass the blood-brain barrier, its clinical performance remains disappointing. This is primarily attributed to the heterogeneous intratumoral environment, particularly the tumor microvasculature. This study investigates the combined convection-enhanced delivery of antiangiogenic drugs and liposomal cytotoxic drugs in a heterogeneous brain tumor environment using a transport-based mathematical model. The patient-specific 3D brain tumor geometry and the tumor’s heterogeneous tissue properties, including microvascular density, porosity and cell density, are extracted from dynamic contrast-enhanced magnetic resonance imaging data. Results show that antiangiogenic drugs can effectively reduce the tumor microvascular density. This change in tissue structure would inhibit the fluid loss from the blood to prevent drug concentration from dilution, and also reduce the drug loss by blood drainage. The comparisons between different dosing regimens demonstrate that the co-infusion of liposomal cytotoxic drugs and antiangiogenic drugs has the advantages of homogenizing drug distribution, increasing drug accumulation, and enlarging the volume where tumor cells can be effectively killed. The delivery outcomes are susceptible to the location of the infusion site. This combination treatment can be improved by infusing drugs at higher microvascular density sites. In contrast, infusion at a site with high cell density would lower the treatment effectiveness of the whole brain tumor. Results obtained from this study can deepen the understanding of this combination therapy and provide a reference for treatment design and optimization that can further improve survival and patient quality of life.