OpenFoam - cryoablation simulation
 

https://ibb.co/4g14DH4

C *∂T/∂t=∂/∂x (Kx * ∂T/∂x)+ ∂/∂y (Kx * ∂T/∂y)+ ∂/∂z (Kx * ∂T/∂z) + Cb
ωb(T)*(Ta-Tt )+Q

Ta is the arterial temperature, Tt is the temperature at time t, Q is the metabolic heat rate of tissue and T is the tissue temperature.

I have to solve this equation to simulate cryoablation of a tumor.
I am confused about which all modules should I combine so as to get the simulation right.
Thank you in advance.