I'm simulating a cold flow analysis (methane-air) on fluent, using species transport model. The model has around 5 millions meshes with tetrahedrons and swept bodies. I use "realizable k-epsilon" turbulence model with non-equilibrium wall function. Though the solution is converged, the results look nonphysical in the combustor. Any body who can give me some tips?