Applying ultra-high dose rates to radiation therapy, otherwise known as FLASH, has been shown to be just as effective while sparing more normal tissue compared to conventional radiation therapy. However, there is a need for a dosimeter that is able to detect such high instantaneous dose, particularly in vivo. To fulfill this need, protoacoustics is introduced, which is an in vivo range verification method with submillimeter accuracy. The purpose of this work is to demonstrate the feasibility of using protoacoustics as a method of in vivo real-time monitoring during FLASH proton therapy and investigating the resulting protoacoustic signal when dose per pulse and pulsewidth are varied through multiple simulation studies. The dose distribution of a proton pencil beam was calculated through a Monte Carlo toolbox, TOPAS. Next, the k-Wave toolbox in MATLAB was used for performing protoacoustic simulations, where the initial proton dose deposition was inputted to model acoustic propagations, which were also used for reconstructions. Simulations involving the manipulation of the dose per pulse and pulsewidth were performed, and the temporal and spatial resolution for protoacoustic reconstructions were investigated as well. A 3D reconstruction was performed with a multiple beam spot profile to investigate the spatial resolution as well as determine the feasibility of 3D imaging with protoacoustics. Our results showed consistent linearity in the increasing dose-per-pulse, even up to rates considered for FLASH. The simulations and reconstructions were performed for a range of pulsewidths from 0.1 to 10μs. The results show the characteristics of the proton beam after convolving the protoacoustic signal with the varying pulsewidths. 3D reconstruction was successfully performed with each beam being distinguishable using an 8cm×8cm planar array. These simulation results show that measurements using protoacoustics has the potential for in vivo dosimetry in FLASH therapy during patient treatments in real time. Through this simulation study, the use of protoacoustics in FLASH therapy was verified and explored through observations of varying parameters, such as the dose per pulse and pulsewidth. 2D and 3D reconstructions were also completed. This study shows the significance of using protoacoustics and provides necessary information, which can further be explored in clinical settings.