Flows within porous media play important roles in many scientific and industrial systems. However, the case wherein such flows become turbulent has not been completely understood, particularly, from a mathematical viewpoint. In this study, $k-\varepsilon$ model (the variable $k$ denotes kinetic energy per unit mass $\varepsilon$ dissipation rate for energy), which widely used usual flow clear f...