The electromagnetic confinement directional solidification (EMCDS) technique is an optimal method for preparing large-size and non-contamination directional TiAl alloy crystals. Despite its advantages, the high temperature gradient inherent to this process induces thermal stress within the ingot, which increases the risk of cracking. To address this challenge, an innovative Integrated Multi-Physics Coupling Model was established to map and study the thermal stress field during EMCDS in this study. It synchronized the computation of the electromagnetic field, temperature field, solute field, flow field, and stress field during the crystal growth of TiAl alloy, and its high accuracy was proved by the micro-indentation experiment. Our analysis reveals that transverse temperature differences are crucial in inducing thermal stresses, and identifies that hot cracks and cold cracks are prone to occur respectively at the area of radial 23R along the sample (X = 23R) and on the sample surface, which aligns with experimental observations impressively. This model can more accurately and efficiently optimize the process parameters of the EMCDS process to avoid cracks and promote its industrial application.