The present study is focused on the development of the RIF (Representative Interactive Flamelet) model which can overcome the shortcomings of conventional approach based on the steady flamelet library. Due to the ability for interactively describing the transient behaviors of local flame structures with CFD solver, the RIF model can effectively account for the detailed mechanisms of $NO_x$ formation including thermal NO path, prompt and nitrous $NO_x$ formation, and reburning process by hydrocarbon radical without any ad-hoc procedure. The flamelet time of RIFs within a stationary turbulent flame may be thought to be Lagrangian flight time. In context with the RIF approach, this study adopts the Eulerian Particle Flamelet Model (EPFM) with mutiple flamelets which can realistically account for the spatial inhomogeneity of scalar dissipation rate. In order to systematically evaluate the capability of Eulerian particle flamelet model to predict the precise flame structure and NO formation in the multi-dimensional elliptic flames, two methanol bluffbody flames with two different injection velocities are chosen as the validation cases. Numerical results suggest that the present EPFM model has the predicative capability to realistically capture the essential features of flame structure and $NO_x$ formation in the bluff-body stabilized flames.