Mohammad Asif


I am solving the adsorption of CO2 taking place in an adsorption bed which is assumed as porous, as shown in the attached picture in my original post. For better understanding, I am sharing my UDF here. The mass balance is looking fine. One more problem is that my program is taking too much time, even though the mesh size is not that big, porosity is 0.4 but I am using default values of Viscous resistance and Inertial resistance. I am adding the adsorption rate of CO2 as a mass source term in the porous region.

I am using linear driving force adsorption model-> dq/dt= k(q^* – q), q^* is calculated by Langmuir isotherm and then adding (1-porosity)*rho*dq/dt as a mass source in porous zone.