A Dual Pressure-Temperature Interpolation Model for Mathematically Exact Adsorption Equilibrium for Performance Simulation of Integrated Absorption-Adsorption Refrigeration System using the MIL-101(Cr)/water Pair
Description
This Data is intended to reproduce the equilibrium plots for silica-gel/water and MIL-101(Cr)/water illustrated in Figure 4 in the manuscript based on the proposed higher-order pressure-temperature interpolation model.
Files
Steps to reproduce
MIL-101/water Isotherm Plotting In order to plot the equilibrium water adsorption isotherms for MIL-101(Cr)/water based on the pressure-temperature interpolation model presented in the study, the following stepwise procedure should be followed: 1. Unzip the folder titled "ISOTHERM_FITTING_DATA_FOR_MIL-101-WATER" 2. Copy the MATLAB subroutine titled "equilibirum_uptake_ads_MIL_101.m" along with the following uploaded files in a single folder: P_P0_ads_70.dat, X_ads_70.dat, coeffs Adsorption Isotherm uptake v.s. relative pressure of MIL-101.dat, diff_X_25_X_70_ads.dat, P_P0_ads_25.dat, P_P0_des_25.dat, X_ads_25.dat 3. Open this folder in MATLAB address bar and type the following commands in command window: x_adsorp_70 = load('P_P0_ads_70.dat'); y = load('X_ads_70.dat'); cs_ads = spline(x_adsorp_70,y); Coeffs_ads = load('Coeffs_ads.dat'); P_P0 = linspace(0.018790850000000,0.953393295000000,10000); x = P_P0; y = load('diff_X_25_X_70_ads.dat'); cs_diff_ads = spline(x,y); x_adsorp = load('P_P0_ads_25.dat'); x_desorp = load('P_P0_des_25.dat'); y = load('X_ads_25.dat'); cs_ads_25 = spline(x_adsorp,y); P_P0_ads = x_adsorp(1); P_P0_des = x_desorp(length(x_desorp)); xx_ads = linspace(P_P0_ads,P_P0_des,10000); P_P0 = 0.05:0.0001:0.9; x_ads_model_prediction = zeros(length(P_P0),1); for i = 1:length(x_ads_model_prediction) x_ads_model_prediction(i) = equilibirum_uptake_ads_MIL_101(T_ads,cs_ads,cs_ads_25,cs_diff_ads,Coeffs_ads,P_P0(i),xx_ads); end plot(P_P0,x_ads_model_prediction) NOTE: 1. T_ads is the temperature in Kelvin at which the adsorption isotherm needs to be plotted 2. The equilibrium desorption isotherms for MIL-101(Cr)/water can be plotted in exactly the same way except the following replacements: P_P0_ads_70.dat by P_P0_des_70.dat, X_ads_70.dat by X_des_70.dat, diff_X_25_X_70_ads.dat by diff_X_25_X_70_des.dat, X_ads_25.dat by X_des_25.dat, coeffs Adsorption Isotherm uptake v.s. relative pressure of MIL-101.dat by coeffs Desorption Isotherm uptake v.s. relative pressure of MIL-101.dat Silica-gel/water Isotherm Plotting In order to plot the equilibrium water adsorption isotherms for silica-gel based on the spline-interpolation model presented in the study, the following stepwise procedure should be followed: 1. Unzip the folder titled "ISOTHERM_FITTING_DATA_FOR_SILICA_GEL_WATER" 2. Copy the MATLAB subroutine titled "equilibirum_uptake_ads_silica_gel.m" along with the following uploaded files in a single folder: Coeffs.dat, X_60.dat 3. Open this folder in MATLAB address bar and type the following commands in command window: Coeffs_ds = load('Coeffs.dat'); y = load('X_60.dat'); P_P0 = linspace(0.065917550328049,0.926910768617581,10000); x_ds = P_P0; cs_ds = spline(P_P0,y); x_ads_model_prediction = zeros(length(P_P0),1); for i = 1:length(x_ads_model_prediction) x_ads_model_prediction(i) = equilibirum_uptake_ads_silica_gel(T_ads,x_ds,P_P0(i),Coeffs_ds,cs_ds); end plot(P_P0,x_ads_model_prediction)