Governing differential equations for natural circulation systems