Mathematica's NDSolve is finding a stiff system very quickly trying to solve a large set of differential equations

Question

I\'m a physics PhD student solving lubrication equations (non-linear PDEs) related to the evaporation of binary droplets in a cylindrical pixel in order to simulate the shap