This work presents an unsteady state theoretical and experimental investigation of natural circulation two phase closed thermosyphon flat plate solar water heaters. The governing equations of the heater storage tank and connecting pipes are presented and generalized in dimensionless form, while the governing equations of the different components of the collector were previously discussed by the author. Also, the author's earlier simulation program of the collector is modified to include the solution of the dimensionless governing equations of the present analysis. For verifying the modified simulation program, a two phase closed thermosyphon solar water heater is designed, constructed and tested at different meteorological conditions, initial storage tank temperatures and different hot water withdrawal load patterns. The comparison between the experimental results and their corresponding simulated ones shows considerable agreement.