Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “capillary”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermodynamic Analysis of Nano-Electrospray Induced Gas Jets

Nano-electrospray (nES) produces a plume of charged liquid droplets which have drag interactions with the surrounding gas after they are emitted from a source capillary. The resulting induced gas flow has often been considered unimportant, but we have demonstrated that a gas microjet of significant velocity up to 10s of m/s can be produced. In this work we present a thermodynamic framework that enables analysis of gas jet generation from electrosprays and introduce the important metrics for such analysis, effectiveness (a measure of momentum transfer from the electrosprayed aerosol to gas) and efficiency (a measure of energy conversion from electrical energy generating the electrospray to gas kinetic energy). This analytical framework is applicable to any two-phase flows consisting of discrete conservative-force-driven particles which exchange momentum with an inert, otherwise quiescent fluid medium to yield a co-flowing two-phase jet. Here, we apply this framework to sprays of water from nano-electrospray emitters to demonstrate that increasing the applied electrical potential difference, increasing liquid mass flowrate, and decreasing droplet size all can increase electrospray induced gas jet strength, but only the latter two do so while also increasing the momentum transfer effectiveness and energy conversion efficiency.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Improved lipid analysis using a 2D-LC-MS system with a novel injection procedure

The aim of this study was to improve analysis of nonpolar lipidomics sample extracts using reversed phase (RP) chromatography. A 4/3/3 (v/v/v) mixture of methanol/methyl tert-butyl ether/chloroform (MeOH/MTBE/CHCl 3 , MMC) was chosen for sample extraction solvent based on its proven extraction capability for several lipid classes. To avoid carry over, loss of analytes and peak distortion the loops and all capillaries of the presented LC system were flushed and filled up with methanol until the analytical column. Here, the choice of methanol was due to its weak elution strength and being infinitely miscible with MMC and several other nonpolar solvents. This allowed injection of a 100 μl sample that was 20 μl nonpolar extraction solvent diluted fivefold with methanol. All lipids of 25 lipid classes were transferred quantitatively to the column head where the online dilution of methanol was carried out with aqueous eluent for focusing the lipid analytes. The weak elution strength of methanol prevented peak distortions. The consecutive reversed phase elution resulted in remarkably narrow peaks (full width at half maximum was 0.07–0.08 min typically) and enhanced sensitivity (limit of detection usually in sub nM region) because of increased sample injection volume and narrow peaks. Calibration and quality control samples made by diluting commercial lipid standards 200–50000 times confirmed the applicability of this approach both for targeted lipid quantification and for untargeted quantitative comparison of lipids from different sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Transparent silica aerogel slabs synthesized from nanoparticle colloidal suspensions at near ambient conditions on omniphobic liquid substrates

This paper presents a novel sol-gel method to synthesize large and thick silica aerogel monoliths at near ambient conditions using a commercial aqueous solution of colloidal silica nanoparticles as building blocks. To achieve slabs with high visible transmittance and low thermal conductivity, the method combines the strategies of (i) synthesizing gels on an omniphobic perfluorocarbon liquid substrate, (ii) aging at temperatures above room temperature, and (iii) performing solvent exchange with a low-surface-tension organic solvent prior to ambient drying. The omniphobic liquid substrates were used to prevent cracking and ensure an optically-smooth surface, while nanoparticle building blocks were small (<10 nm) to limit volumetric light scattering. Gels were aged at temperatures between 25 and 80 °C for up to 21 days to make them stronger and stiffer and to reduce shrinkage and cracking during ambient drying. Ambient drying was achieved by first exchanging water in the gel pores for octane, followed by drying in an octane-rich atmosphere to decrease capillary forces. The synthesized nanoparticle-based silica aerogel monoliths had thicknesses up to 5 mm, diameters up to 10 cm, porosities exceeding 80%, and thermal conductivities as low as 0.08 W m –1 K –1 . Notably, the slabs featured visible transmittance exceeding 75% even for slabs as thick as 5 mm. The as-synthesized aerogel monoliths were exposed to TMCS vapor to induce hydrophobic properties resulting in a water contact angle of 140° that prevented water infiltration into the pores and protected the aerogels from water damage. This simple synthesis route conducted at near ambient conditions produces hydrophobic aerogel monoliths with promising optically transparent and thermally insulating properties that can be adhered to glass panes for window insulation and solar-thermal energy conversion applications.

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

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2024

Nanoconfinement-induced shift in photooxidative degradation pathway of polystyrene

Polymer nanocomposites with high concentrations of nanoparticles (NPs) possess exceptional mechanical, transport, and thermal properties. To enable their widespread use in structural applications and functional coatings, it is crucial to understand how nanoconfinement and the polymer-NP interface influence polymer degradation under various environmental conditions, including prolonged UV exposure. In this study, we investigate the photooxidative degradation of polystyrene (PS)-confined in the interstices of SiO 2 NP films. These nanocomposite films are prepared by the capillary rise infiltration (CaRI) of PS into interstices of SiO 2 NP packings, and subsequently subjected to UV irradiation. Our investigation reveals that PS degradation progresses uniformly across the thickness, with degradation initiating from the center of the NP interstitial pores and extending towards the NP surface. We rationalize this degradation mechanism based on the disparity in the surface energies of PS and the NP surface, as well as slow oxygen diffusion through confined PS. Here, we also demonstrate that packing smaller NPs at a given thickness or thicker packing of a specific NP size facilitates photooxidative degradation, highlighting the critical role of the number of interstitial pores in influencing degradation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Development of a high-temperature inorganic synthetic foam with recycled fly-ash cenospheres for thermal insulation brick manufacturing

In this study, a novel inorganic synthetic solid foam composed of recycled fly-ash cenospheres (FACs) is developed with low-density, good thermal insulation performance, and reasonable mechanical strength. The manufacturing process utilizes low-fraction clay binder materials diluted in water, enabling a self-assembling mechanism assisted by natural capillary effect to form minimal contacts among the FACs. The microstructure and minimal contact formation were confirmed by optical and electron microscope imaging. Both thermal and mechanical properties of the fly-ash cenosphere synthetic solid foam (FAC-SSF) materials were experimentally studied as functions of key material parameters including cenosphere particle size and clay content. A numerical model was developed to predict the effective thermal conductivity of FAC-SSF using a three-step homogenization method which takes into account different thermal percolation paths within the material. The inorganic synthetic foam developed has the potential for thermal insulation applications taking advantage of their relatively simple, low-cost, and scalable manufacturing process, minimal use of binder materials, and relatively high mechanical strength and highly tailorable properties.

42 ENGINEERING↗
DOE OSTI2023

Numerical investigation and experimental validation of water condensation in the gas diffusion layer with different properties

Liquid water in porous media hinders the transport of reactant to catalyst layer, where the electrochemical reactions occur, which in turn, reduces performance and lifetime of a proton exchange membrane fuel cell (PEMFC). In this study, two types of gas diffusion layer (GDL) materials are used to study the two-phase saturation distribution in a PEMFC. The correlations between the effective and anisotropic transport properties and porosity of the two GDLs are solved by the pore scale model. In addition, the anisotropic liquid water permeability and the relationship between saturation pressure and capillary pressure is determined by using the Lattice-Boltzmann method. Fuel cell performance as well as liquid water distribution in the GDL using neutron radiography are obtained for validation. Under wet conditions, the cell with Freudenberg GDL performs better than that using Toray GDL, especially at high current densities. The results of the two-phase model simulation show that the peak water saturation for Toray GDL occurs in the catalyst layer and can reach 40%, while the peak water saturation for Freudenberg GDL occurs inside the GDL directly under the ribs and only reaches 25% saturation. The combined results provide key insights to enable high power density operation of a PEMFC through GDL material optimization.

42 ENGINEERING↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗
DOE OSTI2020

High accuracy in-situ direct gas analysis of Li-ion batteries

Cycling of lithium-ion batteries containing Ni-rich NMC cathodes at high voltage involves intense gas generation. From the safety standpoint, it is critical to understand how different gas species respond to changes of upper cut-off voltages. In this manuscript, we introduce a novel experimental set up for real-time analysis of gas generation in prismatic pouch cells. In a typical experiment, a lithium-ion pouch cell is directly connected to a quadruple mass spectrometer by glass capillary. Pressure difference helps move the generated gases to the mass spectrometer column for full analysis. The gaseous species are probed during both formation cycle and aging cycles in Li-ion pouch cells comprising NMC811 cathode and graphite anode. The gases are generated upon the creation of the solid electrolyte interphase on graphite during the first formation charge. Besides ethylene gas, C 2 H 4 , whose occurrence is related to the decomposition of ethylene carbonate, a minor solvent decomposition happened at a lower voltage is demonstrated to be associated with residue water on the electrode. For aging cycles, three upper cut-off voltages have been selected to investigate the gas evolution behavior upon charging to 4.2, 4.4 and 4.6 V. Higher upper cut-off voltages do not affect the amount of O 2 measured as a part of the generated oxygen could have reacted with electrolytes. Quite similar and low amounts of O 2 remain in the cells upon charging to the selected cut-off voltages. On the other hand, the generation of CO 2 has been found to be very sensitive to upper cut-off voltage. During similar aging cycles, 6167 nmol $g_{NMC}^{-1}$ of CO 2 is generated at 4.6 V, versus 1650 nmol $g_{NMC}^{-1}$ of CO 2 at 4.4 V, versus 91 nmol $g_{NMC}^{-1}$ of CO 2 at 4.2 V.

25 ENERGY STORAGE↗
DOE OSTI2021

Multi-fidelity electrochemical modeling of thermally activated battery cells

Thermally activated batteries undergo a series of coupled physical changes during activation that influence battery performance. These processes include energetic material burning, heat transfer, electrolyte phase change, capillary-driven two-phase porous flow, ion transport, electrochemical reactions, and electrical transport. Several of these processes are strongly coupled and have a significant effect on battery performance, but others have minimal impact or may be suitably represented by reduced-order models. Additionally, assessing the relative importance of these phenomena must be based on comparisons to a high-fidelity model including all known processes. In this work, we first present and demonstrate a high-fidelity, multi-physics model of electrochemical performance. This novel multi-physics model enables predictions of how competing physical processes affect battery performance and provides unique insights into the difficult-to-measure processes that happen during battery activation. We introduce four categories of model fidelity that include different physical simplifications, assumptions, and reduced-order models to decouple or remove costly elements of the simulation. Using this approach, we show an order-of-magnitude reduction in computational cost while preserving all design-relevant quantities of interest within 5 percent. The validity of this approach and these model reductions is demonstrated by comparison between results from the full fidelity model and the different reduced models.

25 ENERGY STORAGE↗
DOE OSTI2021

Non-planar platinum group metal-free fuel cell cathodes for enhanced oxygen transport and water rejection

Proton exchange membrane fuel cells (PEMFC) with cathodes using platinum group metal-free (PGM-free) catalysts could significantly reduce costs, but the lower volumetric oxygen reduction reaction (ORR) activity requires thick electrodes that suffer from liquid water flooding and increased oxygen transport resistance. To address these challenges, we developed a 3D gas diffusion electrode (GDE) architecture to enhance liquid water removal through the diffusion media and reduce cathode saturation. The cathode features a uniform catalyst layer adjacent to the membrane for high ORR activity and then pillars of the hydrophilic catalyst layer that pass through the microporous layer (MPL), providing a low capillary pressure barrier pathway to the carbon fiber paper layer and channel. The non-planar cathode also increases the interfacial area between the catalyst layer and hydrophobic MPL for a greater fraction of the cathode with high O 2 concentration. Our studies included parametric experimental study of the pillar density to identify the optimum pitch between pillars. Our measurements show significant improvements in the mass transport region of the polarization curve with 3D structured electrodes leading to an 8% increase in maximum current density, 19% increase in maximum power density, and 16% increase in current density at 0.67 V with air.

25 ENERGY STORAGE↗
DOE OSTI2022

Multifunctional Superelastic Graphene Aerogels Derived from Ambient-Dried Graphene Oxide/Camphene Emulsions

Graphene oxide (GO) emulsions have been widely adopted to fabricate lightweight graphene aerogels (GAs). However, preparing high-performance GAs from conventional GO emulsions under ambient conditions seems rather challenging because of the structural collapse caused by strong capillary forces during drying. Here, we report a new, low-cost route to prepare high-performance, multifunctional GAs by the ambient drying of GO/camphene emulsions. Benefiting from the liquid-solid and solid-gas phase transitions of camphene, highly porous GO aerogels were obtained under ambient conditions. Moreover, the thermally annealed GO aerogels exhibited excellent mechanical properties, fire resistance, thermal insulation performance, and oil adsorption capability.

ambient drying↗
DOE OSTI2019

Multimode jetting unlocks a trade-off between nanostructure morphology and composition in focused electron beam induced deposition

An ability to accelerate the growth rate while simultaneously controlling shape and composition of deposits is fundamental to expanding the utility of electron beam assisted deposition for “direct write” nanostructure fabrication from gas phase precursors. Here, we report on novel tunable mode of Focused Electron Beam Induced Deposition (FEBID) as enabled by the switchable use of reactive (O2) and inert (Ar) “assist” gases directed at high impingement energy and momentum to a deposition site concurrently with a supply of a precursor in a dual-nozzle FEBID arrangement. Organometallic molecules introduced into a deposition chamber in a free-molecular flow regime through a conventional gas injection nozzle, yielding a broad uniform substrate coverage by the precursor in a quasi-equilibrium adsorbed state. Simultaneously, a supersonic micro-jet from a micron-scale capillary placed in close proximity to a deposition site produces a low Knudsen number continuous flow impingement of energized “assist” gas, which perturbs the precursor state to facilitate different chemical reaction pathways. This process offers modulation of deposition outcomes, e.g., enhanced vertical growth rate (aspect ratio) of nanostructures when using argon micro-jet while locally removing parasitic carbonaceous by-products of organometallic precursor decomposition reaction by e-beam through interactions with oxygen micro-jet. Such in-situ/in operando (during growth and on-demand switchable) modulation of FEBID allows for rapid and precisely controllable deposition for a desired outcome, such as tuning metal-to-carbon composition and maintaining shape fidelity in fabricating high-aspect ratio composite nanostructures, which extends the capabilities of FEBID as a robust multi-dimensional, multi-purpose nanoscale 3D printing tool.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2026

Overview of advanced plasma-facing materials testing for Fusion Pilot Plants at DIII-D

Characterization and testing of advanced plasma-facing materials (PFMs) for Fusion Pilot Plants (FPP) is being conducted at the DIII-D National Fusion Facility through the ongoing two-year FPP Candidate Materials Thrust. Year one tested 17 novel materials utilizing the Divertor Materials Evaluation System (DiMES), with samples analyzed pre- and post-experiment via SEM, EDS, and confocal microscopy. Repeatable reference discharges were developed to ensure uniformity between experiments, including a new strike-point rastering scenario to provide more uniform heat/particle flux across DiMES during ELMing H-mode discharges. Various sample geometries and temperatures were used to achieve FPP-relevant conditions, including samples angled 10° towards the incident plasma flux and pre-heating up to 500 °C. The first exposure of liquid lithium (Li) capillary porous structures in a tokamak demonstrated uniform emission of Li vapor and suppression of Li droplets in H-mode when preheated to 350 °C. Dispersoid-strengthened W with 1 wt% TaC, TiC, and ZrC exposed to H-mode showed cracking and dispersoid ejection for all varieties except TiC, providing a clear down-selection. Ultra-high temperature ceramic materials TiB 2 and ZrB 2 showed minimal degradation under L-mode exposure. Silicon carbide (SiC) fiber composites showed arcing along edges, while CVD SiC remained pristine. Atmospheric plasma-sprayed W and SiC coatings endured H-mode exposure without macroscopic delamination; SiC exhibited granular ejection, while W showed increased outgassing. Additional W-based alloys were stress tested in H-mode, including Ni-based W heavy alloys, W f SiC f /W composites, W multi-principle element alloys, and functionally-graded W/SiC, to varying degrees of success.

DIII-D↗
DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗
DOE OSTI2025

Comparative study between experimental measurements and model predictions of bubble rise velocity in molten LiCl-KCl

This research investigated the shape and rise velocities of gas bubbles (helium, nitrogen, argon, and krypton) in molten LiCl-KCl at 500 °C using a specialized optical visualization apparatus. The apparatus features a custom quartz rectangular prism cell housed in a furnace, enabling the application of the shadowgraph measurement technique. Bubbles with diameters between 0.45 and 3.2 mm were produced via a capillary, and their equivalent diameters and rise velocities were analyzed using high-resolution imaging and centroid tracking algorithms. Important non-dimensional parameters were calculated to characterize behavior. Bubbles with 0.03 < Eo < 1.32, 17 < Re < 718, and Mo order of 10 -11 were generated in the molten LiCl-KCl, revealing that bubbles with diameters less than 1.46 mm remain spherical, whereas those with larger diameters are ellipsoidal. The transition and oscillation onset for bubble shape occurred at larger diameters than those previously observed in air–water systems. Comparative analysis of the gas-molten salt data against existing air–water bubble rise velocity correlations highlighted inaccuracies in predictions, particularly across the different bubble regimes. An existing correlation was calibrated to the molten salt data to improve the bubble velocity predictions. In conclusion, this study contributes essential experimental data for the design and safety evaluation of molten salt reactors and offers insights for the optimization of sparge systems for fission product removal.

Bubble rise velocity↗
DOE OSTI2020

Estimating diffusion coefficients of shale oil, gas, and condensate with nano-confinement effect

Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. This study aims at calculating diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO 2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This paper estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 10 2 to 10 4 times as porosity decreases from 0.1 to 0.03.

04 OIL SHALES AND TAR SANDS↗