A phase angle drift method for loss of mains/grid protection is disclosed. According to one aspect, an accumulated electrical phase angle drift derived from the difference between the current measured local frequency and the estimated frequency using historical data is compared to an angle threshold. An estimated grid frequency may be calculated based on the historical delay, and the window, over which the estimated frequency is calculated. An addition/subtraction of a phase angle offset value is calculated for a half cycle is performed when the frequency difference between the estimated frequency fnest and the measured frequency fn is greater or equal to a first determined value.