  In this paper, we present TLBfind, a GPU code for simulating the
hydrodynamics of droplets along with a dynamic temperature field. TLBfind
hinges on a two-dimensional multi-component lattice Boltzmann (LB) model
simulating a concentrated emulsion with finite-size droplets evolving in a
thermal convective state, just above the transition from conduction to
convection. The droplet concentration of the emulsion system is tenable, and at
the core of the code lies the possibility to measure a large number of physical
observables characterising the flow and droplets. Furthermore, TLBfind includes
a parallel implementation on GPU of the Delaunay triangulation useful for the
detection of droplets' plastic rearrangements, and several types of boundary
conditions, supporting simulations of channels with structured rough walls.
