Histotripsy with ultrasound is an emerging noninvasive therapeutic modality that uses cavitation to precisely destroy diseased soft tissue. Accurate simulations of histotripsy are needed for treatment planning and device design. These simulations are performed in the time-domain, span hundreds of wavelengths, and must handle strong shocks and discontinuities between materials, such as the brain and the skull. The discontinuous Galerkin (DG) method is an outstanding candidate for such simulations. DG methods possess the following qualities: 1) high order accuracy, 2) geometric flexibility, 3) excellent dissipation properties, and 4) excellent scalability on massively parallel machines. The objective of this work is to develop a massively parallel DG method for histotripsy simulations in the brain.