Aluminum production using the Hall-Héroult process requires an intensive electric current. With an increasing demand for aluminum due to the growth in use of raw materials such as aluminum, and economic developments, aluminum producers are striving to reduce their production costs to remain competitive in a difficult market. One way of doing this is to reduce the voltage drop in the Hall-Héroult cells, which yields a remarkable amount of savings in the long term from only a slight optimization. For this reason, a thorough understanding of the phenomena taking place during operation is necessary.In this research work, anode assembly as one of the major components of the HallHéroult cell was modeled using APDL (ANSYS® Parametric Design Language). The newly presented features of the full anode assembly model, which make it a robust one in terms of geometrical modeling, were introduced in detail. A submodel was extracted from the fall model to carry out primary numerical simulations to investigate the ThermoElectro-Mechanical (TEM) phenomena taking place in the stub hole region to where up to 25% of the total voltage drop in the anode assembly, caused by contact resistance at the interfaces, is attributed. Special attention was paid to the good prediction of contact conditions at the cast iron to carbon interface. In comparison to previous research work, a more thorough and precise approach was taken to employ equations used to predict the initial air gap at the cast iron to carbon interface, which has an influential role in controlling contact particularly at lower temperatures (400°-600°