I would couple two different lattices (Navier-Stockes) for velocity and Advection-Diffusion for temperature. I would put a constant temperature around the sphere using a data processor as a boundary condition. The last thing is to couple the two lattices using again a data processor or using the already existing function “latticeToPassive…”