Study on the stability behaviour of two-phase natural circulation systems using a four-equation drift flux model