Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
METHOD AND SYSTEM FOR QUANTIFYING POROUS MATERIAL
Document Type and Number:
WIPO Patent Application WO/2024/036164
Kind Code:
A2
Abstract:
The processes, systems, and methods for quantifying porous material are disclosed. A porous solid adsorbs an adsorbate. The saturated porous solid is heated and an isotherm may be generated such as from measured adsorbate mass loss. The type of isotherm is identified, and subsequent analysis may be performed to calculate specific surface area, pore size distribution, and pore volume. The specific type of subsequent analysis, such as derivation of a BET transform and application of BET theory, depends on the type of isotherm. This allows the porous material to by quantified using a variety of different adsorbates.

Inventors:
STAHLFELD KURT W (US)
BELMONT ERICA L (US)
Application Number:
PCT/US2023/071861
Publication Date:
February 15, 2024
Filing Date:
August 08, 2023
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
UNIV WYOMING (US)
International Classes:
G06T7/40
Attorney, Agent or Firm:
COLEMAN, Kyle S. et al. (US)
Download PDF:
Claims:
Atty Dkt #: UWYO-P0012PCT What is claimed is: 1. A method for quantifying porous material, the method comprising: placing a porous solid in a saturator; adsorbing an adsorbate by the porous solid to form a saturated porous solid; heating the saturated porous solid; generating an isotherm; determining if a BET transform can be calculated utilizing the isotherm; calculating a surface area and BET constant if the BET transform is calculated; and calculating a pore size distribution and surface area if the BET transform is not calculated. 2. The method of claim 1, wherein the heating comprises heating the saturated porous solid to a first temperature until a boiling point of the adsorbate is reached and then heating the saturated porous solid to a second temperature. 3. The method of claim 1, wherein the adsorbate is a liquid. 4. The method of claim 1, wherein the adsorbate is a dissolved solid. 5. The method of claim 1, wherein saturation of the porous solid includes a multilayer formation of the adsorbate on the porous solid. 6. The method of claim 1, further comprising: calculating a surface area of the saturated porous solid. 7. The method of claim 1, further comprising application of Density Functional Theory (DFT) or Dubinin-Radushkevich theory to the calculated isotherm. Atty Dkt #: UWYO-P0012PCT 8. A method for quantifying porous material, the method comprising: placing a porous material in a sealed saturator; saturating the porous material with an adsorbate in the sealed saturator; placing the saturated porous material in a thermogravimetric analysis (TGA) machine and in a controlled chemical composition environment, such as within a sealed pan with a pinhole lid; heating the saturated porous material in the TGA machine, wherein the heating comprises heating the saturated porous material to a first temperature and then heating the saturated porous material to a second temperature; generating an isotherm from the mass loss measured from the heated saturated porous material; and determining if a BET transform may be calculated utilizing the isotherm. 9. The method of claim 8, further comprising: calculating a BET transform; calculating a BET constant; and calculating a specific surface area. 10. The method of claim 8 further comprising: calculating a Kelvin radius distribution; calculating pore size and pore volume distributions; and calculating a specific surface area. 11. The method of claim 8, wherein the isotherm is derived for a liquid adsorbate. 12. The method of claim 8, wherein the isotherm is derived for a dissolved solid adsorbate. Atty Dkt #: UWYO-P0012PCT 13. The method of claim 8, wherein saturation of the porous solid is determined to include a multilayer formation of the adsorbate on the porous solid. 14. The method of claim 8, wherein saturation of the porous solid is determined to include a monolayer formation of the adsorbate on the porous solid. 15. A system for quantifying porous material, the system comprising: a porous solid; an adsorbate, wherein the adsorbate comprises at least one of a gas, a liquid, or a dissolved solid; a saturator, wherein the saturator is configured to saturate the porous solid with the adsorbate; and a thermal analyzer machine configured to heat the saturated porous solid and generate an isotherm; wherein the isotherm determines if a BET transform of the porous solid can be calculated. 16. The system of claim 15, wherein the thermal analyzer machine is configured to heat the saturated porous material at a first temperature until a boiling point of the adsorbate is reached. 17. The system of claim 15, wherein a BET constant of the porous solid is calculated. 18. The system of claim 15, wherein a pore size distribution is calculated utilizing at least the thermal analyzer machine. 19. The system of claim 15, wherein a specific surface area of the porous solid is calculated. 20. The system of claim 15, where the TGA machine maintains an atmosphere surrounding the porous solid containing a known partial pressure of adsorbate.
Description:
Atty Dkt #: UWYO-P0012PCT TITLE: METHOD AND SYSTEM FOR QUANTIFYING POROUS MATERIAL CROSS-REFERENCE TO RELATED APPLICATIONS [0001] This application claims priority under 35 U.S.C. § 119 to provisional application Serial No. 63/396,064 filed on August 8, 2022, which is incorporated by reference in its entirety. FIELD OF THE DISCLOSURE [0002] The present disclosure relates to porous material. More particularly, but not exclusively, the present disclosure relates to a method and system for quantifying porous material. BACKGROUND [0003] Porous solids with nanometer-sized pores and large specific surface areas are a highly important class of materials. Uses of such porous solids include filtration, battery electrodes, catalyst supports, and carbon sequestration. These porous solids are characterized by their specific surface areas, typically > 100 m 2 /g, and pore size distributions. These parameters are typically measured using cryogenic adsorption of nitrogen, argon, or other adsorbates, a physisorption technique often paired with Brunauer- Emmett-Teller (BET) analysis. This method quantifies how a particular solid will interact with cryogenic adsorbates, but this is often a poor predictor of how that solid will interact with other adsorbates, greatly limiting the applicability of these results. Nevertheless, this method remains the standard technique due to lack of other options. Therefore, what is needed is a system and method for quantifying porous material to determine how effective the material will be for various uses. SUMMARY [0004] Therefore, it is a primary object, feature, or advantage of the present disclosure to improve over the state of the art. Atty Dkt #: UWYO-P0012PCT [0005] It is a further object, feature, or advantage of the present disclosure to quantify porous solids using various adsorbates. [0006] It is a still further object, feature, or advantage of the present disclosure to use elevated temperatures to generate isotherms. [0007] Another object, feature, or advantage is to use elevated temperatures to quantify a porous solid. [0008] Yet another object, feature, or advantage is to determine if a BET transform can be calculated for a porous solid. [0009] In one aspect of the present disclosure, a method for quantifying porous material or porous solids is disclosed. The method may include placing a porous solid in a saturator. The method may also include adsorbing an adsorbate by the porous solid to form a saturated porous solid. The method may further include heating the saturated porous solid to desorb the adsorbate. The method may also include generating an isotherm such as from measured adsorbate mass loss. The method may further include determining if a BET transform can be calculated utilizing the isotherm. The method may also include calculating a BET constant if the BET transform is calculated and calculating a surface area and pore size distribution even if the BET transform is not calculated. [0010] In another aspect of the present disclosure, another method for quantifying porous solids is disclosed. The method may include placing a porous material in a sealed saturator. The method may also include saturating the porous material with an adsorbate in the sealed saturator. The method may include placing the saturated porous material in a thermogravimetric analysis (TGA) machine. The method may further include heating the saturated porous material in the TGA machine. The heating step may include heating the saturated porous material to a first temperature and then heating the saturated porous material to a second temperature. The method may further include generating an isotherm such as from mass loss measured from the heated saturated porous material and determining if a BET transform may be calculated utilizing the isotherm. [0011] In another aspect of the present disclosure a system for quantifying porous material is disclosed. The system may include a porous solid. The system may also include an adsorbate. The adsorbate may include at least one of a gas, a liquid, or a dissolved solid. The system may also include a saturator. The saturator may be configured to saturate the Atty Dkt #: UWYO-P0012PCT porous solid with the adsorbate. The system may also include a thermal analyzer machine configured to heat the saturated porous solid and measure mass loss data, from which an isotherm may be generated. The isotherm may determine if a BET transform of the porous solid can be calculated. [0012] One or more of these and/or other objects, features, or advantages of the present disclosure will become apparent from the specification and claims that follow. No single embodiment need provide each and every object, feature, or advantage. Different embodiments may have different objects, features, or advantages. Therefore, the present disclosure is not to be limited to or by any objects, features, or advantages stated herein. BRIEF DESCRIPTION OF THE DRAWINGS [0013] Illustrated embodiments of the disclosure are described in detail below with reference to the attached drawing figures, which are incorporated by reference herein. [0014] FIG.1 is an illustration of the system for quantifying porous material in accordance with one aspect of the present disclosure. [0015] FIG.2 is a chart illustrating isotherms of nitrogen, argon, water, and toluene on a 13x molecular sieve (mol sieve) in accordance with one aspect of the present disclosure. [0016] FIG. 3 is a chart illustrating isotherms of nitrogen, argon, water, and toluene on activated carbon in accordance with one aspect of the present disclosure. [0017] FIG.4 is a chart illustrating BET transforms and linear fits of toluene on activated carbon in accordance with one aspect of the present disclosure. [0018] FIG.5 is a chart illustrating BET transforms and linear fits of water on a 13x mol sieve in accordance with one aspect of the present disclosure. [0019] FIG.6 is a chart illustrating BET transforms and linear fits of toluene on a 13x mol sieve in accordance with one aspect of the present disclosure. [0020] FIG.7 is a chart illustrating isotherms of water, toluene, and nitrogen on silica gel and MCM-41 in accordance with one aspect of the present disclosure. [0021] FIG. 8 is a chart illustrating Kelvin radii of several common fluids at 1 atm in accordance with one aspect of the present disclosure. Atty Dkt #: UWYO-P0012PCT [0022] FIG. 9 is a chart illustrating pore size distributions of silica gel determined from nitrogen, water, and toluene desorption in accordance with one aspect of the present disclosure. [0023] FIG.10 is a chart illustrating pore size distributions of MCM-41 determined from nitrogen, water, and toluene desorption in accordance with one aspect of the present disclosure. [0024] FIG. 11 is a chart illustrating the thermogravimetric desorption method for calculation of porous solid SSA in accordance with one aspect of the present disclosure. [0025] FIG. 12 is a flow chart illustrating a method for quantifying porous material in accordance with one aspect of the present disclosure. [0026] FIG.13 is a flow chart illustrating another method for quantifying porous material in accordance with one aspect of the present disclosure. BRIEF DESCRIPTION OF THE TABLES [0027] Illustrated aspects of the disclosure are described in detail below with reference to attached Tables, which are incorporated by reference herein and where: [0028] Table 1 provides specific surface areas derived from isotherms of nitrogen, argon, water, and toluene on a molecular sieve and activated carbon in accordance with one aspect of the present disclosure. [0029] Table 2 provides BET constants for a molecular sieve and activated carbon with nitrogen, argon, water, and toluene adsorbates in accordance with one aspect of the present disclosure. [0030] Table 3 provides the pore volumes and specific surface areas for silica gel derived from isotherms of water, toluene, and nitrogen in accordance with one aspect of the present disclosure. [0031] Table 4 provides the specific surface areas and pore volumes for MCM-41 derived from isotherms of water, toluene, and nitrogen in accordance with one aspect of the present disclosure. [0032] Table 5 provides the temperature methods used for TGA desorption and analysis methods applied to derived isotherms, categorized by IUPAC isotherm Type. [0033] Table 6 provides the contact angles measured and used for Kelvin analysis. Atty Dkt #: UWYO-P0012PCT DETAILED DESCRIPTION [0034] Porous solids represent a wide range of highly important materials with diverse uses including filtration, catalyst support, drug delivery, carbon removal, hydrogen storage, and electrochemical applications such as battery electrodes. The adsorption of water by porous solids has broad importance in chemical engineering, catalysis, and material science. The physically or chemically adsorbed water on the surface of porous solids may play an important role in governing interface properties, chemical reaction pathways, or catalytic performance. Protein and surfactant adsorption play a fundamental role in the biomaterials field. Biological media, such as blood or serum, coat the biomaterial porous solid with proteins. This protein layer mediates the interaction between the biomaterial and cells. Surfactant adsorption is similar, utilizing surfactant molecules in the place of proteins. [0035] The pores within porous solids can range in size over several orders of magnitude from micropores (<2 nm) to mesopores (2-50 nm) to macropores (>50 nm). Accurate analysis of the specific surface areas (SSAs) and size distribution of pores for porous solids is critical for understanding these materials and their applications. Direct measurement of SSAs and pore sizes is generally not possible due to the small pore size and three- dimensionality of porous solids. Instead, physisorption is the most common method employed for the characterization of these porous solids. Physisorption typically proceeds as follows: a sample of the porous solid material is placed under vacuum and cooled to cryogenic temperature. A small amount of an adsorptive gas is introduced to the solid and the pressure increase is measured, allowing for the amount of gas adsorbed onto the solid to be calculated. Additional gas is incrementally added, with the amount of adsorbed gas measured at each step. This produces a plot of amount of gas adsorbed versus the relative pressure, typically called an isotherm. [0036] Various theories can be applied to an isotherm to calculate the specific surface area and pore size distribution of the solid. The most famous and widely used of these is BET theory, named after the authors who developed it in 1938. BET theory applies a mathematical model to the formation of a multilayer of gas molecules adsorbed onto the surface of a porous solid and calculates a surface area. Pore size distributions can also be Atty Dkt #: UWYO-P0012PCT calculated from physisorption isotherms, commonly done using the Kelvin equation-based Barrett-Joyner-Halenda (BJH) method developed by Barrett et al. in 1951 or newer methods such as non-local density functional theory. [0037] Despite its widespread use, cryogenic physisorption has significant limitations. Physisorption is most commonly conducted using nitrogen at -196 °C and in a pressure range of 10 -7 to 1 atm. The rate of diffusion of adsorptive gasses is slow at these conditions, making it difficult to accurately measure equilibrium isotherms. These kinetic limitations are especially pronounced for microporous materials with pore diameters < 2 nm, and isotherms produced from these materials are frequently of limited value. The choice of adsorptive molecule is also limiting. As historical context, Brauner et al. logically chose nitrogen (N2) as the adsorptive gas for their famous study of iron catalysts for nitrogen fixation, and nitrogen has remained as the most commonly used adsorptive. However, nitrogen has specific interactions with certain surface functional groups and exposed ions due to its quadrupolar nature, which shifts the isotherm to the low-pressure regime and causes both experimental difficulties due to extremely low pressure and analysis issues. [0038] Use of argon at -186 °C is now recommended by the International Union of Pure and Applied Chemistry (IUPAC) instead of nitrogen because argon does not exhibit specific interactions with surface functional groups, thereby overcoming some of the issues incurred with nitrogen, but argon is still used less frequently than nitrogen due to cost. Additionally, argon still incurs the kinetic limitations noted above for nitrogen due to the required cryogenic temperatures. Physisorption is also performed with other gases but, in all cases, results only reflect how that particular gas interacts with a porous material, and this behavior may or may not be applicable to other chemicals of relevance to that material. [0039] The processes, systems, and methods of the present disclosure show that surface areas and pore size distributions of porous solids can be determined by thermogravimetric desorption using only a thermogravimetric analyzer (TGA), and in a high throughput manner. The processes, systems, and methods of the present disclosure show that BET theory can be applied to thermal desorption, and BET constants, which are non- dimensional measures of adsorbate-adsorbent interaction strength, can be calculated for different adsorbate-adsorbent systems. Additionally, the processes, systems, and methods of the present disclosure show that multilayer formation does not occur for an important Atty Dkt #: UWYO-P0012PCT class of porous solids, desiccants, with some adsorbates, in contrast to results suggested by nitrogen physisorption; instead, adsorption occurs exclusively or predominantly by capillary condensation. For these materials, the Kelvin equation can be applied to determine pore size distributions and surface areas. [0040] A novel approach for deriving isotherms using thermogravimetric desorption of adsorbates from porous solids is disclosed. For systems which exhibit multilayer formation, BET theory can then be applied to determine surface areas from the isotherms along with BET constants. For systems where multilayer formation does not occur, BET theory is not applicable and the Kelvin equation can be applied to generate pore size distributions along with surface areas. The novel method disclosed may use elevated temperatures of approximately 100-500 °C, thereby avoiding kinetic limitations associated with cryogenic physisorption. [0041] This method can be conducted with many different adsorbates, or substances that are deposited on the surface of another substance, including those of particular relevance to applications of the porous solids. An adsorbate may include any gas, liquid, or dissolved solid. The atoms, ions, or molecules from the adsorbate adhere to a surface of the adsorbent, creating a film. This unique aspect of the method allows the specific interaction between a porous solid, or adsorbent, and an adsorbate of interest to be analyzed. The novel approach may be demonstrated with water and toluene, or other adsorbates on several porous solids or other adsorbents. This method may work with many relevant adsorbate/adsorbent systems. [0042] Adsorbents may include oxygen containing compounds, carbon-based compounds, metal organic framework compounds, or polymer-based compounds. The adsorbents may be a powder or in the form of spherical pellets, rods, moldings or monoliths. The adsorbent may have a hydrodynamic radius. The adsorbent may have high abrasion resistance, high thermal stability or small pore diameters to allow for a higher exposed surface area or a higher capacity for adsorption. Oxygen containing compounds may be typically hydrophilic or polar. The oxygen-containing compounds may include silica gel, limestone or calcium carbonate, or zeolites. The carbon-based compounds may be hydrophobic or non-polar and may include activated carbon or graphite. The polymer-based compounds may be polar or non-polar depending on the functional groups in the polymer matrix. Atty Dkt #: UWYO-P0012PCT Zeolites and metal-organic framework compounds are adsorbents typically selected for carbon capture and storage. [0043] Fig. 1 illustrates a summary of the thermogravimetric desorption method for calculation of porous solid specific surface area, as well as either a BET constant or a pore size distribution depending on the isotherm type for the tested adsorbent/adsorbate system in accordance with one aspect of the present disclosure. A process schematic 10 for thermogravimetric determination of isotherms, and subsequent analysis of those isotherms to derive specific surface area and either a BET constant or a pore size distribution depending on the appropriate analysis for the tested adsorbent/adsorbate system in accordance is shown. Porous solids 12, having a multitude of pores 14, are first saturated with the adsorbate 16 in a saturator 18. The saturator 18 may be a bench top saturator. The saturator 18 may include a porous solid holder 28, such as a glass jar, containing the porous solid 12 and a test tube or an adsorbate holder 20 filled with the liquid adsorbate 16 to saturate the porous solid 12. The saturated porous solid 22 may be loaded into a sample pan 24. The sample pan or the TGA may have a heating element 26, such as a TGA furnace. The loaded sample pan 24 may be heated while the sample mass loss is measured using a thermogravimetric analyzer (TGA). The volume of liquid adsorbed onto the solid as a function of temperature can be calculated from the TGA mass-temperature trace. This volume-temperature plot may be mapped into a physisorption isotherm. The physisorption isotherm may be defined as adsorbed volume plotted against relative pressure (P/P S ) where P is the partial pressure of the adsorptive gas and PS is the saturation vapor pressure of the adsorptive gas at a given temperature. FIG. 11 provides another summary of the thermogravimetric desorption method for calculation of porous solid specific surface area. [0044] In one aspect of the present disclosure, P is equal to ambient pressure (e.g., 1 atm at sea level and 0.78 atm in Laramie, WY) and PS is a function of temperature which can be readily calculated from tabulated thermodynamics data. The produced isotherm may then be analyzed to determine its Type classification. Types I, II and IV isotherms show large increases in volume adsorbed as P/PS increases from 0 to ~0.2, while Types III and V show minimal, or no volume adsorbed in this region. For Type I, II and IV isotherms, BET theory may be applicable and the BET transform can be applied to calculate a porous solid specific surface area and BET constant. For Type III and V isotherms, the adsorbent- Atty Dkt #: UWYO-P0012PCT adsorbate interaction may be too weak for multilayer formation and BET theory is not valid. For these isotherms, the Kelvin equation can be applied to calculate pore size distribution and specific surface area. [0045] Isotherms derived from thermogravimetric desorption of toluene and water adsorbed on 13x molecular sieve (mol sieve) are shown in Fig. 2. Nitrogen and argon physisorption were conducted using conventional cryogenic techniques and are also shown in Fig.2. Toluene and water isotherms on the mol sieve show Type I behavior with large increases in volume adsorbed from P/PS of 0 to 0.2 followed by a relatively constant slope to P/P S of unity. The nitrogen and argon isotherms show Type II behavior with a sharp increase in volume adsorbed occurring at P/P S less than 0.005, followed by a plateau to P/PS of 0.8 then a slight increase at P/PS greater than 0.9. Slight hysteresis is observed for nitrogen and argon at high P/PS, indicating pore filling by capillary condensation. In contrast, no hysteresis loop can be observed in the thermogravimetric desorption-derived curves as this method only provides insights into the desorption branch of the isotherm. Thermodynamic equilibrium is established on the desorption branch. Water, toluene, and argon all show a maximum adsorbed volume or available pore volume of approximately 0.27 mL/g, while nitrogen shows an available pore volume of 0.35 mL/g, demonstrating that accessible pore volume is dependent on the choice of adsorbate. [0046] Isotherms derived from thermogravimetric desorption of water and toluene and cryogenic physisorption of nitrogen and argon on activated carbon are shown in Fig. 3. Toluene exhibits the same Type I behavior as seen on mol sieve. In contrast, the isotherm of water adsorbed on activated carbon shows essentially no volume adsorbed at P/P S < 0.4 and is best described as a modified Type V isotherm. This behavior can be attributed to weak adsorbent-adsorbate interactions and is typical of water adsorption on hydrophobic microporous adsorbents. Water and toluene isotherms on activated carbon both show available pore volumes of approximately 0.30 mL/g, argon shows 0.44 mL/g, and nitrogen shows 0.62 mL/g. The adsorption characteristics and available pore volumes of activated carbon are clearly dependent on adsorbate identity. Considering that adsorption of aromatic chemicals is a major application of activated carbon, the measurement of pore volume available to relevant chemical species is highly important. The results show that nitrogen Atty Dkt #: UWYO-P0012PCT physisorption overestimates the adsorption capacity of activated carbon for toluene by 106%, while argon overestimates by 47%. BET Analysis [0047] BET transforms are calculated from isotherms using the standard BET method and BET transforms and linear fits for toluene adsorbed on activated carbon, water adsorbed on mol sieve, and toluene adsorbed on mol sieve are shown in Figs.4-6, respectively. Fig. 4 illustrates BET transforms and linear fits of toluene on activated carbon. Fig.5 illustrates BET transforms and linear fits of water on the mol sieve. Fig.6 illustrates BET transforms and linear fits of toluene on the mol sieve. BET transforms show highly linear behavior with best fit lines demonstrating R 2 > 0.99. Porous solid specific surface areas (SSAs) are readily calculated using the slope and intercept of the fitted line, and results are summarized in Table 1 below, and uncertainties reflect standard deviations. Table 1: Specific surface area of isotherms Specific Surface Area (m 2 /g) Nitrogen Argon Water Mol Sieve 5997 5670 5956 ± 27 5651 ± [0048] For mol sieve porous solids, the thermogravimetrically determined water is 595.6 ± 2.7 m 2 /g and toluene BET SSA is 565.1 ± 21.0 m 2 /g. For comparison, cryogenic physisorption gives a nitrogen BET SSA of 599.7 m 2 /g and an argon SSA of 567.0 m 2 /g for mol sieve. All four SSA measurements show good agreement within a range of 35 m 2 /g, or 6% of the mean value. This agreement of the BET SSAs is particularly impressive considering that the corresponding isotherms reveal significantly different adsorbed volumes in the relevant P/P S range of 0.05 to 0.25, as shown in FIG.2. SSAs of activated carbon show much more variance than mol sieve SSAs with nitrogen and argon BET surface areas of 802.0 m 2 /g and 642.3 m 2 /g respectively. Unlike mol sieve, the BET SSAs of activated carbon are clearly dependent upon the choice of adsorbent; toluene on activated carbon demonstrates a surface area of 496.1 ± 5.8 m 2 /g. BET theory is not applicable to water adsorbed on activated carbon as water does not form a multilayer on activated carbon, as shown in FIG.3, so no BET SSA can be calculated for this system. Atty Dkt #: UWYO-P0012PCT [0049] BET SSAs are highly consistent across four different adsorbents for the highly ordered mol sieve, suggesting that these surfaces area could reasonably be interpreted as geometric surface areas. The BET SSAs for activated carbon show much greater variation, indicating that SSAs for activated carbon are instead dependent on the properties of the adsorbent molecule. BET SSAs of activated carbons have been shown to be heavily dependent on the choice of adsorbate molecule. [0050] In addition to the widely used SSA, BET analysis also produces a BET constant. BET constants are generally not reported because they are specific to the adsorbate used, typically nitrogen, which is seldom the chemical of interest to the system being studied. The BET constant is a non-dimensional measure of the strength of the interaction between adsorbate and the adsorbent. BET constants for water and toluene adsorbed on mol sieve are 24.7 ± 1.2 and 10.6 ± 0.9, respectively, as shown in Table 2, and uncertainties reflect standard deviations. The greater value for water may indicate that mol sieve has a higher affinity for water than toluene. Table 2. BET Constants of mol sieve and activated carbon. BET Constants (unitless) [0051] Toluene on activated carbon has a BET constant of 52.0 ± 3.7, demonstrating that activated carbon has a high affinity for toluene. This finding is consistent with and supports the use of activated carbon for toluene adsorption. Water exhibits only weak interactions with the relatively hydrophobic activated carbon, with very little volume adsorbed at P/PS < 0.3, indicating that there is no meaningful multilayer formation within this range and BET theory cannot be applied. The nitrogen BET constant of nitrogen adsorbed on mol sieve is negative, which occurs somewhat frequently for cryogenic physisorption. Negative BET constants violate the assumptions of BET theory and are physically meaningless. BET constants for argon adsorbed on mol sieve and activated carbon are 1699.6 and 950.3 Atty Dkt #: UWYO-P0012PCT respectively. These large values, and the large BET constant of 7,715 for nitrogen on activated carbon, are due to small y-intercept values of the BET linear fits (~4·10 -6 g/mL). Kelvin Analysis [0052] Adsorbent/adsorbate systems that produce Type III or V isotherms, such as water on activated carbon shown in FIG.3, show no evidence of coherent multilayer formation with essentially no volume adsorbed at P/PS < 0.2. Thus, for these systems, adsorption is governed predominantly by capillary condensation. For wetting fluids, liquid is more energetically stable within small capillaries and pores as compared with the bulk liquid state. Sufficiently small capillaries can condense gas into liquid well above the saturation temperature, or below the saturation pressure, of the gas. Adsorption of water on activated carbon is difficult to model accurately due to the large contact angle so adsorption on silica gel and Mobil Composition of Matter No. 41 (MCM-41) are examined here instead. Isotherms of water and toluene adsorbed on MCM-41 and silica gel measured using thermogravimetric desorption are shown in Fig.7, along with nitrogen isotherms measured using cryogenic physisorption. The nitrogen isotherms show standard Type IV behavior with significant volume adsorbed at P/P S < 0.05 indicating multilayer formation. In contrast, toluene adsorbed on MCM-41 shows no appreciable volume adsorbed at P/P S < 0.2 while the remaining three water and toluene systems show no appreciable volume adsorbed at P/PS < 0.5. As no coherent multilayer appears for the water and toluene systems, BET is not applicable and the Kelvin equation may be applied. [0053] The Kelvin equation can be reformulated for constant pressure, thermogravimetric desorption, and Kelvin radii for several common chemicals at 1 atm are shown in Fig.8. The common chemicals may include benzene, toluene, butane, hexane, methanol, and water. Large pores are emptied at low temperatures, just above the boiling point of the adsorbate, and progressively smaller pores empty as temperature increases and P/PS decreases. Kelvin radii show the same general shape across the six fluids indicating that the method applied here for water and toluene can work for other common chemicals including alkanes, aromatics, alcohols, and others. [0054] Applying the Kelvin equation to thermogravimetric desorption data allows the pore size distribution of a porous solid to be calculated. Fig.9 shows the pore size distributions of silica gel, determined using the thermogravimetric desorption technique and water and Atty Dkt #: UWYO-P0012PCT toluene adsorbates. Uncertainties show standard deviations. For comparison, a pore size distribution derived from cryogenic nitrogen physisorption is also shown. All three distributions show the peak pore size falling between 6 nm and 7 nm, suggesting that the thermogravimetric desorption method can determine peak pore size using various adsorbates with comparable accuracy to cryogenic nitrogen physisorption. Total pore volumes and specific surface areas (SSAs) for all three adsorbents are shown in Table 3 below, where nitrogen BET SSA is also shown for comparison. Uncertainties show standard deviation. Table 3. Pore size distribution from nitrogen, water, and toluene on silica gel. SSA (m 2 /g) Pore volume (mL/g) [0055] The nitrogen pore volume for silica gel is 45% greater than the water pore volume, and the toluene pore volume is within experimental error of the water. The importance of these results can be illustrated through consideration of the adsorbent/adsorbate systems. For example, as silica gel is commonly used as a desiccant, water adsorption capacity is an important measurement and nitrogen physisorption overestimates the available pore volume of this adsorbent to water. [0056] The BET nitrogen analysis also overestimates the SSA of silica gel compared to water by 75%, and this difference is due to a greater volume of nitrogen adsorbed and significant water adsorption in pores > 10 nm diameter, as larger pores have a lower surface area-to-volume ratio. The pore size distribution of toluene on silica gel has a longer tail at larger pore diameters than those of water and nitrogen, with significant pore volume in pores larger than 20 nm and a correspondingly lower surface area. Despite this behavior, toluene does accurately predict the peak pore size. Atty Dkt #: UWYO-P0012PCT [0057] FIG. 10 illustrates pore size distributions of MCM-41 as determined using the thermogravimetric desorption technique with water and toluene adsorbates, and the Kelvin equation. For comparison, the BET pore size distribution derived from cryogenic nitrogen physisorption is also shown. Results for MCM-41 show similar trends compared to silica gel: nitrogen shows a narrow pore size distribution with the majority of pore volume falling between 2 nm and 4 nm, and a peak pore size of 3.2 nm. The pore size distribution from water shows the same peak pore size of 3.2 nm, although the distribution is broader with significant specific pore volume between 5 nm and 10 nm. The toluene pore size distribution also shows a broader pore size distribution than for nitrogen, but results show good agreement with peak pore size. [0058] Nitrogen pore volume in MCM-41 is 47% greater than water, as shown in Table 4, consistent with other results showing greater adsorption of nitrogen. Uncertainties show standard deviations. Water shows slightly greater pore volume on MCM-41 compared to toluene, but this difference is within experimental error. SSAs of MCM-41 show greater nitrogen values compared to water and toluene, consistent with the trends observed for silica gel. Table 4. Pore size distribution from nitrogen, water, and toluene on MCM-41. SSA (m 2 /g) Pore volume (mL/g) [0059] Nitrogen isotherms on silica gel, FIG.7, and mol sieve, FIG 3, show that both of these common desiccants form nitrogen multilayers. Thermogravimetric desorption using water shows that mol sieve does not form a coherent multilayer with water, but instead adsorbs water predominantly by capillary condensation. As the main purpose of these Atty Dkt #: UWYO-P0012PCT desiccants is to adsorb water the thermogravimetric method accurately captures the behavior of this common system. Methods Experimental [0060] FIGS. 12 and 13 are flow charts illustrating a method for quantifying porous material. First, the porous solid 12, such as activated carbon, 13x molecular sieve 100-120 mesh, MCM-41, and silica gel 230-400 mesh 60 Å, may be heated to ensure desorption of water and other contaminants (Step 200). For example purposes only, the porous solid 12 may be heated to 200 ^C in a vacuum oven for 24 hours. Next, the porous solid 12 is placed or put in a saturator 18 (Step 202). For example, approximately 0.5 g of solid was placed in the bottom of the porous solid holder 28, for example a glass jar of 7.5 cm diameter by 14 cm height may be used. In other aspects of the present disclosure, the porous solid holder 28 may be larger or smaller. The adsorbate 16 is placed in the adsorbate holder 20. For example, approximately 14 mL of liquid adsorbate may be placed in a glass test tube measuring 16 mm outer diameter and 99 mm length. In other aspects of the present disclosure, the adsorbate holder 20 may be larger or smaller. Once the adsorbate holder 20 is filled with adsorbate 16, the adsorbate holder 20 may be placed in the porous solid holder 28. For example, the test tube may be placed in the jar and an air-tight lid is screwed onto the top of the jar. Next, the porous solid adsorbs the adsorbate (Step 204). The complete saturator 18 may be placed on the bench for a minimum time such as 48 hours to allow the porous solid 12 to adsorb vapor, the minimum time may be higher or lower. Thermogravimetric analysis (TGA) [0061] Next, the saturated porous solid 22 is placed on the sample pan 24 (Step 206; Fig. 13). The sample pan 24 may be sealable, such as by a lid with a pinhole. The sample pan may also be part of a thermal analyzer machine. For example, 4-6 mg of sample may be loaded into an aluminum Tzero Differential Scanning Calorimetry (DSC) pan and pan may be sealed with an aluminum Tzero hermetic lid with 75 μm diameter pinhole. The sample pan 24, or sealed DSC pan, may then be placed on a TGA pan and loaded on to the TGA machine. The TGA pan may be platinum. The balance of the TGA may be tared with an empty sample pan 24 and the TGA pan prior to loading the saturated porous solid 22. In Atty Dkt #: UWYO-P0012PCT some aspects of the present disclosure the hermetically sealed sample pan 24 with pinhole lid allows adsorptive gas to desorb from porous solid sample with evolved gas escaping through the pinhole. The pinhole is sufficiently small to prevent diffusion of ambient gas into the sample pan 24, thus maintaining an environment which can be approximated as pure adsorptive gas within the pan. [0062] In one aspect of the present disclosure, the use of a sealed pan with a pinhole lid may be essential, or necessary. The saturated porous solid 22 is loaded into an aluminum pan that is sealed with a lid, such as one with a 75 μm diameter pinhole. The sealed pan may contain adsorbate gas that is desorbed from the saturated porous solid 22. This may maintain a gas environment around the sample of 100% adsorbate gas, which may be highly advantageous. The sealed pan maintains this environment while the pinhole in the lid allows the adsorbate gas to escape so the pressure does not increase while being small enough to prevent TGA flow gas from diffusing into the pan. In some aspects of the present disclosure the pinhole diameter may range from 20 μm to 120 μm or have a larger or smaller diameter. The diameter size used may be dependent on the mass of the saturated porous solid 22 or the porous solid 12. [0063] A TGA may be used to analyze the samples. The TGA may be calibrated using the Curie temperature method with alumel and nickel. The TGA may be calibrated by other methods. The Curie temperatures of alumel and nickel, 154.2 ^C and 357.0 ^C respectively, correspond to the temperature range of mass loss of the water and toluene, along with other common adsorbents, allowing the TGA to be precisely calibrated in the relevant temperature range. [0064] Next the saturated porous solid 22 may be heated (Step 208). The saturated porous solid 22 may be heated to a first temperature over a first period of time or until the saturated porous solid 22 reaches the boiling point of the adsorbate 16 and then heated to a second temperature over a second period of time or until the porous solid 22 reaches a specific temperature. The boiling point of the adsorbate 16 may be defined as the temperature at which the vapor pressure of the adsorbate is equal to the ambient pressure (e.g., 93 ^C for water at 0.78 atm). For example, the saturated porous solid may be heated at 40 ^C/min to the boiling point of the adsorbate and then 2 ^C/min to 550 ^C. Atty Dkt #: UWYO-P0012PCT [0065] For saturated porous solids 22 displaying Type III isotherms, all mass loss may end by approximately 150 ^C; in other aspects of the present disclosure, all mass loss may occur at a higher or lower temperature. The TGA or other heating element 26 may heat the sample or saturated porous solid 22 at specified rate(s). The first specific rate may start from room temperature to some temperature of interest, such as the boiling point of the adsorbate. The second rate may cover a range of interest, such as the adsorbate boiling point to a maximum temperature. The maximum temperature may range from 200-600 ^C or be higher or lower. There may be a region from room temp to the range of interest where no analysis is conducted. During heating in this region, adsorbate gas can diffuse out of the pinhole and nitrogen can diffuse in, creating an environment within the sample pan 24 that is less than 100% adsorbate. To reduce the time in this region, a fast temperature ramp rate may be used to minimize the amount of time the saturated porous solid 22 spends in the temperature range of this region. The temperature ramp rate used may be dependent on the heating element 26 or on what TGA instrument is used. The temperature ramp rate may be based on the choice of instrument used. In some aspects of the present disclosure the temperature ramp rate may be 40 ^C/min. Reasonable range of temperature ramp rates may include 10 ^C/min to 500 ^C/min and could be higher or lower. [0066] In some aspects of the present disclosure, once the temperature reaches the boiling rate of the adsorbate, the temperature ramp rate may need to be slowed down. This may allow for better resolution. The slowed down temperature ramp rate may be either 2 ^C/min or 0.5 ^C/min but in some aspects, may range between 0.2 ^C/min to 10 ^C/min or be higher or lower. [0067] In some aspects of the present disclosure, the vapor pressure of adsorbate gas in proximity must be non-zero and known. In conventional physisorption methods, the sample container is filled only with the adsorbate gas (typically N2) and the absolute pressure is carefully monitored. In at least one of the methods disclosed in the present disclosure the sample or saturated porous solid 22 may be placed in a sealed pan with a pinhole lid. As the temperature rises adsorbate evaporates, filling the pan. The environment within the pan is thus nearly 100% adsorbate gas and at a pressure equal to the ambient (typically 1atm, or 0.78atm in Laramie WY) thus satisfying that the vapor pressure is non- zero and known. Atty Dkt #: UWYO-P0012PCT [0068] Several TGA instruments exist that can maintain a set humidity level (partial pressure of water) within the sample cell. In some aspects of the present disclosure, the TGA sample cell could be set to maintain a constant, nonzero partial pressure of water with the remainder being inert gas, such as N2, Ar, or other inert gases. Porous solid loaded with water adsorbate, or the saturated porous solid 22, could then be loaded into the TGA without using the sealed pan with pinhole lid. The TGA, while maintaining a constant environment containing set water partial pressure level, could be heated in the same fashion as above, such as at a temperature ramp rate of 2 ^C/min to 550 ^C, to desorb the water from the sample. Analysis would be conducted in the same fashion with the partial pressure of the water vapor being controlled. In this aspect, the heating rates do not have to be fine- tuned as much nor do DSC pans with pinhole lids have to be used. Analogous implementation could be achieved for other, non-water adsorbates if appropriate equipment is available, such as a TGA which can maintain a partial pressure of another adsorbate, other than water, in the sample cell. [0069] The saturated porous solids 22 displaying Type III isotherms may be rerun during the heating step. he rerun may use a lower heat rating, such as 0.5 ^C/min, to allow for better temperature resolution over the mass loss range, as shown in Table 5. Table 5. Temperature methods used for TGA desorption and analysis methods applied to derived isotherms, categorized by IUPAC isotherm Type. Isotherm Type Step 1 Step 2 Analysis Method [0070] Gas, such as a purge gas, may flow through the TGA. The gas may include nitrogen, argon, or any other gas. For example, flow rates of nitrogen purge gas may be 40 mL/min balance and 60 mL/min sample, the standard flows for the TGA. All example tests herein were conducted in Laramie, WY, where the atmospheric pressure is 0.78 atm and all calculations were conducted using this pressure. All TGA tests were performed in triplicate Atty Dkt #: UWYO-P0012PCT and mean values reported. Nitrogen physisorption was conducted using a Micromeritics ASAP 2020. Argon physisorption was performed by MSE analytical services. Calculations [0071] Next, the isotherm may be generated or derived (Step 210). In one aspect of the present disclosure, all calculations were performed in MATLAB. The mass of adsorbate adsorbed onto the porous solid per unit mass of solid at a given temperature, ^^ ^^ ^ ^^^, may be calculated from TGA mass loss traces using Eq.1 ^ ^ ^^ ^ ^^^ ൌ ^^^ ^^^െ ^^ ^^ ^ ^ ^^ [1] where ^^^ ^^^ is the sample mass at temperature ^^ as measured by TGA, and ^^ ^^ is the final sample mass as measured by TGA. All mass loss is assumed to be due to desorption when analyzed microporous and mesoporous materials are synthesized at high temperatures and therefore assumed to be stable at the moderate temperatures (≤ 550 ^C) employed in this technique. Adsorbate volume per unit mass of solid at ^^ may be found by dividing ^^ ^^ ^ ^^^ by the temperature-dependent liquid adsorbate density, ^^^ ^^^, Eq. 2. Density values may be obtained from National Institute of Standard and Technology. ^ ^ ^^ ^ ^ ^ ^ ^ ^ ^^ ^ ^^ ^ ^ ^^ ^^^ [2] [0072] Cryogenic pressure as P is the pressure of adsorbate gas and P 0 is the bulk saturation pressure of the adsorbate. Cryogenic physisorption is typically conducted isothermally and at the saturation temperature of the adsorbate corresponding to atmospheric pressure, making P 0 constant and equal to atmospheric pressure, while P is varied from sub-atmospheric to atmospheric pressure. In contrast, for the thermogravimetric method presented here, P is constant and equal to ambient pressure while the saturation pressure is increased from ambient to approximately 5.0 MPa, depending on the adsorbate and maximum desorption temperature achieved in Atty Dkt #: UWYO-P0012PCT the TGA. The saturation vapor pressure is denoted as PS for the thermogravimetric method to avoid mistaking this value as the constant saturation pressure of the cryogenic physisorption technique, while noting that the physical meaning of P S remains unchanged from P 0 . [0073] Isotherms may be generated from TGA mass loss data by first mapping temperature to P/P S . P/P S, as shown in Eq.3, is a function of temperature only ^^ ൌ ^^ ^ ^ ^^ ^^ ^^^ ^^^ [3] where the temperature at every experimental time step is measured and reported by the TGA. The partial pressure of adsorptive gas, P, within the DSC sample pan 24 is assumed constant and equal to the ambient pressure. The saturation pressure, PS, is found from tabulated thermodynamics data at each analyzed temperature. [0074] Isotherms may then be constructed by plotting ^^ ^^ ^ ^^^ (Eq.2) versus P/P S . Produced isotherms may not extend to P/P S of 0 because this mapping requires the saturation pressure of the adsorbate, which is a function of temperature, and saturation pressure is only defined up to the critical temperature of the adsorbate. These critical temperatures are 319 ^C for toluene and 374 ^C for water, and it is not possible to map to isotherm space above these temperatures, or below the corresponding P/PS, for each adsorbate. Due to the large saturation pressures of toluene and water as temperatures approach the critical temperatures, the unmappable region is relatively small; at a P of 0.078 MPa, toluene could not be mapped at P/PS < 0.025 and water at P/PS < 0.005. [0075] BET transforms may be calculated from isotherms in the conventional manner described by Eq.4 (Step 214). ^ ^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ൌ 1 [4] [0076] The Rouquerol criteria may be applied to determine the appropriate range of P/P S for BET analysis. The normalized adsorption isotherm (NAI) may be calculated by Eq.5. Atty Dkt #: UWYO-P0012PCT ^ ^^ ^^ ^^ ൌ ^^^ ^ 1 െ ^ ^ ^ [5] [0077] The upper may be determined for each adsorbent-adsorbate system by finding the P/PS at which maximum NAI occurred. For the example systems to which BET analysis was applied herein, the upper bounds of P/P S for water/mol sieve, toluene/mol sieve, and toluene/activated carbon were found to be 0.21, 0.27, and 0.19 respectively. The lower limit of the BET range for all systems was selected to be P/PS of 0.05 due to data noise at lower P/PS. [0078] The volume of the BET monolayer, ^^ ^^ ^^ , and the BET constant, ^^, may be obtained using the standard BET equations given by Eq.6 and 7 (Step 216). ^ ^ ^^ ^^ ൌ 1 ^ ^ ^ ^^ [6] ^ ^ ൌ 1 ^ ^^ ^ ^ [7] where ^^ and ^^ are the slope linear fit. BET surface areas may then be found by dividing the monolayer volume, ^^ ^^ ^^ , by the monolayer thickness of the adsorbate. The monolayer thickness for water was taken to be 0.30 nm and the monolayer thickness for toluene was taken to be 0.34 nm in the examples herein. [0079] In some aspects of the present disclosure, pore size distributions may be calculated for saturated porous solids 22 that exhibit multilayer formation and are suitable for BET analysis, such as by applying BJH analysis to the isotherms (Step 220). The calculations applied in steps 222 and 224 below may be used to calculate the pore size distribution for saturated porous solids which do not exhibit multilayer formation. [0080] The Kelvin equation can be used to describe filling and emptying of pores at temperatures above, or pressures below, saturation values based on bulk ambient conditions due to differences in vapor pressure at a curved liquid-vapor interface, such as the surface of adsorbate within a pore. The Kelvin equation can be rearranged to yield Eq. 8. Atty Dkt #: UWYO-P0012PCT ^ ^ ൌ 2 ∙ ^^ ∙ ^^ ∙ cos ^^ ^ ^ ^^ ^ ^^ௌ ^ ^^ ∙ [8] ^^ ^^ where ^^ is adsorbate surface ^^ is the contact angle between liquid adsorbate and the porous solid, ^^ is the ideal gas constant, ^^ is absolute temperature, P S is the saturation pressure of the liquid at a flat interface, and P is the saturation pressure of the liquid over a hemispherical meniscus of radius r. In typical cryogenic physisorption, the temperature of the sample is kept constant (e.g., at -196 ^C for nitrogen physisorption) and the surface tension, molar volume, and flat-interface saturation pressure are approximated as constant over the range of tested pressures. However, when Kelvin pore radius or radii or are calculated (Step 222) isobarically over a temperature range from approximately the adsorbate boiling point to the critical point in the thermogravimetric desorption method, P is constant and equal to the ambient pressure while the values of ^^, ^^, and PS vary. Tabulated values for ^^, ^^, and PS can be obtained from NIST as functions of temperature, allowing Kelvin radii to be determined over the range of temperatures achieved in the TGA during thermogravimetric desorption. Contact angles can be measured using the sessile drop method for systems to which Kelvin analysis was applied, and example results are tabulated in Table 6. Contact angles are not needed for systems to which BET analysis was applied. Table 6. Contact angles measured and used for Kelvin analysis. Silica Gel MCM-41 [0081] Water and toluene have low (≤ 25 ^) contact angles on MCM-41 and silica gel. As the Kelvin equation uses the cosine of the contact angle, the results of the Kelvin equation are relatively insensitive to changes in contact angle at small contact angles, and error resulting from uncertainty in contact angle measurements is relatively low. In contrast, activated carbon is relatively hydrophobic, producing large contact angles with water and making Kelvin equation results sensitive to the contact angle. For example, taking the Atty Dkt #: UWYO-P0012PCT contact angle between activated carbon and water to be 75 ^ produces a surface area of 1,443 m 2 /g, but an uncertainty in contact angle of ± 2 ^ results in the surface area varying from 1,256 m 2 /g to 1,644 m 2 /g. As such, Kelvin analysis is most accurately applied to systems that exhibit low contact angles in order to avoid the uncertainty associated with large contact angles. As an alternative approach for such systems, the surface energies of the solid-vapor and solid-liquid interfaces can be used in lieu of the contact angle for greater accuracy. [0082] In order to calculate a Kelvin pore size distribution from TGA desorption data, approximately 50 temperature points may be chosen between the boiling point of the tested adsorbate and the maximum tested temperature. Average temperatures may then be calculated between adjacent temperature points using Eq.9. ^ ^ ^^ ൌ ^^ ^^ ^ ^^ ^^^1 2 [9] [0083] The loss of adsorbate mass per unit mass of adsorbent between temperature analysis points (∆ ^^ ^^, ^^ ^ may be calculated using Eq.10. ∆ ^^ ^^, ^^ ൌ ^^ ^^, ^^ െ ^^ ^^, ^^^1 [10] [0084] The rate of adsorbate mass loss with temperature increase may be approximated with Eq.11, and a P/P S corresponding to ^^ ^^ ^^ -4 -1 ^ ^ ^^ = 5x10 K may be used as the lower bound for pore size analysis. This lower chosen as it can be consistently applied across different systems and corresponds to >99% of total mass loss. ^^ ^^ ^^ ∆ ^^ ^^, ^^ ^^ ^^ ^^ ^^ [11] [0085] The change in per mass between temperature analysis points (∆ ^^ ^^ ^ may be calculated using Eq.12. ∆ ^^ ∆ ^^ ^^, ^^ ^ ^ Atty Dkt #: UWYO-P0012PCT where ∆ ^^ ^^ is the total pore volume drained during step ^^ per unit mass of adsorbent, and ^^^ ^^ ^ ^ is density of the liquid adsorbate as a function of temperature. The change in Kelvin radius over the temperature step, ∆ ^^ ^ , may be calculated using Eq.13. 2 ∙ ^^^ ^^^^ ∙ ^^^ ^^^^ ∙ cos ^^ 2 ∙ ^^^ ^^^ା^^ ∙ ^^^ ^^^ା^^ ∙ cos ^^ ∆ ^^ ^ ^ ^ ^^ ^^ௌ^ ^^^^ ^ ^^ ∙ ^^ ^^ ^^ ^^ௌ^ ^^^ା^^ ^ ^^ ∙ ^^ [13] [0086] Eq.14. ^ ^ 2 ∙ ^^^ ^^ ^^^ ∙ ^^^ ^^ ^^ ^ ^ ൌ ^ ∙ cos ^^ ^ ^ ^^ ^^ ^^^ ^^ ^^^ ^ ^^ ^^ [14] [0087] Cylindrical in pore radius during step ^^ may then be calculated (Step 224) using Eq.15. ∆ ^^ ^^ ∆ ^^ ^^ 2∆ ^^ ^^ ∆ ^^ ^^ ∙ ^^ ^^ [15] [0088] The total surface surface area over pore radius, and this may be numerically evaluated (Step 226) using the trapezoidal method, Eq.18. ^ ^ ^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^^ ^ ∆ ^^ ^ ^ ^^ 1 ே ^ ^ ^^ ^^ ^ ^ ∆ ^^^ ^ ∆ ^^^ା^ [0090] The disclosure is not to be limited to the particular embodiments described herein. In particular, the disclosure contemplates numerous variations in quantifying porous material. The foregoing description has been presented for purposes of illustration and description. It is not intended to be an exhaustive list or limit any of the disclosure to the precise forms disclosed. It is contemplated that other alternatives or exemplary aspects are considered included in the disclosure. The description is merely examples of embodiments, processes or methods of the disclosure. It is understood that any other modifications, Atty Dkt #: UWYO-P0012PCT substitutions, and/or additions can be made, which are within the intended spirit and scope of the disclosure.