We present numerical simulation of the nitrogen atmospheric pressure plasma jet (APPJ) using three fluid models-namely, laminar model, modified laminar model, and turbulent (k-epsilon) model-coupled with gas-phase reaction kinetics. The spatial profiles of the light emission intensities, gas temperature, and NO density predicted by the turbulent model show a better agreement with the experimental observations, compared with laminar and modified laminar models. We have demonstrated that the turbulent model shows more oxygen entrainment, more mixing with the ambient air, and a lower axial velocity at the downstream. These allow the turbulent model to more precisely capture the APPJ characteristics than the laminar and modified laminar models do.