Peak heat flow occurs at mid-ocean ridges and decreases with the age of the oceanic lithosphere. Several plate models, including the Parsons and Sclater (PSM) model, Global Depth and Heat (GDH1) model and Constant Heat flow Applied on the Bottom Lithospheric Isotherm (CHABLIS) model, have been used to predict heat flow in the ocean lithosphere. The discrepancy between the predicted and measured heat flow in the younger lithosphere (i.e. younger than 55 Myr) influenced by local hydrothermal circulation has been used to estimate hydrothermal heat flux and investigate hydrothermal processes. We can modify the cooling models by substituting the ordinary mass density of lithosphere by fractal density with singularity. This new model provides a modified solution to fit the observed heat flow data used in other models in the literature throughout the age range. This model significantly improves the results for prediction of heat flow that were obtained using the PSM, GDH1 and CHABLIS models. Furthermore, the heat flow model does not exhibit special characteristics around any particular age of lithosphere. This raises a fundamental question about the existence of a "sealing" age and accordingly the hydrothermal flux estimation based on the cooling models.