Zero-dimensional two-stage SOFC stacks dynamic model was developed to investigate the effect of operating parameters on stacks performance. The model resolves spatially thermal and thermo-electrochemical behaviour for electrochemical reactions, Catalytic Partial Oxidation and Steam Reforming processes. Design variables and thermo-electrochemical properties were obtained from in-house-fabricated SOFCs carried out by project partners. The completed SOFCs based Combined Heat and Power, CHP, system model was validated by data