Based on the stability analysis with no linearization and expansion,it is argued that instability in the lattice BGK model is originated from the linearrelaxation hypothesis of collision in the model. The hypothesis stands up only whenthe deviation from the local equilibrium is weak. In this case the computation is abso-lutely stable for real fluids. But for flows of high Reynolds number, this hypothesis isviolated and then instability takes place physically. By performing a transformationa quantified stability criteria is put forward without those approximation. From thecriteria a sufficient condition for stability can be obtained and serve as an estimationof the limited Reynolds number as high as possible.