So it does affect I assume due to density changes? if yes consider the density URF.
Do you have good mesh?
Which Solver are you using?
Regarding source term linearization: you should always try to initialize the source terms otherwise they are handled explicitly and this not a good idea.