#matlab #adsorption #curvefitting
Inspired by : Saikat Roy & Arun S. Moharir (2019). Modeling the Generic Breakthrough Curve for Adsorption Process; https://arxiv.org/pdf/1907.00195
clc;clear;close all;
%parameter
L=0.3; %m
eps=0.4; %-
u=0.005; %m/s
LDF=0.5; %1/s LDF=15*De/Rp^2 (LDF = for simplification)
b=0.1; %m3/kg
qm=0.02; %kg adsd/kg adst
c0=1.5; %kg/m3
rhop=1100; %kg/m3
%time-mesh
t=linspace(0,500,51);
ta=LDF.*t; %ta=t'
%space-mesh
Nz=51;
z=linspace(0,L,Nz);
za=LDF.*z/u; %za=z'
dza=za(2)-za(1); %dimensionless space-difference
%another dimensionless unit
rhopa=(1-eps)/eps*rhop/c0; %rhopa=rhop'
ba=b*c0; %ba=b'
