Heat exchange between connected particles due to conduction.

$$\dot{Q}_{cond} = k A_{b} \kappa_{b} \left( T_{p1} - T_{p2} \right) / L_{b}$$

Integration of heat is done as:

$$\dot{Q}_{p1} = -\dot{Q}_{cond}$$ $$\dot{Q}_{p2} = +\dot{Q}_{cond}$$ $$T_{p1,p2} = T_{p1,p2} + \frac{\dot{Q}_{p1,p2}}{C_{p1,p2} m_{p1,p2}} \Delta t_{sim}$$

Symbol	Parameter name	Description
$\kappa_{b}$		Bond thermal conductivity [W/m/K]
$A_{b}$		Bond cross-cut area [m2]
$C_{p1}$, $C_{p2}$		Particle heat capacities [J/kg/K]
$k$	CONDUCTION_SCALING_FACTOR	Scaling factor for conductivity [-]
$L_{b}$		Bond length [m]
$m_{p1}$, $m_{p2}$		Particle masses [kg]
$\dot{Q}_{p1}$, $\dot{Q}_{p2}$		Particle heat flows [W]
$\dot{Q}_{cond}$		Heat flow due to conduction [W]
$\Delta t_{sim}$		Simulation time step [s]
$T_{p1}$, $T_{p2}$		Particle temperatures [K]