Random nucleation

Dear all,

I am currently involved in a phase change simulation.


I have not yet found the way to deal with random nucleation occurring for instance during the quenching : during the transient heat transfer calculation, I would like to randomly generate at each time step crystal nuclei at a spesific rate of nucleation ( dependant of : temperature, number of nucleation, the activatin energy, many parametres...)

I am aware of the random function but cannot figure out how to use it in this specific context.


Any help would be much appreciated !