In the past two decades, High-Intensity Focused Ultrasound (HIFU) has been actively investigated for inducing tissue ablation for tumours, healing arteries and even skin tightening for inducing facelifts. Being a non-contact and non-invasive procedure, HIFU is an attractive alternative over other procedures like radio-frequency ablations and microwave ablations as these are semi-invasive methods. Computer simulations in this domain can provide us with a fast-track way to test and deduce results virtually before conducting related experiments. In this study, we simulated the properties of bovine liver tissue virtually in a MATLAB-based k-Wave toolbox to track lesion variations with varying HIFU on-time and with a constant amount of heat deposition. Photoacoustic reconstructions were also carried out with a linear and circular array of transducers using the simple Delay-And-Sum (DAS) algorithm to access the real-time imaging capability of this technique for HIFU therapy. All the results obtained are discussed and validated. This study will, in general, aid in deciding the HIFU configuration and operational parameters with respect to the tissue type, lesion size, and lesion location within the tissue.