Partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics in order to develop a prediction model for turbulent flame lift off. The essence of flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of flame holes initially created by local quenching events. The numerical simulation for flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for constant-density fuel-air channel mixing layer to obtain the turbulent flow and mixing fields, from which a time series of two dimensional scalar dissipation rate array is extracted at a fixed virtual flame surface horizontally extending from the end of split plate to the downstream. Then, the Lagrangian simulation of the flame hole random walk mapping projected to the scalar dissipation rate array yields temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. The statistical results exhibit that the chance of partial quenching is strongly influenced by the crossover scalar dissipation rate while almost unaffected by the iteration number of the mapping that can be regarded as a flame-edge speed.