A modification to the f-Chart method has been developed to predict monthly and annual performance of thermosyphon solar domestic hot water systems. Stratification in the storage tank is accounted for through use of a modified collector loss coefficient. The varying flow rate throughout the day and year in a thermosyphon system is accounted for through use of a fixed monthly “equivalent average” flow rate. The “equivalent average” flow rate is that which balances the thermosyphon buoyancy driving force with the frictional losses in the flow circuit on a monthly average basis. Comparison between the annual solar fraction predited by the modified design method and TRNSYS simulations for a wide range of thermosyphon systems shows an RMS error of 2.6 percent.