A metabolic heat-based model was used for estimating the growth of Kluyveromyces marxianus, and the modified Luedeking-Piret kinetic model was used for describing the inulinase production kinetics. For the first time, a relationship was developed to relate inulinase production kinetics directly to metabolic heat generated, which corroborated well with the experimental data (with R values of above 0.9). It also demonstrated the predominantly growth-associated nature of the inulinase production with Luedeking-Piret parameters α and β, having values of 0.75 and 0.033, respectively, in the exponential feeding experiment. MATLAB was used for simulating the inulinase production kinetics which demonstrated the model's utility in performing real-time prediction of inulinase concentration with metabolic heat data as input. To validate the model predictions, a biocalorimetric (Bio RC1e) experiment for inulinase production by K. marxianus was performed. The inulinase concentration (IU/mL) values acquired from the model in were validated with the experimental values and the metabolic heat data. This modeling approach enabled the optimization, monitoring, and control of inulinase production process using the real-time biocalorimetric (Bio RC1e) data. Gas chromatography and mass spectrometry analysis were carried out to study the overflow metabolism taking place in K. marxianus inulinase production.