A simulation method for solar thermosyphon collector