Direct numerical simulations are performed to explore the evolution behaviour of the turbulent/non-turbulent interface (TNTI) in a temporally evolving turbulent plane jet, using the evolution equation for the TNTI surface area. A novel algorithm is used to calculate the surface area of the TNTI and entrainment flux. It is shown that the surface area remains relatively constant, which leads to the mean entrainment velocity being inversely proportional to the square root of time. On average, the effects of the stretching and curvature/viscous terms on the TNTI area roughly counterbalance each other, while the curvature/inviscid term associated with vortex stretching is virtually zero. More specifically, the stretching term contributes to the production of the surface area, while the curvature/viscous term is associated with a destruction in the surface area. The local effect of the curvature/viscous term exhibits high spatial intermittency with small-scale extreme/intense events, whereas the effect of the large-scale stretching term is more continuous. To shed light on the contribution of curvature/viscous term to the evolution of the surface area, we decompose it into three components. The effect of the curvature/normal diffusion term (the curvature/viscous dissipation term) in the bulging regions (the valley regions) mainly contributes to the production of the area. The continuous decrease of the average mean curvature is associated with the production of the bulging regions and the destruction of the valley regions. Finally, although the entrainment velocity is mainly dominated by the normal diffusion effect, all three components related to the viscous effect are indispensable to the production and destruction of the TNTI area. This numerical study contributes to a better understanding of the evolution of the TNTI area. © The Author(s), 2024.