With the increasing demands for groundwater management in industrial and environmental engineering, the research and application of air sparging remediation technology have become crucial for enhancing multiphase flow control. Particularly under complex geological conditions with temperature gradients, understanding and predicting the behavior of multiphase flows are essential for improving the efficiency and safety of remediation technology. Currently, most studies focus on the analysis of multiphase flows under static conditions, with fewer assessments on the impact of dynamically changing temperatures and pressures, thus limiting the accuracy and reliability of models in practical applications. Addressing this research gap, this study proposes a fully transient temperature-pressure field coupled model, along with a corresponding iterative solution algorithm, to accurately simulate the dynamic characteristics of underground multiphase flows under varying temperature gradients. The progress of this research extends our understanding of the mechanisms of multiphase flow in air sparging, supporting technological innovation and environmental protection efforts in related fields. It also offers new solutions for resource development and pollutant management under complex geological conditions.