“…In this work, a complete formulation and implementation of finite-difference method based thermoelectric model is presented. The model contains the system of PDEs for the typical TE module (figure 1): two heat equations (one for the hot side substrate and the other for the cold side substrate), pair of PDEs for both P and N type thermoelectric materials (one for the heat transport, the other for the electric field), set of boundary equations for boundaries A, B, C, and D. The main advantage of the proposed finite-difference model over the previous implementations [15] is that, in this formulation, electric field and heat transport in both P and N type materials are handled separately, and the temperature dependencies of thermal conductivity, electrical conductivity, and Seebeck coefficient of the P/N materials are addressed directly. PDEs, which are timedependent, are discretised using the finite-difference 'explicit' method, which is particularly well-suited in solving heat-related high-speed dynamic events [17].…”