Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pore size distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗
DOE OSTI2024

Mesoscale simulation of granular materials under weak shock compaction–pore size distribution effects

This research established a systematic method to generate various pore-size distributions (PSDs) and studied the effect of PSDs on the shock compaction response of granular materials using two-dimensional mesoscale simulations under identical porosity. Simulations utilized various PSDs for three particle shapes (circle, ellipse, and square). Contacting particle configurations using three PSDs, characterized by spatially uniform distributed pores to heterogeneous distributed pores, and non-contacting particle configurations under a single case of PSD were tested. The PSD of generated particle sets was characterized using coordination number, mean diameter, and bimodality coefficient as statistical metrics. Mesoscale simulations showed that regardless of the conditions of pore distributions, shock compaction of granular materials consistently demonstrates a precursor, shock compaction front, and end. However, the shock compaction velocity of contacting particles was dependent on the PSDs despite the constant initial porosity. The compaction velocity was faster in particle configurations with relatively uniform pore distributions than in heterogeneous pore distributions, which our study demonstrated can be attributed to particle rearrangement during compaction. Circular-shaped particles had high sensitivity in shock compaction response to the various PSDs. Furthermore, a contacting particle configuration tended to propagate the shock compaction wave relatively faster than particles that were in a non-contact configuration. This study established the relative importance of considering PSD as a metric over the coordination number in studies of the shock compaction response of granular materials. Further, insights are provided on the evolving shock substructure to characterize the shock compaction response of granular materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Pore2Chip: All-in-one python tool for soil microstructure analysis and micromodel design

The Pore2Chip Python package is designed to create 2D micromodels using extracted data from 3D X-ray computed tomography (XCT) images. This package helps analyze soil structure and function, allowing for the investigation of hydro-biogeochemical processes that impact mineral extraction and reactivity, oxygen concentrations, and nutrient availability in disturbed or managed soils. Key metrics encompass pore size distributions, pore throat size distributions, and connectivity (pore coordination numbers). The final output is a 2D scalable SVG design representing a core or aggregate. Designs can be fabricated with methods such as laser etching, 3D printing, and photolithography.

lab-on-chip↗
DOE OSTI2026

Impact of Porous Transport Layer Morphology on the Performance of Proton Exchange Membrane Water Electrolyzers with Ultra-Low Iridium Loadings

Reducing Ir loadings in proton exchange membrane water electrolyzer anodes is critical for lowering capital expenses. Loading reduction could be achieved by improving the Ir activity via doping/alloying and/or the development of advanced microstructures. However, the anode porous transport layer (PTL) is a comparatively simple component whose properties also impact Ir utilization. Therefore, well-designed PTLs may also enable reduced Ir loadings. In this work, we survey eight PTLs from various manufacturers to observe their impact on cell performance at low (0.4 mg Ir cm -2 ) and ultralow (0.1 mg Ir cm -2 ) Ir loadings. The PTLs were characterized by their microstructural properties, including porosity, particle size distribution, and pore size distribution. Electrochemical cell performance was correlated to PTL morphology, and it was found that PTLs with lower porosities and smaller particle and pore radii enabled good performance even at ultralow Ir loadings. 1000-h durability testing indicated that using lower porosity PTLs can significantly improve durability behavior. A runaway voltage phenomenon was observed during durability testing of cells with ultralow Ir loadings, which was caused by increases in both anode and cathode overpotentials. Furthermore, we observed that the beginning of test performance of 0.1 mg Ir cm -2 cells correlates to the 1000-h degradation rates of 0.4 mg Ir cm -2 cells, suggesting that for the Ir catalyst used in this work, short-term testing at ultralow loadings can be used as an indicator of long-term degradation at higher loadings.

08 HYDROGEN↗
DOE OSTI2022

Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates

Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .

58 GEOSCIENCES↗
DOE OSTI2023

Functionalized porous organic polymers as uranium nano-traps for efficient uranium extraction

Compositions are provided for efficient uranium extraction, for example from wastewater, seawater, or other water sources. The compositions can include a functionalized porous organic polymer functionalized with one or more uranium binding moieties, e.g. having a plurality of amidoxime or amidrazone groups covalently attached thereto. The compositions can include covalent organic frameworks, porous aromatic frameworks, and various porous organic polymers, especially those having a hierarchical pore size distribution over a range of pore sizes. The compositions can have functional groups such as amidoxime or an amidrazone covalently attached thereto. The hierarchical pore size distribution can be determined based upon at least 60% of the pore sizes in the range of pore sizes having a pore volume of at least 0.01 cm 3 g −1 in the pore size distribution at 77 K. Methods of making the compositions and methods of using the compositions are also provided.

Ma, Shengqian↗
DOE OSTI2022

Multiple experimental studies of pore structure and mineral grain sizes of the Woodford shale in southern Oklahoma, USA

Pore structure study is an important part of unconventional shale reservoir characterization, since the pore system provides the primary petroleum storage space and fluid flow pathways. Previous studies have suggested that the pore structure is related to the total organic carbon (TOC) content, mineral compositions, and the maturity of the organic matter (OM). However, few studies have focused on the mineral grains, the primary grains being deposited but before cementation, which are the building blocks of shale. Eight Woodford Shale outcrop samples from southern Oklahoma were chosen to study the effects of mineral grain size on the pore structure characterization, using multiple and complementary experimental approaches, including laser diffraction, mineralogy, TOC, pyrolysis, liquid immersion porosimetry, mercury intrusion porosimetry, gas physisorption, (ultra) small angle X-ray scattering, scanning electron microscopy, and spontaneous imbibition. The results from different experiments of eight samples show that the Woodford Shale has the mean mineral grain diameters at 3–6 μm, a wide range of porosity at 3–40% and pore diameters at 50–1,000 nm, and various pore connectivity. Grain size variation was probably caused by the sea-level fluctuation during its deposition, which affect the porosity, pore size distribution, and pore connectivity. With decreasing mineral grain sizes, the porosity tends to increase while the pore connectivity worsens. The results also indicate that OM and carbonates in this low-maturity Woodford Shale could block the pores and decrease the porosity. Coupling with the grain size analyses, the control of depositional environment on grain sizes and subsequent effects on pore structure is identified. The pore structure characteristics over a wide pore-diameter range provided by multiple experiments could improve the understanding of storage space and fluid flow in the Woodford Shale to further increase its petroleum production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using ultrasonic attenuation in cortical bone to infer distributions on pore size

Here, in this work we infer the underlying distribution on pore radius in human cortical bone samples using ultrasonic attenuation data. We first discuss how to formulate polydisperse attenuation models using a probabilistic approach and the Waterman Truell model for scattering attenuation. We then compare the Independent Scattering Approximation and the higher-order Waterman Truell models’ forward predictions for total attenuation in polydisperse samples. Following this, we formulate an inverse problem under the Prohorov Metric Framework coupled with variational regularization to stabilize this inverse problem. We then use experimental attenuation data taken from human cadaver samples and solve inverse problems resulting in nonparametric estimates of the probability density function on pore radius. We compare these estimates to the “true” microstructure of the bone samples determined via microCT imaging. We find that our methodology allows us to reliably estimate the underlying microstructure of the bone from attenuation data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Biomass char gasification kinetic rates compared to data, including ash effects

This investigation provides both experimental and theoretical analyses of biomass char gasification kinetics over 1150–1350 °C in the environment of CO 2 , H 2 O and combinations of these reactants. The experimental data include continuous measurements of particle mass, external and internal temperatures, shape, and size with periodic but not continuous measurements of char porosity, pore size distribution, and pore surface area. The model describes heterogeneous char gasification for poplar wood, corn stover and switchgrass in different shapes and industrially relevant sizes (12.5 mm and smaller). Nonlinear least-squares regression produces optimized power-law model parameters for global kinetics that describe biomass char gasification with respect to both CO 2 and H 2 O separately and in combination. A single set of kinetic parameters describes reaction rates for all three fuels. Model simulations agree with measured data at all stages of char conversion even though the observed mass loss rates change significantly with conversion. In conclusion, the change in rate with conversion depends strongly on ash content, and this investigation provides a simple and theoretically based reaction rate expression that includes the effect of inert ash occupying an increasing portion of the reactive surface with increased conversion.

09 BIOMASS FUELS↗
DOE OSTI2024

Thermally Insulating Transparent Barrier (THINNER) coatings on single pane windows

Conventional silica aerogel monoliths can provide remarkable thermal insulation but the presence of large pores (> 30 nm) tends to scatter visible light and render the material opaque or translucent instead of transparent. In addition, they are prone to cracking during synthesis and handling which makes them difficult to integrate in products and in particular in window solutions. This project developed two new solgel synthesis methods using silica precursors or preformed silica nanoparticles and ambient drying to produce mesoporous organo-silica monoliths. The monoliths were (i) thermally insulating, (ii) optically transparent, (iii) flexible, and (iv) hydrophobic. They feature porosity ranging from 50% to 90% with narrow pore size distribution with pore less than 20 nm resulting in excellent optical clarity (haze < 2%) and very low thermal conductivity (< 30 mW/mK). Interestingly, not only porosity but also pore size and mass fractal dimension were found to affect the thermal conductivity of the mesoporous silica. The superior transparency of the monoliths was shown to be attributed to dependent scattering among silica nanoparticles. Flexibility was achieved through trimethylchlorosilane surface modification. Furthermore, process scale-up and integration of the ambigel monoliths into window solutions using optically clear adhesives were also demonstrated in 6”x6” double-pane windows. The aerogel and the window solution were shown to be durable to accelerated aging under UV, moisture, and/or temperature gradient and thermal cycling. Overall, the window solution achieved the technical performance and the cost target of $10/sqft set for the SHIELD program. However, scaling to industry relevant scale (> 10’x10’) remains a challenge due to requirements on the drying process and fume hood size and to the propensity of large aerogel slab to crack during drying.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Hierarchically porous electrospun carbon nanofiber for high-rate capacitive deionization electrodes

Capacitive deionization (CDI) is a promising technology that has gained interest for the desalination of brackish water. Hierarchically porous carbons are commonly used as electrodes for CDI due to their high surface areas and controlled pore size distributions that maximize ion adsorption capacity and rate. Electrospinning is an effective way of generating carbon nanofibers with high inter-fiber macroporosity that can be further modified to improve surface area, total pore volume, and pore size distribution. This work describes the use of sacrificial mesopore formers in tandem with a micropore etching technique to induce hierarchical porosity in electrospun fibers. Mesopores are formed via the dissolution of silica nanoparticles that are introduced into the fibers during the electrospinning step. After mesopore formation, micropores are etched into the resulting surface through KOH impregnation and thermal activation. This sequential technique creates a hierarchical network of pores from the inherent macroporosity of the fiber network, to the mesopores, and finally micropores to simultaneously maximize surface area and accessibility. Micropore formation is optimized to maximize specific surface area while maintaining physical integrity of the fibers. Further, the combination of mesopores and micropores enables fast ion adsorption rates and capacity. Carbon fiber electrodes fabricated in this method achieve specific surface areas exceeding 1400 m 2 g -1 , with pore volumes exceeding 1.0cc g -1 . The pore size distributions are highly controlled, with 80% of total pore volume coming from pores <20nm in radius. In 500 ppm constant voltage CDI tests, these fiber electrodes obtain a salt adsorption capacity of over 14 mg g -1 at a salt adsorption rate of ~4mg g -1 min -1 , showcasing the high capacity matched with high rate of these easily fabricated, inexpensive materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stochastic generation of electrolyzer anode catalyst layers

Here, we introduce a stochastic methodology to reproduce the complex pore structure observed in commercial iridium catalyst layers. This method preserves the α pore (pores smaller than 250 nm) and β pore (pores greater than or equal to 250 nm) regions of the catalyst layer. The morphology of the generated materials was validated by comparing the pore size distributions of generated materials against those obtained from commercial materials imaged using x-ray nano computed tomography. We further demonstrate that the pore size distributions of the generated materials are statistically indistinguishable from the imaged catalyst layers, indicating that the stochastic methodology is capable of accurately reproducing catalyst layer morphology. Pore network modelling was conducted on the generated catalyst materials to simulate single-phase permeability, electrical conductivity, and ionic conductivity, and these properties were found to be within experimentally measured ranges for electrolyzer catalyst layers. Additionally, simulations were performed on the generated materials with varying ionomer and iridium catalyst loadings. As the ionomer loading is added, proton conductivity increases exponentially, which demonstrates the importance of optimizing ionomer loading, considering that these effects will be exacerbated in the hydration and temperature conditions of operating electrolyzers. The stochastic material generation method presented in this work is a powerful tool for the development of novel low loading catalyst layers, where the effect of various structural parameters on electrolyzer performance characteristics can be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Insights into determining pore size properties of ultrafiltration membranes

The selectivity of porous membranes is often characterized using solute rejection tests, where membranes are challenged with dilute aqueous solutions of neutral solutes at operating conditions that minimize concentration polarization and fouling. In single solute tests, a membrane is challenged with one molecular weight (MW) solute at a time from low to high MW. Since single solute methods are time-intensive, mixed solute tests have become more common, where a mixture of several MW solutes challenges a membrane at once. However, the presence of large solutes in a mixture increases the rejection of smaller solutes. Furthermore, there are no universally accepted operating conditions or standard methods used by membrane manufacturers or researchers for the experiments, leading to difficulty in pore size and pore characteristic comparisons. In this paper, commercial ultrafiltration membranes were challenged with single and mixed solute polyethylene glycol (PEG) and dextran aqueous solutions. First, rejection values determined using total organic carbon (TOC) and high-performance liquid chromatography (HPLC) from single solute filtration experiments are compared. Differences in rejection curves obtained by the two techniques are attributed to solute polydispersity. Mixed solute filtration experiments with binary mixtures of solutes showcased co-solute interactions, which increase with both the size and weight percent of large solute in the mixture. Mixed solute filtration experiments at varying operating conditions (i.e., stir speed and flux) were conducted to determine operating conditions that mitigate co-solute interactions. Stir speed had a minimal effect on co-solute interactions. In contrast, low flux conditions can help minimize co-solute interactions, leading to pore size distributions that closely resemble results observed in single solute filtration using narrowly dispersed solutes. Additionally, at low flux conditions, the predicted membrane pore size distributions utilizing mixed solute experiments with PEG and dextran were similar.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A User-Friendly GUI Tool for Automated Microstructural Analysis of Fiber-Reinforced Composites and Porous Structures

Understanding and quantifying microstructural features such as fiber orientation and porosity is critical for predicting the mechanical behavior and performance of fiber-reinforced polymer composites. Traditional manual analysis is time-consuming, subjective, and unsuitable for high-throughput datasets. We present a graphical user interface (GUI) application that automates the analysis of microscopy images to extract key microstructural metrics, including fiber orientation tensors, fiber orientation distribution, porosity and pore size distribution. The app integrates multiple image segmentation techniques including global and local thresholding, clustering, and region-based approaches, offering flexibility for different types of image qualities and features. Users can load microstructural images, select regions of interest and segmentation techniques tailored to their image dataset. It also addresses a critical challenge in fiber orientation analysis: the ambiguities caused by touching, overlapping, or partially cut fibers. It supports autorun examples for standardized workflows, enabling reproducible analysis and facilitating training and benchmarking. This tool significantly reduces manual intervention, enhances consistency, and accelerates data generation for structure–property modeling, process optimization, and digital materials research. The tool is intended for use by materials scientists, engineers, and researchers engaged in composite characterization, quality control, and machine learning-based microstructural studies.

Chawla, Komal [ORNL] (ORCID:0000000190327565)↗
DOE OSTI2022

Reduced diffusion and enhanced retention of multiple radionuclides from pore structure characterization of barrier materials for enhanced repository performance

Fluid flow and chemical transport in porous media are the macroscopic consequences of pore structure, which integrates geometry (e.g., pore size and surface area, pore-size distribution) and topology (e.g., pore connectivity). Low-permeability geological media whose pores are poorly interconnected will exhibit the characteristics of anomalous diffusion and sample size-dependent effective porosity, which will strongly impact long-term net diffusion and retention of radionuclides in geological repository settings involving different host rocks and barrier materials. A suite of innovative and complementary experimental approaches is utilized to study the microscopic pore structure and macroscopic fluid flow & chemical transport for a range of host rocks and barrier materials, in addition to standard clay minerals and reference rocks. With a particular focus on quantifying the presence and magnitude of “isolated” pores for a reduced effective porosity in low-permeability geomedia, the integrated methodologies for basic properties and pore structure characterization of these geomedia include X-ray diffraction, thin section petrography, grain size distribution, water immersion porosimetry after vacuum-pulling for full saturation, mercury intrusion porosimetry, nitrogen physisorption, scanning electron microscopy, X-ray computed tomography, and (ultra-)small angle neutron (X-ray) scattering. In addition, custom-designed gas diffusion, tracer recipe involving a range of anionic and cationic chemicals with subsequent analyses by laser ablation and inductively coupled plasma-mass spectrometry, along with batch sorption, column transport, and imbibition tests were conducted for coupled effects of pore structure and chemical retention/transport. From the perspectives of pore structure in conjunction with multiple and complementary approaches to examining a range of sample sizes under different observational scales, we find that the poor pore connectivity is prevalent in low-permeability media (mudstone and crystalline rock) that is related to geological processes (e.g., compaction, diagenesis and thermal maturation). For example, the deep and organic matter-rich mudstones have a much smaller effective porosity than the total porosity (as a result of poor pore connectivity) and associated diffusion coefficient, and the effective porosity & diffusion coefficients are also dependent upon the sample sizes used in the measurement. Similarly, most of the pore space in the shallow mudstone is also controlled by pore-throat diameters in the 5-50 nm range of intergranular pore types from its fine-grained nature, but with an overall good pore connectivity. However, the nm-sized pore space (physically pore-network architecture) and strong sorption capacities (chemical retention from clay minerals) of both shallow and deep mudstones lead to the synergistic retention of cationic radionuclides and their utilities as effective host rocks and barrier materials. Our unique approaches of studying how the micro-scale pore structure affect macro-scale fluid flow, diffusion & retention, and chemical transport produce improved mechanistic understanding, and realistic quantification, of diffusion and retention of typical radionuclides in a range of generic host rocks and barrier materials (clay/shale, salt, crystalline rock, and tuff), with the overall results leading to scientifically-based understanding of enhanced isolation (from both diffusion and retention) of radionuclides and improved confidence on the long-term performance of geological repository to store high-level radioactive wastes. In addition to the training of 25 undergraduates, graduates, and postdocs of UTA, the scientists (organizations) involved in performing this work (e.g., discussion, sample sharing, and operation of SANS and SAXS instruments) include Ed Matteo, Yifeng Wang, and Kristopher Kuhlman (Sandia National Laboratories), Jens Birkholzer, Liange Zheng, Tim Kneafsey, and Sharon Borglin (Lawrence Berkeley National Laboratory), Mavrik Zavarin (Lawrence Livermore National Laboratory), Yukio Tachi and Yuta Fukatsu (Japan Atomic Energy Agency), Mieke de Craen (Euridice, Belgium), Markus Bleuel (NIST), Wei-Ren Chen, Gergely Nagy, Changwoo Do, William Heller, Larry Anovitz, and Kenneth Littrell (ORNL), as well as Jan Illvsky, Ivan Kuzmenko, Ju-Sang Park and Jon Almers (ANL). Key deliverables include a total of 13 peer-reviewed journal articles (nine published and three under review), 23 presentations at scientific conferences (AAPG, AAPG Southwest Section, AGU, Asian Clay Conference, GSA, GSA South-Central Section, IHLRWM, InterPore, International Conference on Chemistry and Migration Behavior of Actinides and Fission Products in the Geosphere, International Conference on Coupled Processes in Fractured Geological Media: Observation, Modeling and Application), and academic institutions (UTA, New Mexico State University; University of Poitiers, France; University of Helsinki, Finland; Uppsala University, Sweden; Istanbul Technical University, Turkey) and other organizations (Andra, France; Posiva Oy, Finland).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Measuring Multicomponent Adsorption of Tracer Gases on Natural Zeolites

A natural clinoptilolite sample near the Nevada National Security Site was obtained to study adsorption and retardation on gas transport. Of interest is understanding the competition for adsorption sites that may reduce tracer gas adsorption relative to single-component measurements, which may be affected by the multi-scale pore structure of clinoptilolite. Clinoptilolite has three distinct domains of pore size distributions ranging from nanometers to micrometers: micropores with 0.4–0.7 nm diameters, measured on powders by CO2 adsorption at 273 K, representing the zeolite cages; mesopores with 4–200 nm diameters, observed using liquid nitrogen adsorption at 77 K; and macropores with 300–1000 nm diameters, measured by mercury injection on rock chips (~ 100 mesh), likely representing the microfractures. These pore size distributions are consistent with X-ray computed tomography (CT) and focused ion beam scanning electron microscope (FIB-SEM) images, which are used to construct the three-dimensional (3D) pore network to be used in future gas transport modeling. To quantify tracer gas adsorption in this multi-scale pore structure and multicomponent gas species environment, natural zeolite samples initially in equilibrium in air were exposed to a mixture of tracer gases. As the tracer gases diffuse and adsorb in the sample, the remaining tracer gases outside the sample fractionate. Using a quadrupole mass spectrometer to quantify this fractionation, the degree of adsorption of tracer gases in the multicomponent gas environment and multi-scale pore structure is assessed. The major finding is that Kr reaches equilibrium much faster than Xe in the presence of ambient air, which leads to more Kr uptake than Xe over limited exposure periods. When the clinoptilolite chips were exposed to humid air, the adsorption capability decreases significantly for both Xe and Kr with relative humidity (RH) as low as 3%. Both Xe and Kr reaches equilibrium faster at higher RH. The different, unexpected, adsorption behavior for Xe and Kr is due to their kinetic diameters similar to the micropores in clinoptilolite which makes it harder for Xe to access compared to Kr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗