Keeping in mind that evaporation is a surface phenomenon, it will more strongly affect the top of the slurry, while its bottom will be initially less affected. The key idea of the procedure proposed here is to module the speed of CBDs shrinking and density increase (i.e. the solvent removal in our computational workflow) as a function of their position. Particularly, CBD particles at the top of the slurry will shrink faster than the ones in the bottom. In this work, as illustrated in Figures and, the CBD particles are classified in three groups as a function of their z-coordinate along the slurry thickness: the ones in the first third (bottom, referred as CBD1), the ones in the second third (center, referred as CBD2) and the ones in the last third (top, referred as CBD3). Then, we defined the shrinking speed (𝐯 𝑏𝑜𝑡𝑡𝑜𝑚 ) of CBD1 in such a way to remove all the solvent at the end of the simulation. In other terms, and considering that CBD1 particles are the ones with the slower shrinking speed, the latter means assuming that the solvent removal is complete at the end of the drying process. Then, a relative shrinking factor (RSF) is defined for the center and top CBDs as