Linear stability analysis of National Spherical Torus Experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of extended magneto-hydrodynamic (MHD) model in NIMROD. The purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after Licoating of divertor and first wall of the NSTX tokamak. Besides ideal MHD modeling including finite-Larmor radius effect and 2-fluid Hall and electron diamagnetic drift contributions, non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. Unlike earlier conclusion from eigen value code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal 2-fluid limit. After considering the higher plasma resistivity due to higher Z eff , the complete stabilization could be explained. A thorough analysis of both ideal and non-ideal models is presented, after accurately including vacuum-like cold halo region in NIMROD to investigate ELMs.