On the Numerical Simulation of Rapid Pressure Swing Adsorption for Air Separation

Abstract

Model penjerapan buai tekanan deras yang konsisten (RPSA) secara fizikal telah dibangunkan oleh Choong (2000). Model ini boleh diselesaikan secara berangka dengan menjelmakan persamaan kebezaan biasa (ODE). Seterusnya, sistem ODE dikamirkan terhadap masa. Dua kaedah pendiskretan ruang dipertimbangkan, iaitu kaedah penempatan bersama ortogon (OC) dan penempatan bersama ortogon terhadap unsur terhingga (OCFE). Program penyelakuan RPSA disahkan menggunakan model proses supaya dapat dibandingkan keputusan berangka dengan penyelesaian keserupaan yang diperoleh oleh Scott (1991), bagi penekanan dan penyahtekanan udara terhadap lapisan penjerapan separuh tak terhingga. Ralat keabadian jisim digunakan sebagai pengukur kepada kejituan pengiraan berangka. Bagi menjamin hasil pengiraan keabadian jisim yang jitu, semua kamiran masa dijelmakan kepada ODE. Dengan jelmaan ini, kejituan pengiraan berangka hanya bergantung kepada pendiskretan ruang dan had terima yang digunakan dalam algoritma pengamiran ODE. Kesan pendiskretan ruang dan nilai had terima yang digunakan dalam penyelesaian ODE (TOL) terhadap ketepatan keputusan berangka dan masa pemprosesan komputer dikaji. Suatu sekaitan untuk mengganggar TOL dicadangkan. Walaupun kaedah OC memerlukan masa pengkomputeran yang lebih panjang, tetapi keputusan yang dihasilkan lebih jitu berbanding kaedah OCFE. Kaedah OC sesuai digunakan untuk projek masa depan kami seperti membangunkan ciri algoritma bagi proses berkitar. Kata kunci: Penjerapan buai tekanan deras; pemodelan; penyelakuan berangka; pemisahan udara The physically consistent rapid pressure swing adsorption (RPSA) models developed by Choong are solved numerically by spatially discretisating the partial differential equations (PDEs) to a system of ordinary differential equations (ODEs), which are then integrated over time. Two spatial discretisation methods are considered, i.e. the methods of orthogonal collocation (OC) and orthogonal collocation on finite elements (OCFE). The RPSA simulation programs are validated by simplifying the process models to compare the numerical results with similarity solutions obtained by Scott for air pressurisation and depressurisation into a semi–infinite adsorbing bed. The error in conservation of mass is computed as a guide to the numerical accuracy of the calculations. To ensure good accuracy in the calculation of conservation of mass, we transform all the time integrals into ODEs. With the transformation, the accuracy of the numerical calculations depends only on the spatial discretisation and the tolerance used in the ODE integration algorithm. The effect of the spatial discretisation and the value of TOL on the accuracy of the numerical results and the computer processing time is studied. A rule of thumb to estimate the value of TOL is proposed. The method of OC is found to give accuracy comparable to that of the method of OCFE, albeit requiring more computing time. We consider the method of OC sufficient for the purpose of our future work, i.e. to develop novel algorithm features for cyclic processes. Key words: Rapid pressure swing adsorption; modelling; numerical simulation; air separation