In this paper I present a general modelling framework for coupled fluid dynamics and chemistry problems, and apply it to the simulation of a series of complex, homeostatic reaction diffusion systems. The model can incorporate any number of chemical species and reactions. Those chemical species diffuse, react and are advected by fluid flows. I illustrate some characteristic results from the modelling of the Gray Scott reaction diffusion system with thermally resolved reactions. Extending my previous work on ecological dynamics of nonliving structures, I demonstrate that thermal homeostasis of reaction diffusion spots can occur in systems without the use of the porous wall boundary condition that has traditionally been used for the Gray Scott system. I present an initial analysis of the parameter space of this system as well as detailing the mechanism behind the thermal homeostasis.