def test_tplot(self): Qads = np.array([ 116.296, 116.486, 116.634, 116.776, 116.895, 117.021, 117.126, 117.252, ]) Prel = np.array([ 0.301205, 0.351224, 0.400769, 0.450922, 0.50095, 0.551038, 0.601054, 0.65121, ]) Tcheck = np.array([ 5.02005, 5.35198, 5.69661, 6.06841, 6.46997, 6.91209, 7.40569, 7.96934, ]) r = tplot.tplot(Prel, Qads, thickness.HarkinsJura(), tmin=5.0, tmax=8.0, dcf=0.0012800, sacf=1.000, sa=399.8633) np.testing.assert_almost_equal(r.ext_sa, 4.0291, 3) np.testing.assert_almost_equal(r.mv, 0.146955, 4) np.testing.assert_almost_equal(r.ma, 395.8343, 2) for i in range(len(Tcheck)): np.testing.assert_almost_equal(Tcheck[i], r.t_fit[i], 4)
def test_tplot(self): Qads = np.array([116.296, 116.486, 116.634, 116.776, 116.895, 117.021, 117.126, 117.252, ]) Prel = np.array([ 0.301205, 0.351224, 0.400769, 0.450922, 0.50095, 0.551038, 0.601054, 0.65121, ]) Tcheck = np.array([5.02005, 5.35198, 5.69661, 6.06841, 6.46997, 6.91209, 7.40569, 7.96934, ]) r = tplot.tplot( Qads, Prel, thickness.HarkinsJura(), dcf=0.0012800, sacf=1.000, sa=399.8633 ) np.testing.assert_almost_equal( r.ext_sa, 4.0291, 3 ) np.testing.assert_almost_equal( r.mv, 0.146955, 4 ) np.testing.assert_almost_equal( r.ma, 395.8343, 2 ) for i in range(len(Tcheck)) : np.testing.assert_almost_equal( Tcheck[i], r.t[i], 4 )
from micromeritics import tplot, bet, thickness, util, plots from micromeritics import isotherm_examples as ex # Setup the raw data for the calculation. # using example isotherm data from Silica Alumina analyzed with N2. s = ex.silica_alumina() # Compute the BET surface area for this sample # using the range 0.05 to 0.30, and the Cross-sectional area # for Nitrogen 0.162 BET_calc = bet.bet(s.Prel, s.Qads, Pmin = 0.05, Pmax = 0.3 , csa = 0.162) # Do the tplot calculation. TP_calc = tplot.tplot(s.Prel, s.Qads, thickness.HarkinsJura(), tmin = 3.0, tmax = 5.0, dcf = 0.0015468, sacf = 1.0, sa = BET_calc.sa ) plots.plotTPlot( TP_calc ) plots.show() # Show the results of the tplot calculation. print("Micropore Volume: %.4f cm^3/g" % TP_calc.mv) print("Micropore Area: %.4f m^2/g" % TP_calc.ma) print("External surface area: %.4f" % TP_calc.ext_sa) print("Total surface area (BET): %.4f" % TP_calc.tot_sa)
from micromeritics import tplot, bet, thickness, util, plots from micromeritics import isotherm_examples as ex # Setup the raw data for the calculation. # using example isotherm data from Silica Alumina analyzed with N2. s = ex.silica_alumina() # Compute the BET surface area for this sample # using the range 0.05 to 0.30, and the Cross-sectional area # for Nitrogen 0.162 BET_calc = bet.bet(s.Prel, s.Qads, Pmin=0.05, Pmax=0.3, csa=0.162) # Do the tplot calculation. TP_calc = tplot.tplot(s.Prel, s.Qads, thickness.HarkinsJura(), tmin=3.0, tmax=5.0, dcf=0.0015468, sacf=1.0, sa=BET_calc.sa) plots.plotTPlot(TP_calc) plots.show() # Show the results of the tplot calculation. print("Micropore Volume: %.4f cm^3/g" % TP_calc.mv) print("Micropore Area: %.4f m^2/g" % TP_calc.ma) print("External surface area: %.4f" % TP_calc.ext_sa) print("Total surface area (BET): %.4f" % TP_calc.tot_sa)