This thesis is based on a project "Bifacial Photovoltaic Energy Production Analysis" to build a detailed simulation model system accurately simulate bifacial panel performance under real field radiation conditions and deployment configuration, and to predict its corresponding energy yield. To the author’s up-to-date knowledge, the model system is unpreceded among same type simulation software in complexity, details in consideration, ranges of deployment and parameters. The model system can also be used as a platform for more components and variables to be added on, such as adding on more rows of panel arrays to simulate bifacial solar farm scenario; and adding spectral information for more accurate analysis. The system components’ sub-models were carefully chosen based on a broad literature review in related aspects; especially in sky diffuse radiance, ground reflection, and bifacial solar cells. Built in MATLAB© based on mathematical expressions from above said models, the system consists of 5 bifacial panels and their racking as shading objects and the central panel performance is under investigation and has taken consideration of all possible panel azimuth and elevation combinations. Model simplification and resolution are carefully considered so to achieve a good balance in complexity, computation load and output accuracy. Output reliability is confirmed with other people’s work. Furthermore, the model has been fully checked and peer tested. Outputs under different parameter settings are analysed and discussed. Conclusions and recommended future work are provided at the end of the thesis.