Continued Burning

In trying to finish up my thesis, I have been trying to understand (for the last several months) why the gas that I am modeling numerically does not cool down to the constant temperature boundary. The model represents gas burning in a closed spherical vessel 10 cm in diameter. When the gas reacts, it gives off heat and heats all the gas in the domain. Once the heating reaction(s) stops or slows down, all the gas in the domain should cool to the temperature at the boundary. In my case, it doesn't.

At first, when talking with other people, I thought I might have made an error in my equations or the equations were specific to a certain geometry. Just last week, though, I had a professor look over the equations and we realized that the model provides a continuous supply of reactants.

Going back to the original paper that explains the numerical model, I discovered that the authors had mentioned that they had written the equations with a continuous supply of reactants!

It pays to read and reread the directions!!