An efficient computational fluid dynamics (CFD) method for simulating the flow and convective heat transfer process of airflow in a tunnel is required to analyze the freezing and thawing of surrounding rock and to apply the results to the design of the insulation layer for a tunnel located in a cold region. Comparisons of experimental data and CFD results using a standard k − ε turbulence model, a wall function, a thermal function and an adaptive finite element method are presented. Comparison of the results indicated that the proposed model and simulation method are efficient at determining the solid–air interface heat coefficient in a thin and infinitely wide horizontal plate and the hydrodynamic and thermal fields in a 3-D cavity. After demonstrating that the necessary validations are satisfied, this paper presents an analysis of the characteristics of airflow and air–rock heat transfer in a cold-region tunnel.