In this paper, an efficient, numerical procedure is presented to track crack growth modeling without remeshing. The method's key feature is the coupling of a modified quarter-point crack tip element (MQPE) with the level set method (LSM) for crack growth problems. The LSM was used to represent the crack location, including the location of crack tips. The MQPE was used to compute the stress and displacement fields necessary for determining the rate of crack growth. Numerical test cases including various geometrical exceptions (the center-crack plate specimen, the single edge-crack plate specimen, and the double-edge crack plate) demonstrate the accuracy, robustness, and efficiency of the MQPE/LSM coupling. The extrapolation technique was used to estimate numerically the calibration factor for various specimens. This work confirms the feasibility of the MQPE/LSM to model accurately the singularity existing in the vicinity of the cracks. It allows an economic and adequate calculation of the stress intensity factors, which can be introduced into the various criteria of fracture or laws of propagation of the crack. The new method reduces the need for remeshing, and results agree well with reference data.