Search NASA⌕ Search

SEARCH · Search NASA

Results for “SLIP FLOW”

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

Effect of CO 2 on the water slip flow at silica surfaces for nanometer slit pores of talc

In our previous study, Molecular Dynamics (MD) simulated water flow with no dissolved gas in nanometer slit pores of the hydrophobic talc surfaces revealed a slip length of about 0.5 nm, which is close to the size of “water exclusion zone” but much less than experimental slip length. Atomic Force Microscopy (AFM) imaging confirmed the presence of pancake shape nanobubbles at the hydrophobic talc (001) surface (size from tens to hundreds of nanometers). The effect of CO 2 on water slip flow in talc nanopores is reported to reveal the importance of dissolved gas for confined fluid flow. MD simulated 7 nm CO 2 nanobubbles were found to attach and spread at the talc (001) surface, consistent with the AFM observation. Simulated CO 2 saturated water flow in the slit pores of talc (001) surfaces predicted the generation of CO 2 nanobubbles at the end of the slit pore, due to the stabilization of a CO 2 film at the talc (001) surfaces and a water barrier at the edge surfaces of talc. A “critical thickness” of about 1.5 nm was found for CO 2 nanobubbles to be stabilized at the end of the talc (001) surface. With 1.5 nm CO 2 nanobubbles at the talc (001) surface of a 6 nm slit pore, a water slip length of 1.8 nm was determined. Furthermore, the nanometer slit pore simulations with dissolved CO 2 in water revealed an increased slip length. It is expected that an increased slip length value would be found with an increased CO 2 film thickness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation and analysis of slip flow of water at hydrophobic silica surfaces of nanometer slit pores

Water flow in nanopores is controlled fundamentally by pore wettability and interfacial water features. In the research herein reported, Molecular dynamics (MD) was used to simulate water flow in nanometer slit pores of silica surfaces with variable wetting characteristics. Compared to Couette flow, the simulated slit-pores of the talc (001) and octadecyltrichlorosilane (OTS) monolayer surfaces exhibited water flow enhancement and a positive slip length at the interface, whereas water flow velocity at the quartz (001) and amorphous silica surfaces was zero, adhering to the no slip condition. More intensive interactions between the interfacial water and hydrophilic silica surfaces were confirmed by MD simulation and sum-frequency vibrational spectroscopy (SFVS) results. A water exclusion zone of approximately 3 Å was found at the hydrophobic silica surfaces by quantitative MD interfacial water analysis and validated by SFVS. This exclusion zone dimension was similar in magnitude to the slip length (4–5 Å) as determined from simulated water flow in hydrophobic silica slit pores. The simulated slip lengths at the hydrophobic surfaces (0.2–1 nm) were significantly less than the experimental slip lengths (>40 nm) reported in the literature. In the case of the experimental results, this difference appears to be due to nanobubbles expected to be present at hydrophobic surfaces.

42 ENGINEERING↗

Percolation connectivity, pore size, and gas apparent permeability: Network simulations and comparison to experimental data: Connectivity, Pore Size, and Gas Permeability

We modeled single-phase gas flow through porous media using percolation networks. Gas permeability is different from liquid permeability. The latter is only related to the geometry and topology of the pore space, while the former depends on the specific gas considered and varies with gas pressure. As gas pressure decreases, four flow regimes can be distinguished as viscous flow, slip flow, transition flow, and free molecular diffusion. Here we use a published conductance model presumably capable of predicting the flow rate of an arbitrary gas through a cylindrical pipe in the four regimes. We incorporated this model into pipe network simulations. We considered 3-D simple cubic, body-centered cubic, and face-centered cubic lattices, in which we varied the pipe radius distribution and the bond coordination number. Gas flow was simulated at different gas pressures. The simulation results showed that the gas apparent permeability k app obeys an identical scaling law in all three lattices, k app ~ (z-z c ) β , where the exponent β depends on the width of the pipe radius distribution, z is the mean coordination number, and z c its critical value at the percolation threshold. Surprisingly, (z-z c ) had a very weak effect on the ratio of the apparent gas permeability to the absolute liquid permeability, k app /k abs , suggesting that the Klinkenberg gas slippage correction factor is nearly independent of connectivity. We constructed models of k app and k app /k abs based on the observed power law and tested them by comparison with published experimental data on glass beads and other materials.

58 GEOSCIENCES↗

Use of finite volume radiation for predicting the Knudsen minimum in 2D channel flow

In an earlier paper we employed an analogy between surface-to-surface radiation and free-molecular flow to model Knudsen flow through tubes and onto planes. In the current paper we extend the analogy between thermal radiation and molecular flow to model the flow of a gas in a 2D channel across all regimes of rarefaction. To accomplish this, we break down the problem of gaseous flow into three sub-problems (self-diffusion, mass-motion and generation of pressure gradient) and use the finite volume method for modeling radiation through participating media to model the transport in each sub-problem as a radiation problem. We first model molecular self-diffusion in the stationary gas by modeling the transport of the molecular number density through the gas starting from the analytical asymptote for free-molecular flow to the kinetic theory limit of gaseous self-diffusion. We then model the transport of momentum through the gas at unit pressure gradient to predict Poiseuille flow and slip flow in the 2D gas. Lastly, we predict the generation of pressure gradient within the gas due to molecular collisions by modeling the transport of the forces generated due to collisions per unit volume of gas. We then proceed to combine the three radiation problems to predict flow of the gas over the entire Knudsen number regime from free-molecular to transition to continuum flow and successfully capture the Knudsen minimum at Kn ∼ 1.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-dimensional lift-up problem for a rigid porous bed

The present study analytically reinvestigates the two-dimensional lift-up problem for a rigid porous bed that was studied by Mei, Yeung, and Liu [“Lifting of a large object from a porous seabed,” J. Fluid Mech. 152, 203 (1985)]. Mei, Yeung, and Liu proposed a model that treats the bed as a rigid porous medium and performed relevant experiments. In their model, they assumed the gap flow comes from the periphery of the gap, and there is a shear layer in the porous medium; the flow in the gap is described by adhesion approximation [D. J. Acheson, Elementary Fluid Dynamics (Clarendon, Oxford, 1990), pp. 243-245.] and the pore flow by Darcy’s law, and the slip-flow condition proposed by Beavers and Joseph [“Boundary conditions at a naturally permeable wall,” J. Fluid Mech. 30, 197 (1967)] is applied to the bed interface. In this problem, however, the gap flow initially mainly comes from the porous bed, and the shear layer may not exist. Although later the shear effect becomes important, the empirical slip-flow condition might not physically respond to the shear effect, and the existence of the vertical velocity affects the situation so greatly that the slip-flow condition might not be appropriate. In contrast, the present study proposes a more general model for the problem, applying Stokes flow to the gap, the Brinkman equation to the porous medium, and Song and Huang’s [“Laminar poroelastic media flow,” J. Eng. Mech. 126, 358 (2000)] complete interfacial conditions to the bed interface. The exact solution to the problem is found and fits Mei’s experiments well. The breakout phenomenon is examined for different soil beds, mechanics that cannot be illustrated by Mei’s model are revealed, and the theoretical breakout times obtained using Mei’s model and our model are compared. The results show that the proposed model is more compatible with physics and provides results that are more precise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular dynamics simulation of methane transport in confined organic nanopores with high relative roughness

Understanding and characterizing the transport of shale gas (methane) through the nanopores of kerogens are critical for the accurate prediction of shale gas recovery. However, the key factors that regulate shale gas transport through highly roughened nanopores of shale kerogens are not fully understood. Here in this work, methane transport in organic nanopores with a high relative roughness is characterized using equilibrium and non-equilibrium molecular dynamics methods. According to our results, the CH 4 mass flux has a linear relationship with the pressure gradient, consistent with previous studies, while the calculated slip lengths and gas fluxes varied with different roughness geometries in the order of sigmoidal ≥ triangular > rectangular. Surface slip flow can be a major contributor to the overall gas flux, but surprisingly, the relative contribution of surface slip flow is independent of the pressure gradient. In contrast, the contributions of both slip flow and the average gas fluxes vary strongly with pore diameters. Typical contributions of the adsorbed layer to the overall gas flux are in the range 20–40% but vary from as high as 74% in a 4-nm pore to as low as 6% in a 16-nm pore. Compared to smooth nanopores, we find that, in nanopores with realistically high degrees of relative roughness, methane confinement in cavities decouples slip flow from the flow in the pore interior, significantly reducing the overall flux.

03 NATURAL GAS↗

Predicting Fluid Flow Regime, Permeability, and Diffusivity in Mudrocks from Multiscale Pore Characterisation

In geoenergy applications, mudrocks prevent fluids to leak from temporary (H 2 , CH 4 ) or permanent (CO 2 , radioactive waste) storage/disposal sites and serve as a source and reservoir for unconventional oil and gas. Understanding transport properties integrated with dominant fluid flow mechanisms in mudrocks is essential to better predict the performance of mudrocks within these applications. In this study, small-angle neutron scattering (SANS) experiments were conducted on 71 samples from 13 different sets of mudrocks across the globe to capture the pore structure of nearly the full pore size spectrum (2 nm–5 μm). We develop fractal models to predict transport properties (permeability and diffusivity) based on the SANS-derived pore size distributions. The results indicate that transport phenomena in mudrocks are intrinsically pore size-dependent. Depending on hydrostatic pore pressures, transition flow develops in micropores, slip flow in meso- and macropores, and continuum flow in larger macropores. Fluid flow regimes progress towards larger pore sizes during reservoir depletion or smaller pore sizes during fluid storage, so when pressure is decreased or increased, respectively. Capturing the heterogeneity of mudrocks by considering fractal dimension and tortuosity fractal dimension for defined pore size ranges, fractal models integrate apparent permeability with slip flow, Darcy permeability with continuum flow, and gas diffusivity with diffusion flow in the matrix. This new model of pore size-dependent transport and integrated transport properties using fractal models yields a systematic approach that can also inform multiscale multi-physics models to better understand fluid flow and transport phenomena in mudrocks on the reservoir and basin scale.

36 MATERIALS SCIENCE↗

Klinkenberg Effect and Effective Pressure for Gas Permeability of Tight Sandstones

Gas permeability k g of tight sandstones was measured while cycling pore pressure, Pf, and confining pressure, P c , using two different procedures. Firstly, we cycled the pore pressure from a minimum value close to atmospheric pressure to a maximum one slightly lower than P c , while keeping the Pc constant. Secondly, P f and P c were simultaneously varied so as to keep the difference, P c - P f , constant. Significant effects owing to slip flow were observed when the gas pressure was lower than about 5 MPa. Klinkenberg plots showed that, under constant P c , gas permeability kg decreased linearly for high values of 1/P f , reached a minimum and then sharply increased as 1/P f approached zero. Under constant P c - P f , k g again decreased linearly in the high limit of 1/P f and showed additional decrease in the low 1/P f limit. The linear segments in both procedures indicate that k g in tight sandstones was subject to slip flow at low values of P f , while the curved segments revealed the sensitivity of permeability to effective pressure when Pf became comparable in magnitude to confining pressure. The effective pressure coefficient α was estimated and generally appeared to be a well-defined decreasing function of the pressure difference P c - P f . We modified the Walsh permeability model for fissured rocks to include slip-flow effects and obtained satisfactory fits to the experimental data for the samples with lowest porosities and highest pressure-sensitivities. This observation suggests that the intrinsic permeability of these rocks was controlled by inter- and/or intra-granular crack-like pores.

58 GEOSCIENCES↗

Shear-stress-induced structural arrangement of water molecules in nanoscale Couette flow with slipping at wall boundary

This study investigated the structuring of water molecules in a nanoscale Couette flow with the upper plate subjected to lateral forces with various magnitudes and water slipping against a metal wall. It was found that when the upper plate is subjected to a force, the water body deforms into a parallelepiped. Water molecules in the channel are then gradually arranged into lattice positions, creating a layered structure. The structural arrangement of water molecules is caused by the water molecules accommodating themselves to the increase in energy under the application of a lateral force on the moving plate. The ordering arrangement of water molecules increases the rotational degree of freedom, allowing the molecules to increase their Coulomb potential energy through polar rotation that accounts for the energy input through the upper plate. With a force continuously applied to the upper plate, the water molecules in contact with the upper plate move forward until slip between the water and upper plate occurs. The relation between the structural arrangement of water molecules, slip at the wall, and the shear force is studied. The relation between the slip and the locking/unlocking of water molecules to metal atoms is also studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Velocity Measurements in Channel Gas Flows in the Slip Regime by means of Molecular Tagging Velocimetry

Direct measurements of the slip velocity in rarefied gas flows produced by local thermodynamic non-equilibrium at the wall represent crucial information for the validation of existing theoretical and numerical models. In this work, molecular tagging velocimetry (MTV) by direct phosphorescence is applied to argon and helium flows at low pressures in a 1-mm deep channel. MTV has provided accurate measurements of the molecular displacement of the gas at average pressures of the order of 1 kPa. To the best of our knowledge, this work reports the very first flow visualizations of a gas in a confined domain and in the slip flow regime, with Knudsen numbers up to 0.014. MTV is cross-validated with mass flowrate measurements by the constant volume technique. The two diagnostic methods are applied simultaneously, and the measurements in terms of average velocity at the test section are in good agreement. Moreover, preliminary results of the slip velocity at the wall are computed from the MTV data by means of a reconstruction method

74 ATOMIC AND MOLECULAR PHYSICS↗

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗

The influence of interfacial slip on two-phase flow in rough pores

The migration and trapping of supercritical CO 2 (scCO 2 ) in geologic carbon storage is strongly dependent on the geometry and wettability of the pore network in the reservoir rock. During displacement, resident fluids may become trapped in the pits of a rough pore surface forming an immiscible two-phase fluid interface with the invading fluid, allowing apparent slip flow at this interface. We present a two-phase fluid dynamics model, including interfacial tension, to characterize the impact of mineral surface roughness on this slip flow. We show that the slip flow can be cast in more familiar terms as a contact-angle (wettability)-dependent effective permeability to the invading fluid, a nondimensional measurement which relates the interfacial slip to the pore geometry. The analysis shows the surface roughness-induced slip flow can effectively increase or decrease this effective permeability, depending on the wettability and roughness of the mineral surfaces. Configurations of the pore geometry where interfacial slip has a tangible influence on permeability have been identified. The results suggest that for large roughness features, permeability to CO 2 may be enhanced by approximately 30% during drainage, while the permeability to brine during reimbibition may be enhanced or diminished by 60%, depending on the contact angle with the mineral surfaces and degrees of roughness. For smaller roughness features, the changes in permeability through interfacial slip are small. As a result, a much larger range of effective permeabilities are suggested for general fluid pairs and contact angles, including occlusion of the pore by the trapped phase.

54 ENVIRONMENTAL SCIENCES↗

A study on the effect of various design parameters on the natural circulation flow rate of the ex-vessel core catcher cooling system of EU-APR1400

In this paper, a study on the effect of various design parameters such as the channel gap width, heat flux distribution, down-comer pipe size and two-phase flow slip ratio on the natural circulation flow rate is performed based on a physical model for a natural circulation flow along the flow path of the ex-vessel core catcher cooling system of an EU-APR1400, and these effects on the natural circulation flow rate are analyzed and compared with the minimum flow rate required for the safe operation of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Strong Electroosmotic Coupling Dominates Ion Conductance of 1.5 nm Diameter Carbon Nanotube Porins

Extreme confinement in nanometer sized channels can alter fluid and ion transport in significant ways, leading to significant water flow enhancement and unusual ion correlation effects. These effects are especially pronounced in carbon nanotube porins (CNTPs) that combine strong confinement in the inner lumen of carbon nanotubes with the high slip flow enhancement due to smooth hydrophobic pore walls. For this paper, we have studied ion transport and ion selectivity in 1.5 nm diameter CNTPs embedded in lipid membranes using single nanopore measurement setup. Our data show that CNTPs are weakly cation-selective at pH 7.5 and become non-selective at pH 3.0. Ion conductance of CNTPs exhibits an unusual 2/3 power law scaling with the ion concentration at both neutral and acidic pH values. Coupled Navier-Stokes and Poisson-Nernst-Planck (NS-PNP) simulations and atomistic Molecular Dynamics (MD) simulations reveal that this scaling originates from strong coupling between water and ion transport in these channels. These effects could result in development of a next generation of biomimetic membranes and carbon nanotube-based electroosmotic pumps.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fault Slip and Fluid Flow: Seismic Source Analysis to Assess Role of Multiple Slip Patches in Fault Permeability

The relationship between fault reactivation, microearthquakes (MEQs), and permeability evolution during fluid injection plays a critical role in energy harvesting and waste disposal. Recent studies have demonstrated the possibility of predicting fault permeability using cumulative seismic moments of MEQs quantitatively. To understand the underlying physical processes, we conduct fault reactivation experiments using Utah FORGE granitoid and analyze acoustic emission (AE) signals generated during stepwise increases in fluid injection pressure. Frequency analysis of thousands of calibrated AE signals reveals that fault reactivation produces multiple AE source patches with millimeter-scale radii—smaller than the sample fault radius. The cumulative area of the reactivated patches covers the fault multiple times over (∼10x–50x area) for each pressure step. These findings provide mechanistic insight that measured permeability enhancement is not driven by a single large slip event, but by the sequential and interacting activation of multiple slip patches that create a continuous flow pathway.

Nurshal, M. E. M. [Pennsylvania State University, ↗

Concurrent coupling of atomistic simulation and mesoscopic hydrodynamics for flows over soft multi-functional surfaces

We develop an e?cient parallel multiscale method that bridges the atomistic and mesoscale regimes, from nanometers to microns and beyond, via concurrent coupling of atomistic simulation and mesoscopic dynamics. In particular, we combine an all-atom molecular dynamics (MD) description for specific atomistic details in the vicinity of the functional surface with a dissipative particle dynamics (DPD) approach that captures mesoscopic hydrodynamics in the domain away from the functional surface. In order to achieve a seamless transition in dynamic properties we endow the MD simulation with a DPD thermostat, which is validated against experimental results by modeling water at di?erent temperatures. We then validate the MD-DPD coupling method for transient Couette and Poiseuille flows, demonstrating that the concurrent MD-DPD coupling can resolve accurately the continuumbased analytical solutions. Subsequently, we simulate shear flows over grafted polydimethylsiloxane (PDMS) surfaces (polymer brushes) for various grafting densities, and investigate the slip flow as a function of the shear stress. We verify that a ‘‘universal’’ power law exists for the slip length, in agreement with published results. Having validated the MD-DPD coupling method, we simulate timedependent flows past an endothelial glycocalyx layer (EGL) in a microchannel. Coupled simulation results elucidate the dynamics of the EGL changing from an equilibrium state to a compressed state under shear by aligning the molecular structures along the shear direction. MD-DPD simulation results agree well with results of a single MD simulation, but with the former more than two orders of magnitude faster than the latter for system sizes above one micron.

Wang, Yuying↗

Gas sorption and non-Darcy flow in shale reservoirs

Gas sorption and non-Darcy flow are two important issues for shale gas reservoirs. The sorption consists of dissolution and adsorption. Dissolved gas and adsorbed gas are different. The former is dissolved in the shale matrix, while the latter is concentrated near the solid walls of pores. In this paper, the Langmuir equation is used to describe adsorption and Henry’s law is used to describe dissolution. The K coefficient in Henry’s law of 0.052 mmol/(MPa g TOC) is obtained by matching experimental data. The amount of dissolved gas increases linearly when pressure increases. Using only the Langmuir equation without considering dissolution can lead to a significant underestimation of the amount of sorbed gas in shales. For non-Darcy gas flow, the apparent permeability model for free gas is established by combining slip flow and Knudsen flow. For adsorbed gas, the surface diffusion effect is also considered in this model. Here, the surface diffusion coefficient is suggested to be of the same scale as the gas self-diffusion coefficient, and the corresponding effective permeability is derived. When

02 PETROLEUM↗

Electrostatic gating of ion transport in carbon nanotube porins: A modeling study

Carbon nanotube porins (CNTPs) are biomimetic membrane channels that demonstrate excellent biocompatibility and unique water and ion transport properties. Gating transport in CNTPs with external voltage could increase control over ion flow and selectivity. Herein, we used continuum modeling to probe the parameters that enable and further affect CNTP gating efficiency, including the size and composition of the supporting lipid membrane, slip flow in the carbon nanotube, and the intrinsic electronic properties of the nanotube. Our results show that the optimal gated CNTP device consists of a semiconducting CNTP inserted into a small membrane patch containing an internally conductive layer. Moreover, we demonstrate that the ionic transport modulated by gate voltages is controlled by the charge distribution along the CNTP under the external gate electric potential. The theoretical understanding developed in this study offers valuable guidance for the design of gated CNTP devices for nanofluidic studies, novel biomimetic membranes, and cellular interfaces in the future.

36 MATERIALS SCIENCE↗