A new set of computational techniques has been devised to directly couple the local, transient flow conditions and the local, transient conditions of particles and other bodies interacting with the flow. The direct coupling method developed may be used for particle laden flow with large particle sizes, high particle concentrations, and frequent particle interactions, where the traditional Eulerian-Lagrangian methods may not be appropriate. In the technique, the particles/solids are discretized in the computational domain. The particles move within the fluid, and the local, internal temperature distributions within the particles is numerically solved for simultaneously with the flow. Thus, the ability to access local flow and temperature variation throughout the system, including within the particles and other solid bodies, facilitates the study of thermally dependent phenomena in a manner not previously achieved. Fluid induced forces and fluid-particle heat exchange can be directly calculated and models can be implemented to simulate the particle-particle/solid body mechanical and thermal interactions. In this work, the capabilities of this method are demonstrated through the study of the effects of temperature dependent material properties, adhesion forces, and phase changes in mainly the compressible, turbulent flow regime. However, the technique can be extended to other phenomena such as chemically reacting flows, erosion and fouling, and other flow regimes. The expanded capabilities to study temperature dependent phenomena in particle laden flow offered by the method developed can provide the information needed for more effective and efficient utilization and control of particle laden flows.