The turbulent counter-flow flame (TCF) has proven to be a useful benchmark to study turbulence-chemistry interactions, however, the widely observed bulk flow fluctuations and their influence on the flame stability remain unclear. In the present work, premixed TCFs are studied numerically using a Large Eddy Simulation (LES) method. A transported probability density function (pdf) approach is adopted to simulate the sub-grid scale (sgs) turbulence-chemistry interactions. A solution to the joint sgs-pdf evolution equation for each of the relative scalars is obtained by the stochastic fields method. The chemistry is represented using a simplified chemical reaction mechanism containing 15 reaction steps and 19 species. This work compares results with two meshing strategies, with the domain inside nozzles included and excluded respectively. A conditional statistical approach is applied to filter out the large scale motions of the flame. With the use of digital turbulence, the velocity field in the flame region is well reproduced. The processes of local extinction and re-ignition are successfully captured and analysed together with the strain rate field, and local extinctions are found correlated to the turbulent structures in the reactant stream. The predicted probability of localised extinction is in good agreement with the measurements, and the influence of flame stoichiometry are also successfully reproduced. Overall, the current results serve to demonstrate the capability of the LES-pdf method in the study of the premixed opposed jet turbulent flames.