Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pore pressure”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Porous membrane enabled mass spectrometry characterization of microfluidic devices

A system for sampling a liquid includes a sample fluid conduit including a membrane having pores. The membrane prevents the passage of the sample liquid through the pores at a first pressure of the sample liquid in the sample fluid conduit. A surface sampling capture probe has a distal end. The capture probe includes a solvent supply conduit and a solvent exhaust conduit. A solvent composition flowing at the distal end of the capture probe establishes a liquid junction with the membrane and establishes a second pressure within the liquid junction at the membrane. The second pressure is lower than the first pressure. Sample liquid will be drawn through the pores of the membrane by the second pressure at the liquid junction. A method for sampling a liquid and for performing chemical analysis on a liquid are also disclosed.

Cahill, John F.↗
DOE OSTI2024

Porous membrane enabled mass spectrometry characterization of microfluidic devices

A system for sampling a liquid includes a sample fluid conduit including a membrane having pores. The membrane prevents the passage of the sample liquid through the pores at a first pressure of the sample liquid in the sample fluid conduit. A surface sampling capture probe has a distal end. The capture probe includes a solvent supply conduit and a solvent exhaust conduit. A solvent composition flowing at the distal end of the capture probe establishes a liquid junction with the membrane and establishes a second pressure within the liquid junction at the membrane. The second pressure is lower than the first pressure. Sample liquid will be drawn through the pores of the membrane by the second pressure at the liquid junction. A method for sampling a liquid and for performing chemical analysis on a liquid are also disclosed.

Cahill, John F.↗
DOE OSTI2021

Deformation of pores in response to uniaxial and hydrostatic stress cycling in Marcellus Shale: Implications for gas recovery

One of the main challenges during gas production from shale reservoirs is low recovery rate. One contributing factor to this outcome is an insufficient understanding of pore systems, especially pore behaviour following changes in reservoir conditions or resulting from gas production practices. Because the pressure in the producing well can be controlled, understanding the effects of pressure variation on the pore size distribution and methane trapping is necessary to help design optimal conditions to improve the gas recovery rate. This work is the first systematic study of sub-millimeter pore deformation in shale caused by uniaxial and hydrostatic stress up to 100 MPa. Overmature samples from the Middle Devonian Marcellus Shale were analyzed using neutron scattering (SANS and USANS) techniques to interpret the response of nanopores to stress cycling of magnitude and duration compatible with the hydraulic fracturing treatments. Experiments reported here are performed at a series of uniaxial pressures up to 100 MPa and at hydrostatic pressures of deuterated methane 0 and 50 MPa. Since at the original depth of the shale samples’ burial of 2184 m the hydrostatic pressure is approximately 27 MPa and the lithostatic pressure is about 55 MPa, the experimental conditions reasonably well simulate the reservoir pressure regime. Our SANS and USANS results show that different pore sizes are affected by uniaxial stress in different ways. Specifically, in the pore size range from 1 nm to 800 nm, a decrease of pore density with pressure is observed, with the most depleted being mesopores of about 100 nm in diameter. The observed decrease is likely related to deformation of kerogen, followed by a loss of pore nano-volume, as well as methane trapped in the micropores. For pores larger than 5 μm, USANS data suggest that the negative trend is reversed at above 74 MPa, and the number density of large macropores may increase with increased stress even above the original value. The increased number of macropores at high pressure may create new interconnected conduits for gas migration, resulting in a better recovery rate. Another important finding of this study is an irreversible rearrangement of pore size distribution taking place after pressure cycling. Furthermore, this irreversible reorganization of pore size distribution should be taken into account during management of well production to maximize recovery rate.

03 NATURAL GAS↗
DOE OSTI2020

Microstructure Experiments-Enabled MARMOT Simulations of SiC/SiC-based Accident Tolerant Nuclear Fuel System

We have undertaken an experimental-computational project that addresses a few key technology gaps associated with the use of SiC/SiC composites for light water reactor fuel cladding. Two principal endeavors in this project are to assess the irradiation-induced microstructural changes and swelling in SiC/SiC composites, and to characterize and model the porous oxide surface layer (when exposed to steam) along SiC recession that is detrimental to the clad integrity under accident conditions. The project tasks are: (i) ion irradiation and characterization (University of Tennessee, Knoxville), (ii) mechanical tests and analysis (University of South Carolina), (iii) steam exposure tests (Oak Ridge National Laboratory), (iv) electron microscopy and spectroscopic characterization (NC State University), (v) x-ray microscopy and reconstruction (University of South Carolina and NC State University), (vi) phase field modeling and simulations (NC State University and Idaho National Laboratory). The effects of 10 MeV Au ion irradiation at 350°C on the microstructure evolution in SiC/SiC composites are investigated at doses up to 400 displacements per atom (dpa) at the University of Tennessee, Knoxville. Atomic force microscopy and optical profilometry reveal irradiation induced axial and radial shrinkage of the fibers for doses greater than 10 dpa. Based on detailed electron microscopy characterization, the primary cause of the fiber shrinkage is attributed to irradiation-induced loss of carbon packets. Additionally, the multilayer PyC interface is observed to portray high resistance to irradiation damage. The mechanical response and failure mechanisms of un-irradiated samples is also assessed through loading tests and X-ray imaging at the University of South Carolina. Steam exposure tests are performed at the Oak Ridge National Laboratory in a facility that represents a reactor pressure vessel under a loss-of-coolant-accident scenario. The samples that are analyzed methodically in this report are tested for 32/31 hours at 1200°C with a velocity of 0.25 cm/s for the pressures: 0.1 MPa, 0.45 MPa, 0.92 MPa and 1.38 MPa. Scanning/transmission electron microscopy analysis conducted at the NC State University (NCSU) shows that the oxide layer thickness increases with the steam pressure. While the oxide layer is crystalline (α- cristobalite) for the pressures 0.45 MPa, 0.92 MPa and 1.38 MPa, the layer is amorphous at 0.1 MPa. Results from Raman spectroscopy have also confirmed the formation of α-cristobalite phase of SiO₂. The abrupt increase in the integrated Raman intensity ratio between steam pressures 0.1 and 0.45 MPa suggests the onset of accelerated crystallization. Non-destructive three dimensional X-ray microscopy/tomography (XCT) and computational image processing techniques are employed by the University of South Carolina to probe the porous oxidation features on SiC samples at varying pressures. Interestingly, most pores are observed to be located away from the surface as well as the oxide-SiC interface. The average oxide layer thickness assessed from the XCT analysis is seen to be in excellent agreement with the values determined through scanning electron microscopy analysis for the highest pressures where the oxide layer is relatively more uniform. A phase-field model developed by the NC State University and Idaho National Laboratory for simulating oxidation of SiC by steam captures the paralinear kinetics of SiC oxidation with a high degree of fidelity. Results from quasi-one-dimensional and two dimensional simulations show that the pores with oxidizing species lead to a higher volatilization rate. These simulations indicate that the enhanced apparent volatilization rates at higher pressures observed in experiments can be rationalized by the increased volatilization from the pores. The project team has also successfully developed the capability to import images from experiments for realistic evolution of the oxidizing microstructure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Lattice expansion and ligand twist during CO2 adsorption in flexible Cu bipyridine metal–organic frameworks

Flexible metal–organic frameworks (MOFs) can show exceptional selectivity and capacity for adsorption of CO 2 . The incorporation of CO 2 into flexible MOFs that have Cu 2+ coordination centers and organic pillar ligands is accompanied by a distortion of the framework lattice arising from chemical interactions between these components and CO 2 molecules. CO 2 adsorption yields a reproducible lattice expansion that is enabled by the rotation of the pillar ligands. The structures of Cu 2 (pzdc) 2 (bpy) and Cu 2 (pzdc) 2 (bpe), CPL-2 and CPL-5, were evaluated using in situ synchrotron X-ray powder diffraction at room temperature at CO 2 gas pressures up to 50 atm. The structural parameters exhibit hysteresis between pressurization and depressurization. The pore volume within CPL-2 and CPL-5 increases at elevated CO 2 pressure due to a combination of the pillar ligand rotation and the overall expansion of the lattice. Volumetric CO 2 adsorption measurements up to 50 atm reveal adsorption behavior consistent with the structural results, including a rapid uptake of CO 2 at low pressure, saturation above 20 atm, and hysteresis evident as a retention of CO 2 during depressurization. Finally, a significantly greater CO 2 uptake is observed in CPL-5 in comparison with predictions based on CO 2 pressure-induced expansion of the pore volume available for adsorption, indicating that the flexibility of the CPL structures is a key factor in enhancing adsorption capacity.

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

Surface ionization waves propagating over non-planar substrates: wavy surfaces, cut-pores and droplets

Atmospheric pressure plasmas intersecting with dielectric surfaces will often transition into surface ionization waves (SIWs). Several applications of these discharges are purposely configured to be SIWs. During propagation of an SIW over a dielectric surface, the plasma charges the surface while responding to changes in geometrical and electrical material properties. This is particularly important for non-planar surfaces where polarization of the dielectric results in local electric field enhancement. In this paper, we discuss results from computational investigations of negative and positive SIWs propagating over nonplanar dielectrics in three configurations—wavy surfaces, cuts through porous materials and water droplets on flat surfaces. Here we found that negative SIWs are particularly sensitive to the electric field enhancement that occurs at the crests of non-planar surfaces. The local increase in ionization rates by the electric field enhancement can result in the SIW detaching from the surface, which produces non-uniform plasma exposure of the surface. Positive SIWs tend to adhere to the surface to a greater degree. These trends indicate that treatment of pathogen containing droplets on surfaces may be best performed by positive SIWs. The same principles apply to the surfaces cut through pores. Buried pores with small openings to the SIW may be filled by plasma by either flow of plasma into the pore (large opening) or initiated by photoionization (small opening), depending on the size of the opening compared to the Debye length.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effects of Applied Interfacial Pressure on Li-Metal Cycling Performance and Morphology in 4 M LiFSI in DME

Lithium-metal anodes can theoretically enable 10x higher gravimetric capacity than conventional graphite anodes. However, Li-metal anode cycling has proven difficult due to porous and dendritic morphologies, extensive parasitic solid electrolyte interphase reactions, and formation of dead Li. We systematically investigate the effects of applied interfacial pressure on Li-metal anode cycling performance and morphology in the recently developed and highly efficient 4 M lithium bis(fluorosulfonyl)imide in 1,2-dimethoxyethane electrolyte. We present cycling, morphology, and impedance data at a current density of 0.5 mA/cm 2 and a capacity of 2 mAh/cm 2 at applied interfacial pressures of 0, 0.01, 0.1, 1, and 10 MPa. Cryo-focused ion beam milling and cryo-scanning electron microscopy imaging in cross section reveal that increasing the applied pressure during Li deposition from 0 to 10 MPa leads to greater than a fivefold reduction in thickness (and therefore volume) of the deposited Li. This suggests that pressure during cycling can have a profound impact on the practical volumetric energy density for Li-metal anodes. A “goldilocks zone” of cell performance is observed at intermediate pressures of 0.1–1 MPa. Increasing pressure from 0 to 1 MPa generally improves cell-to-cell reproducibility, cycling stability, and Coulombic efficiency. However, the highest pressure (10 MPa) results in high cell overpotential and evidence of soft short circuits, which likely result from transport limitations associated with increased pressure causing local pore closure in the separator. All cells exhibit at least some signs of cycling instability after 50 cycles when cycled to 2 mAh/cm 2 with thin 50 μm Li counter electrodes, though instability decreases with increasing pressure. In contrast, cells cycled to only 1 mAh/cm 2 perform well for 50 cycles, indicating that capacity plays an important role in cycling stability.

25 ENERGY STORAGE↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
DOE OSTI2020

Implications of post-processing induced microstructural changes on the deformation and fracture response of additively manufactured Ti–6Al–4V

Ti–6Al–4V parts fabricated via electron beam melting (EBM) powder bed fusion are often subjected to hot isostatic pressing below the β-transus temperature to mitigate defects. During which the pressure aids in pore closure and the thermal exposure results in coarsening of α phase while retaining the columnar prior-β phase grain morphology present in the as-fabricated condition. The same post-processing treatment can also be carried out above β-transus temperature which is an effective way to modify both the grain morphology and the associated α textures. Here, the objective of this work is to correlate thermally induced microstructural changes to the deformation and fracture response of the EBM processed Ti–6Al–4V. To this end, we have carried out both sub-transus and super-transus heat-treatments of the as-fabricated material. The mechanical response of the as-fabricated and all heat-treated materials are characterized by in-situ tensile tests under a high-resolution digital optical microscope. This enabled us to capture large-scale panoramic images of the deforming microstructure, and overcome the trade-off between the image resolution and the field of view during in-situ experiments. The series of images captured throughout the imposed deformation are subsequently used to perform microstructure-based digital image correlation to measure microstructural-scale strains over a large area. The results of the in-situ tests together with detailed fractographic analyses are then used to elucidate how the heterogeneous deformation spanning over multiple microstructural length-scales, for example, at the scale of lamellae, colonies of lamellae and grain boundaries, affect the overall deformation and fracture response of the post-processed materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ethane-CO2 Mixture Adsorption in Silicalite: Influence of Tortuosity and Connectivity of Pores on Selectivity

Selective adsorption using nanoporous materials is an efficient strategy for separating gas mixtures. In a nanoporous material, pores can exist in different shapes and can have different degrees of inter-connectivity. In recent studies, both pore connectivity and tortuosity have been found to affect the adsorption and dynamical properties of ethane and CO2 in silicalite differently. Here, using Monte Carlo simulations, we investigate if these two attributes can affect the selective adsorption of one component from a mixture of ethane and CO2 in silicalite. For this, the adsorption of an equimolar mixture of ethane and CO2 is simulated in 12 models of silicalite—SnZm (n, m = 0, 1, 2, 3 or 4; with n and m denoting, respectively, the fraction (out of 4) of straight and zigzag channels of silicalite that are available for adsorption)—differing in degrees of pore connectivity and tortuosity. The adsorption selectivity in this system is found to exhibit a reversal with the adsorption dominated by ethane at low pressures (below ~1 atm) and by CO2 at higher pressures (above ~10 atm). Pore connectivity is found to suppress the selective adsorption of CO2 at higher pressures and also shifts the selectivity reversal to higher pressures. The selectivity reversal results from a competition between the polarizability-affected adsorption at lower pressures and efficient packing at higher pressures. The efficient packing of CO2 is a compounded effect resulting from the larger effective pore volume available for CO2 due to its stronger interaction with the pore surface and smaller molecular volume. CO2 molecules show a preference to adsorb in non-tortuous pores, and this preference is found to be stronger in the presence of ethane. The effects of pore connectivity and tortuosity elucidated here should be applicable to a wide range of natural and engineered nanoporous materials, and this knowledge could be used to identify materials with better capability for separating and storing CO2 based on their pore attributes.

Gautam, Siddharth (ORCID:0000000314435382)↗
DOE OSTI2024

Integrated experimental studies of pore structure and fluid uptake in the Bossier Shale in eastern Texas, USA

Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.

58 GEOSCIENCES↗
DOE OSTI2026

Pore-Scale Study on the Positive Feedback Between Stress and Porosity Caused by Pressure Solution in Porous Media

Pressure solution is an important process in the evolution of sedimentary rocks, which provide storage space for most of our petroleum resources. It directly influences the generation, migration, and storage of petroleum fluids in subsurface sedimentary rocks. Here, in this paper, we develop a pore-scale, mechanochemical model to demonstrate a possible positive feedback between the local porosity and pore surface stress, in which a higher local porosity causes a higher local pore surface stress, thus enhancing pressure solution and consequently further increasing the local porosity. Pore surface stress represents stress on a solid grain adjacent to a pore. Specifically, the pore-scale, mechanochemical model directly simulates the stress distribution over solid and pore surfaces using a finite element model. The dissolution of solids at the solid-pore interfaces under a far-from-equilibrium condition is simulated using a first-order kinetics model that accounts for the local stress distribution. The updated pore geometry, caused by pore surface dissolution, is then used in the stress simulation in the next numerical iteration. Two types of porous media, the Oriskany sandstone and an artificial porous medium with spherical pores, were tested in the mechanochemical simulation. The positive stress-porosity feedback during pressure solution was observed in both samples. In addition, the model quantitatively illustrated the distribution of local mineral dissolution rates on all pore surfaces, as well as its relation to the effective mineral dissolution rate of the entire sample. Based on the comparison between the two porous media, the local mineral dissolution was regulated by pore space distribution, geometry, and coalescence during pressure solution. This work is the first that uses direct, pore-scale numerical simulation to demonstrate the positive stress-porosity feedback during pressure solution, which has the potential to advance the understanding of the mechanical-chemical (MC) coupling in many geological processes that are relevant to subsurface energy systems, such as the recovery of petroleum hydrocarbons and geothermal energy.

CT scanning↗
DOE OSTI2025

Non-equilibrium simulations of hydraulic permeation: Role of mechanical boundary conditions in dense membranes

Understanding the mechanisms of water transport in reverse osmosis membranes is critical for improving membrane performance and guiding material design. While classical models describe transport as either solution-diffusion (SD)—involving concentration-driven diffusion through a homogeneous medium—or pore-flow (PF)—involving pressure-driven convection through percolated water channels—their applicability to crosslinked polyamide membranes remains debated. Here, using non-equilibrium molecular dynamics simulations, we investigate the impact of mechanical support conditions on pressure-driven water transport in polyamide membranes across varying crosslink densities and pressure differentials (1000–5000 bar). Two support conditions are considered: graphene-restrained, representing experimentally relevant supported membranes, and freeze-restrained, mimicking a self-supported structure. In graphene-restrained systems, water concentration gradients and constant pressure profiles emerge, consistent with SD theory and incompatible with PF assumptions due to the absence of percolated pores and sub-nanometer voids. In contrast, freeze-restrained systems display uniform water concentration and linearly decreasing pressure at 1000 bar, and exhibit compressibility-induced water gradients and partial percolation at 5000 bar, resembling PF-like behavior. However, the underlying assumptions of PF theory—continuous solvent pathways and pressure transmission through water-filled pores—are not met under most conditions. Our results demonstrate that accurate modeling of reverse osmosis membranes must incorporate realistic mechanical boundary conditions to distinguish between transport mechanisms. For dense polyamide membranes supported by porous substrates, graphene-restrained simulations best reflect experimental setups and support the SD model as the dominant mechanism of water permeation.

molecular dynamics↗
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
DOE OSTI2024

Molecular‐Level Insight into the Chlorofluorocarbons Adsorption by Defective Covalent Organic Polymers

Abstract Halocarbons have important industrial applications, but because of their contribution to global warming and the fact that they can cause ozone depletion, they are considered highly toxic. Hence, the techniques that can capture and recover the used halocarbons with energy‐efficient methods have been recently received greater attention. In this contribution, we report the capture of dichlorodifluoromethane (R12), which has high global warming and ozone depletion potential, using covalent organic polymers (COPs). The defect‐engineered COPs were synthesized and demonstrated outstanding sorption capacities, ~226 wt % of R12 combined with linear‐shaped adsorption isotherms. We further identified the plausible microscopic adsorption mechanism of the investigated COPs via grand canonical Monte Carlo simulations applied to non‐defective and a collection of atomistic models of the defective COPs. The modeling work suggests that significant R12 adsorption performance is attributed to a gradual increment of porosities due to isolated/interconnected micro‐/meso‐pore channels and the change of the long‐range ordering of both COPs. The successive hierarchical‐pore‐filling mechanism promotes R12 molecular adsorption via moderate van der Waals adsorbate‐adsorbent interactions in the micropores of both COPs at low pressure followed by adsorbate‐adsorbate interactions in the extra‐voids created at moderate to high pressure ranges. This continuous pore‐filling mechanism makes defective COPs as promising sorbents for halocarbon adsorption.

Shen, Jian↗
DOE OSTI2021

Molecular dynamics simulations of energy accommodation between gases and polymers for ultra-low thermal conductivity insulation

Determining the energy accommodation between gases and solids is essential to developing porous thermal insulation materials with ultra-low effective thermal conductivity that reduce energy use, greenhouse gas emissions, and fossil fuel consumption. The energy accommodation coefficients of most gases, however, have been rarely studied, especially with respect to solids that have relatively high thermal resistivity, e.g., polymers. Here, by using all-atom nonequilibrium molecular dynamics simulations, we reveal the accommodation coefficients of He, Ar, N 2 , and O 2 with polymers, mainly polystyrene. We find that their values are around 0.51, 0.72, 0.79, and 0.90, respectively, suggesting a critical reexamination of the commonly used theoretical maximum value of 1. We have also conducted experiments and validated the value for air, which is about 0.81. Such a change in accommodation coefficients can lead to a reduction of about 70%, 50%, 35%, and 20% in the thermal conductivity of He, Ar, N 2 , and O 2 gases in nano pores (below 100 nm) or at low pressures (below 1 millibar). With these new accommodation coefficients, we find that in a 10 nm pore with ambient pressure at 300 K, the gas thermal conductivity of He, Ar, N 2 , and O 2 in porous polystyrene can be as low as 9.7 × 10 –4 , 3.4 × 10 –4 , 7.3 × 10 –4 , and 8.5 × 10 –4 W·m –1 ·K –1 , respectively, which are two to three orders of magnitude lower than their bulk values, promising higher thermal resistivity of insulation materials. This work reveals the fundamental energy exchange between gases and polymers, providing important guidance for designing high-performance thermal insulation materials for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structure and dynamics of ethane confined in silica nanopores in the presence of CO 2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C 2 H 6 ) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C 2 H 6 and CO 2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C 2 H 6 , C 2 H 6 :CO 2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C 2 H 6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C 2 H 6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO 2 molecules are more strongly attracted to the pore surface compared to C 2 H 6 . The C 2 H 6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO 2 molecules play in enhancing C 2 H 6 diffusion (“molecular lubrication”) at high CO 2 :C 2 H 6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗