Photopolymerization processes play an integral role in the manufacturing and performance of modern inkjet printing technologies, enabling curing of materials under light exposure. This study presents a baseline mathematical model for photopolymerization in inkjet printers, formulated as a multi-physics problem encompassing chemical reactions, diffusion, and electromagnetic phenomena. The model is constructed within a thermodynamical framework, ensuring consistency with the first and second laws of thermodynamics. The entropy of mixing is incorporated to describe species interactions during polymerization, while conservation laws for mass and energy guide the derivation of species fluxes. The framework captures the interplay between light absorption, radical formation, and monomer diffusion by coupling electromagnetic field equations with chemical kinetics. Numerical simulations, implemented via the finite element method, solve the nonlinear coupled equations. Results illustrate the spatial and temporal evolution of critical variables, including concentrations of monomers, radicals, photoinitiators, and temperatures. Furthermore, the impact of varying light intensity profiles and material properties on curing efficiency and final polymer structure is analysed, providing insights into process optimization. This baseline model establishes a robust foundation for future work, including the development of an inverse problem to further optimize photopolymerization processes. The presented framework applies not only to inkjet printing but also to other photopolymerization-based manufacturing processes. This study aims to advance the understanding and control of photopolymerization in complex engineering systems by bridging fundamental research and industrial applications.