Drag is a common topic in everyday activities and academic research. Drag crisis, which associates with the transition of boundary layers from laminar to turbulent, is a phenomenon that can gain advantages when increase velocity to the critical value for blunt bodies. Meanwhile, the appearance of laminar separation bubble (LSB) is thought to be the cause of drag crisis. The mechanisms, and structure as well as the influencing factors of laminar separation bubble have also been comprehensively studied. The paper focuses on the numerical approach to simulate three - dimensional (3D) teardrop model for the study of aerodynamic characteristic and LSB mechanisms. Results from the open-source code OpenFOAM are compared to that of the experiments. The 3D simulation results are exploited to analyze the velocity vector field, pressure distribution, turbulent flow field ... in the critical Reynolds regime, in which the drag crisis occurs, from Re = 4×10^4 to Re = 7×10^4. At the same time, the appearance of laminar separation bubble in the critical Reynolds regime is also analyzed. The results highlight the advantages of the transition k-ω SST turbulence model compared to the k-ω SST turbulence model in simulating the drag crisis phenomenon and the footprint of bubble formation.