Heat Conduction Model

The Heat Conduction contact model calculates the heat flux based on the relative temperatures and the particle overlap.

This model calculates the heat flux between Particle-Particle and Particle-Geometry contacts. This model must be used with Update Temperature particle body force to complete the calculation. The Temperature Update particle body force allows the inclusion of an external heat source and is where each particle temperature is calculated based on the heat flux calculated in the contact model and a given external heat flux.

The calculated heat flux is applied to particles only. You can assign the Geometry as a heat source. The Geometry temperature does not change whereas the particle temperature will be updated based on the calculated heat flux.

For dilute phase simulations, convective heat transfer is dominant and conduction between the particles or wall can be neglected. However, for dense phase, contacts between particles are significant such that conductive heat transfer must be taken into account. A single phase DEM simulation on heat transfer in granular flow in rotating vessels provides a simple approach in modeling inter-particle heat transfer. This model is based on the work of (H&A) Chaudhuri (Chaudhuri, Muzzio, and Tomassone 2006).

The heat flux between the particles is defined as:

Qp1p2=hcΔTp1p2

Where the contact area is incorporated in the heat transfer coefficient c and is defined as:
hc=4kp1kp2kp1+kp2[3FNr*4E*]13
 Where FN is the normal force, r* the geometric mean of the particles radii from the Hertz’s elastic contact theory and E* is the effective Young’s modulus for the two particles. The bracketed term on the RHS of the equation models indicates the contact area between particles.
Note: You must use this contact model along with the Temperature Update Body Force model.