Thermal battery modeling is important for further battery development and optimization. The temperature strongly influences the performance and aging behavior. In the cell stack, electrochemical processes take place resulting in a large amount of heat release, which, in turn, affects the temperature distribution. Therefore, the main focus is on the cell stack, the most complex structure inside the cell. In particular, the discontinuous and anisotropic material properties represent a major challenge for simulations due to the layering. This work proposes self-developed methods, based on the Finite Volume Method and the Finite Element Method, taking on these challenges. First, for both methods the functionality is verified and numerical convergence is validated. These, and also classical methods, are compared based on test problems with a known analytical solution in view of numerical errors as well as computing time. It if found that their accuracy and efficiency depends strongly on the specific problem, which makes their numerical investigation necessary and inevitable. Second, the methods are evaluated on a specific battery problem. ... mehrTheir results are plausible and correspond to the physical phenomena.