“…PSA for hydrogen or air separation, the partition parameter K c is relatively small, and the order of time-derivative slope, O(1= ), is typically less than 10. Consequently, algorithms without rigorous consideration on the diagonal dominance may give a reasonable numerical simulation (Farooq et al, 1993;Chou and Huang, 1994;Yang et al, 1997;Mendes et al, 2001;Webley and He, 2000). However, it should be noted that even for PSA processes with weak adsorbates the coexistence of convection and adsorption in a di erence equation provides matrices with stronger diagonal dominance and consequently leads to more robust algorithms.…”