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)
