We propose a model that explains the fermion mass hierarchy by the Froggatt-Nielsen mechanism with a discrete Z F N symmetry. As a concrete model, we study a supersymmetric model with a single flavon coupled to the minimal supersymmetric Standard Model. Flavon develops a TeV scale vacuum expectation value for realizing flavor hierarchy, an appropriate µ-term and the electroweak scale, hence the model has a low cutoff scale. We demonstrate how the flavon is successfully stabilized together with the Higgs bosons in the model. The discrete flavor symmetry Z F N controls not only the Standard Model fermion masses, but also the Higgs potential and a mass of the Higgsino which is a good candidate for dark matter. The hierarchy in the Higgs-flavon sector is determined in order to make the model anomaly-free and realize a stable electroweak vacuum. We show that this model can explain the fermion mass hierarchy, realistic Higgs-flavon potential and thermally produced dark matter at the same time. We discuss flavor violating processes induced by the light flavon which would be detected in future experiments. N sym-1 Recently, it is proposed that a flavon can be identified as the QCD axion [3-8] 2 See, e.g., Refs. [35,36] for recent models with a flavor modular symmetry and also for such interactions in string models.3 Models with a combination of Higgs doublets H u H d as a flavon are studied in Refs. [41][42][43].