Hi, I'm working on modeling a multiphase flow with reactions occurring in the gaseous phase only. I have a well converged solution for the flow with inter-phase mass transfer between the liquid flow and the gaseous flow which I am using as my initial conditions. Yet when I introduce the reactions, the solution oscillates about a constant value. I've attempted refining the mesh, but that didn't provide me with any better results.

I've run other simulations with reactions without problem before, but this case is causing me a lot of grief. If anyone has any insights in dealing with this type of system I'd appreciate it.

Did you change the time step also?