This paper presents the detection and interpretation of numerical power density anomalies in the full-core PWR reactor model, developed at the level of RPV. The study shows that in the Monte Carlo neutron transport coupled with burnup modeling, power density anomalies may occur and should be carefully investigated. The power density anomalies originate from the initial core asymmetry due to different designs of top and bottom reflectors and non-uniform initial axial neutron flux distribution. The power density anomalies increase with increasing fuel burnup due to spatial changes in 135Xe concentrations and following power density redistribution along the reactor core. A method for anomaly detection for the semi-symmetric reactor core is proposed.